来源论文: https://arxiv.org/abs/2603.20705v1 生成时间: Mar 24, 2026 12:20

0. 执行摘要

无限层镍酸盐作为一类备受关注的高温超导候选材料，其电子结构和超导机制的理解仍面临诸多挑战。近期角分辨光电子能谱（ARPES）实验揭示了镍酸盐中存在一个独特的间隙s轨道，且Ni-3d$_{x^2-y^2}$轨道的kz方向色散被显著压制，与密度泛函理论（DFT）的预测存在明显差异。本研究利用大规模行列式量子蒙特卡洛（DQMC）模拟，在一个扩展的三轨道Emery模型（d-p模型）中纳入了具有三维色散的间隙s轨道，构建了d-p-s模型，旨在显式地处理强电子关联效应和多轨道物理。

研究结果令人瞩目：首先，间隙s轨道引起的电子口袋在强关联作用下显著缩小，但在20%空穴掺杂下依然存在，其尺寸与ARPES实验观测结果高度吻合。其次，Ni-3d$_{x^2-y^2}$轨道的色散因强关联而被强烈重整化，表现出微弱的kz依赖性，这与ARPES测量结果一致，纠正了DFT对kz色散高估的问题。此外，与传统的d-p模型相比，d-p-s模型中的短程反铁磁关联显著增强。这些发现共同强调了强关联和多轨道效应在塑造无限层镍酸盐低能电子结构和多体关联中的关键作用，并为理解其潜在的非常规超导机制提供了新的视角。本研究不仅成功地将理论预测与实验观测结果相匹配，还为未来探索此类复杂量子材料提供了重要的理论基础和计算方法学支持。

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

1.1 核心科学问题

无限层镍酸盐 (infinite-layer nickelates, ILN) 是一类在结构上与铜氧化物高温超导体相似的材料，因此引起了广泛的关注，被认为是潜在的高温超导材料家族。然而，与铜氧化物中普遍认为CuO₂平面是其低能物理的核心不同，镍酸盐展现出更为复杂的轨道结构，例如稀土5d和4f轨道的影响。这种复杂性导致了镍酸盐中出现了一系列独特的物理现象，如近藤效应 (Kondo effect)、洪特耦合 (Hund’s coupling) 以及Ni-d轨道的自掺杂 (self-doping)。对这些复杂现象的理解对于阐明镍酸盐的超导配对对称性以及与铜氧化物的异同至关重要。

近年来，一系列ARPES实验为镍酸盐的电子结构提供了直接证据，这些实验尤其关注两个核心特征：

  1. 间隙s轨道的存在与特性： ARPES实验证实了间隙s轨道在一定掺杂浓度下依然存在，并对费米面产生贡献。这与一些早期DFT预测的s轨道完全处于费米能级以下的情况不符。如何准确描述这个间隙s轨道的行为，尤其是在强关联背景下，是理解镍酸盐电子结构的关键。
  2. Ni-3d$_{x^2-y^2}$轨道的kz色散与Lifshitz相变： ARPES观测到Ni-3d$_{x^2-y^2}$轨道具有有限的kz色散，但其幅度远小于许多DFT计算的预测。更重要的是，在Γ和Z点（即kz=0和kz=π处），并未观测到Lifshitz相变。这表明Ni-3d$_{x^2-y^2}$轨道的kz依赖性比DFT预测的弱得多。一些DFT研究通过调整s轨道的在位能量试图消除Γ点费米口袋的s轨道特征，但dx2-y2轨道的显著kz色散仍然难以避免。如何通过理论模型，在强关联的框架下，调和这些与DFT预测不符的实验观测，是当前面临的核心科学问题。

1.2 理论基础

为了解决上述问题，本研究选择了一个基于强电子关联的多轨道模型。传统的d-p模型（Emery模型）是描述铜氧化物强关联物理的成功范式，它考虑了Cu 3d$_{x^2-y^2}$轨道与O 2p$_{x/y}$轨道之间的杂化和关联。鉴于镍酸盐与铜氧化物在结构上的相似性，以及强关联效应在高温超导中扮演的核心角色，将Emery模型作为研究镍酸盐的理论起点是自然的选择。然而，为了更真实地反映镍酸盐的复杂性，需要对其进行扩展。

多轨道d-p-s模型： 本文的核心理论模型是在标准d-p模型基础上，引入了一个额外的间隙s轨道，从而构建了一个四轨道d-p-s模型。这个s轨道被认为是位于NiO₂平面之间的空间间隙中，并且具有三维色散。引入s轨道的动机直接来源于ARPES实验观测。该模型试图捕捉镍酸盐中Ni 3d$_{x^2-y^2}$、O 2p$_{x/y}$以及间隙s轨道之间的相互作用。重要的是，该模型考虑了强电子关联效应，特别是Ni 3d$_{x^2-y^2}$轨道上的在位库仑排斥（Udd），因为d轨道的局域性更强，关联效应通常更显著。

1.3 技术难点

在理论和计算上，研究具有强电子关联的多轨道模型面临诸多挑战：

  1. 强关联效应的精确处理： 强关联系统无法通过微扰理论有效处理。电子-电子相互作用导致的强关联效应是非微扰的，需要发展和应用能够精确处理这些效应的数值方法。DQMC就是其中一种，但它本身也存在挑战。
  2. 费米子符号问题 (Fermion Sign Problem)： 这是DQMC方法在模拟费米子系统时面临的最主要障碍。在某些参数区域（尤其是低温度、高掺杂或多轨道系统），蒙特卡洛采样的权重可能出现负值，导致统计平均的方差急剧增大，计算结果变得不可靠，甚至无法得到有效解。这严重限制了DQMC在低温和/或大掺杂区域的适用性，也限制了可模拟的系统尺寸。
  3. 多轨道和三维色散的复杂性： 引入额外的s轨道和三维色散（kz依赖性）会显著增加哈密顿量的复杂性，导致更大的矩阵计算和更高的计算成本。这使得在合理的时间内模拟足够大的系统尺寸（以捕捉kz色散的细节）变得困难。
  4. 解析延拓的挑战： DQMC模拟在虚时域进行，而物理量的实际观测（如谱函数）在实频域。将虚时数据转换为实频数据是一个数学上不适定 (ill-posed) 的问题，需要使用如最大熵方法 (Maximum Entropy Method, MEM) 等复杂的数值技术进行解析延拓。MEM本身也有其局限性，结果的精度可能受限于输入数据的质量和先验信息的选择。
  5. 有限尺寸效应： 任何数值模拟都只能在有限的系统尺寸上进行。为了从模拟结果中推断出宏观材料的物理性质，需要确保有限尺寸效应足够小，或者能够通过系统性地增加尺寸来推断无限大系统的行为。对于具有kz色散的系统，需要至少两层的周期性边界条件，而捕捉细致的kz依赖性则需要更大的三维尺寸，这进一步加剧了计算负担。

