Fermi-Hubbard模型中自旋纠缠与伪能隙起源的量子关联：基于冷原子模拟与DΓA理论的深度解析

来源论文: https://arxiv.org/abs/2605.31240v1 生成时间: Jun 01, 2026 14:18

0. 执行摘要

二维费米-哈巴德（Fermi-Hubbard）模型是现代凝聚态物理中理解强关联电子系统（如铜氧化物高温超导体）的基石模型。然而，由于强关联带来的非微扰特性和量子蒙特卡洛（QMC）中的符号问题，该模型在低空穴掺杂和低温区间的微观物理图景一直充满争议。其中，**伪能隙（Pseudogap, PG）**相的起源究竟是源于经典的超前磁涨落，还是量子力学自旋单态纠缠（如安德森提出的共振价键理论 RVB）的涌现，是凝聚态物理的核心未解之谜。

本项工作（由 Frederic Bippus, Thomas Chalopin, Karsten Held 等人合作完成）发表于2026年，通过**冷原子量子模拟实验（基于 $^6\text{Li}$ 原子二维光晶格）与先进的图规级非局域关联理论——动力学顶点近似（Dynamical Vertex Approximation, $\text{D}\Gamma\text{A}$）**相结合，首次从量子信息和纠缠熵的微观视角对这一问题给出了决定性的回答。研究表明：

1. 在超选择规则（Superselection Rules, SSR）约束下的自旋单态纠缠，精确地在伪能隙温度 $T^*$ 以下涌现。
2. 这种可操作的量子纠缠严格局限于最近邻（Nearest-Neighbor, NN）格点，在次近邻（Next-Nearest-Neighbor, NNN）及更远距离处完全消失。
3. 与之相反，经典的磁关联（如互信息和自旋相关函数）呈现更长程的衰减特征（相关长度约 $2.5$ 个晶格常数）。

这一发现有力地否定了完全基于经典自旋涨落的伪能隙理论，证实了伪能隙的微观图景必须包含短程量子自旋单态纠缠的形成，为高温超导配对机制的建模提供了极其严格的微观约束。

1. 核心科学问题、理论基础与技术细节

1.1 核心科学问题：伪能隙的微观物理本质

在铜氧化物超导体中，当温度低于 $T^*$ 且处于欠掺杂区时，费米面上的特定区域（反节点附近）会展现出谱权重的压制，即形成“伪能隙”。学术界对此主要存在两类微观解释：

• 经典/半经典涨落理论：认为伪能隙是由长程的、经典的自旋涨落（或电荷/条纹相涨落）热激发引起的，在这些理论中，系统可以通过经典的玻尔兹曼权重混合来描述，不依赖于量子相干相叠加。
• 量子纠缠/价键配对理论：以安德森的共振价键（Resonating Valence Bond, RVB）理论为代表，认为伪能隙本质上是局部自旋单态（Singlet）量子纠缠的预配对（Pre-pairing）状态。在 $T^*$ 处，电子自旋开始形成两体纠缠的单态（价键），但由于相位涨落，此时还没有建立全局的超导相位凝聚。

由于传统的实验探测手段（如中子散射、静态磁易受率）无法有效区分经典的自旋热涨落与纯粹的量子纠缠，上述争论持续了数十年。本工作引入了**量子纠缠度量（纠缠负度）**作为探针，从根本上过滤掉了经典关联，从而直接检验上述两种图景的对错。

1.2 理论基础：超选择规则（SSR）约束下的纠缠负度

在强关联费米子系统中，直接计算两格点间的降低密度矩阵（Reduced Density Matrix, RDM）会得到一个 $16 \times 16$ 的复矩阵。然而，由于物理测量必须遵守局部电荷守恒（Charge Conservation）和费米子奇偶性守恒（Parity Superselection Rule），这被称为超选择规则（SSR）。

