来源论文: https://arxiv.org/abs/2604.17812v1 生成时间: Apr 21, 2026 12:37

0. 执行摘要

近年来，镍氧化物（Nickelates）超导体的发现为高温超导研究开辟了除铜氧化物之外的新战场。特别是 La$_3$Ni$_2$O$_7$ 和 La$_4$Ni$_3O_{10}$ 等多层 Ruddlesden-Popper 相体系，在压力下展现出的超导特性引起了物理学界的极大关注。 La$_4$Ni$_3O_{10}$ 作为一种三层结构体系，其 Ni 离子的平均价态为 $d^{7.33}$，这在电子填充上非常接近 1/3 填充（即每个位点平均拥有 2/3 个电子）。

本研究由大阪大学的 Yushi Yamada 等人完成，核心目标是通过一维（1D）的三路梯子模型（Three-leg ladder）来模拟三层镍氧化物的物理行为。研究者采用了量子多体计算领域最为精确的算法之一——密度矩阵重整化群（DMRG）方法，针对三路 $t-J$ 模型和 Hubbard 模型展开了系统性调查。研究发现，当系统处于 1/3 填充且存在强层间耦合时，会形成一个自旋间隙态（Spin-gapped state）。关键结论在于：向该态掺杂空穴会诱导自旋关联指数级衰减，而配对关联则呈现幂律衰减，这正是超导（或超导涨落）的典型信号。 相反，电子掺杂则不具备这种特性。这一不对称性为寻找更高临界温度的镍基超导体提供了重要的理论指导。

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

1.1 核心科学问题：三层体系中的配对起源

在 La$_4$Ni$_3O_{10}$ 中，Ni 的 $d_{3z^2-r^2}$ 轨道通过顶角氧原子在层间形成强耦合。由于该轨道的网络在 1/3 填充附近运行，其多体关联态是否能支撑库珀对（Cooper pairs）的形成？如果能，这种配对是来自层内（链内）还是层间（链间）？为什么空穴掺杂和电子掺杂的效果完全不同？这些问题构成了本工作的科学核心。

1.2 理论基础：从 $t-J$ 到梯子模型

为了处理强关联系统，研究者使用了 $t-J$ 模型。该模型是 Hubbard 模型在强排斥极限（$U \gg t$）下的有效模型，通过剔除双占据能级，保留了动能项（$t$）和交换项（$J$）。

三路梯子模型的哈密顿量定义如下：

$$\hat{H} = - t_{\parallel} \sum_{j,l,\sigma} (\tilde{c}_{j,l,\sigma}^{\dagger} \tilde{c}_{j+1,l,\sigma} + H.c.) - t_{\perp} \sum_{j,\langle l,l' \rangle, \sigma} (\tilde{c}_{j,l,\sigma}^{\dagger} \tilde{c}_{j,l',\sigma} + H.c.)$$

