来源论文: https://arxiv.org/abs/2201.05495 生成时间: Mar 01, 2026 02:07

0. 执行摘要

自 2019 年在无限层镍氧化物 NdNiO2 中发现超导电性以来，科学界对其超导机制是否与铜氧化物（Cuprates）相同存在激烈争论。作为同构且同电子配置（$d^9$）的体系，NdNiO2 被视为研究高温超导的关键参照。然而，由于镍的电负性与铜不同，其电荷转移能级和关联效应表现出独特特征。

本研究采用最先进的多参照量子化学方法（Multireference Quantum Chemistry Methods），通过对电子移除准粒子态（即空穴掺杂态）的 $ab initio$ 波函数分析，得出了几个关键结论：

  1. 空穴定位差异：在 NdNiO2 中，掺杂空穴主要定位于 Ni $3d_{x^2-y^2}$ 轨道，形成闭壳层单态（Closed-shell singlet），而非铜氧化物中典型的 Zhang-Rice 单态。
  2. 三带模型失效：传统的、在铜氧化物中极其成功的“三带 Hubbard 模型”仅能描述 NdNiO2 波函数的 60% 左右，而对铜氧化物的描述力高达 90%。
  3. 径向关联效应：镍氧化物中存在极强的动态径向关联（Dynamic radial-type correlations），这使得 $d^8$ 多重态效应更加显著，空穴的分布呈现出更强的三维特征。

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

核心科学问题

长期以来，强关联超导体系的描述依赖于简化的有效模型。对于铜氧化物，Zhang-Rice 单态（ZR singlet）是理解其低能物理的核心。科学界的核心疑问是：在镍氧化物中，掺杂空穴究竟是像铜氧化物那样主要分布在氧的 $2p$ 轨道并与金属离子耦合，还是更多地停留在金属 Ni 位点上？这种微观层面的波函数构成直接决定了超导配对的对称性与相互作用强度。

理论基础：从 ZSA 框架到多参照波函数

根据 Zaanen-Sawatzky-Allen (ZSA) 分类，铜氧化物通常被归类为电荷转移型绝缘体（Charge-transfer insulators），其电荷转移能 $\Delta$ 小于库仑排斥能 $U$。而早期的 DFT 计算表明，NdNiO2 更接近 Mott-Hubbard 极限。这意味着其基态波函数中，$d^8 L$（L 代表氧孔穴）的权重远低于铜氧化物。为了在不引入经验参数的前提下定量分析这一点，必须使用能够平衡描述静态和动态电子关联的波函数方法。

技术难点：强关联与大活性空间的挑战

  1. 静态关联（Static Correlation）：在强关联体系中，单行列式（如 HF 或 DFT）无法描述电子在近简并轨道间的重分布。需要使用完整活性空间自洽场（CASSCF）方法。
  2. 动态关联（Dynamic Correlation）：电子间的瞬时库仑避让对能量精度至关重要，这通常需要二阶微扰（CASPT2/NEVPT2）或组态相互作用（CI）来补偿。
  3. 计算规模：为了包含必要的 Ni $3d, 4d$ 以及 O $2p$ 轨道，活性空间（Active Space）会迅速膨胀。传统的 CASSCF 算法由于受限于行列式数量，难以处理超过 18 个电子在 18 个轨道（CAS(18,18)）以上的体系。

方法细节：嵌入式簇模型与 FCIQMC

本研究采用了嵌入式簇模型（Embedded Cluster Approach）。研究人员选取了包含五个 $NiO_4$（或 $CuO_4$）单元的量子簇，并将其置于拟合的 Madelung 势场（由点电荷阵列构成）中，以精确模拟晶格环境的静电效应。

为了克服活性空间的限制，研究引入了以下高级算法：

  • FCIQMC（全构型相互作用量子蒙特卡洛）：通过随机采样 Slater 行列式空间来求解接近 FCI 极限的基态波函数，支持高达 CAS(24,30) 的超大空间。
  • DMRG（密度矩阵重整化群）：利用矩阵乘积态（MPS）表示波函数，有效处理一维链状或准二维体系的纠缠。在本工作中，DMRG 被用于求解大型活性空间下的 eigenvalue 问题。
  • NEVPT2（n-电子价态微扰理论）：相比 CASPT2，NEVPT2 解决了侵入态（intruder states）问题，能够更稳健地提供动态关联修正。