在这些局部限制下，Alice 和 Bob 无法通过局部操作和经典通信（LOCC）来实现不同电荷或奇偶性扇区之间的量子相干干涉。因此，只有与 SSR 兼容的密度矩阵分量所承载的纠缠，才是操作上可利用的（Operationally Accessible）。本工作通过形式化投影，构造了 SSR 限制下的两体密度矩阵：

其中 $P_{\mathcal{N}_i}$ 是将第 $i$ 格点的局部电子数投影到特定扇区 $\mathcal{N}_i \in \{0, 1, 2\}$ 的投影算符。经过这一投影，密度矩阵被块对角化，非对角项连接不同局部粒子数扇区的元素全部归零（见 FIG. 1(c) 中的灰色区域）。

为了量化这一矩阵中的纠缠，研究采用了**对数负度（Logarithmic Negativity, $N$）**作为测度，它基于 Peres-Horodecki 判据（部分转置偏正性判据, PPT）：

其中 $\text{T}_j$ 代表对格点 $j$ 的局部密度矩阵进行部分转置，$\left\| \cdot \right\|_{\text{tr}}$ 为迹范数（Trace Norm）。当且仅当 $N > 0$ 时，表明系统存在无法用经典热混合描述的非平庸纠缠。

特别地，在自旋 SU(2) 对称性下，对于局限于两格点（单粒子占据）的 $2 \times 2$ 局域自旋子空间（对应 FIG. 1(c) 中的红色块），其纠缠判定可进一步化简。其部分转置后的特征值为：

由于对角元 $\rho_{7,7} \ge 0$ 恒成立，发生自旋纠缠的充分必要条件为：

根据算符重构，其中对角项 $\rho_{7,7}$ 对应空穴-双空穴背景下的非纠缠组分，而非对角项 $\rho_{8,9}$ 代表自旋翻转相干项，可重构为等时自旋相关函数：

这表明在冷原子气体显微镜中，只要高精度测得局部电子数关联和等时自旋-自旋关联，即可无偏地重构出这一纠缠测度。

1.3 技术难点与方法细节：DΓA 理论框架

普通的动力学平均场理论（DMFT）仅考虑自旋和电荷的局域动力学涨落，忽略了非局域的关联效应，因此无法正确描述具有强空间相干性的伪能隙。本研究采用了更高级的动力学顶点近似（Dynamical Vertex Approximation, $\text{D}\Gamma\text{A}$）。

DΓA 核心物理步骤：

1. 自洽 DMFT 求解：首先通过局域杂质求解器（如连续时间量子蒙特卡洛算法 CT-HYB，基于 w2dynamics 软件）自洽求解局域格林函数 $G_{\text{loc}}(i\omega_n)$ 和自能 $\Sigma_{\text{loc}}(i\omega_n)$。