1.4 方法细节

本研究主要采用行列式量子蒙特卡洛 (DQMC) 方法，并辅以最大熵方法 (MEM) 进行谱函数的解析延拓。

1.4.1 哈密顿量与裸色散 (Hamiltonian and bare dispersion)

研究从一个包含Ni 3d$_{x^2-y^2}$、O 2p$_{x/y}$和间隙s轨道的四轨道哈密顿量开始：

$\hat{H}_0 = \sum_{\mathbf{k}, \sigma} \Psi^{\dagger}_{\mathbf{k}, \sigma} H(\mathbf{k}) \Psi_{\mathbf{k}, \sigma}$

其中，$\Psi^{\dagger}_{\mathbf{k}, \sigma} = (d^{\dagger}_{\mathbf{k}, \sigma}, p_{x,\mathbf{k}, \sigma}^{\dagger}, p_{y,\mathbf{k}, \sigma}^{\dagger}, s^{\dagger}_{\mathbf{k}, \sigma})$ 是在动量空间中的四轨道算符。哈密顿量矩阵 $H(\mathbf{k})$ 描述了各轨道之间的在位能量和杂化项：

$H(\mathbf{k}) = \begin{pmatrix} \varepsilon_d - \mu & -2it_{pd}S_x & 2it_{pd}S_y & 0 \\ -2it_{pd}S_x C_z & \varepsilon_p - \mu & 4t_{pp}S_xS_y & -4it_{ps}S_x C_z \\ 2it_{pd}S_y C_z & 4t_{pp}S_xS_y & \varepsilon_p - \mu & -4it_{ps}S_y C_z \\ 0 & -4it_{ps}S_x C_z & -4it_{ps}S_y C_z & \varepsilon_s(\mathbf{k}) - \mu \end{pmatrix}$

这里，$\varepsilon_d, \varepsilon_p, \varepsilon_s$ 是相应轨道的在位能量，$\mu$ 是化学势。$S_x = \sin(k_x/2)$, $S_y = \sin(k_y/2)$, $C_z = \cos(k_z/2)$。各项杂化参数如下：

  • $t_{pd}$: Ni 3d$_{x^2-y^2}$与周围氧p轨道之间的最近邻杂化。
  • $t_{pp}$: NiO₂平面内最近邻氧p轨道之间的杂化。
  • $t_{ps}$: 间隙s轨道与邻近氧p轨道之间的杂化。由于对称性限制，d轨道与s轨道之间没有直接杂化。

间隙s轨道的色散 $\varepsilon_s(\mathbf{k})$ 具有以下形式，体现了其三维特性：

$\varepsilon_s(\mathbf{k}) = \varepsilon_s + 2t^{ss}_{100}(\cos k_x + \cos k_y) + 4t^{ss}_{110}\cos k_x \cos k_y + 2t^{ss}_{001}\cos k_z + 4t^{ss}_{101}(\cos k_x + \cos k_y)\cos k_z + 8t^{ss}_{111}\cos k_x \cos k_y \cos k_z$

其中，$t^{ss}_{ijk}$ 是s轨道之间沿着不同晶向的最近邻跳跃项。模型中所有参数的数值（如表1所示）均根据无限层镍酸盐的实验观测和DFT计算进行选择，以确保裸色散能捕捉到镍酸盐的典型特征。特别地，Ep=5.0 eV 的选择旨在解释镍酸盐中相对较大的电荷转移能量，这与铜氧化物有所不同。较大的电荷转移能量意味着掺杂的空穴不太可能占据氧轨道，从而支持了忽略p轨道上在位库仑排斥（Upp）的假设，这在实践中也有助于缓解QMC符号问题。

对这个 Bloch 哈密顿量 $H(\mathbf{k})$ 进行对角化，可以得到系统的四个非相互作用能带，这些能带构成了后续包含强相互作用计算的起点。裸色散分析表明，Ni 3d$_{x^2-y^2}$轨道在X点和R点之间表现出显著的kz色散，而间隙s轨道在Γ和A点形成电子口袋。选择这些参数是为了与前人研究中镍酸盐的能带结构计算保持一致。

强关联项： 除了上述非相互作用哈密顿量，模型还包含了d轨道上的在位库仑排斥 $U_{dd}$。由于其局域性，d轨道上的关联效应通常被认为是最强的。

$H_{int} = U_{dd} \sum_i n_{d,i,\uparrow} n_{d,i,\downarrow}$

本研究中，$U_{dd}$ 的值设定为6.0 $t_{pd}$，这是一个典型的强关联强度值。而p轨道上的在位库仑排斥 $U_{pp}$ 被忽略，这是基于较大的电荷转移能量 $\varepsilon_p - \varepsilon_d$，使得掺杂的空穴主要在d轨道上，或者说p轨道的关联效应相对较弱。

1.4.2 行列式量子蒙特卡洛 (DQMC)