2. 关键 benchmark 体系，计算数据，性能数据

Benchmark 体系：NdNiO2 vs CaCuO2

研究选择了无限层镍氧化物 NdNiO2 及其同构的铜氧化物参照物 CaCuO2。这种对比实验可以抵消掉晶体结构（正方平面配位）带来的共性，直接凸显金属离子（Ni+ vs Cu2+）性质的差异。

关键计算数据分析

(1) 多重态能级分布（Table II）

通过 CASSCF+CASPT2 计算得到的电子移除态能级：

  • 基态对称性：两者均为 $^1A_{1g}$。但在 NdNiO2 中，$^1A_{1g}$ 到 $^3B_{1g}$（三线态）的激发能约为 1.88 eV，而 CaCuO2 中为 2.50 eV。这说明镍氧化物中单态-三线态的劈裂较小。
  • 稳定化能：引入动态关联（CASPT2）后，NdNiO2 的 $^1A_{1g}$ 态稳定了 0.53 eV，而 CaCuO2 仅稳定了 0.24 eV。这表明镍氧化物中的动态关联远比铜氧化物活跃。

(2) 波函数权重分析（Table III）

这是本论文最核心的数据。在 $d^8$（空穴掺杂）基态波函数中：

  • NdNiO2 (CAS-1)：其主要的构型是 $|d^0_{b1} p^{\uparrow\downarrow}_{b1}\rangle$，权重占 51.87%。这意味着空穴主要留在 Ni 位点，且氧轨道保持闭壳层。ZR 单态构型权重仅约 14%。
  • CaCuO2 (CAS-1)：ZR 单态构型（$|d^\uparrow p^\downarrow\rangle$）权重占 67.58%。这证实了铜氧化物空穴主要分布在氧上的物理图景。
  • 动态关联的影响：在扩大活性空间（CAS-2，引入 $3d$ 轨道径向关联）后，NdNiO2 的主要构型权重下降至 42.40%，进一步验证了三带模型无法完全捕捉其物理特性。

(3) 电子密度差异分析（Figure 3）

电子密度差图（$\rho_{d^8} - \rho_{d^9}$）直观显示：

  • 在 $xy$ 平面内，NdNiO2 的空穴密度紧紧包裹在 Ni 核周围，且半径较大；CaCuO2 的空穴则均匀分布在四个 O 离子和中央 Cu 离子上。
  • 在 $xz$ 平面内，NdNiO2 的空穴在 $z$ 方向上延伸达 2 Å，表现出显著的三维性；而 CaCuO2 的空穴被限制在 1 Å 以内。这与 Ni $3d_{z^2}$ 轨道的参与有关。

3. 代码实现细节，复现指南，软件包及开源链接

软件包清单

  1. OpenMolcas (v21.11)：主要量子化学计算平台。负责积分计算、CASSCF 轨道优化以及 CASPT2/NEVPT2 微扰修正。
  2. QCMaquis：高性能 DMRG 模块，作为 OpenMolcas 的插件使用，负责大活性空间的配置相互作用求解。
  3. Multiwfn：由卢天开发的多功能波函数分析软件，用于提取自然轨道密度和生成密度差图。
  4. Jmol：用于轨道展示和晶体结构可视化。

复现指南

若要复现本文结果，需遵循以下步骤：

第一步：建立嵌入式模型

  1. 根据实验晶格参数生成五个 $NiO_4$ 单元。相邻的 Ni 离子需使用全离子势（TIPs）或大核有效核心势（ECP）代替，以保持配位环境。
  2. 使用 MOLCAS 内置的 GATEWAY 模块输入基组（推荐 ANO-L 质量）。
  3. 拟合点电荷阵列：计算 Madelung 常数，通过调整外层点电荷电荷量，使量子簇内的电势梯度与无限晶体一致。

