来源论文: https://arxiv.org/abs/2604.21533v1 生成时间: Apr 24, 2026 10:09

0. 执行摘要

自2023年发现高压下La3Ni2O7（2222相）具有近80 K的高温超导电性（SC）以来，镍基超导体成为凝聚态物理的前沿焦点。然而，近期实验发现具有交替层状结构的1313相La3Ni2O7（由单层La2NiO4和三层La4Ni3O10组成）在压力下的最高Tc仅约3.6 K。本项工作系统地采用密度泛函理论结合动力学平均场理论（DFT+DMFT）以及随机相近似（RPA）方法，研究了该混合相的电子性质与配对对称性。研究表明，1313相中的单层（SL）子系统呈现莫特（Mott）绝缘体特征，而三层（TL）子系统保持金属态并主导超导电性。Tc的显著抑制归因于两个核心因素：一是TL子系统相对于体相La4Ni3O10发生了显著的空穴掺杂，导致配对强度减弱；二是SL子系统作为“弱链路”构成了S-N-S约瑟夫森结，限制了层间相位干涉，从而抑制了全局Tc。这一结论明确了2222相才是La3Ni2O7家族高温超导的真正载体。

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

1.1 核心科学问题

镍基超导体的核心谜题之一是结构与Tc的关联。La3Ni2O7存在两种主要的结构变体：

  1. 2222相（Bilayer）：双层NiO6平面结构，展示出~80 K的高Tc。
  2. 1313相（Hybrid）：由单层（SL）和三层（TL）交替堆叠。实验观察到其Tc仅为3.6 K，且伴随丝状超导特征。

科学问题：为什么化学组分完全相同、且都包含NiO2平面的1313相，其Tc会比2222相低一个数量级以上？这种混合结构内部是否存在轨道选择性的强关联效应？其超导配对对称性是否发生了改变？

1.2 理论基础：混合Ruddlesden-Popper（RP）相的电子结构

1313相La3Ni2O7在结构上是La2NiO4（n=1 RP相）和La4Ni3O10（n=3 RP相）的混合体。在~13 GPa的高压下，其对称性从Cmmm转变为P4/mmm。理论研究必须考虑Ni的$e_g$轨道（$d_{x^2-y^2}$和$d_{z^2}$）在不同层环境下的关联演化。SL层的Ni原子处于二维极限，易于发生Mott物理；而TL层的Ni原子通过顶角氧形成三层网络，具有更强的金属性。

1.3 技术难点

  • 层依赖的关联处理：传统的DFT无法处理Ni-d轨道上的强库仑排斥。必须引入DMFT来捕捉非扰动的多体效应，特别是轨道选择性的相变。
  • 混合结构的跨层迁移：SL和TL子系统之间存在电荷重新分布（Self-doping），这改变了各自的有效填充数，对超导配对至关重要。
  • 动态磁化率的计算：在强关联背景下，计算动量空间依赖的自旋磁化率并求解线性化Eliashberg方程，其数值稳定性要求极高。

1.4 方法细节

  1. DFT+DMFT架构：
    • 使用投影局部轨道（Wannier函数）构建低能Hamiltonian，包含Ni-3d和O-2p轨道。
    • 局部关联部分使用连续时间量子蒙特卡洛（CT-HYB）杂质求解器。库仑相互作用参数设为$U=5.0$ eV, $J_H=0.8$ eV。
    • 计算自能函数$\\Sigma(i\\omega_n)$和谱函数$A(\mathbf{k}, \\omega)$，以确定准粒子重整化因子$Z$。
  2. RPA分析：
    • 基于重整化后的Tight-Binding参数（包含准粒子权重$Z$缩放后的跃迁能$t_{eff} = Zt$）。
    • 构建多轨道哈伯德模型。自旋磁化率矩阵表达为： $$\\chi^{(s)}(\mathbf{q}, i\\omega) = [I - \\chi^{(0)}(\mathbf{q}, i\\omega)U^{(s)}]^{-1} \\chi^{(0)}(\mathbf{q}, i\\omega)$$
    • 求解配对矩阵方程 $V \\phi = \\lambda \\phi$，提取最大特征值$\\lambda$对应的对称性。

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

2.1 SL子系统的近绝缘行为

通过DMFT计算的自能函数$Im\\Sigma(i\\omega_n)$显示，SL层的$d_{z^2}$轨道在零能附近具有显著的截距，表明其处于Mott绝缘态的边缘。$d_{x^2-y^2}$轨道虽然具有一定的谱权重，但重整化因子极小，表现为“坏金属（Bad Metal）”。

  • 数据点：SL层$d_{z^2}$轨道的局域谱函数在费米能级处几乎为零。这意味着SL层无法提供有效的超导载流子，其物理角色转变为介电隔离层或约瑟夫森弱链路。

2.2 TL子系统的电荷转移与重整化

计算发现TL子系统发生了显著的空穴掺杂。相对于块体La4Ni3O10的电子填充，1313相TL层的总占据数减少了约0.17。

  • Table I 关键数据对比：
    • 1313 TL-inner: $n(d_{z2})=1.040, n(d_{x2-y2})=0.986$
    • 1313 TL-outer: $n(d_{z2})=1.032, n(d_{x2-y2})=1.055$
    • Bulk La4Ni3O10-outer: $n(d_{z2})=1.088, n(d_{x2-y2})=1.060$
    • 结论：TL子系统的空穴掺杂水平提高，导致费米面拓扑结构改变。

