来源论文: https://arxiv.org/abs/2603.11037v1 生成时间: Mar 11, 2026 23:14

在光晶格中实现 Emery 模型：铜氧化物与镍氧化物高温超导量子模拟的新范式

0. 执行摘要

高温超导体的微观起源，特别是铜氧化物（Cuprates）和镍氧化物（Nickelates）的物理机制，一直是凝聚态物理中最为核心且极具挑战性的课题。尽管单能带费米-哈伯德模型（Single-band Fermi-Hubbard model）在描述某些特性上取得了成功，但越来越多的证据表明，要捕获高温超导的完整图景，必须回归到其更本质的多能带描述——Emery 模型（三能带哈伯德模型）。

本文基于 Hannah Lange 等人的最新研究，深度解析了一种在超冷原子光晶格平台上实现三能带 Emery 模型的可行方案。该方案巧妙利用了 Lieb 晶格几何结构，并通过干涉激光束的偏振调节实现了关键的电荷转移能（Charge-transfer energy, $\Delta_{pd}$）的可调性。通过大规模密度矩阵重整化群（DMRG）数值模拟，研究证明了该平台能够有效模拟从铜氧化物到镍氧化物的不同参数区间，并观测到了张-赖斯单态（Zhang-Rice singlets）的形成以及空穴/电子掺杂在磁性关联上的显著不对称性。这一突破为在量子模拟器上超越经典计算极限、探索真实材料超导机制奠定了坚实基础。

1. 核心科学问题，理论基础，技术难点与方法细节

1.1 核心科学问题：为何必须转向 Emery 模型？

在过去的四十多年里，物理学家试图将铜氧化物层的物理简化为单能带模型。该模型假设物理学主要由铜原子的 $d_{x^2-y^2}$ 轨道主导。然而，铜氧化物的基本单元实际上是由铜（Cu）和氧（O）组成的 $\text{CuO}_2$ 平面。氧的 $2p$ 轨道与铜的 $3d$ 轨道之间存在强烈的杂化。

单能带模型忽略了以下关键因素：

  1. 电荷转移性质：铜氧化物属于电荷转移绝缘体（Charge-transfer insulators），而非纯粹的 Mott 绝缘体。这意味着掺杂的空穴主要进入氧的 $p$ 轨道，而非铜的 $d$ 轨道。
  2. 张-赖斯单态：空穴在氧轨道上与铜旋转产生的局域单态，其动力学可能包含单能带模型无法完全捕捉的下一近邻跳跃项（$t'$）或密度辅助跳跃项。
  3. 镍氧化物的差异：新发现的无限层镍氧化物超导体具有更大的电荷转移能 $\Delta_{pd}$，使得其物理行为介于哈伯德极限和 Emery 极限之间。单能带模型难以统一描述这两类材料。

1.2 理论基础：三能带 Emery 模型哈密顿量

论文研究的核心是定义在 Lieb 晶格上的 Emery 模型。其空穴图象下的哈密顿量为：

