来源论文: https://arxiv.org/abs/2605.10101v1 生成时间: May 17, 2026 10:00

0. 执行摘要

自 2023 年在双层镍氧化物 $La_3Ni_2O_7$ 中发现近 80 K 的高压超导性以来，该材料已成为非常规超导研究的核心。然而，密度泛函理论（DFT）预测的“三费米面”图像与角分辨光电子能谱（ARPES）实验观测到的“$\gamma$ 坑（pocket）缺失”现象之间存在显著矛盾。本项研究通过结合基于时间演化变分原理（TDVP）的团簇微扰理论（CPT）以及大规模密度矩阵重整化群（DMRG）计算，系统地研究了现实的双层两轨道 Hubbard 模型。研究发现，电子关联驱动了显著的轨道选择性重整化：随着关联强度增加，$d_{z^2}$ 衍生的 $\gamma$ 带下沉至费米能级以下，而 $d_{x^2-y^2}$ 轨道在费米能级附近保留了费米弧结构。更为重要的是，计算表明超导配对通道经历了从 $d_{z^2}$ 主导到 $d_{x^2-y^2}$ 主导的跨越式转变，但始终保持着层间单态 $s_{\pm}$ 配对对称性。这一发现协调了实验与理论的冲突，指明了关联效应在镍基超导机理中的决定性作用。

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

1.1 核心科学问题：费米面的“消失”之谜

$La_3Ni_2O_7$ 的电子结构由 Ni 的两个 $e_g$ 轨道（$d_{x^2-y^2}$ 和 $d_{z^2}$）主导。在 DFT 框架下，这两个轨道通过成键与反成键组合形成 $\alpha, \beta, \gamma$ 三个费米面。其中 $\gamma$ 坑被认为是驱动 $s_{\pm}$ 配对的关键。然而，多项 ARPES 实验显示 $\gamma$ 带实际上位于费米能级下方约 70 meV 处，且能带发生了剧烈的重整化。核心问题在于：这种由关联引起的费米面重构是否会破坏超导？如果超导依然存在，其配对通道和对称性发生了怎样的演化？

1.2 理论基础：双层两轨道 Hubbard 模型

为了描述该体系，研究采用了包含 $d_{x^2-y^2}$ 和 $d_{z^2}$ 轨道的双层 Hubbard 模型。哈密顿量由非相互作用项 $H_0$ 和相互作用项 $H_U$ 组成：

$$H_0 = \sum_{ijlm\alpha\beta\sigma} t_{ij,lm}^{\alpha\beta} c_{il\alpha\sigma}^{\dagger} c_{jm\beta\sigma} + \sum_{il\alpha\sigma} (\epsilon_{\alpha} - \mu) n_{il\alpha\sigma}$$

其中，跃迁参数 $t$ 采用了真实材料的紧束缚参数。相互作用项 $H_U$ 包含了经典的库仑排斥 $U$、层间/轨道间排斥 $U'$ 以及 Hund 耦合 $J_H$：

