Beyond the conventional Emery model: crucial role of long-range hopping for cuprate superconductivity

利用动力学顶点近似，本研究证明，尽管传统的埃默里模型能够捕捉铜氧化物超导性的定性特征，但引入超出标准三个参数的长程跃迁参数对于获得定量准确的相图和正确的d波序参量至关重要。

要点

想象一下，你正试图烤出一个完美的纸杯蛋糕，但你不是在和面粉与糖打交道，而是在试图理解某些材料（称为铜氧化物）如何毫无电阻地传导电流。这种现象被称为超导性。

几十年来，科学家们一直使用一种特定的“配方”来模拟这些材料，即Emery 模型。可以将这个模型想象成一座城市的简化地图。在这座城市里，有铜制的“房子”和氧制的“房子”。电子（即“人们”）在房子之间跳跃。

这张地图的传统“配方”只允许人们跳到他们的紧邻邻居（隔壁的铜房或氧房）。这就像在说：“你只能走到你正隔壁的那栋房子。”

旧地图的问题

这篇论文的作者是 Eric Jacob 和 Karsten Held 领导的团队，他们决定使用一种非常强大的计算机模拟（称为“动力学顶点近似”方法）来测试这张旧地图。他们发现，旧地图遗漏了一个至关重要的要素。

在现实世界中，人们不仅仅走到隔壁；他们也可以走到隔两栋的房子，甚至如果路径畅通，可以跳过几栋房子。在物理学中，这些被称为长程跳跃。

当科学家们使用旧有的、受限的地图（仅允许紧邻跳跃）时，模拟未能产生正确类型的超导性，尤其是在材料被“掺杂”（混入额外的电子或空穴）到实际铜氧化物通常表现最佳的水平时。这就像试图只用一半的配料来烤蛋糕；结果根本无法正常膨胀。

新发现

团队意识到，为了得到“完美的蛋糕”（正确的超导行为），他们需要在地图中加入长程跳跃。他们必须让电子跳跃到更远的房子，而不仅仅是紧邻的邻居。

当他们加入这些额外的跳跃时，发生了以下情况：

1. “穹顶”出现了：在超导研究中，科学家们在图表中寻找一个“穹顶”形状。这个穹顶显示了超导性表现最佳的条件范围。旧地图产生了一个微小、狭窄的穹顶，与现实不符。而带有长程跳跃的新地图，产生了一个巨大、健康的穹顶，其外观与真实铜氧化物材料中观察到的超导行为完全一致。
2. “有序”变得合理：旧地图为电子配对方式（称为“序参量”）创造了一种奇怪、凹凸不平的模式。这就像一场舞蹈，舞伴们以一种奇怪的节奏互相踩脚。新地图则创造了一种平滑、经典的"d 波”舞蹈模式，这正是科学家在这些材料中预期看到的。

为何重要（根据论文观点）

该论文指出，长期以来，科学家们一直使用一种物理学的“简化”版本，它在粗略估算时或许还行，但在需要精确数值时就会失效。

• 旧方法：就像用一张手绘的城市草图来规划地铁系统。它能把握大体概念，但列车会相撞，因为地图遗漏了长长的隧道。
• 新方法：就像使用一种高科技 GPS，它考虑了所有可能的路线，包括长距离路线。这使得模拟能够精确预测超导性发生的时间和地点。

核心结论

作者们得出结论：如果你想准确描述这些超导材料的工作原理，你就必须包含电子的长距离跳跃。忽略它们会导致关于材料何时变为超导态以及如何行为的错误预测。他们并没有发明一种新的超导体或新的医疗设备；他们只是修正了我们用来理解现有材料的数学“地图”，表明旧地图遗漏了至关重要的道路。

技术摘要

技术摘要：超越传统 Emery 模型