2. 计算局域两粒子不可约顶点函数 $\Gamma_{\text{loc}}$：这是 DΓA 的精髓。$\Gamma_{\text{loc}}(\nu, \nu', \omega)$ 在低能下具有高度局域性，是连接局域波动与空间涨落的桥梁。

3. 求解 Bethe-Salpeter 方程（梯子近似）：利用 $\Gamma_{\text{loc}}$ 作为构建块，在自旋（Spin）和电荷（Charge）通道中分别引入空间涨落（梯子图级联），从而获得非局域的两粒子易受率 $\chi(\mathbf{q}, i\Omega_m)$：

   $$\chi_{\text{sp/ch}}(\mathbf{q}, i\Omega_m) = \left[ \chi_{0}^{-1}(\mathbf{q}, i\Omega_m) - \Gamma_{\text{loc, sp/ch}} \right]^{-1}$$

4. 重构自能 $\Sigma(\mathbf{k}, i\omega_n)$：非局域涨落（特别是反铁磁自旋涨落）通过 Schwinger-Dyson 方程反馈到单粒子自能中：

   $$\Sigma(\mathbf{k}, i\omega_n) - \Sigma_{\text{loc}}(i\omega_n) \propto \sum_{\mathbf{q}} V_{\text{eff}}(\mathbf{q}) G(\mathbf{k}+\mathbf{q})$$

5. 解析延拓（Analytic Continuation）：通过最大熵方法（Maximum Entropy Method, MEM）中的 “$\chi$-to-kink” 方案将 Matsubara 频率自能延拓至实频 $\Sigma(\mathbf{k}, \omega)$，进而获得谱函数 $A(\mathbf{k}, \omega)$，判定反节点处的能隙行为。

通过以上复杂的场论计算，研究者得以在热力学极限下定量模拟空间涨落，同时完整保留了量子多体纠缠信息。

2. 关键 Benchmark 体系与计算所得数据

2.1 模拟体系参数设置

• 模型：二维正方晶格 Fermi-Hubbard 模型（公式 4）。
• 相互作用强度：$U/t = 6.5$（处于中等偏强关联区，非常接近典型铜氧化物的物理区间）。
• 化学势/掺杂度：电子掺杂度 $\delta \in [-0.2, 0.2]$，其中 $\delta = 0$ 代表半满（Half-filling），负值代表空穴掺杂。
• 温度区间：$T/t \in [0.14, 0.6]$。在冷原子实验中，通过调控原子保持时间实现不同的热化温度。

2.2 伪能隙相边界的确定

利用 DΓA 实频谱函数 $A(\mathbf{k}, \omega)$（见 FIG. S1），研究者系统地观测了布里渊区中**节点（Node, $\mathbf{k} \approx (\pi/2, \pi/2)$）与反节点（Antinode, $\mathbf{k} \approx (0, \pi)$）**处的能谱演化：

• 在高温区，反节点和节点均展现出单峰金属行为，即典型的费米液体或准粒子状态。
• 随着温度降至临界温度 $T^*$ 以下，节点处依然保持高谱权重，但反节点处的谱权重出现明显分裂或压制，形成了局部能隙——即伪能隙。
• DΓA 确定的 $T^*(\delta)$ 边界与冷原子实验中测得的等时自旋易受率极大值温度 $\Theta(\delta)$ 高度吻合。

2.3 核心计算与实验数据对比

最近邻 SSR 纠缠负度 $N_{NN}$

在确定了伪能隙边界 $T^*$ 之后，研究团队绘制了纠缠相图（FIG. 1(a)）：

• 常规相（$T > T^*$）：无论是实验（白色菱形）还是 DΓA 计算（背景色梯度），纠缠负度 $N_{NN}$ 均严格为 0（黑色区域）。这意味着在金属和高温区，系统虽然存在自旋关联，但没有任何可操作的自旋相干纠缠，表现为纯经典的自旋随机热混合状态。
• 伪能隙相（$T < T^*$）：当且仅当温度降到 $T^*(\delta)$ 曲线以下时，$N_{NN}$ 陡然转为正值（红色和黄色区域）。实验测得的极大纠缠负度为 $N_{NN} = 0.0087(14)$，与 DΓA 理论预测高度相符。这一关键结果极其明确地建立了量子纠缠的出现与伪能隙 onset 之间的瞬时对应关系。

单态-三重态比例 $p_{\text{singlet}}/p_{\text{triplet}}$

为了从微观上剖析自旋纠缠的驱动力，论文对比了单态与三重态的概率比（FIG. 2）：

• 根据自旋单态投影算符，在 SU(2) 对称下，纠缠判据可重构为：

  $$\frac{p_{\text{singlet}}}{p_{\text{triplet}}} > 1$$

• 最近邻 (NN)：在 $T > T^*$ 时，该比例小于 $1.0$；在进入伪能隙区时，最近邻比例迅速跨越 $1.0$ 的纠缠阈值，在半满低温处达到最大值（~1.2）。

• 次近邻 (NNN)：对于次近邻（FIG. 2(b)），无论是实验还是理论，该比例在所有温度和掺杂下都严格小于 0.5。这表明次近邻格点之间连最低限度的纠缠门槛都无法达到，从而从实验和数值两方面证明了自旋单态纠缠是极度局域化的（仅限于一阶邻居）。

经典关联与量子纠缠的长度尺度对比（FIG. 3 & FIG. S3）

研究者进一步对比了表征总关联（包含经典和量子部分）的 SSR 互信息 $I$ 与自旋纠缠负度 $N$ 的空间特征：

• 纠缠负度：相关长度 $\xi_N = 0$，因为其在 $d \ge 2$ 处瞬间截断为零。
• 自旋相关函数 $\langle S_i^z S_j^z \rangle$：磁相关长度 $\xi_{\text{spin}}$ 在进入伪能隙时显著增加，在最低温下可达 $2.5$ 个格点。
• 互信息 $I$：呈现短程但显然超越最近邻的特征，其空间相关长度 $\xi_I \approx 1.5$ 格点。

这组数据提供了极为深刻的物理洞察：伪能隙状态下的量子相干（纠缠）是局域化在格点对上的，而导致经典宏观易受率增加的磁涨落则是长程的。

3. 代码实现细节与复现指南

为了便于科研人员复现该论文的 DΓA 计算结果，以下提供基于开源软件库的计算流指南。

3.1 核心开源软件包及链接

1. w2dynamics：基于连续时间量子蒙特卡洛（CT-QMC）的强关联多体解算器，用于计算 DMFT 自洽场及局域两粒子不可约顶点 $\Gamma$。
   • Github Repository: w2dynamics/w2dynamics
2. DGApy：用于执行动力学顶点近似（DΓA）梯子方程组求解的 Python 框架。
   • Github Repository (参考类似框架): TRIQS/dga 或同类 DΓA 社区代码
3. ana_cont：用于将松原频率（Matsubara Frequency）格林函数或自能解析延拓至实频的 Python 工具包。
   • Github Repository: krivenko/ana_cont

3.2 逐步复现流程（Step-by-Step Guide）

步骤 1：DMFT 基础计算与顶点提取

利用 w2dynamics 在正方晶格上求解单带 Fermi-Hubbard 模型（$U=6.5t$，$T=0.2t$，掺杂由化学势 $\mu$ 控制）：

# w2dynamics 配置文件 w2d.ini 关键参数示例
[general]
beta = 5.0          # 对应 T/t = 0.2
U = 6.5
mu = -0.5

[atoms]
atom1 = [0]

[cthyb]
measure_g2 = True   # 必须开启此项以测量两粒子格林函数 G2
N_iw = 1024         # 频率切片数
N_iw2 = 64          # 两粒子频率切片数

运行 w2dynamics 收集输出的 HDF5 文件，并提取两粒子局域顶点 $\Gamma_{\text{loc}}(\nu, \nu', \omega)$。

步骤 2：运行 DΓA 梯子方程求解器

导入局域顶点和格林函数至 DΓA 框架（例如基于 DGApy）：

import numpy as np
from dgap_solver import solve_dga

# 读取 DMFT 关联量
dmft_data = "w2dynamics_output.h5"

# DGA 参数设置
k_grid = (16, 16, 1) # 16x16 的二维晶格

# 执行 DGA 梯子计算，返回非局域格林函数与自能
Sigma_k_iw, G_k_iw = solve_dga(dmft_data, k_grid=k_grid, channels=['spin', 'charge'])

步骤 3：解析延拓寻找 $T^*$

利用 ana_cont 库对自能 $\Sigma(\mathbf{k}, i\omega_n)$ 进行实频延拓：

from ana_cont.solvers import MaxEnt

# 定义 Matsubara 频率和待延拓的自能数据
im_freq = np.array([ (2*n+1)*np.pi/beta for n in range(N_iw) ])
raw_data = Sigma_k_iw[antinode_idx, :] # 反节点处的自能

# 初始化最大熵解算器
me = MaxEnt(im_freq, raw_data, kernel_mode='fermionic')
result = me.solve(alpha_mesh=np.logspace(-3, 3, 100))

# 获得实频自能与谱函数
omega_real = result.grids['real']
A_antinode = -1.0/np.pi * np.imag(G_k_real)

如果 $A(\omega)$ 在 $\omega = 0$ 处展示出双峰结构且中间有凹陷，则代表系统已进入伪能隙区。

步骤 4：构建两格点降低密度矩阵与纠缠分析

根据论文公式（12）和（14），利用 DΓA 计算所得的两粒子格林函数，逆向重构基矢下的密度矩阵。对于两体系统（格点 $i$ 和 $j$），其密度矩阵块：

# 依据 SU(2) 对称性，重构 SSR 密度矩阵中的关键对角元与非对角元
rho_7_7 = measure_density_correlator(G2_k_iw) # 对应空穴关联
rho_8_9 = 2.0 * spin_correlator_DGA[d=1]     # 最近邻等时自旋关联

# 计算部分转置特征值
lambda_minus = rho_7_7 - np.abs(rho_8_9)

# 判定纠缠负度
if lambda_minus < 0:
    N_NN = np.log2(1.0 - 2.0 * lambda_minus)
    print(f"最近邻存在量子纠缠！N_NN = {N_NN:.5f}")
else:
    N_NN = 0.0
    print("系统处于无纠缠的经典关联状态。")

—%n%n## 4. 关键引用文献及局限性批判评论

4.1 关键引用文献及其在本工作中的作用

1. Anderson, P. W. (1973/1987) [1, 2]：提出了共振价键（RVB）理论。本工作最核心的科学论点（自旋单态作为伪能隙的微观基石）正是安德森 RVB 理论的直接微观表现。
2. Toschi, A. et al. (2007) [32] 与 Rohringer, G. et al. (2018) [33]：奠定了动力学顶点近似（DΓA）的场论数学框架，使超越局域近似的强关联计算成为可能。
3. Chalopin, T. et al. (2026) PNAS [43]：实验观测了 Fermi-Hubbard 体系中自旋-电荷尺度关联在伪能隙 onset 处的涌现。本工作在此基础之上，进一步将该实验数据提取并转换为量子信息学中的“纠缠负度”。
4. Peres, A. (1996) [80] 与 Horodecki, M. et al. (1996) [48]：定义了偏转置矩阵负特征值与纠缠的等价性（PPT 判据），提供了判定自旋纠缠的严格数学工具。

4.2 局限性批判评论

尽管本工作在凝聚态物理和量子信息的交叉领域取得了里程碑式的进展，但从量子化学和高精度多体计算的视角来看，仍存在以下值得深入探讨的局限性：

1. 次近邻跃迁 $t'$ 的缺失

• 问题：本研究在计算和实验中均设定次近邻跃迁振幅 $t' = 0$。然而，在实际的铜氧化物超导体中，$t'/t$ 通常在 $-0.15$ 到 $-0.35$ 之间。$t'$ 会显著破坏电子-空穴对称性，移动范霍夫奇点（Van Hove Singularity），并对超导配对区间和电荷条纹相产生决定性影响。
• 评述：虽然作者宣称对于较小的 $t'$，物理图景“定性上一致”（第7页），但没有给出定量支撑。增加 $t'$ 可能会显著增强局域电荷流，从而引入非自旋扇区的非平庸纠缠，使基于纯自旋单态（公式 10）的简明判定失效。

2. 最近邻纠缠度量极小（$N \sim 10^{-3}$）

• 问题：在伪能隙区间，实验和理论测得的 $N_{NN}$ 处于 $0.001 \sim 0.009$ 的极低水平。与理想的纯自旋 singlet（$N=1$）相比，这一纠缠度极度稀释。
• 评述：如此微弱的纠缠信号在实际量子器件中很难被有效利用，其作为“操作上可利用的纠缠（operationally accessible entanglement）”在实际的量子通信和计算协议中是否真的具备鲁棒性，仍然存疑。这也反映出系统在有限温度下，热噪声对量子相干性的压制极其严重。

3. 梯子近似（Ladder Approximation）对多粒子散射的简化

• 问题：本工作所采用的 DΓA 计算基于通道限制的梯子图求和，未能完整求解自洽的 Parquet 方程组。
• 评述：在伪能隙区，自旋和电荷通道之间存在强烈交叉耦合（Cross-channel coupling）。单纯依赖自旋通道的梯子求和，可能会过高估计反铁磁涨落对自能的贡献，从而将能隙开裂的临界温度 $T^*$ 稍微向高温区推移。未来有必要引入重正化群（fRG）或全自洽 Parquet DΓA 来进行基准校正。

5. 补充扩展：深入解析量子纠缠与高超机理

为了更全面地理解本项工作对量子物理和材料科学的深远影响，以下对部分核心物理机制进行补充探讨。

5.1 为什么是超选择规则（SSR）？从全费米子纠缠谈起

在凝聚态物理中，如果忽略超选择规则（SSR），直接计算所谓的“全费米子负度”（Fermionic Negativity, $N^F$），会得到一个完全不同的图像。在论文的补充信息 S4（第30-31页，FIG. S4）中，作者计算了 $N^F$：

• 全费米子纠缠 $N^F$ 的行为：在伪能隙边界 $T^*$ 以外的金属区，最近邻的 $N^F$ 依然显著大于零，甚至在零相互作用（$U=0$）的自由费米气体中也是正值。
• 物理机制差异：$N^F$ 包含了由于费米子**相干跃迁（Hopping）**导致的非平庸相位相干。这种纠缠是由自由电子在晶格间的平移对称性运动（波包扩展）带来的，它是金属输运的本征属性，与强关联多体物理中的伪能隙无关。更重要的是，在局部粒子数守恒的限制下，这种非 SSR 纠缠根本无法通过局部的物理测量加以提取或操作。
• SSR 的功绩：SSR 的引入就像一个完美的**“物理过滤器”**。它将由于电荷流跃迁带来的常规相干背景无情滤除，只保留了纯粹的、由于多体自旋交换作用（Exchange Coupling $J \sim 4t^2/U$）诱导的量子自旋单态相干。这也是为什么一旦应用 SSR，金属区的纠缠信号瞬间归零，而伪能隙相的信号则显露无遗。

5.2 安德森 RVB 理论的局部证实与整体悬决

本工作对于安德森著名的 RVB 理论给出了极其重要的局部微观实证，但同时也指出了其局限：

1. 局部的证实：RVB 理论预言伪能隙是由局域自旋单态（Valence bonds）组成的。本工作在 $T < T^*$ 下观测到的 $p_{\text{singlet}}/p_{\text{triplet}} > 1$ 以及纠缠只局限于 $d=1$ 格点，无可辩驳地证实了在伪能隙相中，电子确实以最近邻自旋单态价键的形式锁在了一起。这为“局部单态”的图像提供了最直接的量子信息学证据。
2. 整体的悬决：经典的 RVB 理论认为这些价键在晶格间是“共振（Resonating）”的，即价键在空间中存在宏观的量子叠加。然而，正如作者在第8页所坦承：“RVB 的共振特性无法从本工作测量和计算的静态降低密度矩阵中检测出来。”

要证实这种真正的“共振”特性，必须走向两体以上的、甚至多体全局纠缠熵（如拓扑纠缠熵、纠缠谱）的测量。这是下一步冷原子量子模拟和高精度数值计算必须攻克的高地。

5.3 展望：对未来量子化学与材料设计的启示

对于从事强关联量子化学（如过渡金属催化剂设计、重费米子体系计算）的学者，本工作的交叉视角带来了全新启示：

• 关联性的新指标：在传统的量子化学中，我们常用 D1/D2 诊断或轨道自然占据数来判定静态/动态电子关联。本工作证明了结合 SSR 的纠缠负度是一个极其灵敏的物理量，能够精确指涉出体系发生特定相变的关节点。这一指标未来可以引入到多组态自洽场（CASSCF）或密度矩阵重正化群（DMRG）的后处理中，用于评估过渡金属多核簇的自旋纠缠结构，指导精准磁性设计。