来源论文: https://arxiv.org/abs/2606.00223v1 生成时间: Jun 07, 2026 12:46

执行摘要

自2019年无限层镍氧化物 $R\text{NiO}_2$（$R$ = 稀土元素，如 La, Nd, Pr）中发现非常规超导电性以来，凝聚态物理学界掀起了继铜氧化物和铁基超导体之后的又一轮高温超导研究热潮。由于其具有与铜氧化物高度相似的二维 $\text{Ni-O}$ 二氧化镍平面以及名义上的 $3d^9$ 电子组态，镍氧化物一直被视为理解高温超导机制的“完美参照物”。然而，镍氧化物独特的自掺杂效应（源于稀土 $5d$ 能带穿过费米面）使其低能物理图像显著复杂化。

本研究采用最先进的第一性原理计算方法，将**密度泛函理论（DFT）与动力学平均场理论（DMFT）相结合，首次系统、全面地对比研究了无限层镍氧化物 $\text{LaNiO}_2$ 在电子掺杂（$n$型）与空穴掺杂（$p$型）**两个极限维度下的物理演化规律。计算结果揭示了极强的不对称性：

自掺杂效应的不对称抑制与促进：空穴掺杂强烈抑制了稀土 $R(5d)$ 的自掺杂效应，使得费米面处的 $\Gamma$ 点电子口袋迅速消失；而电子掺杂则显著增大了 $R(5d)$ 电子口袋的尺寸，但令人意外的是，增加的电子并未有效填充 $\text{Ni-}d_{x^2-y^2}$ 能带，而是大量流入了稀土 $5d$ 轨道与间隙状态，导致 $\text{Ni-}d_{x^2-y^2}$ 轨道依然处于“事实上的空穴掺杂”状态。
磁基态的稳健性差异：反铁磁（AFM）序在空穴掺杂下被迅速压制，这与铜氧化物中观察到的现象高度一致；然而，在电子掺杂侧，反铁磁性作为基态显示出极高的稳健性，甚至比未掺杂母体更为稳定。这一发现源自自旋极化态密度（DOS）中独特的非线性载流子填充机制。
强关联效应的单带主导性：尽管在电子结构上表现为多带特征，动力学关联效应分析（通过自能函数与有效质量增强因子 $m^*/m$）表明，低能物理和强关联特征几乎完全由 $\text{Ni-}d_{x^2-y^2}$ 单个轨道主导，而 $\text{Ni-}d_{z^2}$ 轨道的关联极弱，且几乎不受掺杂影响。

本研究勾画出了镍氧化物完整的电子与空穴掺杂双侧相图轮廓，极大地丰富了人们对过渡金属氧化物强关联物性的认识，并为实验上探索 $n$ 型镍氧化物超导体以及拓展超导相图边界提供了决定性的理论支撑。

一、核心科学问题、理论基础与技术方法

1.1 核心科学问题：自掺杂与掺杂不对称性

在强关联电子体系中，过渡金属氧化物的物理特性对其过渡金属离子的 $d$ 电子数（即填充度）极度敏感。对于名义上具有 $d^9$ 组态的无限层镍氧化物 $R\text{NiO}_2$，最核心的科学争论之一在于其究竟是类铜氧化物的单带 Mott 物理，还是类似于铁基超导体的多带物理。

在未掺杂的母体 $\text{LaNiO}_2$ 中，费米面附近不仅存在主要的 $\text{Ni-}d_{x^2-y^2}$ 能带，还存在穿过费米面的 $\text{La-}5d$ 能带。这些 $\text{La-}5d$ 能带在布里渊区中心（$\Gamma$ 点）和顶点（$A$ 点）分别形成了电子口袋。这一物理现象被称为**“自掺杂效应”（Self-doping effect）**，它扮演了电荷库（Charge reservoir）的角色，自发地将部分电子从 $\text{Ni-}3d$ 轨道转移至 $\text{La-}5d$ 轨道，使得 $\text{Ni-}d_{x^2-y^2}$ 轨道偏离了半满的 $d^9$ 状态，显现出约 $d^{8.8}$ 左右的有效组态。这一机制使得镍氧化物的母体呈现金属星，而非铜氧化物中典型的共线反铁磁 Mott 绝缘体体行为。

