来源论文: https://arxiv.org/abs/2603.13604v1 生成时间: Mar 21, 2026 17:42

双层镍氧化物超导关联的层级结构：从 $d_{z^2}$ 轨道配对到 $d_{x^2-y^2}$ 关联重布的深度解析

0. 执行摘要

近期，双层镍氧化物 $La_3Ni_2O_7$ 在高压下发现约 80 K 的超导转变温度（$T_c$），这使其成为继铜氧化物之后又一极具潜力的高温超导体系。然而，其超导配对的微观起源——究竟是类铜氧化物的 $d_{x^2-y^2}$ 轨道主导，还是双层结构特有的 $d_{z^2}$ 轨道主导——一直是理论界争论的焦点。本文解析的最新研究通过变分蒙特卡洛（VMC）方法，对双层二轨道 Hubbard 模型进行了非微扰数值模拟。研究发现，$La_3Ni_2O_7$ 展现出一种独特的“层级配对结构”：初级配对相互作用起源于 $d_{z^2}$ 轨道的成键-反键分裂（bonding-antibonding splitting），而轨道杂化（orbital hybridization）则将超导关联重新分布到 $d_{x^2-y^2}$ 通道。这种机制解释了为何两个轨道通道表现出量级相当的长程超导关联，并证明了 $s_{\pm}$ 配对对称性对费米面拓扑变化的鲁棒性。这一发现为理解多层强关联超导体提供了一个全新的统一框架。

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

1.1 核心科学问题：配对的“第一推动力”是什么？

在 $La_3Ni_2O_7$ 中，Ni 离子的价态接近 $d^{7.5}$，这意味着 $d_{z^2}$ 和 $d_{x^2-y^2}$ 轨道都参与了低能物理过程。目前的理论争议主要集中在：

  1. 类铜氧化物方案：超导电性主要由 $d_{x^2-y^2}$ 轨道驱动，类似于单层铜氧化物的物理。
  2. 双层成键方案：配对由 $d_{z^2}$ 轨道在层间的强跳迁（$t_{\perp}$）引起的成键-反键分裂驱动。
  3. Hund 耦合方案：认为 $d_{z^2}$ 电子形成局域自旋单态，通过强 Hund 耦合诱导 $d_{x^2-y^2}$ 电子配对。

本研究试图回答：在考虑真实轨道杂化的情况下，超导关联如何从微观相互作用演化为宏观序参数？

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

研究采用的模型包含了 Ni 的 $d_{z^2}$ 和 $d_{x^2-y^2}$ 轨道，其哈密顿量由动能项（$H_{kin}$）和相互作用项（$H_{int}$）组成：

  • 动能项：包含了层内跳迁、层间跳迁（特别是 $d_{z^2}$ 轨道的强层间耦合 $t_{\perp}$）以及轨道间的杂化项。
  • 相互作用项：包括层内同轨道库仑排斥 $U$、异轨道排斥 $U'$、Hund 耦合 $J$ 以及配对跃迁项 $J'$。为了满足旋转不变性，设定 $U' = U - 2J$，$J = J'$。

关键参数是 $\Delta E/t_{\perp}$，其中 $\Delta E$ 是 $d_{x^2-y^2}$ 与 $d_{z^2}$ 的能级差。这决定了体系的“轨道层级”。

1.3 技术难点：非微扰处理与费米面敏感性

传统的 RPA（随机相位近似）或 FLEX（涨落交换近似）方法高度依赖于费米面上的嵌套（nesting）向量。但在 $La_3Ni_2O_7$ 中，费米面拓扑随压力和掺杂变化剧烈（例如 $\gamma$ 片层是否消失）。

  • 难点 1：如何处理强关联效应，避免弱耦合理论在接近莫特绝缘体边缘时的失效？
  • 难点 2：如何区分“变分参数中的配对强度”与“物理上的长程关联函数”？
  • 难点 3：如何精确描述轨道杂化对超导序参数在动量空间分布的影响？

1.4 方法细节：变分蒙特卡洛（VMC）

研究采用了 Gutzwiller-Jastrow 型试探波函数：