$$V_i = \sum_{\alpha} U n_{i\alpha\uparrow} n_{i\alpha\downarrow} + \sum_{\alpha < \beta, \sigma\sigma'} (U' - \frac{J_H}{2}) n_{i\alpha\sigma} n_{i\beta\sigma'} - \sum_{\alpha < \beta} 2J_H S_{i\alpha} \cdot S_{i\beta} + \sum_{\alpha \neq \beta} J_H c_{i\alpha\uparrow}^{\dagger} c_{i\alpha\downarrow}^{\dagger} c_{i\beta\downarrow} c_{i\beta\uparrow}$$

1.3 技术难点：强关联下的谱函数计算

处理这种具有多个轨道和层间耦合的复杂体系，传统方法如平均场论无法捕捉轨道选择性的 Mott 物理。技术难点在于如何在较大的团簇上精确模拟动力学关联函数。传统的精确对角化（ED）受限于尺寸，而量子蒙特卡罗（QMC）在处理 Hund 耦合项时面临严重的正负号问题。此外，计算超导配对相关函数需要极高的数值精度以区分不同的对称性通道。

1.4 方法细节：TDVP-CPT 与 DMRG 的结合

  1. TDVP-CPT：研究者利用 TDVP 演化矩阵乘积态（MPS），在 16 个物理位点（2x2x2 团簇）上计算格林函数 $G_{ij}(\omega)$。相比于传统的 Krylov 子空间方法，TDVP 能提供更长的相干演化时间，从而获得更高的频率分辨率。通过 CPT 公式将团簇格林函数扩展到全布里渊区： $$G^{-1}(\mathbf{k}, \omega) = G_{cluster}^{-1}(\omega) - V(\mathbf{k})$$
  2. DMRG：为了确定基态配对相关函数 $\Phi(r)$，在 $2 \times L_y \times L_x$（$L_x=24$）的条带几何上进行大规模 DMRG 计算。研究者保留了高达 15000 的键维（bond dimension $D$），以确保配对相关函数的代数衰减行为（即准长程序）能够被清晰观测到。

2. 关键 Benchmark 体系与数据解析

2.1 轨道选择性能带重整化（图 2）

计算展示了随着 $U$ 从 1.5 eV 增加到 6.0 eV 的演化：

  • 弱耦合 (U=1.5, 2.2 eV)：$\gamma$ 坑依然存在，此时 $d_{z^2}$ 轨道的相干性较强，态密度（DOS）在费米面附近有显著贡献。
  • 中强耦合 (U=3.7 eV)：这是一个关键临界点。$\gamma$ 带被推至费米面以下，此时 $\gamma$ 坑消失，系统发生 Lifshitz 相变。这与 ARPES 观测到的薄膜数据高度吻合。
  • 强耦合 (U=6.0 eV)：$\alpha$ 带和 $\beta$ 带在 $\Gamma-X$ 方向开启拟能隙，导致费米面重构成由 $d_{x^2-y^2}$ 主导的“费米弧（Fermi arcs）”。

2.2 超导配对相关函数 $K_{sc}$（图 3）

研究定义了配对相关函数 $\Phi(r) \sim r^{-K_{sc}}$。根据 Luttinger 液体理论，$K_{sc} < 2$ 意味着配对易感性发散。数据解析如下：

  • 在 $U=1.5$ eV 时，$K_{sc}^{(dz2)} = 1.41$，远小于 $K_{sc}^{(dx2-y2)} = 1.94$。这表明在 $\gamma$ 坑存在时，$d_{z^2}$ 轨道主导超导。
  • 在 $U=6.0$ eV 时，情况发生逆转：$K_{sc}^{(dx2-y2)} = 1.38$，而 $K_{sc}^{(dz2)} = 1.99$。这定量地证明了：即使 $\gamma$ 坑消失，超导依然存在，但配对动力学切换到了 $d_{x^2-y^2}$ 轨道。

2.3 相互作用参数的影响（图 4 & 表 S1/S2）

通过控制变量法，研究发现：

  1. 轨道间杂化 $t^{xz}$：在弱耦合下是驱动超导的关键，移除 $t^{xz}$ 会导致配对相关函数指数衰减。
  2. Hund 耦合 $J_H$：在强耦合下增强了层间反铁磁交换，从而稳定了 $d_{x^2-y^2}$ 的配对。若关闭 $J_H$，系统倾向于形成电荷密度波（CDW）而非超导。

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

3.1 软件包与开发环境

该工作主要基于 Julia 语言的高性能张量网络库构建：

  • MPSKit.jl：用于执行 TDVP 时间演化。它支持非阿贝尔对称性的张量算子，是处理多轨道 Hubbard 模型的利器。
  • TensorKit.jl：作为底层张量库，处理 $U(1)_{charge} \times SU(2)_{spin}$ 对称性，极大减少了希尔伯特空间的有效维度。

3.2 核心算法步骤：TDVP-CPT

  1. 基态寻找：使用 DMRG 寻找 $2\times 2\times 2$ 团簇的基态 $|0\rangle$，截断误差 $\epsilon \lesssim 1 \times 10^{-6}$。
  2. 算子作用：作用产生算子 $c^{\dagger}_{il\alpha\sigma}|0\rangle$ 得到初始态。
  3. 时间演化：利用 TDVP 将状态演化至 $t_{max} = 30$。采用二阶积分方案，步长 $\Delta t = 0.05$。演化过程中动态调整 MPS 键维，保持截断误差在 $10^{-5}$ 量级。
  4. 傅里叶变换：引入高斯展宽因子 $e^{-\eta^2 t^2}$ ($\eta=0.08$) 进行积分，获得单粒子谱函数。

3.3 复现难点：大规模 DMRG

复现 $L_x=24$ 的配对相关函数需要极高的计算资源。建议：

  • 并行化策略：利用多线程或分布式内存并行处理张量收缩。
  • 收敛判定：配对相关函数对键维 $D$ 极其敏感。在 $L_y=2$ 的体系中，$D$ 需达到 15000 以上才能观察到 $K_{sc}$ 的稳定值。
  • 开源链接：读者可参考 MPSKit.jl 官方文档 进行基础 MPS 算子的搭建。

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

4.1 关键引用

  • [1] Sun et al., Nature (2023)：首次报道 80 K 超导。本文的基础实验背景。
  • [4] Luo et al., PRL (2023)：提出了双层两轨道紧束缚模型。本文的模型参数来源。
  • [33, 35] Wang et al. / Sun et al. (2025)：ARPES 实验发现 $\gamma$ 坑下沉。本文直接回应的实验争议。
  • [51] Sénéchal et al., PRB (2002)：CPT 理论的奠基性工作。

4.2 局限性评论

尽管本工作通过关联效应优雅地解决了能带冲突，但仍存在以下局限：

  1. 尺寸效应：虽然 16 位点团簇在 CPT 中已算庞大，但对于动量空间的分辨率仍有限，特别是对于费米弧细节的刻画可能存在模糊。
  2. 氧空位与非化学计量比：实验中的 $La_3Ni_2O_{7-\delta}$ 往往存在氧空位，这会改变掺杂水平。本文固定在 $n=3/8$ 的平均填充，忽略了局域无序的影响。
  3. 层间耦合的动态性：模型假设 $t_{\perp}$ 是常数，但在高压下，Ni-O-Ni 键角的变化可能引入更复杂的层间依赖。

5. 补充解析：镍基与铜基超导的异同

5.1 轨道角色的互换

在铜氧化物中，单一的 $d_{x^2-y^2}$ 轨道包揽了费米面和配对。但在 $La_3Ni_2O_7$ 中，存在“分工与协作”：

  • $d_{z^2}$ 轨道：提供强烈的层间反铁磁背景，像“胶水”一样提供配对势。
  • $d_{x^2-y^2}$ 轨道：提供巡游性，负责在费米面上形成相干凝聚。 这种“轨道选择性”是镍基超导区别于铜基的核心特征。

5.2 关联驱动的 $s_{\pm}$ 对称性稳定性

本文最重要的结论之一是：尽管主导轨道变了，但配对对称性保持为 $s_{\pm}$。这是因为双层结构自带的镜像对称性决定了能带的“成键/反成键”性质。只要配对函数在成键带（$\alpha, \gamma$）和反成键带（$\beta$）上异号，系统就能有效利用层间排斥能降低总能量。这种鲁棒性解释了为什么即使材料生长条件（薄膜 vs 块材）导致费米面拓扑不同，超导迹象依然广泛存在。

5.3 对未来实验的启示

研究建议，未来的 ARPES 实验应重点关注 $d_{x^2-y^2}$ 轨道衍生的费米弧随压力的演化。如果能观测到本文预测的拟能隙行为，将有力地支持关联驱动的超导机制。同时，中子散射实验应寻找与 $J_H$ 和层间交换相关的特征共振峰，以验证配对“胶水”的来源。