本研究聚焦于以下核心科学问题：

当我们在体系中引入额外的空穴或电子时，自掺杂效应将如何演化？
多带特征（$\text{La-}5d$ 与 $\text{Ni-}3d$ 的共存）是否会导致电子和空穴掺杂之间出现物理性质的不对称？
在有/无晶格弛豫的情形下，强关联效应（局域库仑排斥、洪特耦合）在双侧掺杂相图中的动态响应是什么？

1.2 理论物理基础：多轨道 Hubbard 模型与动态关联

为了准确描述上述物理过程，单粒子近似下的独立粒子模型（如常规的 LDA/GGA 泛函）往往会因为无法正确处理强的电子局域库仑排斥而失效。必须诉诸于多轨道 Hubbard 模型：

$$ H = \sum_{k,\sigma, a, b} \epsilon_{ab}(k) c^{\dagger}_{k a \sigma} c_{k b \sigma} + H_{int} $$

其中，第一项表示非局域的紧束缚跃迁项（能带结构），第二项 $H_{int}$ 表示局域的多体库仑相互作用。对于 $\text{Ni-}e_g$ 轨道（包含 $d_{x^2-y^2}$ 和 $d_{z^2}$），$H_{int}$ 采用标准的 Kanamori 参数化形式：

$$ H_{int} = U \sum_{m} n_{m\uparrow} n_{m\downarrow} + (U' - J) \sum_{m < m'} n_{m} n_{m'} - 2J \sum_{m < m'} \mathbf{S}_m \cdot \mathbf{S}_{m'} + J' \sum_{m < m'} d^{\dagger}_{m\uparrow} d^{\dagger}_{m\downarrow} d_{m'\downarrow} d_{m'\uparrow} $$

此处，$U$ 为同轨道库仑排斥，$U' = U - 2J$ 为异轨道库仑排斥，$J$ 为洪特耦合（Hund’s coupling），$J'$ 为对跃迁振幅。局域电子的动力学行为通过动力学自能 $\Sigma_{ab}(i\omega_n)$ 来捕获，这是常规 DFT+U 计算（提供静态的、无频散的势修正）所无法企及的。

1.3 技术难点与挑战

自掺杂口袋的精细刻画：$\text{La-}5d$ 电子口袋尺寸较小，其与 $\text{Ni-}d_{z^2}$ 轨道存在杂化，但由于空间群对称性的限制，它与主要的 $\text{Ni-}d_{x^2-y^2}$ 能带无法杂化。准确模拟这一选择性杂化要求基组和投影区间的选择必须非常精确，任何轨道投影窗口的偏差都会导致物理图像失真。
虚拟晶格近似（VCA）与超胞（Supercell）计算的均衡：化学掺杂（如用 Sr 替代 La 引入空穴，用 Ce 或 Th 替代 La 引入电子）会带来双重效应——既改变了体系的载流子浓度，又由于离子半径的不同改变了局部晶格常数及配位多面体体积。如何在计算中既能模拟连续的掺杂浓度演化，又不会由于晶格常数简化导致局域失真，是一大技术难点。
双重计数（Double Counting）修正的敏感性：在 DFT+DMFT 中，由于 DFT 已经部分包含了静电 Hartley 相互作用和平均场近似下的交换关联能，必须通过双重计数项 $\Sigma_{dc}$ 对其进行扣除。在具有多带自掺杂的镍氧化物中，$\Sigma_{dc}$ 的选择对 $\text{Ni-}d$ 和 $\text{La-}5d$ 能带之间的相对能量位置（电荷转移能）有极大影响，本项研究中精细采用了 Held 公式来规避这一敏感性。

1.4 DFT+DMFT 联合计算框架细节

本工作构建了极其严谨的双层级计算流程：

第一层级（结构与非关联电子结构）：首先利用投影缀加波（PAW）方法（基于 VASP 软件包）对 $12.5\%$ 和 $25\%$ 的掺杂超胞进行完全晶格弛豫。同时利用虚拟晶格近似（VCA）插值出任意掺杂浓度下的平衡晶格参数。随后，利用全电势线性缀加平面波加局部轨道方法（FP-LAPW+lo，基于 WIEN2k 代码）在非磁性状态下求解 Kohn-Sham 方程，获得高精度的非磁性 Kohn-Sham 轨道能谱作为 DMFT 的输入基组。
第二层级（局域强关联多体物理）：利用 TRIQS 软件包及 TRIQS/DFTtools 工具箱，在 $\pm 10\text{ eV}$ 的超大能量窗口内构造包含 $\text{Ni-}d_{x^2-y^2}$ 和 $\text{Ni-}d_{z^2}$ 的局域 Wannier 函数。使用连续时间量子蒙特卡洛（CT-HYB）杂化膨胀不纯物解协器对多轨道杂化问题进行自洽求解。通过勒让德多项式（Legendre polynomials，40阶）对虚时格林函数 $G(\tau)$ 进行拟合以抑制蒙特卡洛统计噪声。计算在反温度 $\beta = 40\text{ eV}^{-1}$（相当于约室温 $290\text{ K}$）下自洽收敛，收敛标准为格林函数自洽差值小于 $10^{-5}$。

二、关键测试体系、计算所得数据与物理图像剖析

2.1 结构弛豫与晶格演化

在进行一切电子结构探讨前，首先必须明确掺杂对无限层四方结构（空间群 $P4/mmm$）物理常数的影响。基于 VASP PBE 泛函优化得到的 LaNiO2 基态晶格常数为：$a = 3.97\text{ \AA}$， $c = 3.37\text{ \AA}$。对于两种极端的完全掺杂相，其晶格常数演化规律为：

空穴掺杂端（以 $\text{SrNiO}_2$ 为插值终点）：由于 $\text{Sr}^{2+}$ 的引入及其伴随的空穴化，晶格体积收缩，晶格常数变为 $a = 3.99\text{ \AA}$， $c = 3.40\text{ \AA}$。
电子掺杂端（以 $\text{CeNiO}_2$ 为插值终点）：由于高价态 $\text{Ce}^{4+}$ 的引入及其引起的电子注入，晶格常数缩减为 $a = 3.92\text{ \AA}$， $c = 3.35\text{ \AA}$。这一晶格常数的线性插值构成了 VCA 方法的基础，同时也作为输入用于约束超胞计算。

2.2 费米面拓扑演化与载流子分布不对称性

图 2 直观展示了在 $0.2\text{ e}^-$（电子掺杂）和 $0.2\text{ h}^+$（空穴掺杂）下非磁性能带结构及费米面的演化。这是本工作最引人瞩目的成果之一：

2.2.1 空穴掺杂（$p$-type, 0.2 $\text{h}^+$）下的演化

在未掺杂母体中，由于自掺杂效应，穿过费米面的能带有：主要贡献的 $\text{Ni-}d_{x^2-y^2}$ 能带，以及在布里渊区 $\Gamma$ 点形成的 $\text{La-}d_{z^2}$ 电子口袋，以及在 $A$ 点形成的 $\text{La-}d_{xy}$ 电子口袋。当引入 $0.2$ 空穴时（如图 2(c) 所示）：

$\Gamma$ 点处的 $\text{La-}d_{z^2}$ 电子口袋彻底消失。它的能量位置被整体推高到费米能级 $E_F$ 之上。
$A$ 点处的 $\text{La-}d_{xy}$ 电子口袋也显著收缩，且伴随其轨道向更高能量方向平移。
这一现象表明：空穴掺杂强力压制了自掺杂效应，体系的低能能带结构变得高度“单一化”（即向铜氧化物的纯单带 $d_{x^2-y^2}$ 图像靠拢）。

2.2.2 电子掺杂（$n$-type, 0.2 $\text{e}^-$）下的演化

相反地，当注入 $0.2$ 电子时（如图 2(a) 所示）：

$\Gamma$ 点处的 $\text{La-}d_{z^2}$ 电子口袋和 $A$ 点处的 $\text{La-}d_{xy}$ 电子口袋急剧扩大。
然而，关键的物理细节在于电荷的重新分配（图 2(g)）。注入的电子并未简单填充到名义上的阻挫带 $\text{Ni-}d_{x^2-y^2}$ 中。事实上，由于体积收缩和自能重新分布的竞争，$\text{Ni-}d_{x^2-y^2}$ 轨道的实际填充数变化极小。新增的电荷大量涌入了稀土 $\text{La-}d$ 能带以及位于格点间隙的间隙 $s$ 态（Interstitial states）。
这表明：电子掺杂非但没有实现有效的 $\text{Ni-}3d$ 轨道的电子充填，反而极大地扩充了“电荷库”（即稀土 $5d$ 带和间隙态）的尺寸，维持了 $\text{Ni-}d_{x^2-y^2}$ 轨道偏离半满的强关联不对称状态。

下表（引自原文 Table I）给出了通过 DMFT 提取的在不同掺杂水平下，$\text{Ni-}e_g$ 两轨道的有效质量增强因子 $m^*/m$（反映了强关联导致的电子重整化效应）：

Doping Level (added $e^-$)	$m^*/m$ ($d_{x^2-y^2}$)	$m^*/m$ ($d_{z^2}$)
-0.20 (0.20 $\text{h}^+$)	2.75	1.21
0.00 (Undoped)	2.91	1.23
0.20 (0.20 $\text{e}^-$)	2.51	1.24

2.2.3 关联效应的不对称减弱

从上表数据中可以得出几个至关重要的多体物理结论：

$\text{Ni-}d_{x^2-y^2}$ 是绝对主导的强关联轨道：在所有掺杂浓度下，$d_{x^2-y^2}$ 的质量增强都在 $2.5$ 倍以上，最高在母体处达到 $2.91$。而 $d_{z^2}$ 轨道的质量重整化极小（$\sim 1.2$），行为接近弱关联的独立粒子，这直接证明了低能区强关联物理的单轨道特征。
掺杂引起关联性的不对称减弱：无论是空穴掺杂还是电子掺杂，都会减弱 $d_{x^2-y^2}$ 的关联度（$m^*/m$ 均下降）。但是，电子掺杂侧的减弱幅度（由 $2.91$ 降至 $2.51$，降幅达 $14\%$）显著大于空穴掺杂侧（由 $2.91$ 降至 $2.75$，降幅仅为 $5\%$）。这对应着自能函数 $\text{Im}\Sigma(i\omega_n)$ 斜率在低频区的改变（见图 10），表明电子掺杂侧由于费米面拓扑结构的急剧变化和多带杂化的增强，局域关联被更加有效地屏蔽。

2.3 磁基态稳定性与能量差异分析

为了厘清关联效应如何决定体系在宏观上的磁学行为，作者利用带有库仑势修正的 $\text{DFT}+U$ 框架模拟了三种最具代表性的磁有序构型：铁磁序（FM）、C型反铁磁序（C-AFM，面内反铁磁、面间铁磁）以及G型反铁磁序（G-AFM，共线反铁磁）。以非磁性（NM）状态的能量作为参照基准（$E_{NM} = 0$），计算得到的能量差 $\Delta E = E_{magnetic} - E_{NM}$ 随掺杂浓度的变化曲线见图 3(a)：

未掺杂母体 ($x=0$)：G-AFM 和 C-AFM 的能量极度接近且均远低于 NM 态（$\Delta E \approx -350\text{ meV/f.u.}$），其中 C-AFM 微弱占优，构成磁基态。
空穴掺杂端 ($x > 0$)：随着空穴浓度增加，$\Delta E$ 急剧向 $0$ 靠拢。在 $x=0.2$ 处，反铁磁态与 NM 态的能量差已被极大地压缩，磁矩（图 3(b)）由于空穴的稀释效应而发生坍塌。这说明空穴掺杂能够极其迅速地压制磁性长程序，使其进入超导温区所需的无序、强自旋涨落涨状态。
电子掺杂端 ($x < 0$)：相反，在电子掺杂侧，G-AFM 和 C-AFM 展现出了惊人的稳定性。随着电子浓度的增加，其与非磁性态的能量差不仅没有减小，反而进一步扩大（在 $-0.4$ 掺杂下，$\Delta E$ 甚至突破了 $-500\text{ meV/f.u.}$）。同时，镍离子的局域磁矩（$\mu_{\text{Ni}}$）在电子掺杂侧不降反升，在 $U=5\text{ eV}$ 下稳定保持在 $\sim 1.0\ \mu_B$ 附近（图 3(b)）。

2.3.1 磁性不对称演化的物理根源

为什么电子掺杂会导致磁基态如此稳健？这一现象的物理本质可以通过图 3(c) 展示的**自旋极化态密度（SP-PDOS）**来完美阐释：

在 C-AFM 状态下，$\text{Ni-}3d$ 轨道的多数自旋通道（Majority spin）几乎被完全填满，而少数自旋通道（Minority spin）则在费米面附近留有很大的空腔。
当注入额外的电子时：根据能量极小原理，新增的电子将专一性地填充少数自旋通道（Minority spin component）。由于少数自旋通道被电子填充，其对应的能带向费米面下方移动，但这并不会导致未配对电子自旋数的急剧减少，反而通过增强的局域交换分裂使得局域磁矩和反铁磁自旋对易作用得以巩固，这就维持并强化了反铁磁长程序。
当抽离电子（即空穴掺杂）时：空穴化并不会排空少数自旋通道，而是直接从几乎饱和的多数自旋通道中剥离电子。这导致多数自旋通道的电子态密度在费米面处急剧剥落，直接削弱了局域自旋交换分裂能，使得磁有序在一瞬间崩溃。

这种磁性响应的不对称性不仅是镍氧化物独特的自掺杂能带结构的必然结果，也对实验上通过电子掺杂调控超导电性提供了重要的启示：在电子侧，强大的反铁磁序将难以被涨落打破，超导电性可能会受到强烈压制，除非能够通过极强的门电压或特定衬底应力破坏其磁序对称性。

三、代码实现细节、算法复现指南与开源工具链

为了确保科研人员能够无缝复现本研究中报道的全部 DFT+DMFT 关联计算成果，本节将深入拆解计算流程中的每一步实现细节，并提供完整的开源 Repo 工具链及参数映射表。

3.1 理论计算工具链清单

结构优化与 VCA 构建：VASP (Version 5.4.4+) VASP 官网
高精度非磁性自洽场计算：WIEN2k (Version 19.1+) WIEN2k 官网
多体关联自洽求解器及投影工具：TRIQS (Version 3.1.x) 结合 TRIQS/DFTtools 以及基于连续时间强耦合杂化膨胀算法的杂化不纯物求解器 TRIQS/cthyb。
- TRIQS 核心库：GitHub - TRIQS/triqs
- DFTtools 投影接口：GitHub - TRIQS/dfttools
- cthyb 求解器：GitHub - TRIQS/cthyb

3.2 复现指南与关键输入参数设定

步骤 1：VASP 超胞弛豫与 VCA 模型建立

超胞设置：构建 $2 \times 2 \times 2$ 的三维 $\text{LaNiO}_2$ 晶格超胞（包含 32 个原子）。在 $12.5\%$ 的掺杂浓度下，使用 1 个 Sr 原子（空穴端）或 Ce/Th 原子（电子端）替换其中的 1 个 La 原子。
K点采样网格：对于超胞，使用 $6 \times 6 \times 6$ 均匀 Monkhorst-Pack 网格；对于单胞 VCA 晶格，采用 $12 \times 12 \times 12$ 密集网格。
能量截止（ENCUT）：统一设置为 $520\text{ eV}$，以保证电子和晶格应力在弛豫过程中完全收敛。
收敛判据：电子自洽步收敛至 $10^{-6}\text{ eV}$，离子力收敛标准为每个原子受力小于 $0.01\text{ eV/\AA}$。

步骤 2：WIEN2k 非磁性能带计算

势能泛函：选择标准的广义梯度近似泛函 GGA-PBE。
关键参数 $R_{MT} \times K_{max}$：设置为 $6.0$（这也是高精度关联计算中通常推荐的参数平衡点）。
Muffin-Tin 半径选择：
- $R_{MT}(\text{La}) = 2.50\text{ a.u.}$
- $R_{MT}(\text{Ni}) = 2.10\text{ a.u.}$
- $R_{MT}(\text{O}) = 1.62\text{ a.u.}$
网格密度：使用 $8 \times 8 \times 10$ 的 K 空间网格进行全布里渊区积分，生成后续 Wannier 投影所需的能谱数据文件（即 .struct, .vector, .energy 文件）。

步骤 3：TRIQS/DFTtools 局域投影与格林函数构造

在这一阶段，我们首先利用 python 脚本调用 DFTtools 的 Wien2kConverter 模块，将 WIEN2k 的能带输出转换为 HDF5 存档格式，并在指定的投影窗口内提取 Wannier 函数。

from triqs_dft_tools.converters.wien2k import Wien2kConverter

# 初始化转换器，指定 WIEN2k 的 case 名字
Converter = Wien2kConverter(filename="LaNiO2")

# 执行数据转换，投影窗口设定在费米面上下 10 eV 以内
Converter.convert_to_hdf5(
    dft_interface = 'wien2k',
    file_analysis = True,
    energy_window = [-10.0, 10.0]  # eV
)

投影窗口将提取名义上杂化最强的 $\text{Ni-}3d$ 轨道。在 HDF5 文件中，将定义一个包含 $d_{x^2-y^2}$ 和 $d_{z^2}$（属于 $e_g$ 流形）的局域关联块（block）。

步骤 4：TRIQS/cthyb 自洽 DMFT 求解

以下为自洽 DMFT 循环的核心 Python 伪代码，展示了 Held 公式双重计数（Equation 1）以及不纯物求解器的核心调用参数设置：

from triqs.gf import *
from triqs_dft_tools.sumk import SumkDFT
from triqs_cthyb import Solver
import numpy as np

# 设定多体相互作用参数
U_physical = 5.0  # eV
J_hund = 0.8      # eV
beta = 40.0       # 对应室温 290K

# 计算 Held 公式定义的双重计数项 (Equation 1 in text)
def calculate_held_dc(U, J, n_density, d_orbitals=2):
    # d 轨道总数，本研究中 eg 轨道 d = 2
    numerator = U + (d_orbitals - 1) * (U - 2*J) + (d_orbitals - 1) * (U - 3*J)
    denominator = 2 * d_orbitals - 1
    dc_sig = (numerator / denominator) * (n_density - 0.5)
    return dc_sig

# 实例化 SumKDFT 类用于晶格自洽整合
SK = SumkDFT(hdf_file="LaNiO2.h5", suite="wien2k")

# 初始化 CTHYB 不纯物求解器
S = Solver(beta=beta, gf_struct=SK.gf_struct_solver)

# 自洽循环控制
n_loop = 40
for iteration in range(n_loop):
    
    # 1. 晶格自洽步骤，计算非自洽格林函数 G0
    SK.symm_sort_gf()
    
    # 估算当前的格点局域电荷密度，用于 Held 双重计数修正
    current_density = SK.get_m_fractional_occupations()
    Sigma_dc = calculate_held_dc(U_physical, J_hund, current_density)
    
    # 注入双重计数自能项
    SK.set_dc(Sigma_dc)
    
    # 计算不纯物有效杂化格林函数
    G0_iw = SK.extract_G0()
    
    # 2. 调用 QMC 求解器在虚频（Matsubara 轴）上求解不纯物自能
    S.G0_iw << G0_iw
    S.solve(
        h_int = SK.construct_kanamori_hamiltonian(U=U_physical, J=J_hund),
        n_cycles = 1000000,          # 蒙特卡洛采样步数，保障高精度
        length_cycle = 40,
        n_warmup_cycles = 10000,
        random_name = "",
        max_time = -1                # 不限制最大计算时长，直到采样完毕
    )
    
    # 3. 将不纯物解得的自能拟合并反馈回 SumKDFT，进入下一步循环
    SK.put_Sigma(S.Sigma_iw)
    
    # 判断是否达到收敛阈值
    if iteration > 0 and np.abs(SK.density_diff) < 1e-5:
        print(f"DMFT self-consistent loop converged at iteration {iteration}!")
        break

在自洽计算完成后，可利用 Padé 近似或最大熵方法（MaxEnt）将虚频自能 $\Sigma(i\omega_n)$ 解析延拓到实频轴上，绘制出如图 4 所示的高分辨率、动量解析的谱函数 $A(k, \omega)$，从而直观读出有效带宽的收缩和能带弯曲。

四、关键引用文献与局限性批判

4.1 核心参考文献解析

本研究立足于以下数篇奠基性、关键性的前沿文献之上：

[1] Li, D. et al., Nature 572, 624 (2019): 该文献是镍氧化物高温超导领域的“大爆炸”工作，首次报道了在 Sr 掺杂的无限层 $Nd_{0.8}Sr_{0.2}NiO_2$ 薄膜中观察到超导电性。本项工作将其掺杂范围拓展到几乎未被探索的电子掺杂区，完成了完整的双侧掺杂物理图谱绘制。
[34] Armitage, N. P., Fournier, P. & Greene, R. L., Rev. Mod. Phys. 82, 2421 (2010): 铜氧化物超导体中电子-空穴掺杂不对称性的经典综述。本工作通过与该文献中铜氧化物的经典相图对比（见图 1），尖锐地指出：尽管在化学掺杂手段上看似对称，镍氧化物由于其独有的、扮演电荷储备库角色的 $R(5d)$ “旁观者能带（Spectator bands）”存在，导致其物理不对称性比铜氧化物更加显著。
[12] Kitatani, M. et al., npj Quantum Mater. 5, 59 (2020): 该理论计算证实了自掺杂效应对镍氧化物关联程度的自发调节作用。本工作在此基础之上，进一步利用 DMFT 动态自能技术精确量化了这种自掺杂在双侧外加化学掺杂下的动态消长机制。

4.2 本文理论工作的局限性与潜在挑战

尽管本研究在 DFT+DMFT 结合多带图像上取得了令人瞩目的理论进展，其模型及算法中仍不可避免地存在若干物理上的“理想化假设”和局限性：

1. 顶角氧缺陷（Apical Oxygen Vacancies）的简化

在实验制备无限层 $R\text{NiO}_2$ 的过程中，需要通过拓扑还原法（通常使用 $\text{CaH}_2$ 作为还原剂）强制将前驱体钙钛矿相 $R\text{NiO}_3$ 中的顶角氧（Apical oxygen）脱除。这一脱氧过程在实际样品中极难做到百分之百的完美。实验样品中普遍存在局域的顶角氧残留或额外的氧空位缺陷。残留的顶角氧会在局部引入强的晶格场畸变，使得原本呈纯平面的 $\text{Ni-O}_2$ 层产生局域的三维化倾向（使得本应退简并的 $d_{z^2}$ 轨道能量发生剧烈波动）。本研究假设的是完美的无限层结构，未将局域氧缺陷导致的化学无序纳入考量。

2. 氢杂质嵌入（Hydrogen Intercalation）效应的缺失

近年来的多项尖端实验和理论分析表明，在强化学还原剂 $\text{CaH}_2$ 的作用下，还原产物中往往会伴生显着的氢嵌入，形成 $R\text{NiO}_2\text{H}_x$ 间隙化合物。嵌入的 $\text{H}^-$ 离子通常占据顶角氧空位，具有 $-1$ 价，这会极大地改变周围镍离子的名义价态（进一步空穴化 Ni 轨道），并显著改变费米面附近的电子拓扑结构。本文的理论模型并未对氢原子的局域共存进行建模。

3. 三维晶格应力与衬底耦合的忽略

超导镍氧化物样品目前几乎无一例外地都以超薄薄膜的形式生长在 $\text{SrTiO}_3$（STO）衬底上。由于 STO 衬底晶格常数（$\approx 3.905\text{ \AA}$）与 $\text{LaNiO}_2$ 面内晶格常数存在晶格失配，薄膜会受到显着的双轴相干拉伸应力，导致晶格长宽比 $c/a$ 显著偏离体相。这种衬底应变不仅会调控能带的带宽，还会强烈改变 $\text{Ni-}d_{z^2}$ 与稀土 $5d$ 能带之间的杂化强度。本研究采用 VCA 和体相超胞进行计算，无法完全还原超薄膜样品中复杂的界面和应变物理。

五、补充：对超导机制与实验可观测性的深度探讨

为了将本项理论研究更紧密地与最前沿的实验观测及未来超导应用相联系，本节将对该理论计算结果在真实材料中的可观测性以及未来的实验探索路径进行深度拓展。

5.1 超导配对对称性的多轨道竞争

本工作最重要的多体结论之一是：即便在电子掺杂侧存在庞大的稀土 $5d$ 电子口袋，DMFT 有效质量和自能函数依然有力地支持单带（以 $\text{Ni-}d_{x^2-y^2}$ 为核心）的低能关联物理。这一结论对于界定超导配对对称性具有极强的方向性约束：

如果镍氧化物的强关联核心完全集中在 $d_{x^2-y^2}$ 轨道上，那么其最自然的超导配对对称性应该与铜氧化物高度一致，即呈现典型的单带 $d_{x^2-y^2}$ 波对称性（$d$-wave pairing），配对主要由局域的反铁磁自旋涨落所媒介。
然而，正如本工作所揭示，电子掺杂会使 $\text{La-}5d$ 与间隙态极度膨胀，这不排除在大电子掺杂极限下，$\text{La-}5d$ 电子口袋可能会与 $\text{Ni-}3d$ 发生实质性的多带超导竞争，导致超导间隙（Superconducting gap）展现出多带特征（例如具有符号改变的 $s_{\pm}$ 对称性），类似于铁基超导体。这表明在电子掺杂侧探索超导，极有可能开辟出一种介于铜氧化物与铁基超导体之间的新型非常规超导机制形态。

5.2 实验实现 $n$ 型掺杂镍氧化物的可行路径

由于无限层镍氧化物前驱体钙钛矿相 $R\text{NiO}_3$ 本身在热力学上处于亚稳态，制备难度已经极高。若要在 La 晶格位引入高价态的阳离子（如 $\text{Ce}^{4+}$ 或 $\text{Th}^{4+}$）来实现电子掺杂，实验上必须面临极其严苛的物理化学挑战：

多相共存与热力学不稳定性：Ce 或 Th 极易形成高度稳定的二氧化物（如 $\text{CeO}_2$, $\text{ThO}_2$）次相，这会导致在薄膜生长过程中出现严重的相分离。这要求在真空薄膜制备（如脉冲激光沉积，PLD）过程中，必须对氧分压、衬底温度以及激光能量密度进行极其精细的微调，以保证多组分固溶体的单相生长。
拓扑还原工艺的温和性控制：由于 $\text{Ce}^{4+}$ 的还原势与 $\text{Ni}^{3+}$ 存在巨大差异，在利用 $\text{CaH}_2$ 脱氧时，如果还原温度过高或时间过长，极易导致 Ce 被过度还原，进而引起三维晶格框架的全面坍塌。因此，实验上必须发展出更为温和、精细控制时间的“非破坏性”拓扑还原工艺。
非平衡态离子液体调控（Ionic Liquid Gating）：如果化学固溶体掺杂在材料制备上难以逾越热力学障碍，另一种极其有前景的替代路径是利用离子液体静电场调控技术。通过在纯 $\text{LaNiO}_2$ 薄膜表面施加极高的正向栅极偏压，可以在界面处诱导极高的静电累积电子。这种“无化学无序”的纯静电掺杂手段，能够直接绕开化学多相共存的难题，使实验人员在单一器件上即可横跨空穴到电子的双侧掺杂极限，从而直接对本工作所预测的电子-空穴不对称性进行最直接、最干净的实验检验。