第二步：轨道优化与活性空间选择

  1. 执行单行列式 RHF 计算作为起始。
  2. 定义 CAS 空间：
    • CAS-1: 包含 Ni/Cu $3d_{x^2-y^2}$, $3d_{z^2}$，对称匹配的 O $2p$ ZR-like 轨道，以及对应的双壳层（double-shell）$4d$ 轨道。共 18 电子 / 24 轨道。
    • CAS-2: 在 CAS-1 基础上增加 Ni/Cu $3d_{xy}, 3d_{xz,yz}$ 等其余 $d$ 轨道。共 24 电子 / 30 轨道。
  3. 使用 DMRG-CASSCF 进行轨道自洽优化。

第三步：关联能修正与熵分析

  1. 在优化后的轨道基础上，运行 NEVPT2 计算以获得动态关联能。
  2. 调用 QCMaquis 计算轨道纠缠熵（Orbital Entanglement Entropy）。通过分析单轨道熵 $s(1)_i$ 和互信息 $I_{i,j}$ 来量化 Ni $3d$ 与 $4d$ 之间的径向关联强度。

4. 关键引用文献与评论

关键参考文献

  1. Li et al., Nature 572, 624 (2019)：镍氧化物超导的实验发现论文，奠定了领域基础。
  2. Zhang & Rice, Phys. Rev. B 37, 3759 (1988)：定义了铜氧化物中的 ZR 单态，是本文对比的理论基石。
  3. Katukuri et al., Phys. Rev. B 102, 241112 (2020)：作者的前期工作，详细讨论了 NdNiO2 的亲本态（Parent state）磁性相互作用。
  4. Roos et al., J. Phys. Chem. A 109, 6575 (2005)：ANO-L 基组的开发，保证了重元素计算的精度。

局限性评论

尽管本工作在波函数层面达到了极高的精度，但仍存在以下局限：

  • Nd 4f 电子的处理：本研究为了计算效率，将 Nd 的 $4f$ 电子冻结在核心（Large-core ECP）。然而，近期有研究指出 Nd $5d$ 与 Ni $3d$ 的杂化（即所谓的“电子库”效应）对费米面拓扑有重要影响。忽略 $4f$ 的极化效应可能会轻微改变 Ni 位点的局部屏蔽环境。
  • 簇尺寸限制：虽然 5 单元簇能捕捉短程关联，但超导是一种宏观量子现象，长程自旋涨落和动量空间中的拓扑特性无法在小簇模型中完全体现。
  • 动态筛选（Dynamic Screening）：计算中虽然考虑了点电荷环境，但环境对空穴的极化响应（Polarizable Embedding）并未动态处理，这可能导致激发起能量略微偏高。

5. 补充：量子信息视角下的电子关联

本论文的一个亮点是引入了轨道纠缠熵的概念。在传统的固态物理视角下，我们习惯讨论轨道杂化强度；而在量子化学视角下，纠缠熵提供了更量化的指标。

径向关联（Radial Correlation）的物理图景

论文中提到的“径向关联”或“双壳层效应（Double-shell effect）”，本质上是描述电子在原子核周围“呼吸”的一种动态过程。当 Ni 离子从 $d^9$ 变为 $d^8$ 时，剩余电子感受到的有效核电荷增加，电子云会发生收缩。量子化学通过引入占位很低的 $4d$ 辅助轨道来描述这种径向收缩（Reverse Breathing）。

Figure 4 清楚地展示了：在 NdNiO2 中，Ni $3d$ 轨道与 $4d$ 轨道之间的互信息（图中蓝色粗线）显著强于 CaCuO2。这解释了为什么三带模型失效：三带模型（包含 $d_{x^2-y^2}, p_x, p_y$）本质上是一个“能量空间”极其受限的模型，它假设电子的径向波函数是静态的。但在镍氧化物中，这种“呼吸”效应非常剧烈，如果不包含 $4d$ 轨道，就无法正确描述空穴移除后的多重态能级，进而无法准确预测配对势。

对未来超导研究的启示

本项研究实际上在暗示：镍氧化物超导可能不需要完全遵循铜氧化物的“ ZR 单态 + 单带 Hubbard 模型”路径。由于 $d^8$ 组态在镍中更具三维性且局部关联更强，我们可能需要考虑多轨道超导模型，甚至必须考虑 Nd $5d$ 轨道提供的金属化背景。这一发现为寻找更高转变温度的镍基超导材料提供了理论导向：即如何通过应力或化学取代来调控 Ni $d$ 轨道的径向分布，从而增强或抑制特定的关联通道。