DQMC是一种强大的数值精确多体方法，能够在有限温度下无偏地评估系统各种物理性质。它在强关联模型中得到了广泛应用，尤其在理解铜氧化物的物理方面取得了巨大成功。DQMC的主要优势在于它以非微扰的方式处理电子-电子相互作用，不依赖于任何失控的近似，除了统计误差和有限尺寸效应之外，可以被认为是数值精确的方法。在此框架下，可以直接计算各种热力学量和关联函数，包括单粒子格林函数、自旋和电荷关联以及配对磁化率等。

DQMC的基本思想是：

  1. 虚时演化： 将系统的配分函数 $Z = \text{Tr}[e^{-\beta H}]$ 离散化为一系列短的虚时切片，即 $e^{-\beta H} = (e^{-\Delta\tau H})^L$，其中 $\beta = L\Delta\tau$ 是逆温度。
  2. Hubbard-Stratonovich 变换： 为了处理非对角的两体相互作用项（如 $U_{dd} n_{d,\uparrow} n_{d,\downarrow}$），引入辅助场，将两体相互作用转化为与单体哈密顿量耦合的辅助场。这样，每个虚时切片上的哈密顿量都变成了费米子的单体问题，可以通过行列式有效地计算。
  3. 蒙特卡洛采样： 辅助场的配置通过蒙特卡洛方法进行采样，每个配置的权重由费米子行列式决定。通过对大量辅助场配置的采样，可以计算出物理量的平均值。
  4. 格林函数计算： DQMC可以直接计算虚时单粒子格林函数 $G_{\alpha\beta}(i, j, \tau)$，这是谱函数解析延拓的输入数据。

尽管DQMC功能强大，但其臭名昭著的符号问题 (sign problem) 限制了可访问的参数空间，使得本文中的许多结果都仅限于相对较高的温度（例如 $\beta t_{pd} = 10.0$）和有限的晶格尺寸。为了平衡计算复杂性（也受限于QMC符号问题）和包含重整化能带结构的kz依赖性，研究采用了多种晶格尺寸，例如8×8×2（用于能带结构）和4×4×2（用于更耗时的自旋关联函数）。对于更高kz分辨率的计算，也采用了4×4×4的晶格。

1.4.3 最大熵方法 (Maximum Entropy Method, MEM)

DQMC模拟在虚时 ($\tau$) 域生成格林函数 $G_{\alpha\beta}(\mathbf{k}, \tau)$。为了获得实频 ($\omega$) 域的谱函数 $A_{\alpha}(\mathbf{k}, \omega)$，需要进行解析延拓。这是一个典型的弗雷德霍姆第一类积分方程问题：

$G_{\alpha\alpha}(\tau) = \int_{-\infty}^{\infty} d\omega \frac{e^{-\tau\omega}}{1 + e^{-\beta\omega}} A_{\alpha}(\omega)$

MEM是一种常用的、鲁棒的解析延拓方法，它通过最大化信息熵来寻找最符合虚时数据且最“平滑”的实频谱函数。MEM的核心在于选择合适的先验信息（默认模型），并在拟合数据精度的同时，避免产生过多的结构。它能够有效处理数据的噪声和不完整性，从虚时格林函数中提取出可靠的实频谱函数。本文利用轨道分辨的局域态密度（LDOS）$A_{\alpha}(\omega)$ 和动量分辨的谱函数 $A_{\alpha}(\mathbf{k}, \omega)$ 来表征关联电子态的色散特性，这些结果可以直接与ARPES实验进行比较。

通过以上理论模型和计算方法，研究团队得以在强关联框架下，对无限层镍酸盐的电子结构和磁性关联进行深入探测，并成功解释了ARPES实验中的关键观测结果。

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

本研究旨在通过DQMC模拟，将理论预测与无限层镍酸盐的ARPES实验观测结果相匹配。为此，研究聚焦于两种具有代表性的轨道占据情况，以近似对应实验中的未掺杂和20%空穴掺杂情况。

基准体系与掺杂情况：

  1. 未掺杂情况 (Parent compound LaNiO$_2$): 对应于d-p平面上的空穴密度 $n_{dp} \approx 1.104$，以及s轨道上的电子密度 $n_s \approx 0.053$。这里的 $n_{dp} = 1.0$ 对应于传统d-p模型的半填充情况，$n_{dp} > 1.0$ 表示空穴掺杂。s轨道电子密度 $n_s \approx 0.05$ 维持在一个较小的数值，与ARPES实验中s轨道特征较弱的观测结果相符。
  2. 20%空穴掺杂情况 (20% Ca-doped case): 对应于d-p平面上的空穴密度 $n_{dp} \approx 1.259$，以及s轨道上的电子密度 $n_s \approx 0.049$。这一掺杂水平在实验中观察到超导性出现。

为了兼顾计算复杂性（受限于QMC符号问题）和kz依赖性的描述，研究采用了不同尺寸的晶格：8×8×2（用于能带结构）和较小的4×4×2（用于更耗时的自旋关联函数）。对于更高kz分辨率的分析，也采用了4×4×4晶格，尽管这伴随着更严重的符号问题，导致逆温度略低（$\beta=8.0$ 而非 $\beta t_{pd}=10.0$）。

2.1 裸色散与强关联下的重整化

裸色散： 图2展示了未引入相互作用时的d$_{x^2-y^2}$和s轨道能带结构，投影在高对称路径上。D$_{x^2-y^2}$轨道在X和R点之间表现出显著的kz色散，而间隙s轨道在Γ和A点形成电子口袋。虚线表示的费米能级对应于未掺杂和20%空穴掺杂情况。这些裸色散是理解关联效应影响的基准。

