来源论文: https://arxiv.org/abs/2604.07185v1 生成时间: Apr 09, 2026 12:48

0. 执行摘要

自 2023 年在压力下发现 $T_c \sim 80\text{ K}$ 的双层镍氧化物 $\text{La}_3\text{Ni}_2\text{O}_2$ 以来，镍基超导体的物理机制成为凝聚态物理领域最受关注的焦点之一。特别是对于薄膜体系，如何在常压下实现并调控高 $T_c$ 超导电性具有重要的应用价值。本文所解析的论文《Perpendicular electric field induced $s^{\pm}$-wave to $d$-wave superconducting transition in thin film $\text{La}_3\text{Ni}_2\text{O}_7$》由苏州大学江密教授团队与北京理工大学范仰教授等合作完成。该工作采用动力学集群近似（DCA）结合连续时间辅助场（CT-AUX）量子蒙特卡洛（QMC）方法，系统研究了垂直电场对非对称双层两轨道 Hubbard 模型的影响。

核心发现包括：在无电场时，体系以 $d_{z^2}$ 轨道主导的 $s^{\pm}$ 波配对为主；随着垂直电场的增加，层间势能差导致电子在层间重新分布，抑制了 $s^{\pm}$ 波，并在中间强度电场下诱导出了由 $d_{x^2-y^2}$ 轨道主导的 $d$ 波配对，且其 $T_c$ 随电场呈现出明显的“圆顶”（dome-like）行为。这一发现不仅为实验上通过电场调控镍基超导薄膜提供了理论支撑，也深化了对多轨道关联体系中超导对称性竞争的理解。

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

核心科学问题

镍氧化物超导体的微观机制目前存在两种主流观点：一种认为 $d_{z^2}$ 轨道是超导的起源，电子通过强层间自旋反铁磁涨落配对，并通过轨道杂化获得相干性；另一种则认为 $d_{x^2-y^2}$ 轨道在超导中起主导作用，洪德耦合（Hund’s coupling）在轨道间传递磁关联。本文试图回答：在薄膜体系中，外部垂直电场能否作为一种有效的调控手段，改变体系的轨道占据和有效关联强度，从而在不同配对对称性之间切换？

理论基础：非对称双层两轨道 Hubbard 模型

为了描述 $\text{La}_3\text{Ni}_2\text{O}_7$ 的物理，研究者构建了一个包含 $d_{x^2-y^2}$ (简写为 $x$) 和 $d_{z^2}$ (简写为 $z$) 两个轨道的双层模型。其哈密顿量由单体项 $H_0$ 和相互作用项 $H_U$ 组成：

  1. 单体项 ($H_0$)：包含层内最近邻、次近邻及次次近邻跳符 $t_1, t_2, t_3$。特别地，包含了层间 $d_{z^2}$ 轨道的跳符 $t^\perp_{zz}$ 和层间/层内轨道杂化 $t_{xz}$。垂直电场的影响通过在底层引入一个额外的势能场 $\epsilon_b$ 来模拟，即 $E = (\epsilon_b - \epsilon_t)/e$，其中顶层势能 $\epsilon_t$ 设为 0。
  2. 相互作用项 ($H_U$)：采用标准的多轨道 Hubbard 形式，包含层内同轨道库仑排斥 $U$、层内异轨道库仑排斥 $U'$ 以及洪德耦合 $J$。参数设定遵循常规经验值：$U=4.0\text{ eV}$, $J=U/5=0.8\text{ eV}$，$U'=U-2J=2.4\text{ eV}$。

技术难点：强关联与费米子负符号问题

研究多轨道 Hubbard 模型最大的技术障碍在于其巨大的希尔伯特空间以及量子蒙特卡洛计算中的负符号问题（Sign Problem）。特别是在掺杂和低温下，平均符号值会呈指数级下降，导致计算无法收敛。为了解决这一问题，作者采用了动力学集群近似（DCA）。