2.3 费米面（FS）与配对对称性

重整化后的TL子系统具有五个费米口袋：四个电子口袋（$\\alpha, \\beta, \\beta', \\epsilon$）和一个空穴口袋（$\\gamma$）。

  • ** nesting 效应**：自旋磁化率$\\chi^s(\mathbf{q})$在波矢$\\mathbf{Q}$处出现强峰，该波矢连接了$\\epsilon$和$\\gamma$口袋。
  • 配对特征：RPA计算给出最大的配对特征值$\\lambda$对应于$s^{\\pm}$波对称性。Gap函数在$\\epsilon$和$\\gamma$口袋上具有相反的符号，这与块体La4Ni3O10的物理机制一致。
  • Tc抑制量化：由于空穴掺杂使得$\\epsilon$和$\\gamma$口袋的大小失配（Size mismatch），导致 nesting 强度减弱。计算得到的$\\lambda$值在相同相互作用强度下显著低于La4Ni3O10，解释了Tc从30 K降至个位数的内因。

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

3.1 软件栈推荐

  1. DFT阶段：使用 WIEN2k 或 VASP 进行高压结构下的全电子计算。建议使用PBE泛函，并在13 GPa下进行晶格弛豫。
  2. Wannier投影：使用 Wannier90 软件。投影目标应包含Ni-3d（5个轨道）和O-2p（3个轨道）。对于TL子系统，需构建一个6轨道的低能模型（3层 $\\times$ 2个 $e_g$ 轨道）。
  3. DMFT阶段：推荐使用基于 TRIQS 框架的 DFTTools。求解器可选择 TRIQS/CTHYB。
    • 参数设置：温度设为 $\\beta=40 eV^{-1}$（约 290 K），混合步数不少于 40 步以确保自能收敛。
  4. RPA阶段：由于本研究涉及多轨道（6轨道），传统的单轨道代码不再适用。需编写自定义的Julia或Python脚本实现多轨道RPA公式。核心步骤是处理 $36 \\times 36$ 的磁化率张量。

3.2 关键参数复现（基于Table II）

为了复现1313相TL层的能带，需在Tight-Binding模型中填入以下跃迁参数：

  • $t_{x,O}^{[1,0]} = -0.220$ eV (外层 $d_{x^2-y^2}$ 最近邻跃迁)
  • $t_{z,IO}^{[0,0]} = -0.223$ eV (层间 $d_{z^2}$ 跃迁，显著低于非关联值的-0.650)
  • $t_z^{TL} = 1.12 \\times 10^{-4}$ eV (三层间的有效耦合，极弱)

3.3 开源仓库推荐

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

4.1 关键引用文献

  1. Sun et al., Nature 621, 493 (2023)：首次发现La3Ni2O7双层超导。这是所有后续工作的起点。
  2. Zhu et al., Nature 631, 531 (2024)：三层La4Ni3O10超导的发现，提供了重要的对比基准。
  3. Zhang et al., PRB 110, L180501 (2024)：讨论了空穴掺杂对La4Ni3O10超导的压制作用，本项工作延伸了这一观点。
  4. Nie et al., Nature 10.1038/s41586-026-10352-7 (2026)：薄膜实验揭示1313相不显超导，与本工作互为佐证。

4.2 工作局限性评论

  • SL子系统的简化：本工作将SL层视为完全不贡献配对的绝缘层。然而，在极高压下，SL层可能存在微弱的金属性，其对层间关联的动态反馈在RPA中被忽略了。
  • 约瑟夫森耦合的唯象处理：公式(4)中的约瑟夫森耦合强度$\\eta$是基于跃迁能平方的估计，缺乏基于Ginzburg-Landau理论的微观推导。更精确的模拟需要考虑由于热涨落导致的相位滑移效应。
  • 轨道占据数的不确定性：DFT+DMFT对占据数的计算依赖于投影窗口的选择，这可能导致约0.05-0.1个电子的系统误差，从而影响费米面嵌套（nesting）的精确度。

5. 其他必要补充：镍基超导的设计准则

通过对1313相的研究，我们可以提炼出设计高温镍基超导体的几条黄金准则：

  1. 避免“死层”干扰：1313相证明，在超导层（TL）之间引入Mott绝缘层（SL）虽然在结构上维持了RP相，但功能上构成了约瑟夫森屏障。为了获得高Tc，应追求超导层的连续性或具有极强穿透深度的耦合层。
  2. 电荷补偿机制：在混合相中，电荷倾向于从SL层流向TL层或反之。这种自掺杂效应往往会将体系推离最佳掺杂点。设计时需通过化学取代（如Sr取代La）来抵消这种效应。
  3. Ni-O-Ni键角：本工作再次验证了180°键角的重要性。在1313相中，尽管SL层存在畸变，但TL层内部在高压下趋于180°，这是形成强$\\sigma$-成键态和$s^{\\pm}$波配对的前提。
  4. 维度协同：高温超导需要准二维性来提高态密度，但也需要三维相干来建立全局超导态。1313相失败在三维相干性过弱（$\\eta$极小），这为未来寻找新镍基体系指明了方向：寻找那些具有适度层间跃迁能、能有效平衡局部配对与全局相位的结构。