局域态密度 (LDOS) 重整化 (图3)： 图3比较了未掺杂和20%空穴掺杂情况下的轨道分辨LDOS。除了围绕-4 eV的宽O-2p能带外，研究主要关注费米能级附近的特征。

  • s轨道： 两种掺杂情况下，小的 $n_s \approx 0.05$ 反映为s轨道谱中围绕-1 eV的一个小峰。随着掺杂，这个小峰从大约-1 eV向-0.5 eV移动，与s电子口袋的缩小以及ARPES实验观测结果一致。这表明，即使在有强关联的情况下，s轨道的电子密度仍然很小，其对费米面的贡献也相应减小，尤其是在Γ点。
  • d轨道： 与之前铜氧化物的DQMC模拟相比，d轨道在费米能级附近的准粒子峰变得明显更尖锐。这反映了无限层镍酸盐中d$_{x^2-y^2}$轨道更强的局域性。掺杂后（图3b），上Hubbard带（UHB）的谱权重向低能区域转移，同时与d轨道杂化的p轨道峰也向低能移动。这与铜氧化物在空穴掺杂下的行为相似。
  • 总而言之： 强关联使得s轨道的电子占据率降低，d轨道准粒子峰更尖锐，并且谱重分布随着掺杂发生显著变化，这些都指向了强关联在镍酸盐电子结构中的核心作用。

2.2 s轨道动量分辨谱函数 (图4)

图4展示了沿高对称路径的间隙s轨道的动量分辨谱函数 $A_s(\mathbf{k}, \omega)$。白色虚线表示非相互作用下的色散。

  • Γ点电子口袋： 无论掺杂与否，一旦考虑强关联，s轨道在Γ点的电子口袋被重整化到费米能级以上。这意味着Γ点附近的s轨道特征在ARPES实验中是缺失的，与实际观测吻合。
  • A点电子口袋： A点处的电子口袋在强关联作用下显著缩小。更重要的是，即使在20%空穴掺杂下，这个口袋也并未完全消失，这与ARPES实验的观测结果高度一致。20%的掺杂水平恰好对应于镍酸盐中出现超导的区域，这暗示了s轨道在超导中的潜在作用。
  • X点附近的展宽： 谱函数在X点附近显示出较大的展宽，这暗示了该区域s轨道与d$_{x^2-y^2}$轨道之间存在间接影响，这种影响可能通过p轨道介导。

这些结果表明，强关联效应不仅改变了s轨道的在位能量，还深刻地重塑了其动量空间中的分布，使其与实验观测更加一致，尤其是在费米面附近的特征。

2.3 d$_{x^2-y^2}$轨道kz色散的抑制 (图5与图6)

强电子关联的另一个重要后果是d轨道能带kz色散的强烈抑制。在8×8×2晶格的DQMC模拟中，强烈的 $U_{dd} = 6.0$ eV 重整化了费米能级附近d轨道的带宽，导致其低能行为偏离了DFT结果。

  • 未掺杂情况 (图5a)： 沿高对称路径（kz=0和kz=π平面），d$_{x^2-y^2}$轨道色散并未显示出显著差异，除了Γ和Z点附近。虽然这里的空穴口袋比ARPES观测到的要小，但整体上与实验结果一致。
  • 20%空穴掺杂情况 (图5b)： 即使在实验相关的20%空穴掺杂水平下，色散仍未显示出Lifshitz相变的迹象。UHB的大部分谱权重转移到低能区域，使得其谱函数更接近掺杂载流子存在的非相互作用色散。图5清晰地提供了数值证据，表明d$_{x^2-y^2}$轨道的kz依赖性比DFT计算预测的要弱得多，且没有出现Lifshitz相变。
  • 高kz分辨率晶格 (图6)： 为了排除有限尺寸效应的影响，研究还对4×4×4晶格（具有更高kz分辨率）进行了模拟。图6显示，非相互作用下的d$_{x^2-y^2}$轨道在X、X’、R点之间存在明显的kz差异。然而，实际的谱函数，无论未掺杂（图6a）还是掺杂（图6b），在固定kz的任何高对称路径上仍然表现出非常微弱的kz色散。这进一步证实了强关联能够有效地抑制d$_{x^2-y^2}$轨道的kz依赖性，使其行为更趋近于二维。

2.4 反铁磁自旋关联的增强 (图7)

为了获取更多关于镍酸盐物理的信息，研究还探讨了动量分辨的零频自旋磁化率 $\chi_{\alpha_1\alpha_2}(\mathbf{q}, i\nu=0)$，这是一个理解多轨道系统行为的重要两体量。

  • d-p-s模型与d-p模型比较： 图7展示了d$_{x^2-y^2}$轨道静态自旋磁化率 $\chi_{dd}(\mathbf{q}, i\nu=0)$ 随空穴掺杂 $\delta$ 的变化。与d-p模型（虚线）相比，d-p-s模型（实线）在面内 $\mathbf{q}=(\pi,\pi,0)$ 反铁磁（AFM）自旋-自旋关联方面表现出更强的特性。同时，铁磁关联在 $\mathbf{q}=0$ 时略有减弱，而在 $\mathbf{q}=(\pi,0)$ 时几乎没有变化。这些观测结果表明，引入s轨道能够增强系统的反铁磁涨落。
  • 局域磁矩增强： 图7的插图显示了d$_{x^2-y^2}$轨道的局域磁矩 $(m^2)_d$，它也随着s轨道的引入而增强。这表明s轨道确实使d$_{x^2-y^2}$轨道更加局域化，从而促进了彼此之间的AFM关联。
  • kz依赖性： 图7还直接比较了d-p-s模型中静态自旋磁化率的kz依赖性。注意 $\mathbf{q}=(0,0,0)$ 和 $\mathbf{q}=(0,0,\pi)$ 之间的曲线几乎重叠，同样 $\mathbf{q}=(\pi,0,0)$ 和 $\mathbf{q}=(\pi,0,\pi)$ 之间也几乎重叠。更重要的是，立方 $\mathbf{q}=(\pi,\pi,\pi)$ AFM自旋关联几乎与面内绿色实线相同。所有这些特征进一步证实了d-p-s模型中d$_{x^2-y^2}$轨道准二维的本质。