方法细节：DCA 与 Bethe-Salpeter 方程

  1. DCA 框架：DCA 将无限大格点问题映射到一个嵌入在平均场环境中的有限集群上。相比于传统的单位置动力学平均场理论（DMFT），DCA 保留了集群内的空间相关性，这对于研究非 $s$ 波配对（如 $d$ 波）至关重要。本工作使用了 $N_c = 8$ (即 $4 \times 2$) 的集群，能够捕捉 $(π,0)$ 和 $(0,π)$ 等关键波矢点。
  2. CT-AUX 求解器：集群问题通过连续时间辅助场量子蒙特卡洛方法求解，这是一种无偏的数值方法，能够精确处理局域和短程关联。
  3. 超导性质的提取：通过解粒子-粒子通道的 Bethe-Salpeter 方程 (BSE)： $-\frac{T}{N_c} \sum_{K'} \Gamma^{pp}(K, K') \bar{\chi}^{pp}_0(K') \phi_\alpha(K') = \lambda_\alpha(T) \phi_\alpha(K')$ 其中 $\Gamma^{pp}$ 是不可约顶角函数，$\bar{\chi}^{pp}_0$ 是粗粒化的裸磁化率。当特征值 $\lambda_\alpha(T)$ 随温度降低接近 1 时，对应的配对对称性 $\alpha$ 发生超导转变。外推 $\lambda(T)=1$ 处的温度即为 $T_c$。

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

体系配置

计算针对三种典型的载流子浓度进行：

  • 未掺杂 (Undoped): $n=3.0$（每个 Ni 原子平均 1.5 个电子）。
  • 空穴掺杂 (Hole-doped): $n=2.8$。
  • 电子掺杂 (Electron-doped): $n=3.2$。

关键计算数据分析

1. 配对对称性的竞争 (Fig. 2 & Fig. 4)

  • $E = 0$ 时：在未掺杂和空穴掺杂情况下，$s^{\pm}$ 波对称性的 $T_c$ 最高。其特征向量（Fig. 3）显示配对主要发生在层间的 $d_{z^2}$ 轨道上，这与“层间自旋涨落驱动配对”的理论吻合。
  • 电场诱导转变：随着 $E$ 增加到 $0.5\text{ V}$ 以上，$s^{\pm}$ 波的特征值迅速下降，而 $d$ 波（由 $d_{x^2-y^2}$ 驱动）的特征值开始上升。在 $E=0.5\text{ V}$ 附近发生了对称性交叉。
  • d 波“圆顶”行为：$d$ 波的 $T_c$ 并非随电场单调增加，而是在 $E \sim 0.5 - 1.0\text{ V}$ 达到最大值，随后在高场强下下降。这表明电场不仅提供了层间反转对称性的破缺，还优化了有效填充。

2. 轨道占据与电子重新分布 (Fig. 5)

计算发现，垂直电场 $E$ 驱动电子从底层流向顶层。

  • 对于 $d_{x^2-y^2}$ 轨道，顶层的填充数 $n^t_x$ 随着 $E$ 增加而增加。数据表明，当 $n^t_x$ 接近 0.68 左右时，$d$ 波 $T_c$ 达到最高。这与单层 Hubbard 模型中 $d$ 波最强点在空穴掺杂区域的物理图像相呼应。
  • 层间密度差 $\Delta n_z$ 和 $\Delta n_x$ 的增大是导致 $s^{\pm}$ 波被抑制的主要原因，因为层间轨道能量失配破坏了层间配对的相干性。

3. 掺杂效应的影响

  • 空穴掺杂 ($n=2.8$)：展现了最高的总体 $T_c$，且 $d$ 波配对在较大电场范围内保持稳定。这解释了实验上通过 Sr 掺杂（空穴掺杂）调控超导性的潜力。
  • 电子掺杂 ($n=3.2$)：$d$ 波配对基本消失，且 $s^{\pm}$ 波 $T_c$ 随电场变化较小。这反映了电子掺杂会使费米面结构远离 $d$ 波配对的不稳定性区域。

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

核心软件栈：DCA++

虽然论文未明确指明具体的代码二进制文件，但根据作者团队（如江密教授、Maier 等）的学术背景，该工作通常基于 DCA++ 框架实现。DCA++ 是一个高性能的 C++/CUDA 代码库，专门用于在大规模并行机上解决关联电子模型。

计算复现指南

  1. 格点参数设置：
    • 参考 Table I 设置跳符：$t^x_1=-0.445, t^z_1=-0.131$ 等（单位 eV）。
    • 注意设置 $t^\perp_{zz} = -0.503$。本研究为了简化模型，忽略了 $t^\perp_{xx}$。
  2. 电场引入： 在自能计算中，底层轨道的局域势能项添加 $\epsilon_b$，其范围从 $0.0$ 到 $1.5\text{ V}$。
  3. QMC 采样策略：
    • 使用 CT-AUX 算法。由于存在多轨道和双层结构，测量的内能和格林函数需要分轨道进行存储。
    • 松原频率（Matsubara frequency）截断建议至少取 32-64 个点以保证低能物理的准确性。
  4. BSE 求解流程：
    • 第一步：运行 QMC 获得收敛的集群自能 $\Sigma(\mathbf{K}, i\omega_n)$。
    • 第二步：计算不可约顶角 $\Gamma$。这需要测量四点关联函数，计算量巨大。
    • 第三步：构建 BSE 矩阵。矩阵维度为 $(2 \text{ orbitals} \times 2 \text{ layers} \times N_c \times N_{freq})^2$。对于 $N_c=8$，$N_{freq}=32$，矩阵大小约为 $2048 \times 2048$。

相关开源资源

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

关键引用文献

  1. Sun et al., Nature 621, 493 (2023): 镍基超导体的发现之作，确定了高压下 $T_c \sim 80\text{ K}$ 的实验事实。
  2. Maier et al., Rev. Mod. Phys. 77, 1027 (2005): DCA 方法的基础性综述，详细阐述了集群近似的理论框架。
  3. Shao et al., Nature Communications (2026) (预印本或预期发表): 本文多次提及该工作关于电场诱导超导的 SBMF 分析，本文的 DCA 结果是对其在强关联限下的重要验证。
  4. Ushio et al., arXiv:2506.20497: 关于薄膜 $\text{La}_3\text{Ni}_2\text{O}_7$ 结构的理论构建。

局限性评论

尽管本文在数值精确性上具有显著优势，但仍存在以下局限性：

  1. 集群尺寸限制：由于负符号问题的限制，$N_c=8$ 仍然是一个相对较小的集群。这可能导致它对某些长程超导涨落的捕捉不够充分，或者导致 $T_c$ 的外推值略高于真实热力学极限。
  2. 温度标度问题：计算所得的 $T_c$ 约为 $0.04 - 0.08\text{ eV}$（相当于数百至上千开尔文），远高于实验观测到的 $40-80\text{ K}$。这是因为 QMC 在极低温下无法运行，研究者通常在较高温度下进行计算并外推。这在量级上是可以理解的，但对于定量预测实验电场临界值仍有挑战。
  3. 模型简化：忽略了层间 $d_{x^2-y^2}$ 的直接跳符 $t^\perp_{xx}$。虽然该值通常较小，但在精确竞争能标时可能起到微调作用。
  4. 轨道空间完整性：未包含氧的 $p$ 轨道。虽然有效 Hubbard 模型捕捉了低能物理，但镍基体系中 $p-d$ 杂化的电荷转移特性在某些压力区间可能很重要。

5. 补充：物理直觉与未来研究方向

物理直觉：为何电场能驱动对称性转变？

从物理图像上看，$s^{\pm}$ 波依赖于层间 $d_{z^2}$ 轨道的“配对粘合剂”（由 $t^\perp_{zz}$ 介导的反铁磁交换）。当垂直电场产生时，两层 Ni 轨道的能级不再简并。这就像是在一个联通器两端施加压力差，电子被迫向一侧集中。这种层间对称性的破缺直接削弱了层间相干配对的能力。

与此同时，电场将顶层（Top layer）的填充推向了一个特定的“黄金区域”。对于 $d_{x^2-y^2}$ 轨道，其本身就具有很强的内禀 $d$ 波配对不稳定性。电场就像是一个手柄，通过调节局部化学势，使得顶层的 $d_{x^2-y^2}$ 轨道进入了类似空穴掺杂铜氧化物的超导最强区域。因此，我们看到了 $s^{\pm}$ 波的谢幕与 $d$ 波的登场。

跨学科意义：镍基与铜基超导的桥梁

这项研究的一个深远意义在于它建立了镍基超导体与铜基超导体之间的逻辑关联。镍氧化物原本被认为是“不完美的铜氧化物”，因为其多轨道特性复杂化了物理图像。然而，本文展示了通过电场这种简单的外部激励，我们可以将多轨道的镍基体系“转化”为一种更接近铜基（单轨道 $d$ 波主导）的物理状态。这为统一高温超导理论提供了新的视角。

未来展望

  1. 实验验证：利用离子液体栅极（Ionic Liquid Gating）或铁电层薄膜异质结来产生巨大的垂直电场，观察 $\text{La}_3\text{Ni}_2\text{O}_7$ 薄膜的超导能隙对称性变化。
  2. 自旋涨落能谱：未来可以使用 DCA 计算电场下的动态自旋易感性 $\chi(\mathbf{q}, \omega)$，以确认 $d$ 波配对的具体磁性起源。
  3. 三层及以上体系：推广到 $\text{La}_4\text{Ni}_3\text{O}_{10}$ 等体系，研究电场在更复杂层状结构中的非线性调控效应。

关键词： 强关联电子系统、数值重整化群、高 $T_c$ 机理、层间相互作用、非费米液体行为。