$$ |\Psi\rangle = P_G^{(2)} P_{Jc} P_{Js} |\Phi\rangle $$
  1. $|\Phi\rangle$：有效 BdG 哈密顿量的基态，允许 $s_{\pm}$ 或 $d$ 波配对。其中变分参数 $\tilde{\Delta}_{zz}$ 和 $\tilde{\Delta}_{xx}$ 分别代表 $d_{z^2}$ 和 $d_{x^2-y^2}$ 的层间配对场。
  2. $P_G^{(2)}$：二轨道 Gutzwiller 因子，调节 16 种局部电子配置的权重，捕获强关联带来的电荷涨落抑制。
  3. $P_{Jc}, P_{Js}$：电荷和自旋 Jastrow 因子，捕获长程关联。
  4. 优化算法：使用随机重组（Stochastic Reconfiguration, SR）方法最小化能量，获得最优变分参数。

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

2.1 模拟体系配置

  • 格点规模：主要在 $2 \times 2 \times 20 \times 20$ 的晶格上运行（共 800 个位点，1600 个轨道）。
  • 电子密度：填充度设为 $n = 1.5$（每个 Ni 原子平均 1.5 个电子），对应实验中的高压超导相。
  • 物理参数：设定 $U/t_{\perp} = 8.0, J/U = 0.15$。$t_{\perp}$ 设定为能量单位（约 0.659 eV）。

2.2 核心数据解析：电荷重组与轨道占据

研究发现，随着 $\Delta E/t_{\perp}$ 的增加，体系发生了明显的电荷重分配（Fig. 2a）：

  • 在非相互作用情况下，$n_{z^2}$ 随能级偏移线性变化。
  • 在相互作用下，$n_{z^2}$ 被“锁定”在接近半填充的状态。这种电荷钉扎效应是强关联体系的典型特征，为超导配对提供了稳定的背景。

2.3 变分参数与物理关联的背离（核心发现）

这是该研究最深刻的数据点：

  • 变分参数（Fig. 3）：优化后的 $\tilde{\Delta}_{zz}$ 在 $\Delta E/t_{\perp} > 0.65$ 时显著增大，而 $\tilde{\Delta}_{xx}$ 始终趋于 0。这说明从变分能的角度来看，配对动力学几乎完全由 $d_{z^2}$ 驱动。
  • 关联函数（Fig. 4）：然而，计算长程超导关联函数 $P_{zz}$ 和 $P_{xx}$ 时发现，尽管 $d_{x^2-y^2}$ 通道没有变分配对场，其物理关联 $P_{xx}$ 却与 $P_{zz}$ 量级相当，甚至在某些区域更大。
  • 结论：这证明了超导关联并非直接在 $d_{x^2-y^2}$ 轨道产生，而是通过轨道杂化从 $d_{z^2}$ “流向”了 $d_{x^2-y^2}$。

2.4 动量空间能隙结构

  • $s_{\pm}$ 对称性：在 $\gamma$ 和 $\delta$ 带上表现出几乎各向同性的能隙，且符号相反。
  • 轨道杂化节点：在 $\alpha$ 和 $\beta$ 带上，能隙表现出强烈的动量依赖性，在杂化禁戒的 $k_x = \pm k_y$ 线上消失（Fig. 6c, d）。这一特征完美解释了为何该体系表现出 $s_{\pm}$ 却又具有某种“类 d 波”的动量分布形态。

3. 代码实现细节，复现指南与开源工具

3.1 软件包选择

该研究主要基于**变分蒙特卡洛（VMC）**框架。虽然作者未直接给出 repo 链接，但根据论文描述及领域惯例，复现此类计算通常使用以下开源或半开源软件包：

  1. mVMC (Many-variable Variational Monte Carlo): 由东京大学开发的开源 VMC 软件包，能够处理多轨道 Hubbard 模型及复杂的 Jastrow 因子。 GitHub Link
  2. Variational Monte Carlo (VMC) Codes by Sorella: 参考文献 [70] 提到的随机重组算法（SR）通常与此算法框架相关。

3.2 复现指南

  1. Step 1: 构建单粒子哈密顿量：
    • 提取 DFT 计算所得的 $d_{z^2}$ 和 $d_{x^2-y^2}$ 的跳迁参数。
    • 包含层间耦合 $t_\perp$ 和轨道杂化项 $V_{xz}(k)$。
    • 引入双重计数校正（Double-counting correction），这是处理多轨道模型时确保能级结构正确的关键步。
  2. Step 2: 构造试探波函数：
    • 定义 BdG 准粒子基态，设置初试 $\Delta_{zz} (s_{\pm})$ 参数。
    • 施加 Gutzwiller 算符以抑制双占据。
  3. Step 3: 优化过程：
    • 使用 SR 算法优化。建议先在小格点（如 $L=8$）上稳定参数，再扩展至 $L=20$。
    • 注意 $n=1.5$ 填充度的化学势 $\mu$ 调整。
  4. Step 4: 计算关联函数：
    • 在优化后的波函数上，采样计算 $P_{mm}(\mathbf{r})$。必须确保 $\mathbf{r}$ 取到格点的最大距离（$L/2$）以验证超导的长程有序性。