性能数据总结：

  • 计算复杂性： 尽管DQMC在多轨道和三维色散模型上的计算成本很高，但通过采用不同的晶格尺寸和在一定温度范围内进行模拟，成功获得了有物理意义的结果。
  • 符号问题： 符号问题是主要限制，它决定了可访问的逆温度（$\beta t_{pd}=10.0$ 或 $\beta=8.0$）和晶格尺寸（8×8×2，4×4×4）。这使得在极低温度或更大掺杂下研究超导配对等现象成为挑战。
  • 数据量： 大规模模拟生成了大量的虚时格林函数数据，并通过MEM进行处理，以生成高分辨率的谱函数。

这些计算数据和性能指标共同展示了DQMC在处理复杂强关联多轨道系统方面的能力，同时也突出了其在面对符号问题时的局限性。

本研究使用的行列式量子蒙特卡洛（DQMC）方法是一个复杂的数值技术，虽然论文中并未直接提供具体的代码实现或开源仓库链接，但可以根据其描述和标准DQMC算法的原理推断其实现细节和复现所需步骤。鉴于DQMC在多体物理领域的广泛应用，许多通用或特定模型的DQMC代码库都已开源或可供学术界使用。

3.1 代码实现细节推断

DQMC算法的实现通常包括以下核心模块：

  1. 哈密顿量构建 (Hamiltonian Construction)：

    • 输入参数： 读入跳跃参数（$t_{pd}, t_{pp}, t_{ps}, t^{ss}_{ijk}$）、在位能量（$\varepsilon_d, \varepsilon_p, \varepsilon_s$）、化学势（$\mu$）和库仑相互作用强度（$U_{dd}$）。
    • 动量空间矩阵 $H(\mathbf{k})$： 根据论文中的公式（3）和（5），构建在布里渊区各个动量点 $\mathbf{k}$ 上的哈密顿量矩阵 $H(\mathbf{k})$。这需要计算 $S_x, S_y, C_z$ 等几何因子。
    • 实空间哈密顿量： 为了DQMC的虚时演化，通常需要构建实空间的跳跃矩阵 $T_{ij}$。这可以通过对 $H(\mathbf{k})$ 进行傅里叶逆变换获得，或者直接从实空间最近邻和次近邻跳跃项构建。

  2. 虚时演化与Hubbard-Stratonovich (HS) 变换：

    • Trotter 分解： 将总的逆温度 $\beta$ 分成 $L$ 个虚时切片 $\Delta\tau = \beta/L$。每个切片上的演化算符 $e^{-\Delta\tau H}$ 通过Trotter分解近似处理，将相互作用项和动能项分离。
    • HS 变换： 对相互作用项 $U_{dd} n_{d,\uparrow} n_{d,\downarrow}$ 引入辅助场。对于 $U_{dd} > 0$，常用的方法是实HS变换，将四费米子项转化为二费米子项与一个玻色子辅助场耦合的形式。这个辅助场在每个晶格点、每个虚时切片上都有一个值。
    • 单粒子格林函数： 在每个辅助场配置下，系统成为一个非相互作用的费米子系统，其格林函数可以解析地表示为矩阵指数 $B(\Delta\tau, \{h_i\})$ 的乘积 $C = B_L B_{L-1} \cdots B_1$。其中 $B_l$ 包含了动能项和与辅助场耦合的单体势能项。

  3. 蒙特卡洛采样循环：

    • 初始化： 随机初始化辅助场配置。
    • 热化： 运行一定数量的蒙特卡洛步，使系统达到平衡态。
    • 采样与测量： 在热化后，继续运行蒙特卡洛步并收集测量数据。在每次采样步中，随机选择一个晶格点和虚时切片上的辅助场值进行翻转或更新。
    • 接受/拒绝准则： 根据Metropolis算法，计算接受新辅助场配置的概率，该概率依赖于费米子行列式 $P = \det(I + C_{\text{new}}) / \det(I + C_{\text{old}})$。如果概率为负（符号问题），则处理方式取决于具体的DQMC变体（例如，通常将其视为零概率接受）。
    • 稳定性： 由于矩阵乘积的数值不稳定性，通常需要采用如QR分解或奇异值分解 (SVD) 等技术，在每个虚时切片或周期性地对格林函数矩阵进行稳定化。

  4. 物理量计算：

    • 虚时格林函数： 在每次采样步中，根据稳定化后的 $B$ 矩阵和 $C$ 矩阵，计算不同轨道组合（d-d, p-p, s-s, d-p等）的实空间和动量空间虚时格林函数 $G_{\alpha\beta}(\mathbf{k}, \tau)$。
    • 自旋关联函数： 计算零频自旋关联 $\chi_{\alpha_1\alpha_2}(\mathbf{q}, i\nu=0)$。这需要计算实空间中的自旋-自旋关联函数，然后进行傅里叶变换到动量空间。
    • 局域磁矩： 测量每个晶格点上的局域磁矩 $m^2_d = \langle (n_{d,\uparrow} - n_{d,\downarrow})^2 \rangle$。

  5. 解析延拓 (Analytic Continuation)：

    • 输入： 虚时格林函数 $G_{\alpha\alpha}(\tau)$。
    • 方法： 本研究使用最大熵方法 (MEM)。MEM的实现通常需要一个线性代数库（如LAPACK），一个优化库（用于最大化熵函数），以及一个默认模型（通常为常数或高斯展宽的能带）。
    • 输出： 实频谱函数 $A_{\alpha}(\omega)$ 或 $A_{\alpha}(\mathbf{k}, \omega)$。

3.2 复现指南