问题陈述
Emery 模型是一个包含铜 $d_{x^2-y^2}$ 轨道和面内氧 $p_x, p_y$ 轨道的三带模型，是描述铜氧化物超导体的标准框架。然而，该模型的传统应用通常将跃迁参数限制为最近邻铜 - 氧 ($t_{pd}$) 以及最近邻和次近邻氧 - 氧 ($t_{pp}, t'_{pp}$) 相互作用。先前利用这种“传统”参数集的研究难以定量复现铜氧化物的超导相图，往往依赖于受限于有限尺寸效应的小团簇尺寸，或依赖于可能缺乏必要材料真实性的单带 Hubbard 模型。一个关键的未决问题是：传统的跃迁参数是否足以描述 $d$ 波超导穹顶和正确的序参量对称性，特别是在较高的空穴掺杂浓度下。

方法论
作者采用动力学顶点近似 (D$\Gamma$A) 来研究 Emery 模型。作为动力学平均场理论 (DMFT) 的图解扩展，D$\Gamma$A 允许在精细动量网格（高达 $130 \times 130$）上处理非局域关联，有效地模拟无限系统，并避免了在大团簇上进行量子蒙特卡洛模拟时伴随的符号问题。

本研究利用了一个特定的材料实现：$\text{CaCuO}_2$。非相互作用哈密顿量源自使用广义梯度近似 (GGA) 的密度泛函理论 (DFT) 计算，并投影到最大局域化 Wannier 函数上。相互作用强度 $U = 6.5$ eV 通过约束随机相位近似 (cRPA) 计算得出。作者比较了三种不同的模型：

1. 传统 Emery 模型：仅包含 $t_{pd}$、$t_{pp}$ 和 $t'_{pp}$。
2. 共价 Emery 模型：源自先前文献中为 $\text{Bi}_2\text{Sr}_2\text{CaCu}_2\text{O}_8$ (Bi-2212) 推导的传统参数集。
3. 全 d-p 模型：包含源自 DFT-Wannier 投影的所有跃迁元素，延伸至最近邻和次近邻之外。

通过在 D$\Gamma$A 的阶梯形式中求解线性化 Eliashberg 方程来识别超导不稳定性，并提取超导本征值 $\lambda_{SC}$。$\lambda_{SC} \to 1$ 的值指示临界温度 $T_c$。

主要贡献与结果

• 传统跃迁的不足：作者证明，当传统 Emery 模型向下折叠为有效单带 Hubbard 模型时，会显著低估有效跃迁参数。具体而言，向下折叠模型中的次近邻跃迁 $t'$ 过小（例如，对于 $\text{CaCuO}_2$，$|t'/t| \approx 0.07$），与真实值（$|t'/t| \approx 0.15 - 0.25$）相比存在差距。这种缺陷源于传统模型忽略了全 d-p 哈密顿量中物理存在的长程跃迁路径。
• 能带结构修正：在全 d-p 模型中包含长程跃迁参数，正确地将范霍夫奇点 (van Hove singularity) 移至较低能量，使能带结构与 DFT-Wannier 结果一致。相比之下，传统模型将范霍夫奇点置于过于接近费米能级的位置。
• 超导相图：
  • 传统模型：无法在空穴掺杂 $\delta \geq 0.15$（最佳掺杂范围）时显示出超导不稳定性。它产生一个非常狭窄的超导穹顶（9%–14% 掺杂），并预测出一个具有反节点处谱权重抑制的强调制 $d$ 波序参量。这种调制归因于范霍夫奇点位置不当导致的伪能隙过早形成。
  • 全 d-p 模型：成功复现了与铜氧化物现象学一致的宽超导穹顶，跨度约为 7% 至 22% 的空穴掺杂。在 $\delta = 0.15$ 时，观察到了清晰的超导不稳定性。由此产生的 $d$ 波序参量是常规的，没有受限模型中看到的虚假调制。
• 长程跃迁的作用：该研究确立，长程跃迁参数（超越三个传统参数）在定量上是获得正确相图和序参量所必需的。这些参数的量级与 Hubbard 模型中的 $t'$ 项相当，而后者已知对超导性至关重要。

意义
该论文声称，虽然传统 Emery 模型在定性上是有用的，但在描述低掺杂以外的铜氧化物超导性时，其在定量上是不足的。作者认为，将计算限制在三个传统跃迁参数会导致关于超导穹顶范围和序参量性质的错误预测。通过利用源自第一性原理的全 d-p 模型和 D$\Gamma$A 方法，这项工作提供了一种更贴近材料真实的铜氧化物描述，自然地包含了氧轨道和长程跃迁效应。这种方法解决了关于范霍夫奇点位置及其导致的伪能隙行为的差异，为更准确地理解高温超导性提供了一条途径，而无需诉诸经验参数拟合。