3.3 计算性能要求

  • 由于涉及多变量优化（SR 矩阵反转），计算量随单胞内轨道数 $N_{orb}^3$ 扩展。
  • $L=20$ 的体系建议使用至少 512 核以上的并行计算集群，单次优化迭代可能需要数小时。

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

4.1 关键引用

  1. [1] Sun, H. et al. (Nature 2023)：镍氧化物超导的发现实验，奠定了 80 K 转变温度的基础。
  2. [18] Nakata, M., Kuroki, K. et al. (Phys. Rev. B 2017)：本文作者之一在实验发现前就预测了 $La_3Ni_2O_7$ 可能是 $s_{\pm}$ 超导体。
  3. [53] Kuroki, K. et al. (Phys. Rev. B 2002)：提出了双层 Hubbard 模型中通过成键-反键分裂产生超导的普适机制。
  4. [38] Ryee, S. et al. (Phys. Rev. Lett. 2024)：使用 CDMFT 方法研究了费米面拓扑，本文与其结论进行了重要对比。

4.2 局限性评论

尽管本工作在解释配对起源上非常出色，但仍存在以下局限：

  • 模型简化：仅考虑了 Ni 3d 轨道。虽然有效的两轨道模型捕获了主要物理，但氧 2p 轨道在强关联下的电荷转移效应（Charge Transfer character）被隐式地包含在了参数 $U$ 中。这在研究镍氧化物这类电荷转移能较小的体系时可能存在系统误差。
  • 三维效应：虽然研究自称是“双层”，但实际材料是准三维的，层间的 c 轴方向相干性可能进一步提升 $T_c$。VMC 目前主要在有限尺寸的 2D 晶格上运行，对层间长程有序的描述仍有限。
  • 变分偏置（Variational Bias）：虽然 VMC 比 RPA 更可靠，但它始终受限于试探波函数的形式。如果存在更复杂的非共线自旋态或交织序（intertwined orders），简单的 Gutzwiller-Jastrow 波函数可能无法完全捕捉。

5. 补充讨论：镍氧化物 vs 铜氧化物 vs 铁基超导体

为了更全面地理解这项工作，我们需要将其置于高温超导的全景图中：

5.1 “混合血统”的超导体

$La_3Ni_2O_7$ 表现出一种奇特的混合特征：

  • 像铜氧化物的地方：它具有 Ni-O 平面，主要的超导关联在 $d_{x^2-y^2}$ 轨道上看起来非常“类铜”。
  • 像铁基的地方：它是多轨道的，具有 $s_{\pm}$ 对称性，且超导性对费米面空穴/电子口袋的细节相对不敏感。
  • 独特之处：它利用了垂直方向的“二聚体”（dimer）能量分裂。这种“垂直配对”是铜氧化物（单层或弱耦合多层）所不具备的。

5.2 轨道杂化的关键作用

本研究最精彩的逻辑在于：如果没有轨道杂化，$d_{z^2}$ 将是孤立的强配对通道，而 $d_{x^2-y^2}$ 将是死寂的（不配对）。杂化项 $V_{xz}(k)$ 充当了“导管”，将 $d_{z^2}$ 产生的“配对流”导入到了在层内运动更快的 $d_{x^2-y^2}$ 电子中。这解释了为什么实验上观察到的 $T_c$ 如此之高——它结合了 $d_{z^2}$ 的强相互作用能标和 $d_{x^2-y^2}$ 的高迁移率。

5.3 压力调控物理的本质

根据论文的 $\Delta E/t_{\perp}$ 分析，高压的作用可能不仅在于减小晶格常数，更在于通过改变八面体旋转或压缩，精准地调节这两个轨道的相对能级偏移，使体系进入“轨道钉扎”的理想超导区间。这一结论为通过化学掺杂（如用 Pr 取代 La）或应力工程来调节超导性提供了理论指南。

5.4 对未来研究的启示

这项工作提示我们，寻找下一个高温超导体的方向不应局限于“寻找类铜的单带体系”，而应寻找具有“强配对源轨道”和“高效传输轨道”的杂化体系。这种“配对与输运分工”的层级结构，可能是打破 $T_c$ 瓶颈的关键路径。