复现本研究结果需要一个功能完善的DQMC代码，并熟悉多轨道模型和三维晶格的设置。以下是复现的关键步骤：

  1. 选择DQMC框架： 可以选择从头开始编写DQMC代码，但这通常耗时且容易出错。更实际的方法是使用现有的开源DQMC库或框架，并对其进行修改以适应本研究的特定哈密顿量。
    • 现有框架： 例如，ALF (Algorithms for Lattice Fermions) 是一个用Fortran 90编写的通用DQMC软件包，支持各种哈密顿量和相互作用。另一个是基于C++的 QMC-BLAS，提供了QMC的核心数值库。还有一些基于Python的框架，如 PyQMC，虽然主要用于变分蒙特卡洛，但其数值工具链也有借鉴意义。
  2. 模型哈密顿量实现：
    • 跳跃矩阵： 在选定的DQMC框架中，需要精确实现论文表I中的所有跳跃参数 $t_{pd}, t_{pp}, t_{ps}$ 以及所有s-s跳跃项 $t^{ss}_{ijk}$。这些参数定义了实空间中的非相互作用哈密顿量。
    • 在位能量与化学势： 设置正确的在位能量 $\varepsilon_d, \varepsilon_p, \varepsilon_s$ 和化学势 $\mu$。化学势通常需要迭代调整，以达到目标电子或空穴掺杂密度。
    • 相互作用项： 实现d轨道上的在位库仑排斥 $U_{dd}$。由于忽略了 $U_{pp}$，这简化了HS变换的复杂性。
  3. DQMC参数设置：
    • 晶格尺寸： 根据论文设置8×8×2或4×4×4的晶格。
    • 逆温度： 设置 $\beta t_{pd} = 10.0$ (或 $\beta=8.0$ 用于4×4×4晶格)。确定虚时切片数 $L$ 以保证 $\Delta\tau$ 足够小（通常 $\Delta\tau U < 1$ 或 $\Delta\tau t < 1$）。
    • 测量频率： 设定热化步数和测量步数，以获得足够的统计样本。论文中没有给出具体数字，但通常需要数千到数万热化步和数十万测量步。
    • 稳定性参数： 设置矩阵稳定化（如QR分解）的频率。
  4. 数据分析与可视化：
    • 解析延拓： 获取虚时格林函数数据后，需要使用MEM程序进行解析延拓。开源的MEM实现包括 MaxEnt 或一些数值库中的MEM模块。
    • 谱函数可视化： 使用Python (matplotlib, numpy)、Julia (Plots.jl, DrWatson.jl) 或Matlab等工具对动量分辨谱函数、LDOS和自旋关联进行可视化，并与论文图2-7进行比较。

论文中并未明确提及使用了哪些具体的开源软件包或提供了代码库链接。这是学术论文的常见做法，因为重点在于方法和结果，而非代码实现本身。然而，要复现这样的研究，通常会依赖以下类型的工具和库：

  • DQMC框架：

    • ALF (Algorithms for Lattice Fermions): 一个广泛使用的开源DQMC软件包，支持多种哈密顿量和晶格几何，用Fortran 90编写。其模块化设计允许用户添加自定义哈密顿量。ALF的官方网站通常提供其最新版本和文档：https://alf.physik.uni-wuerzburg.de/
    • PyQMC: 虽然主要侧重于变分蒙特卡洛，但其底层数值工具和Python接口可能对DQMC的快速原型开发有所帮助。其仓库可能在GitHub上：https://github.com/WagnerGroup/pyqmc
    • QuSpin: 一个Python库，用于精确对角化和时间演化，虽然不是DQMC，但对于小尺寸系统的验证和概念理解很有用。仓库：https://qspin.github.io/QuSpin/

  • 数值库：

    • BLAS (Basic Linear Algebra Subprograms) / LAPACK (Linear Algebra PACKage): 高效的线性代数例程库，是DQMC中矩阵运算（乘法、求逆、分解）的基石。这些通常是操作系统自带或通过Intel MKL、OpenBLAS等优化库提供。
    • Boost (C++ Libraries): 如果DQMC代码是用C++编写的，Boost库提供了许多通用工具，如多线程、数学函数等。

  • 解析延拓 (MEM) 工具：

    • MaxEnt Code: 有一些由研究团体开发的MaxEnt代码，例如由Scalettar小组维护的MaxEnt代码，常用于QMC数据的解析延拓。具体链接可能需要通过学术搜索获得。
    • Numerical Recipes: 书中包含MEM的算法描述，但需要自行实现。
    • SciPy/NumPy (Python): Python生态系统提供了丰富的数值计算和优化工具，可以用来实现MEM。

  • 数据可视化：

    • Matplotlib (Python): 强大的绘图库，用于生成论文中类似的2D和伪彩色图。
    • Gnuplot: 另一个常用的命令行绘图工具。

重要提示： 鉴于多轨道和三维系统的复杂性，以及符号问题的存在，复现本研究需要深厚的DQMC理论知识和丰富的编程经验。如果没有直接可用的代码，从头开始实现并验证其正确性将是一个巨大的工程。

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键引用文献

本研究建立在大量先前工作的基础上，并引用了多个关键文献来支持其动机、方法和结果的解释。以下是一些重要的引用及其贡献：

  1. 铜氧化物超导的理论基础 (Refs. 1-4)： 这些文献建立了铜氧化物高温超导的理论框架，特别是Anderson的共振价键理论（Ref. 1）、Zhang-Rice单线态（Ref. 2）和Emery模型的提出（Ref. 3）。这些工作构成了本研究中d-p模型的基础，并提供了比较镍酸盐与铜氧化物的背景。
  2. 无限层镍酸盐的复杂轨道结构 (Refs. 5-10)： 这些文献强调了镍酸盐复杂的轨道物理，包括稀土5d/4f轨道、近藤效应、洪特耦合和自掺杂等。这些研究表明镍酸盐的物理远不止简单的单带模型，需要多轨道模型来捕捉其丰富性。例如，Ref. 6和Ref. 9详细讨论了Ni-d轨道与巡游电子的杂化，以及DFT对镍酸盐电子结构的预测。
  3. ARPES实验证据 (Refs. 18-19)： 这是本研究的核心实验驱动力。Ref. 18和Ref. 19是最近的ARPES实验文章，它们证实了间隙s轨道的存在和其在费米能级附近的贡献，同时观测到d$_{x^2-y^2}$轨道kz色散的弱化和Lifshitz相变的缺失。本研究的许多结果正是为了与这些实验观测结果保持一致。
  4. DQMC方法学 (Refs. 27-32, 37-39)： 这些文献详细介绍了DQMC算法及其在强关联系统中的应用，包括其理论基础、实现细节和在 Hubbard 模型及 Emery 模型中的成功案例。它们为本研究中DQMC模拟的实施提供了方法学支持。特别是Ref. 30、31、34讨论了在三轨道Emery模型中应用DQMC的经验。
  5. 最大熵方法 (MEM) (Refs. 40-43)： 这些文献描述了MEM的理论和应用，以及如何将虚时数据解析延拓到实频域，这对于从DQMC结果中提取谱函数至关重要。
  6. 镍酸盐中的kz色散与DFT预测 (Refs. 6, 9, 20-23)： 这些文献讨论了DFT对镍酸盐能带结构的预测，特别是关于d$_{x^2-y^2}$轨道显著的kz色散和Γ点s轨道电子口袋的特征。本研究的许多发现（如kz色散的抑制）与这些DFT预测形成对比，并试图解释与ARPES的差异。
  7. 镍酸盐中的磁性关联 (Ref. 44)： Ref. 44是近期关于无限层镍酸盐中局部磁交换作用起源的核磁共振（NMR）测量，它报道了比传统理论预测更强的自旋关联，为本研究中反铁磁关联增强的发现提供了实验背景。