$$\hat{H} = - t_{pd} \sum_{\langle ij \rangle, \sigma} (\hat{d}_{i\sigma}^\dagger \hat{p}_{j\sigma} + \text{h.c.}) - \sum_{\langle \langle jj' \rangle \rangle, \sigma} t_{jj'}^{pp} (\hat{p}_{j\sigma}^\dagger \hat{p}_{j'\sigma} + \text{h.c.}) + \Delta_{pd} \sum_{j, \sigma} \hat{p}_{j\sigma}^\dagger \hat{p}_{j\sigma} + \sum_{\nu \in \{p,d\}} U_\nu \sum_i \hat{n}_{i\uparrow}^\nu \hat{n}_{i\downarrow}^\nu$$

其中：

  • $\hat{d}_{i\sigma}$ 和 $\hat{p}_{j\sigma}$ 分别是铜（d 节点）和氧（p 节点）上的空穴湮灭算符。
  • $t_{pd}$ 是 $d$ 与 $p$ 轨道间的跳跃项。
  • $t_{jj'}^{pp}$ 包含氧原子间的近邻（$t_{pp}$）和次近邻（$t'_{pp}$）跳跃。
  • $\Delta_{pd} = \epsilon_p - \epsilon_d$ 是电荷转移能，控制着空穴在 $p$ 轨道和 $d$ 轨道间的分布。
  • $U_d$ 和 $U_p$ 分别是 $d$ 和 $p$ 轨道的原位 Hubbard 排斥能。

1.3 技术难点：如何在实验中构造受控的 $\Delta_{pd}$？

传统的超冷原子实验通常实现的是简单的哈伯德模型。要在光晶格中实现 Emery 模型，面临两大挑战：

  1. 几何结构的复杂性：需要构造一个 Lieb 晶格，即在正方形晶格的边缘中点添加额外的节点。
  2. 能量偏移的可调性：必须能够精确且稳定地调节 $p$ 轨道和 $d$ 轨道之间的能量差 $\Delta_{pd}$。此前的方法如使用数字微镜器件（DMD）投影势场往往面临相位稳定性或空间分辨率的问题。

1.4 方法细节：偏振可调的干涉光晶格方案

作者提出了一种创新的实验方案，利用“领结型”（bow-tie）几何结构的单光束被动相位稳定光晶格。其核心在于控制激光在两个垂直路径上的偏振关系：

  • 势场构造：光轴沿 $x$ 方向的激光偏振设为 $V$（垂直）。在穿过原子云并经过半波片（HWP）后，其在回程（$y$ 方向）的偏振被旋转一个角度 $\theta$。这产生的势场为： $$V(x,y)/V_0 = -\frac{1}{2} [\cos^2(k_Lx) + \cos^2(k_Ly)] - \sin\theta \cos(k_Lx) \cos(k_Ly)$$
  • 偏振调节 $\Delta_{pd}$：通过调节角度 $\theta$，可以连续改变 $d$ 轨道（位于格点中心）和 $p$ 轨道（位于连接中心）之间的能量差。当 $\theta=0$ 时，对应标准的正方形晶格；增加 $\theta$ 则会产生所需的偏移量。
  • Lieb 晶格的实现：利用 DMD 投影一组排斥性的势场阵列，有选择性地抬高每四个“铜”位点中三个的能量，从而在剩余位点上留下完美的 Lieb 几何结构。这种方案结合了被动相位稳定性和 DMD 的灵活性。

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

为了验证实验方案的可行性，作者使用了基于矩阵乘积态（MPS）的 DMRG 方法 对 $12 \times 2$ 单元格（共 72 个位点）的柱状几何体系进行了深入计算。计算针对铜氧化物制度（$\Delta_{pd}/t_{pd} = 3.5$）和镍氧化物制度（$\Delta_{pd}/t_{pd} = 9.0$）展开。

2.1 关键计算参数

  • 相互作用强度：$U_d/t_{pd} = 8.0, U_p/t_{pd} = 3.0$。
  • 跳跃项：$t_{pp}/t_{pd} = 0.1, t'_{pp}/t_{pd} = -0.2$。这些值通过拟合实际光晶格能带结构获得，非常接近真实材料参数。
  • 掺杂水平：从电子掺杂（$\delta < 0$）到空穴掺杂（$\delta > 0$），其中 $\delta = N_h / (L_x L_y)$。

2.2 核心物理观测数据

A. 自旋关联函数的非对称性（Figure 3a）

  • 铜位点关联 ($C^S_{dd}$)：在两种材料制度下，近邻铜位点间的反铁磁（AFM）关联都随掺杂浓度的增加而减弱。这与单能带哈伯德模型的结果一致。
  • 铜-氧关联 ($C^S_{dp}$)：这是 Emery 模型的独特观测项。在电子掺杂侧，$C^S_{dp}$ 几乎为零。然而在空穴掺杂侧（特别是铜氧化物制度下），$C^S_{dp}$ 随 $\delta$ 线性增强且呈负值（AFM 关联）。这有力证明了空穴主要占据 $p$ 轨道并与相邻的 $d$ 轨道自旋形成 AFM 关联。

B. 张-赖斯单态（ZR Singlet）的形成（Figure 3b）

  • 键占据数 ($n_{bond}$)：作者定义了占据的 $d-p-d$ 键。结果显示，在空穴掺杂的铜氧化物中，$n_{bond} \approx \delta/2$。这正是 ZR 单态图象的直接体现：每一个注入的空穴都会在一个 $d$ 位点周围的对称氧轨道上形成一个单态。
  • 磁化强度 ($M_{bond}$)：在 ZR 单态区域，$M_{bond}$ 稳定在 $0.75$ 左右。由于单态贡献为 0，剩余的一个 $d$ 位点自旋贡献 $0.75$（即 $S(S+1)$，其中 $S=1/2$），这再次从微观层面验证了 ZR 单态的稳固性。

C. 镍氧化物与铜氧化物的跨越

  • 随 $\Delta_{pd}$ 增加（从 3.5 变为 9.0），ZR 单态的形成受到明显抑制。镍氧化物显示出更强的 Mott 绝缘体特征，空穴开始更多地进入 $d$ 轨道，而非仅仅局限于 $p$ 轨道形成的 ZR 单态中。这一发现解释了为何镍氧化物的物理性质与铜氧化物存在细微但关键的差异。

3. 代码实现细节，复现指南与开源链接

该研究的数值模拟部分具有极高的技术含量，以下是针对科研人员的复现指南。

3.1 核心软件包：SyTen

  • 名称：SyTen (Symmetry-protected Tensor networks)
  • 主要开发者：C. Hubig, F. Lachenmaier 等 (LMU Munich/Max Planck)
  • 链接：https://www.syten.eu/
  • 特点：该工具包专门针对具有非阿贝尔对称性（如 SU(2)）或多个 U(1) 对称性的强关联体系进行了优化，是目前最高效的张量网络库之一。

3.2 DMRG 实现细节

  • 对称性利用：计算显式地保持了 $U(1) \times U(1)$ 对称性，分别对应总粒子数守恒和总自旋 $S_z$ 守恒。对于偶数（奇数）粒子数，计算在 $S_z = 0 (0.5)$ 扇区进行。
  • 收敛策略：
    • 采用多阶段扫频。初始阶段截断维度 $\chi$ 从 256 开始，逐步增加至 4096。
    • 最后阶段使用 two-site DMRG 算法（2S-DMRG）以确保跳出局部能量极小值。
    • 截断误差：保持在 $10^{-8}$ 以下，确保了长程关联的准确性。
  • 抽样技术：利用 MPS 的 Metropolis 抽样生成 2000 个“快照”（Snapshots），模拟实验中量子气体显微镜（Quantum Gas Microscope）捕获的原位图像。

3.3 实验参数复现（针对超冷原子平台）

  • 原子种类：建议使用 $^6\text{Li}$ 原子，因为其具有优异的 Feshbach 共振可调性。
  • 激光波长：$\lambda_L = 1064$ nm。
  • 反冲能量：$E_r = h \times 7.3$ kHz。
  • 关键拟合结果：对于铜氧化物区域（$\theta = 0.025$ rad），轨道跳跃能 $t_{pd} \approx 0.31 E_r$。这保证了实验能够在当前的亚纳开尔文（sub-nK）温度范围内观察到 AFM 关联。

4. 关键引用文献与局限性评论

4.1 关键引用文献

  1. [2] V. J. Emery, Phys. Rev. Lett. 58, 2794 (1987)：Emery 模型的开创性论文，提出了三能带描述。
  2. [3] F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 (1988)：提出了张-赖斯单态概念，建立了三能带到单能带的映射基础。
  3. [25] M. Lebrat et al., arXiv:2404.17555 (2025)：在实验中首次实现多能带 Lieb 晶格，为本方案提供了硬件基础。
  4. [19] M. Xu et al., Nature 642, 909 (2025)：展示了目前超冷原子体系能达到的最低温度（约 $0.05 T/t$），证明了观测 AFM 关联的可能性。

4.2 局限性评论

尽管该方案非常精妙，但在实际执行中仍存在若干挑战：

  1. 温度限制：尽管计算表明关联在当前实验温度下是可见的，但要真正模拟“超导”相（即观测到长程的对关联），目前的超冷原子制冷技术仍需进一步突破。当前的关联仍处于短程或中程阶段。
  2. 次近邻项的固定性：在光晶格中，$t_{pp}$ 和 $t'_{pp}$ 的比例很大程度上由 Wannier 函数的形状决定，很难独立于 $t_{pd}$ 进行大幅度调节。虽然本文证明了当前参数足以捕获核心物理，但对于某些极端材料体系，这种局限性可能会显现。
  3. 有限尺寸效应：DMRG 计算是在 $12 \times 2$ 的小系统上完成的。尽管这足以研究局域的 ZR 单态形成，但在研究长程序（如 Stripe 相或超导序）时，需要更大规模的二维模拟，这超出了经典计算机的负荷，也正是量子模拟器的用武之地。
  4. 三体项缺失：Emery 模型中可能存在的密度辅助跳跃项（Density-assisted tunneling）在当前势场方案中尚未被显式包含，这在极高精度模拟中可能是一个缺失环节。

5. 补充解析：从 Emery 模型到 Altermagnetism（分交磁性）

除了超导机制外，本文提出的 Lieb 晶格 Emery 模型还意外地为研究新兴的**分交磁性（Altermagnetism）**打开了大门。

5.1 什么是 Altermagnetism？

分交磁性是除了铁磁和反铁磁之外的第三种磁序。它在实空间没有净磁化强度（类似于 AFM），但在动量空间具有自旋分裂的能带（类似于 FM）。这种性质由晶格对称性与自旋序的特殊耦合产生。

5.2 在本方案中的体现

在 Lieb 晶格中，由于 $d$ 轨道和 $p$ 轨道占据不同的对称性中心，当系统处于反铁磁相时，调节 $\Delta_{pd}$ 会破坏原本的某种空间-自旋反转对称性。作者在文中暗示（参考引用 [34]），这种设置可以用来直接模拟和观测分交磁性。这是该量子模拟方案的一个重要“副产品”，极大地扩展了其研究范畴。

5.3 规范变换（Gauge Transformation）的深度意义

在 Appendix A 中，作者详细讨论了规范变换。在真实的 $\text{CuO}_2$ 平面，由于轨道对称性，$d$ 和 $p$ 之间的跳跃系数符号是交替的。为了在光晶格中简化实现，作者采用了一种规范变换将所有跳跃项转为负值。这在密度关联测量上没有影响，但会改变动量空间的准动量（$k \to k + (\pi, \pi)$）。这种处理在数值计算中极其常见，但在解释实验动量空间快照（如 ARPES 类测量）时必须格外小心。

5.4 总结与展望

这项工作的核心价值在于将“材料特异性”引入了量子模拟。过去我们讨论哈伯德模型，往往是在讨论一种抽象的物理。而现在，通过调节 $\Delta_{pd}$，我们可以明确地说：“现在我们在模拟铜氧化物”，或者“现在我们在模拟镍氧化物”。这种从“通用模型”到“专用材料模拟”的转变，标志着量子模拟技术进入了 2.0 时代。未来，结合更高分辨率的自旋/密度空间关联测量，我们有望在光晶格中真正解开高温超导的百年谜题。