$$+ J_{\parallel} \sum_{j,l} (\mathbf{S}_{j,l} \cdot \mathbf{S}_{j+1,l} - \frac{1}{4} n_{j,l} n_{j+1,l}) + J_{\perp} \sum_{j,\langle l,l' \rangle} (\mathbf{S}_{j,l} \cdot \mathbf{S}_{j,l'} - \frac{1}{4} n_{j,l} n_{j,l'})$$

其中 $t_{\parallel}$ 和 $J_{\parallel}$ 代表链内参数，$t_{\perp}$ 和 $J_{\perp}$ 代表链间参数。在镍氧化物的背景下，层间耦合通常强于层内，因此 $t_{\perp} > t_{\parallel}$ 是本研究关注的物理区间。

1.3 技术难点：多路梯子的 DMRG 收敛性

虽然 DMRG 在处理单路链（1D chain）时几乎是完美的，但随着“梯子”路数的增加（从两路到三路），系统纠缠熵（Entanglement Entropy）会迅速增长。对于三路梯子，基态的局部结构更加复杂，特别是当 $t_{\perp}$ 很大时，单胞内会形成分子轨道状的单态。为了保证结果的可靠性，必须保持极高的截断维数（Bond Dimension）。

1.4 方法细节：DMRG 计算设置

  • 边界条件：采用开边界条件（OBC），这是 DMRG 的标准做法。
  • 晶格尺寸：主要计算在 $L_x = 80$ 的长度上进行，总位点数为 $3 \times 80 = 240$。
  • 截断维数 $m$：最高达到了 $m=5000$，确保截断误差控制在 $10^{-6}$ 以下。
  • 能量单位：以 $t_{\parallel}$ 为单位，设定 $J_{\parallel} = t_{\parallel}/3$，并根据 $J_{\perp}/J_{\parallel} = (t_{\perp}/t_{\parallel})^2$ 调整层间交换常数。

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

2.1 1/3 填充态的定性特性

在 $n=2/3$（即 1/3 填充）时，当 $t_{\perp}/t_{\parallel} = 2.5$（强耦合极限），系统进入 C0S0 态（电荷有能隙，自旋有能隙）。通过对三位点单元的精确解（Exact Diagonalization of 3-site unit cell），研究者发现基态呈现出一种“层间单态”的组合，类似于两路梯子中的共振价键（RVB）状态。如图 2 所示，层间自旋关联 $F_{12}$ 和 $F_{13}$ 呈现强负相关（约 -0.25），表明形成了稳定的局域自旋单态。

2.2 空穴掺杂下的幂律配对关联

这是本论文最重大的发现。当从 1/3 填充减少电子（空穴掺杂 $\delta_h$），自旋关联 $|F_l(r)|$ 依然保持指数级衰减（图 3a, 3b），这意味着自旋间隙依然存在。然而，配对关联函数 $P_{ll'}(r)$ 则完全不同：

  • 在对数-对数坐标系下，配对关联呈现线性下降（图 3c, 3d），这标志着幂律衰减（Power-law decay）。
  • 拟合公式为 $P(r) \sim r^{-K_{SC}}$，其中 $K_{SC}$ 在 $0.9$ 到 $1.3$ 之间波动。
  • 根据 Luttinger 液体理论，当 $K_{SC} < 2$ 时，配对关联占据主导地位，预示着超导趋势。

2.3 电子掺杂的失败

向 1/3 填充增加电子（$n > 2/3$）时，自旋关联的衰减变慢，甚至在长程表现出震荡趋势，而配对关联 $P(r)$ 极小且衰减极快（图 4c, 4d）。这说明电子掺杂倾向于破坏已有的自旋单态结构，而不是利用它来形成库珀对。

2.4 Hubbard 模型与 $t-J$ 模型的对比

在第 IV 节中，研究者对比了 Hubbard 模型在 $U/t_{\parallel}=10$ 下的结果（图 5）。结果显示，Hubbard 模型的自旋间隙 $\Delta_S$ 显著小于 $t-J$ 模型。这说明双占据的物理排斥在稳定自旋间隙及促进配对方面起到了关键作用。在 $t-J$ 模型中被严格禁止的双占据，在 Hubbard 模型中以高能激发的形式存在，削弱了层间单态的稳定性。

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

3.1 软件包：ITensor

本项研究的所有数值演示均使用开源张量网络库 ITensor（可能是 C++ 版本或 Julia 版本）完成。ITensor 是目前处理 1D 和准 1D 强关联系统最流行、最高效的工具之一。

3.2 复现步骤建议

  1. 定义晶格（Lattice）： 构建一个 ThreeLegLadder 的 SiteSet。在 ITensor 中，可以通过自定义单胞（Unit Cell）来实现。每个单胞包含 3 个 Electron 位点。

  2. 构建哈密顿量（MPO）： 使用 AutoMPO 构建方程 (1) 中的四项。注意，$t-J$ 模型需要特殊的位点定义（tJ site set），它会自动处理禁止双占据的约束。

  3. 执行计算：

    • 设置 Sweeps 策略：起始 bond dimension 从 100 开始，逐步增加到 5000。
    • 设定截断阈值（Cutoff）：通常设为 $1e-10$ 或更小。
    • 噪音调节：在早期 sweep 中加入 noise 以防止陷入局部极小值。

  4. 测量关联函数：

    • 基态收敛后，测量 $\langle S_{j,l} \cdot S_{j+r,l'} \rangle$ 和 $\langle \Delta_{j,ll'}^{\dagger} \Delta_{j+r,ll'} \rangle$。
    • 其中配对算符 $\Delta$ 需要在 MPO 框架外通过 MPS 的轨道收缩来实现，或者构建专门的测量 MPO。

3.3 关键参数参考值

  • $L_x = 80$
  • $t_{\parallel} = 1.0, t_{\perp} = 2.5$
  • $J_{\parallel} = 0.333, J_{\perp} = 2.083$
  • 电子数 $N$ 设为 $160$ (对应 $n=2/3$) 以及 $160-2, 160-4$ 等用于掺杂分析。

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

4.1 关键引用文献

  1. Sun et al., Nature 621, 493 (2023)： La$_3$Ni$_2$O$_7$ 超导电性的首次报道，是整个镍基超导研究的热潮起点。
  2. White and Scalapino, Phys. Rev. B 57, 3031 (1998)： 建立了利用 DMRG 研究 $t-J$ 梯子模型配对特性的标准范式。
  3. Sakakibara et al., Phys. Rev. B 109, 144511 (2024)： 针对三层镍氧化物 La$_4$Ni$_3O_{10}$ 的理论分析，强调了 $d_{z^2}$ 轨道的重要性。
  4. Dagotto et al., Phys. Rev. B 45, 5744 (1992)： 二路梯子模型中的超导性研究基础。

4.2 工作局限性评论

尽管该研究在数值上非常严谨，但仍存在以下局限性：

  • 维度限制：三路梯子本质上还是 1D 系统。根据 Mermin-Wagner 定理，1D 系统在有限温度下不存在长程超导序，只有 0 温下的准长程关联。将其外推至 3D 块体材料需要谨慎，必须考虑层间相干耦合。
  • 单轨道近似：论文主要关注了类 $d_{z^2}$ 的有效单轨道模型。实际上，镍氧化物是典型的多轨道系统（$d_{x^2-y2}$ 和 $d_{z^2}$ 轨道同时存在），轨道间的杂化和电荷转移（Charge Transfer）可能对 $T_c$ 有显著影响。
  • 弱耦合对比：论文提到其结果与弱耦合理论（RG 方法）预言的 C1S0 态一致性有限。在弱耦合下，电子掺杂和空穴掺杂是对称的，而强关联下的 DMRG 结果展现了剧烈的不对称性。这种“强弱耦合之争”暗示了镍氧化物可能处于极强的关联区。

5. 其他补充：物理直觉与未来展望

5.1 物理直觉：为什么是空穴掺杂？

一个直观的物理图景是：在 1/3 填充时，每个三位点单元被两个电子占据形成一个强力的自旋单态（Singlet）。当你抽出一个电子（空穴掺杂）时，原本锁死的单态获得了流动的空间，这些空穴以“空穴对”的形式移动，从而保留了单态的配对特性并赋予其相干性。而当你增加电子时，由于 Pauli 不相容原理和强库仑排斥，新增电子会强行拆散原有的单态，导致系统进入一种类似杂乱自旋的态，不利于配对。

5.2 对实验的启示

目前的实验合成 La$_4$Ni$_3O_{10}$ 主要还是通过高压手段。本研究暗示，如果实验上能够通过化学替代（例如用 Sr 替代 La）实现精确的空穴掺杂控制，或许能在更低的压力甚至常压下稳定超导相。

5.3 镍氧化物 vs. 铜氧化物

在铜氧化物中，超导主要发生在半填充（Half-filling）附近。而镍氧化物（特别是三层体系）展示了在 1/3 填充附近的奇特物理。这表明过渡金属氧化物的超导机制具有极大的多样性。三路梯子模型不仅是一个计算基准，更是连接理想模型与复杂真实材料的“理论桥梁”。

5.4 未来研究方向

  • 含时 DMRG：研究掺杂后的动力学特性，如光激发下的配对增强。
  • 多轨道扩展：在三路梯子基础上引入 $d_{x^2-y^2}$ 轨道，构建双轨道三路梯子（Two-orbital three-leg ladder）。
  • 纠缠谱分析：通过计算纠缠谱（Entanglement Spectrum）来进一步确认系统的拓扑性质和边缘态，判定其是否属于特殊的 Luttinger 液体相。

作者注：这项工作再次证明了 DMRG 在解决具体材料物理问题时的强大生命力。在量子计算真正成熟之前，这种基于张量网络的经典模拟依然是人类理解强关联量子世界的“最强显微镜”。