这些引用共同构建了本研究的科学背景和技术框架，使其能够在一个坚实的基础上探索镍酸盐的复杂物理。

4.2 对这项工作局限性的评论

尽管本研究通过DQMC模拟取得了显著的进展，成功地解释了ARPES实验中的关键观测结果，并为理解无限层镍酸盐的强关联物理提供了深入洞察，但任何数值模拟都存在固有的局限性。对这些局限性的认识对于正确评估研究结果的适用范围和指导未来工作至关重要。

  1. 费米子符号问题 (Fermion Sign Problem)： 这是DQMC方法最主要的挑战，也是本研究中最大的局限。符号问题限制了：
    • 可访问的温度： 无法达到极低温度。论文中使用的 $\beta t_{pd} = 10.0$ 相对较高，这意味着研究仍处于一个相对较高的温度区域，可能无法完全捕捉到与超导配对直接相关的极低温物理。超导通常发生在极低温度下，此时量子涨落占据主导地位，而符号问题恰好限制了DQMC在这些区域的有效性。
    • 可访问的掺杂范围： 无法模拟更大范围的掺杂，特别是低空穴掺杂区域或电子掺杂区域可能更受符号问题的影响。这限制了对整个超导相图的探索。
    • 可访问的晶格尺寸： 为了缓解符号问题，模拟通常被迫在较小的晶格尺寸上进行（如8x8x2，4x4x4）。虽然研究者已经尝试增加kz分辨率，但仍然无法探索更高分辨率或更大尺寸的系统，这可能导致有限尺寸效应影响对精确kz色散的捕捉。
  2. 有限尺寸效应： 即使采用了8x8x2和4x4x4等尺寸，对于捕捉长程关联、精确描述三维色散（尤其是细微的Lifshitz相变）以及在更大的布里渊区中分辨特征，仍然可能存在不足。例如，图6中对kz色散的分析，研究者也承认受限于QMC符号问题，无法探索更高分辨率的更大晶格。
  3. 近似处理与参数选择：
    • 忽略p轨道上的U$_{pp}$： 论文中基于较大的电荷转移能量而忽略了p轨道上的库仑相互作用 $U_{pp}$。虽然有其合理性，但这仍是一个近似。在某些情况下，p轨道上的关联效应也可能不可忽略，并对最终结果产生影响。
    • 参数的唯象性： 模型中的跳跃积分和在位能量参数是唯象的，旨在与实验和DFT结果定性匹配，而非从头算得到。这种参数的选择可能会影响模型的普适性和预测能力。
  4. MEM的固有挑战： 最大熵方法是一个不适定问题，其结果的精度和可靠性可能受限于输入数据（虚时格林函数）的噪声水平、数据量以及所选默认模型的影响。虽然MEM是目前最先进的解析延拓方法之一，但仍然存在一定的数值不确定性。
  5. 未直接计算超导配对强度： 尽管研究发现反铁磁自旋关联增强，并指出这与非常规超导中的配对机制密切相关，但并未直接计算或证明d$_{x^2-y^2}$轨道相关的d波配对强度。论文中仅提到“我们的模拟确实显示出d$_{x^2-y^2}$轨道相关的d波配对强度增强的一些迹象（未显示）”，这有待未来更专门的计算来证实。
  6. 静态自旋磁化率： 本研究主要关注零频（静态）自旋磁化率。然而，动态自旋涨落（在有限频率下）在介导超导配对中可能扮演更关键的角色。静态结果提供了重要信息，但对动态涨落的探索会更全面。

总的来说，本研究在强关联多轨道模型研究领域迈出了重要一步，但其结论的强度和普适性仍受制于DQMC方法的固有局限性。未来的工作需要克服这些挑战，例如通过发展新的QMC算法或结合其他多体方法（如动力学平均场理论DMFT、变分蒙特卡洛VMC等），以探索更低温度、更大掺杂和更大尺寸的系统，从而更全面地理解无限层镍酸盐的超导机制。

5. 其他必要的补充

5.1 本研究的深远影响与贡献

这项工作在无限层镍酸盐的研究领域具有多方面的深远影响和重要贡献：

  1. 弥合理论与实验鸿沟： 本研究最显著的贡献在于成功地将强关联理论模型（d-p-s模型）的预测与ARPES实验观测结果相匹配。长期以来，DFT计算在预测镍酸盐的s轨道电子口袋大小和d$_{x^2-y^2}$轨道kz色散方面存在偏差。本研究明确展示了强电子关联如何显著重整化这些特征，使理论与实验高度吻合。这为理解这类复杂材料的真实电子结构提供了坚实的理论支撑。
  2. 揭示间隙s轨道的双重角色： 传统观念可能将间隙s轨道视为一个独立的电子口袋。本研究通过详细分析其重整化行为和对d轨道物理的影响，揭示了s轨道并非仅仅是费米面上的一个独立特征，它在强关联背景下其行为被深刻改变，且其存在能够增强d轨道的局域性和反铁磁关联。这种间接但关键的作用对于理解镍酸盐的整体电子性质和关联效应至关重要。
  3. 强调强关联效应的普适性： 镍酸盐与铜氧化物在结构上的相似性促使人们将其视为铜氧化物的潜在“表亲”。本研究通过DQMC模拟，再次验证了强关联效应在这些过渡金属氧化物体系中的普适性和主导作用。它表明，要准确描述这些材料的物理，必须超越单粒子近似（如DFT），显式地处理多体关联。
  4. 提供非常规超导的新线索： 研究发现，d-p-s模型相较于d-p模型显著增强了短程反铁磁自旋关联。鉴于反铁磁涨落被广泛认为是许多非常规超导体中配对胶水（pairing glue）的起源，这一发现为镍酸盐的超导机制提供了重要的启示。它暗示了间隙s轨道可能通过影响d轨道的局域性和磁性涨落，间接促进了超导配对。
  5. 推动多轨道模型的发展： 镍酸盐的复杂性要求使用多轨道模型。本研究成功扩展了经典的Emery模型，纳入了间隙s轨道并处理了其三维色散，为未来研究其他具有复杂轨道结构的量子材料提供了范例。它证明了在现有框架下，通过引入实验驱动的新轨道，能够更准确地描述材料的物理。
  6. 挑战传统DFT预测： 本研究的结果直接挑战了部分DFT计算对镍酸盐能带结构的预测，特别是对d$_{x^2-y^2}$轨道kz色散和Lifshitz相变的描述。这强调了在关联电子体系中，仅依赖DFT可能不足以捕捉所有关键物理，数值多体方法是不可或缺的。

5.2 与铜氧化物的比较和独特之处

无限层镍酸盐常被称为“铜氧化物表亲”，但本研究揭示了它们之间的关键相似点和独特差异：

相似之处：

  • 准二维d$_{x^2-y^2}$轨道物理： 与铜氧化物类似，镍酸盐中的d$_{x^2-y^2}$轨道在强关联作用下表现出显著的准二维特征，其kz色散被强烈抑制。这可能是其作为高温超导体的共同基础。
  • 反铁磁涨落： 两种材料都存在强大的反铁磁自旋涨落，这些涨落被认为是介导非常规超导配对的可能机制。本研究中d-p-s模型下的AFM关联增强，进一步强化了这一观点。
  • UHB谱重转移： 掺杂导致UHB谱权重向低能转移，这在铜氧化物中也是一个典型的强关联现象。

独特之处：

  • 间隙s轨道的存在： 这是镍酸盐独特的特征，铜氧化物中没有与之直接对应的轨道。本研究证实了s轨道在镍酸盐中扮演着不可或缺的角色，其行为在强关联下被重塑，并间接影响d轨道的物理。
  • 更大的电荷转移能量： 镍酸盐通常具有比铜氧化物更大的电荷转移能量（在Emery模型中表现为 $\varepsilon_p - \varepsilon_d$）。本研究中 $\varepsilon_p = 5.0$ eV 的选择反映了这一点，并导致了忽略 $U_{pp}$ 的合理性。这意味着镍酸盐中的空穴可能更倾向于占据d轨道而非p轨道。
  • 更复杂的轨道结构： 除了间隙s轨道，镍酸盐还可能涉及稀土元素的5d或4f轨道，这些进一步增加了其复杂性，可能导致近藤效应等独特现象，而这些在标准铜氧化物模型中不常见。
  • d轨道更强的局域性： 研究指出d轨道的准粒子峰在镍酸盐中更加尖锐，这反映了d轨道在镍酸盐中可能具有更强的局域性，导致更强的关联效应。

5.3 未来研究方向与展望

基于本研究的发现和局限性，未来的研究可以沿着以下几个方向深入探索：

  1. 克服符号问题，实现极低温模拟： 符号问题是DQMC的主要瓶颈。未来的工作可以探索新的DQMC算法变体（例如，纠缠态DQMC、复平面DQMC）或其他无符号问题的方法（如动力学平均场理论+量子蒙特卡洛DMFT+QMC、变分蒙特卡洛VMC），以在更低的温度和更广的掺杂范围内进行模拟。这将使得对超导配对机制的直接探测成为可能。
  2. 直接计算超导配对强度： 本研究仅暗示了AFM涨落对超导的促进作用。未来的研究应直接计算不同配对通道（如d波）的超导配对磁化率或序参量，以量化间隙s轨道对超导配对强度的具体影响。
  3. 探索更复杂的轨道间相互作用： 虽然本研究忽略了 $U_{pp}$，但未来的工作可以评估 $U_{pp}$ 和其他轨道间库仑相互作用（如 $U_{pd}$）在镍酸盐中的作用，以构建更全面的哈密顿量。
  4. 动态自旋涨落与配对机制： 除了静态自旋关联，动态自旋涨落对于理解自旋涨落介导的超导配对至关重要。未来的研究可以尝试计算动态自旋磁化率，尽管这在MEM解析延拓中更具挑战性。
  5. 与更广泛的实验数据对比： 将模拟结果与更多实验技术（如中子散射、共振X射线散射）获得的数据进行对比，以验证模型的普适性和预测能力。
  6. 掺杂相图的探索： 进一步探索不同掺杂浓度下间隙s轨道和d$_{x^2-y^2}$轨道的演变，构建更完整的电子和磁性相图，并与实验中的超导相图进行对照。
  7. 有限尺寸效应的系统性研究： 通过更先进的计算资源和算法，系统性地研究更大尺寸晶格下的物理行为，以更精确地捕捉三维色散和长程关联。

总而言之，本研究为无限层镍酸盐的复杂物理打开了一个新的窗口，强调了强关联和间隙s轨道的不可或缺性。未来的工作将在这些基础上继续深入，有望最终揭示这类材料的超导奥秘。