来源论文: https://arxiv.org/abs/2604.08085v1 生成时间: Apr 10, 2026 12:11

执行摘要

铜氧化物高温超导体是强关联电子物理领域最活跃的研究方向之一。理解其非传统超导性和奇异金属态的关键在于构建并求解合适的低能格点模型。本文深入探讨了两种最常用的模型：单轨道Hubbard模型和三轨道Emery模型。通过结合密度泛函理论（DFT）、瓦尼尔函数（Wannier functions）和动力学平均场理论（DMFT）并使用数值重整化群（NRG）杂质求解器，作者系统地比较了这两种模型在铜氧化物超导体La$_2$CuO$_4$基体系中的光谱和直流（dc）输运特性。研究发现，虽然两种模型在定性上呈现相似的物理图景，但在定量上存在显著差异。Hubbard模型在重现实验数据时需要比传统估计值更大的相互作用强度U，而Emery模型在常规U值下表现良好。此外，本文识别了费米能级附近局部谱重、Lifshitz相变临界掺杂和有效质量等对U值敏感的特性，这些可作为未来实验估算有效耦合常数的潜在途径。值得注意的是，作者还对Planckian耗散理论中实验报告的有效散射率计算中存在的混淆进行了批判性分析，指出回旋有效质量和输运有效质量的混用可能导致对散射率的严重误判。本研究为未来理论和实验工作提供了重要指导，以期更好地理解铜氧化物中的强关联物理。

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

核心科学问题

理解铜氧化物超导体的输运性质是强关联电子物理学中的核心挑战之一。最常见的途径是定义和求解一个低能格点模型，但目前尚不清楚究竟哪个“最小模型”才能捕获所有相关机制并提供定量预测。主要的不确定性在于选择模型中应包含的轨道自由度，以及如何准确定义有效耦合常数。这篇论文的核心目标正是通过系统比较单轨道Hubbard模型和三轨道Emery模型，来解答以下关键问题：

  1. 哪个模型能更好地描述La$_2$CuO$_4$基铜氧化物？ 这涉及到对两种模型光谱和输运特性的全面比较，以评估它们与实验数据的契合程度。
  2. 如何确定每种模型中正确的有效耦合常数U？ 由于U值对模型的物理行为至关重要，但其精确值难以确定，因此需要寻找对U敏感的物理量，并通过与实验对比来约束U的取值范围。
  3. 两种模型在物理上是否等价，或者说，Hubbard模型能否作为Emery模型在低能范围的有效单带描述？ 论文通过比较两者计算出的各种物理量，探讨它们之间的异同。

理论基础

本研究建立在以下几个关键的理论支柱之上：

  • 低能格点模型 (Low-Energy Lattice Models): 这是研究强关联电子系统的标准方法，通过简化实际材料的复杂性，聚焦于少数几个关键自由度和相互作用。论文主要关注两种模型：

    • Hubbard 模型 (Hubbard Model): 描述电子在晶格中跳跃并仅在同一格点上存在局域库仑排斥的单带模型。其哈密顿量通常包含动能项和局域相互作用项U。它能定性描述许多非常规超导体的正常态输运性质和光谱特性，例如铜氧化物中的节点-反节点二分性和赝能隙。
    • Emery 模型 (Emery Model): 一种三轨道模型，通常包含铜d轨道和氧p轨道，并考虑了它们之间的跳跃和相互作用。由于其更精细的轨道自由度，Emery模型被认为能更准确地描述铜氧化物中的电荷转移性质，这对于理解铜氧化物超导至关重要。

  • 密度泛函理论 (Density Functional Theory, DFT): 用于计算材料的电子基态性质，包括能带结构。DFT是一种第一性原理方法，为后续的格点模型参数化提供了基础。通过DFT计算得到的能带结构，可以进一步“下折叠”得到低能格点模型的紧束缚（tight-binding, TB）参数。

  • 瓦尼尔函数 (Wannier Functions): 一种局域化的轨道基组，用于将DFT计算得到的非相互作用能带结构转换成紧束缚模型的参数。通过极大局域化瓦尼尔函数（MLWF）方法，可以从DFT能带结构中提取Hubbard和Emery模型所需的跳跃积分和晶格内能。

  • 动力学平均场理论 (Dynamical Mean Field Theory, DMFT): 一种处理强关联电子系统的非微扰方法。DMFT将格点模型映射到一个有效的安德森杂质模型（Anderson Impurity Model），通过自洽地求解该杂质模型来获得局域自能（self-energy）。DMFT的核心近似是自能只依赖于频率而不依赖于动量，即 Σk(ω) ≈ Σ(ω)。

  • 数值重整化群 (Numerical Renormalization Group, NRG): DMFT中求解安德森杂质模型的常用数值方法之一。NRG的优势在于它直接在实频率轴上计算结果，避免了解析延拓这一在处理数值数据时通常遇到的不稳定且定义不清的问题。尽管NRG引入了额外的近似，但已被证明在许多工作中与精确量子蒙特卡洛杂质求解器等方法进行了基准测试，误差较小。

  • Kubo 泡法 (Kubo Bubble Formalism) 用于直流电导率计算: 描述系统线性响应的通用理论框架。在DMFT中，由于自能的局域性，电导率的顶点修正（vertex corrections）消失，从而使得Kubo泡法可以直接用于计算直流电导率，无需额外近似。这使得DMFT能够直接预测输运性质，并与实验进行比较。

  • Peierls 代换 (Peierls Substitution): 一种在格点模型中引入规范势（或电场）的方法，用于推导电流算符和速度矩阵。这种方法假设轨道是高度局域化的，这在Emery模型中可能比Hubbard模型更合理。

  • Lifshitz 相变 (Lifshitz Transition): 费米面拓扑结构发生变化的相变。在铜氧化物中，随着掺杂水平的变化，费米面可能从空穴型变为电子型，这可以通过角分辨光电子能谱（ARPES）实验观察到，并且与材料的输运性质和比热等宏观物理量密切相关。

  • 有效质量 (Effective Mass): 描述电子在晶格中有效行为的物理量。论文区分了两种重要的有效质量：

    • 回旋有效质量 (Cyclotron Mass, m*c): 出现在de Haas-van Alphen效应和Sommerfeld比热表达式中，描述费米面附近电子的状态密度和费米面的面积对能量的响应。它可以从动量分辨光谱函数中提取。
    • 输运有效质量 (Transport Mass, m*tr): 出现在Drude电导率公式中，描述电子在电场作用下的响应。它可以通过应力张量（kinetic-energy operator）的期望值计算，通常与能带结构曲率有关。

  • Planckian 耗散理论 (Planckian Dissipation Hypothesis): 一种理论假设，认为在某些条件下（例如量子临界点附近），散射率1/τ会饱和到由Planck常数定义的线性温度依赖关系 akBT/ħ，其中a≈1。这与奇异金属态中的线性电阻率密切相关，并被多项实验数据所支持。

技术难点

本研究所面临和解决的技术难点包括：

  1. 模型选择与参数化: 确定“最小模型”本身就是挑战。此外，如何从第一性原理计算（DFT）中可靠地“下折叠”出低能格点模型的参数（如跳跃积分t和相互作用U）是一个复杂问题。特别是，有效库仑相互作用U的计算存在较大不确定性，文献中U的估计范围宽泛（3-14 eV），且不同下折叠方法（如cRPA）和不同模型（单带 vs. 三带）会给出不同的U值。
  2. 双重计数问题 (Double Counting): 在DFT+DMFT方法中，由于DFT在平均场水平上已经部分包含了电子间相互作用，而DMFT又引入了强关联效应，因此需要一个精确的“双重计数”修正，以避免重复计算相互作用能。对于Emery模型，d轨道的Hartree修正Edc的精确值并不为人所知，这引入了额外的估计。
  3. DMFT自洽循环的数值稳定性与收敛性: DMFT需要迭代求解自洽方程。在计算局域格林函数时，布里渊区（BZ）上的动量积分需要高精度，因为光谱函数ImGk,dd(ω)在许多情况下会非常尖锐，这使得积分收敛困难。作者采用cubepy自适应积分库来解决这一问题。
  4. 实频率计算的挑战: 动力学关联函数（如自能和格林函数）通常在虚频率轴上计算，然后通过解析延拓转换到实频率轴。解析延拓是一个病态问题，容易引入数值噪声。NRG作为杂质求解器，直接在实频率轴上工作，从而避免了这一挑战，但NRG本身也存在近似性。
  5. 无序效应 (Disorder Effects) 的处理: 实际的铜氧化物材料通常存在显著的无序，导致残余电阻率和谱线展宽。然而，文中使用的格点模型是“干净”的，不包含无序。为了与实验进行比较，作者通过引入频率无关的虚自能Γ来 phenomenologically 模拟谱线展宽，并通过Matthiessen定律引入残余电阻率。这种处理方式的准确性仍有待商榷，因为无序对光谱和输运的影响可能不是简单的叠加。
  6. 有效质量定义的混淆: 论文发现，在现有实验文献中，回旋有效质量和输运有效质量经常被互换使用，尤其是在Planckian耗散的背景下。然而，这两种质量在物理上是不同的，特别是在Lifshitz相变附近，它们的值可能相差很大。本研究明确区分并计算了这两种质量，揭示了这种混淆可能导致对散射率的严重误判。
  7. DMFT的局限性: 单格点DMFT假设自能是局域的，这忽略了空间关联效应。尽管对于输运性质，顶点修正的消失使得其结果在一定程度上可靠，但在更广泛的强关联物理中，空间关联可能扮演重要角色。此外，单格点DMFT在有限掺杂下不产生量子临界点。

方法细节

本研究采用了一套综合性的计算方法，将第一性原理计算与强关联格点模型求解相结合：

  1. DFT能带结构计算 (DFT Band Structure Calculation):

    • 软件: 使用Quantum Espresso软件包 [145, 146] 进行DFT计算。
    • 晶体结构: 针对La$_2$CuO$_4$的两种晶体结构（T结构和T’结构）进行计算。T结构适用于空穴掺杂的LSCO和LBCO，T’结构适用于电子掺杂的LCCO。
    • 参数: 采用SG15 ONCV PBE赝势 [148]，使用原始晶胞，6x6x6 k点网格，动能截止能量80 Ry，电荷密度截止能量400 Ry，高斯展宽0.02 Ry。
    • 结果: 获得了轨道投影的DFT能带结构，显示了Cu 3d和O 2p轨道的贡献 (图2)。

  2. 下折叠到格点模型参数 (Downfolding to Lattice Model Parameters):

    • 方法: 利用极大局域化瓦尼尔函数（MLWF）方法 [103-105]，通过REBPACK [105, 151] 或Wannier90 [152] 软件将DFT能带结构下折叠到格点模型。
    • 轨道选择: 对于Emery模型，选取Cu 3d$_{x^2-y^2}$轨道和相应的O 2p轨道。对于Hubbard模型，通过单带瓦尼尔化，将Emery模型中最高能量带（α=2）的色散关系作为Hubbard模型的色散关系，以方便比较。
    • 紧束缚参数: 获得了Emery模型的オンサイト能（Ed, Ep）和跳跃积分（tdp, tpp, t’pp）。T结构和T’结构的具体参数列于表I。为了简化模型，截断了所有长程跳跃。

  3. DMFT自洽求解 (DMFT Self-Consistent Solution):

    • 框架: 使用DMFT方法求解Hubbard和Emery模型，其中将格点模型映射到一个有效的安德森杂质模型。
    • 杂质求解器: 采用数值重整化群（NRG）方法 [106-108] 求解杂质模型，以获得实频率的局域自能Σ(ω)。
    • 自洽循环: 具体步骤如下 (详见III.B节):
      1. 初始猜测自能Σ(ω)。
      2. 计算局域格林函数Gloc,dd(ω)。对于Emery模型，涉及BZ上的3x3矩阵积分 (方程7, 8)，使用cubepy自适应积分库 [156]。对于Hubbard模型，为标量积分 (方程10)。
      3. 计算新的杂化函数Δ(ω)。
      4. 通过NRG求解杂质问题，计算杂质自能Σimp。
      5. 近似格点自能Σ(ω) ≈ Σimp，返回步骤2，直到收敛。
    • 双重计数: 对于Emery模型，d轨道上的オンサイト能Ed需要进行双重计数修正Edc (取3.12 eV)。

  4. 掺杂定义与相图扫描 (Doping Definition and Phase Diagram Scans):

    • 掺杂: 通过计算每个格点每个自旋的总占据数 (n_σ) 来确定掺杂水平 δ = 2(1/2 - <n_σ>) (Hubbard) 或 δ = 2(1/2 - (n_d)) (Emery) (方程19-22)。正的 δ 对应空穴掺杂，负的 δ 对应电子掺杂。
    • 扫描: 在固定化学势和温度网格上求解模型，然后通过插值获得固定掺杂水平下的物理量，例如dc电阻率对温度的依赖关系。

  5. 物理量计算:

    • 局域光谱函数 (Local Spectral Function, A(ω)): 通过动量积分获得 (方程23)。
    • 动量分辨光谱函数 (Momentum-Resolved Spectral Function, Ak(ω)): 用于与ARPES实验进行比较。为了与实验匹配，计算结果会额外引入一个频率无关的虚自能Γ进行展宽，以模拟真实样品中的无序效应 (图10)。
    • 直流电导率 (dc Resistivity, ρdc): 使用Kubo泡法计算 (方程13, 17)。由于DMFT中自能的局域性，顶点修正消失。电流算符的推导通过Peierls代换完成 (详见III.C节和附录C)。
    • Lifshitz 相变临界掺杂 (δLT): 通过费米面拓扑结构的变化（由零频率光谱函数A_k(ω=0)的极大值线定义）来识别 (图16, 17)。
    • 有效质量 (Effective Mass):
      • 回旋有效质量 m*c: 从光谱函数中提取费米动量kF和重整化费米速度vF(θ)计算得到 (方程28，附录A)。
      • 输运有效质量 m*tr: 从Drude公式和应力张量 (K) 计算得到 (方程A11，附录A)。
    • Planckian 耗散因子 (α): 基于散射率与Planckian极限的比值计算 (方程29)。

通过上述方法，研究人员能够对Hubbard和Emery模型进行详细的定量分析，并与各种实验数据进行直接比较，从而解决铜氧化物超导中的核心科学问题。

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

本研究的核心在于通过对La$_2$CuO$_4$基铜氧化物这一关键体系的深入研究，系统地比较了Hubbard和Emery模型的预测。La$_2$CuO$_4$是铜氧化物超导体的母体化合物，通过掺杂（Sr、Ba或Ce）可以产生空穴掺杂或电子掺杂的超导电性。论文选择了以下三种代表性化合物作为benchmark体系：

  1. Sr掺杂La$_2$CuO$_4$ (LSCO): 典型的空穴掺杂体系，呈现T结构晶格。广泛用于研究高温超导体的赝能隙、超导态和奇异金属态。
  2. Ba掺杂La$_2$CuO$_4$ (LBCO): 也是空穴掺杂体系，同样具有T结构晶格。在某些掺杂浓度下展现出复杂的结构相变和电荷序，但在此研究中主要关注其与LSCO相似的输运性质。
  3. Ce掺杂La$_2$CuO$_4$ (LCCO): 典型的电子掺杂体系，呈现T’结构晶格。与空穴掺杂体系相比，电子掺杂铜氧化物的物理性质存在显著差异，为模型提供了不同的测试场景。

计算所得数据与关键发现

论文生成了大量的计算数据，并与多项实验测量结果进行了详细对比。以下是主要发现和数据点：

1. 能带结构与参数化

  • DFT和Wannier化能带结构 (图2): 左侧展示了Cu 3d$_{x^2-y^2}$、O 2p$_{x}$和O 2p$_{y}$轨道投影的T和T’结构DFT能带。右侧显示了经过Wannier化后的1带和3带紧束缚能带结构，它们与DFT结果吻合良好。这证明了下折叠方法的有效性。
  • 紧束缚参数 (表I): T和T’结构Emery模型的关键紧束缚参数（オンサイト能差 ε_d - ε_p、跳跃积分 t_pd、t_pp 和 t’_pp）被精确列出。这些参数是模型量化分析的基础。
  • Hubbard模型裸态态密度 (图4): 展示了两种Hubbard模型（对应T和T’结构）的裸态态密度，它们在DMFT求解中作为输入。

2. 局域光谱函数 (A(ω)) 与能隙

  • Hubbard vs Emery 光谱比较 (图6): 对于T结构（LSCO/LBCO），在不同掺杂和温度下，Hubbard (U=5eV) 和Emery (U=8eV) 模型的光谱函数被比较。两者都显示出较低的Hubbard带 (LHB) 或Zhang-Rice单线态带 (ZRSB) 和较高的Hubbard带 (UHB)。
    • 能隙 (Δgap) 依赖性 (图7a): Hubbard模型中的Mott能隙随掺杂变化比Emery模型中的电荷转移能隙更显著。在低掺杂时Mott能隙小于电荷转移能隙，高掺杂时则相反。最佳能隙匹配发生在Hubbard U ≈ 5eV。
    • 能带权重 (图7b): Hubbard模型中LHB和UHB的权重之和为1。Emery模型中ZRSB和UHB的权重之和小于1。Hubbard模型的UHB权重小于Emery模型的UHB，而LHB权重高于ZRSB。
  • 占据态谱重 (Pwindow) 依赖性:
    • 掺杂依赖性 (图7c): Pwindow（费米能级以下-0.12eV到0eV的谱重积分）随掺杂增加而增加。在UHubbard = 8eV和UEmery = 8eV时，两种模型的掺杂依赖性更加相似，表明U值在模型间的可比性。
    • U依赖性 (图8, 12): Hubbard模型的Pwindow(δ)在特定掺杂δmax处达到最大值。δmax随U的增加而从约20%线性增加到40%（对于T结构，图8），对于T’结构（LCCO），最大值发生在更低的掺杂水平（图12）。这一“特征掺杂”δmax对U值高度敏感，可能为未来光电子能谱实验提供估算U的方法。
    • 准粒子峰宽 (图9): 费米能级附近准粒子峰的宽度强烈依赖于U。随着U从4eV增加到16eV，最大值和最小值之间的距离从约0.45eV减小到0.15eV，表明准粒子寿命增加。

3. 动量分辨光谱函数与ARPES对比

  • 与ARPES对比 (图10): Hubbard (U=10eV) 和Emery (U=8eV) 模型在节点处计算得到的动量分辨光谱函数（在费米能级）与LSCO的ARPES实验数据进行了比较。
    • 展宽效应: 理论计算预测的准粒子峰比实验观察到的更锐利。通过额外引入约0.5eV的虚自能Γ进行展宽，可以与实验数据获得良好吻合。这暗示真实样品中存在显著的无序效应。
    • U值适配: Hubbard模型在UHubbard = 10eV时与ARPES实验结果最佳匹配。Emery模型在UEmery = 8eV时表现良好，但仍需展宽。

4. 直流电阻率 (ρdc)

  • Hubbard vs Emery 电阻率 (图13, 14):
    • 定性相似性: 两种模型的电阻率随掺杂和温度变化的定性行为非常相似，几乎难以区分。
    • 定量差异: Emery模型通常表现出更高的电阻率。当Hubbard模型的UHubbard ≈ 1.2 UEmery 时，两种模型的电阻率才能很好地匹配 (图14)。然而，此时光谱函数上的能隙匹配不再最佳。
  • 与实验对比 (图15):
    • LSCO/LBCO (T结构): Emery模型 (UEmery = 8eV) 的结果与LSCO和LBCO的实验数据非常吻合。Hubbard模型需要UHubbard ≈ 10eV才能获得 excellent agreement。
    • LCCO (T’结构): Hubbard模型 (UHubbard ≈ 8eV) 在低于150K时与LCCO实验数据吻合良好，但在更高温度下实验曲线表现出更强的弯曲，理论模型未能捕获。Emery模型 (UEmery = 8eV) 则高估了LCCO的电阻率。
    • 残余电阻率: 为了与实验对比，理论曲线通过Matthiessen定律向上平移了一个与掺杂相关的常数项 ρ0，以考虑真实样品中的无序导致的残余电阻率。

5. Lifshitz 相变

  • 费米面拓扑变化 (图16): 在Emery模型中，费米面拓扑随掺杂从空穴型变为电子型，清晰地展示了Lifshitz相变。例如，在LSCO中，δLT ≈ 21%。
  • δLT 的U依赖性 (图17): 临界掺杂δLT随U的增加而向下移动。当U值与电阻率实验数据吻合良好时，δLT的预测也与LSCO实验值（17%-22%）一致。这表明模型中空穴数量（δ）与实验中的掺杂浓度（x）一致，支持了格点模型方法的合理性。

6. 有效质量与Planckian耗散

  • 回旋有效质量 m*c (图18): mc强烈依赖于U和掺杂。Hubbard模型在UHubbard = 6-12eV时与LSCO实验数据（Legros等人的曲线）吻合最好。Emery模型 (UEmery = 8eV) 也与实验吻合良好，且对温度依赖性较弱。然而，对于LCCO，Emery模型显著高估了mc，而Hubbard模型需要在UHubbard = 4eV时才能匹配，这与电阻率分析结果不一致。
  • 输运有效质量 m*tr (图21): mtr比mc受到的重整化程度小。在δ = 26%时，mc可达20m_e，而mtr仅为2-3m_e，两者相差约5倍。这突出了两种有效质量的本质区别，尤其是在Lifshitz相变附近。
  • Planckian 耗散因子 α (图19): 如果按照Legros等人的方法（使用mc而非mtr）计算，Hubbard模型在UHubbard = 10eV时可以再现实验报告的α ≈ 1。但如果使用正确的输运有效质量m*tr，计算出的α值将约为5，远高于Planckian极限。

7. 轨道分辨与本征带分辨光谱函数 (附录B, 图22)

  • Emery模型中的轨道贡献: ZRSB能带具有铜d轨道和氧p轨道大致相等的混合特征，而UHB主要具有d轨道特征。这表明Hubbard模型中的“单一轨道”并非纯粹的铜d轨道，而是d和p轨道的混合。
  • 本征带分辨光谱: 费米能级附近的能带几乎是纯α=2特征（对应最高能量本征带），这暗示一个基于Eα=2,k色散的单带模型可能足以描述Emery模型在低能区的物理。

性能数据

论文没有详细给出具体的计算性能数据，但提到了几个关键点：

  • 计算复杂性: DMFT方法本身是计算密集型的，特别是布里渊区积分需要自适应算法来处理尖锐的谱函数，例如cubepy库 [156]。
  • NRG的优势: NRG求解器在实频率轴上工作，避免了解析延拓的计算难题，这在处理动态响应函数时尤其重要。
  • 并行化: 未明确提及大规模并行计算的细节，但DMFT自洽循环的性质使其在一定程度上可并行化（例如BZ积分）。
  • 数据量: 作者生成了“extensive dataset”，这表明进行了大量的计算以扫描掺杂-温度相图以及U值。

总的来说，这些计算结果和性能相关信息为我们提供了对Hubbard和Emery模型在铜氧化物中行为的深刻理解，并强调了在模型选择和参数化过程中需要注意的关键定量差异。

本研究的计算流程融合了第一性原理计算与强关联多体理论，其核心实现依赖于一系列成熟的软件包以及作者团队可能开发的专用代码。以下将详细介绍涉及的软件、关键算法细节以及复现该研究的指南。

所用的软件包

  1. 密度泛函理论 (DFT) 计算:

    • Quantum Espresso [145, 146]: 一个广泛使用的开源DFT软件包，用于计算La$_2$CuO$_4$T和T’结构的电子能带结构。作者使用其pw.x模块进行平面波基组的DFT计算。
    • 赝势库: 采用了SG15 ONCV PBE [148] 赝势。

  2. 瓦尼尔函数 (Wannier Functions) 生成:

    • REBPACK [105, 151] 和 Wannier90 [152]: 这两个工具通常结合DFT软件包（如Quantum Espresso）使用，用于从DFT能带结构中提取极大局域化瓦尼尔函数（MLWF）和相应的紧束缚（TB）参数。REBPACK是一个专门用于生成瓦尼尔函数和下折叠能带的Python库，而Wannier90是生成MLWF的Fortran程序。论文指出，他们检查了Wannier90的结果与REBPACK一致。

  3. DMFT 杂质求解器:

    • 数值重整化群 (Numerical Renormalization Group, NRG) [106-108]: 这是一个由作者团队（如R. Žitko）长期开发和维护的NRG求解器实现。NRG能够直接在实频率轴上高精度地求解安德森杂质模型，为DMFT循环提供杂质自能。

  4. 数值积分库:

    • cubepy [156]: 一个Python的自适应积分库，特别适用于处理布里渊区（BZ）上的高振荡或尖锐函数积分。该库实现了Genz-Malik自适应正交算法 [157] 和Berntsen的误差估计 [158]，对于DMFT中计算格林函数所需的BZ积分至关重要。

代码实现细节

虽然论文没有提供具体的代码库链接，但详细阐述了核心算法和公式，为复现提供了坚实基础：

  1. 格点模型哈密顿量的构建:

    • Emery 模型 (方程1-3): 哈密顿量在轨道基组 l = (d, p_x, p_y) 下表示为一个3x3矩阵。其紧束缚参数 (ε_d, ε_p, t_pd, t_pp, t’_pp) 来自Wannier化结果 (表I)。相互作用仅限于d轨道的局域库仑相互作用U。
    • Hubbard 模型 (方程5): 使用Emery模型中最高能量带 (α=2) 的色散关系作为其非相互作用动能项。相互作用仅为单轨道上的局域U。

  2. DMFT 自洽循环:

    • 格林函数计算:
      • Emery 模型 (方程7-8): 局域d轨道格林函数G_loc,dd(ω)通过对BZ上的动量k积分Gk,dd(ω)得到。其中，Gk(ω)是格点哈密顿量h_k和自能Σ(ω)的逆矩阵运算。自能Σ(ω)在轨道空间中为对角矩阵，仅d轨道有非零分量 (方程9)。
      • Hubbard 模型 (方程10): 局域格林函数是标量形式，通过对裸态态密度ρ₀(ε)的频率积分得到。
    • 杂化函数更新 (Δ(ω)): 根据G_loc,dd(ω)和上一步自能Σ_dd(ω)计算得到 (详见III.B节步骤3)。
    • 杂质模型求解: 将Δ(ω)作为输入，通过NRG求解器获得新的杂质自能Σ_imp。
    • 自能更新: 格点自能Σ(ω)近似为杂质自能Σ_imp。

  3. 直流电导率计算 (Kubo Bubble):

    • 公式 (方程13, 17): 基于Kubo公式的dmft实现，涉及到对格林函数和速度算符的积分。由于DMFT的局域自能近似，顶点修正为零，简化了计算。
    • 速度算符 (方程14, 15，附录C): 通过对引入规范势的哈密顿量求导，根据Peierls代换得到。Emery模型中的速度矩阵推导更为复杂，在附录C中详述。

  4. 掺杂与谱函数:

    • 掺杂定义 (方程19-22): 严格定义了Hubbard和Emery模型中掺杂 δ 的计算方法。
    • 谱函数 (方程23-27，附录B): 局域A(ω)、动量分辨A_k(ω)以及轨道分辨和本征带分辨A_l(ω), A_α(ω) 的具体计算公式。

  5. 有效质量与 Planckian 耗散:

    • 回旋有效质量 (方程A1-A3): 通过费米面面积对费米能级的导数或费米速度计算。
    • 输运有效质量 (方程A6-A12): 基于Drude电导率公式和应力张量（K）计算。应力张量K可以通过动量分辨谱函数得到 (方程A12)。
    • Planckian 耗散因子 (方程29): 将计算出的输运有效质量m*tr和电阻率斜率A_ρ 代入公式。

复现指南

要复现本研究，研究人员需要具备DFT、Wannier化、DMFT和NRG的背景知识，并可能需要访问高性能计算资源。以下是大致的复现步骤：

  1. DFT能带结构计算:

    • 下载并安装Quantum Espresso软件包。
    • 获取La$_2$CuO$_4$T和T’结构的晶体结构文件（如CIF格式）。可以从Materials Project数据库 [143, 144] 获取，论文提供了ID：mp-19735 (T-structure) 和 mp-1077929 (T’-structure)。
    • 根据论文III.A节的参数（赝势、k点、截止能量、展宽等），运行Quantum Espresso进行SCF和NSCF计算，生成用于Wannier化的文件。

  2. Wannier化生成紧束缚参数:

    • 安装REBPACK或Wannier90。
    • 使用DFT计算结果作为输入，按照论文第二节和第三节的描述，生成Cu 3d$_{x^2-y^2}$和O 2p轨道（Emery模型）或单带（Hubbard模型）的MLWF。
    • 从MLWF计算中提取紧束缚哈密顿量的参数（オンサイト能和跳跃积分）。论文中的表I提供了关键参数，可以直接用于模型构建。

  3. DMFT 求解器准备:

    • Hubbard 模型: 构建一个Hubbard模型求解器。由于Hubbard模型是单带的，其DMFT自洽循环相对简单。
    • Emery 模型: 构建一个三带DMFT求解器，处理3x3矩阵运算。需要实现论文中方程7和8的BZ积分，以及方程9的自能矩阵形式。由于作者使用的NRG求解器可能是内部开发的，复现者可能需要使用自己的NRG实现或选择其他开源DMFT杂质求解器（如TRIQS库中包含的）。
    • 双重计数: 确保在Emery模型中正确实现d轨道的双重计数修正Edc。

  4. DMFT 自洽计算:

    • 根据论文III.D节，定义掺杂和化学势网格。
    • 在不同U值、掺杂和温度下，运行DMFT自洽循环。确保收敛准则（例如，格林函数或自能的变化）得到满足。

  5. 物理量计算与分析:

    • 光谱函数: 从收敛的DMFT结果中计算局域和动量分辨光谱函数 (A(ω), Ak(ω))。
    • 直流电阻率: 实施Kubo泡法计算，包括速度算符和格林函数积分。
    • Lifshitz 相变: 分析A_k(ω=0)以识别费米面拓扑变化。
    • 有效质量: 实施回旋有效质量和输运有效质量的计算，注意区分两者。这需要对光谱函数和应力张量进行后处理。
    • 数据可视化: 使用matplotlib或其他绘图库重现论文中的图表，以便与论文结果进行直接比较。

开源 Repo 链接

虽然核心DMFT+NRG代码似乎是作者团队内部实现，但上述开源工具足以构建一个类似的计算框架，并有望复现论文的关键发现。在复现过程中，对数值稳定性、积分精度和DMFT收敛性的细致处理至关重要。

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

关键引用文献

本研究建立在大量先前工作的基础上，并引用了诸多领域内的经典和最新文献，以下是其中一些特别值得关注的类别：

  1. 非常规超导体与强关联电子物理学综述 [1-3]: 这些综述奠定了研究铜氧化物超导体的宏观物理背景，强调了低能格点模型的重要性。
  2. Hubbard 模型在铜氧化物和相关体系中的应用 [4-18, 25-39]: 大量工作使用Hubbard模型解释铜氧化物的正常态输运、光谱特性、赝能隙和超导态，构成了本研究中Hubbard模型分析的基石。特别是 [35, 109, 153] 详细介绍了DMFT方法及其在Hubbard模型上的应用。
  3. Emery 模型在铜氧化物中的应用 [19-24, 40-53]: 这些文献突出了Emery模型在更精细描述铜氧化物电荷转移和轨道自由度方面的优势，解释了其与Hubbard模型的区别和互补性。
  4. DFT+DMFT 和下折叠方法 [56-60, 66, 67, 102]: 这些工作提供了从第一性原理（DFT）计算中获取格点模型参数的方法，并处理了DFT与DMFT结合中的双重计数问题，为本研究的模型构建提供了理论和技术支持。
  5. cRPA (constrained Random Phase Approximation) 估算 U 值 [27, 61, 63-81, 83]: cRPA是计算有效库仑相互作用U的常用方法，论文讨论了其在不同材料中U值范围宽泛的现象，并指出了直接将cRPA的频率依赖耦合函数转换为瞬时U值的挑战。
  6. Wannier 函数方法 [103-105, 151, 152]: 这些文献详细介绍了极大局域化瓦尼尔函数（MLWF）的理论和应用，它是从DFT能带结构提取紧束缚参数的关键工具。
  7. 数值重整化群 (NRG) [106-108, 154, 155]: 提供了NRG方法的基础理论和在DMFT中作为杂质求解器的应用细节，强调了其在实频率计算方面的优势。
  8. Kubo 公式与电导率计算 [10, 109, 160, 161, 187]: 这些文献奠定了计算电导率的理论基础，特别是在DMFT框架下，局域自能导致顶点修正消失，简化了Kubo公式的应用。
  9. 实验数据对比 [119, 122-124, 127-129, 137]: 论文直接引用了LSCO、LBCO和LCCO的dc电阻率、ARPES光谱、Lifshitz相变和有效质量的实验测量数据，作为理论模型验证的基准。
  10. Planckian 耗散理论和奇异金属 [126, 130-138, 169-172]: 这些文献讨论了奇异金属态中线性电阻率和Planckian耗散的现象，论文对此进行了深入分析和批判性评估，特别是对有效质量定义混淆的纠正。
  11. 无序效应 [119-121, 161-163]: 讨论了真实材料中无序对输运和光谱特性的影响，为论文中 phenomenologically 引入残余电阻率和谱线展宽提供了背景。

对这项工作局限性的评论

尽管本研究对Hubbard和Emery模型进行了全面的比较和深入分析，并取得了显著成果，但仍存在一些局限性，值得在未来的工作中加以改进或考虑：

  1. 单格点DMFT的局限性:

    • 忽略空间关联: 单格点DMFT的核心近似是自能只与频率有关，而与动量无关（Σk(ω) ≈ Σ(ω)）。这意味着它忽略了格点间的空间关联效应。尽管论文指出在局域自能近似下，Kubo公式的顶点修正消失，使得输运计算相对可靠，但对于其他物理量（如非局域序参数、量子临界行为），空间关联可能至关重要。未来可考虑使用簇DMFT（cluster DMFT）或非局域DMFT来捕获短程空间关联。
    • 不产生量子临界点: 单格点DMFT在有限掺杂下通常不产生量子临界点，这可能限制了其对奇异金属态中更深层次物理机制的描述能力。

  2. 无序效应的现象学处理:

    • Matthiessen定律与谱线展宽: 论文通过Matthiessen定律（引入常数残余电阻率ρ₀）和频率无关的虚自能Γ（模拟谱线展宽）来 phenomenologically 地处理无序效应。然而，无序对材料物理性质的影响可能更为复杂，例如，不同类型的无序可能对光谱线宽和输运散射率产生不同甚至不相称的影响。此外，无序强度可能随掺杂和温度变化，简单的常数修正可能不够精确。更先进的无序处理方法，如相平均DMFT，可能需要引入。

  3. 有效耦合常数U的确定:

    • U值的模型依赖性和物理量依赖性: 论文发现，Hubbard模型在不同U值下可以很好地匹配不同的实验数据（例如，电阻率需要U=10eV，ARPES需要U=10eV，回旋有效质量需要U=6-12eV），而Emery模型在常规U=8eV下表现良好。这表明没有一个单一的U值或单一模型能完全描述所有物理现象。这种“最佳U值”的依赖性使得模型的选择和参数化仍然具有一定程度的模糊性。
    • U值过高的问题: Hubbard模型在重现实验数据时需要的U值（例如10eV）远高于cRPA等方法通常预测的4-5eV。这可能暗示Hubbard模型本身作为“最小模型”的局限性，或者下折叠过程、掺杂定义等方面存在未被完全理解的因素。

  4. 下折叠过程的近似性:

    • 长程跳跃截断: 为了简化模型，论文截断了长程跳跃积分。这会引入对能带结构的轻微修改。更精确的下折叠可能需要保留更多长程项，甚至包括非局域相互作用或频率依赖的有效跳跃。
    • 双重计数修正: Emery模型中d轨道的Hartree修正Edc的值（3.12 eV）是根据文献估计的，这可能不是材料特定的精确值。Edc的选择会影响模型的物理特性。

  5. Peierls 代换的有效性:

    • 局域化假设: Peierls代换依赖于轨道高度局域化的假设，这在Emery模型（轨道更局域）中可能比Hubbard模型（能带色散更广）更合理。未来工作可以探索不依赖Peierls代换的电流算符推导方法。

  6. 温度范围的限制:

    • DMFT，特别是与NRG结合时，在极低温度下可能存在局限性，尽管NRG在实频率计算方面有优势。论文中的一些估计（如Lifshitz相变和有效质量）是在中间温度下进行的，认为这些温度下的估计“可能最有意义”，暗示了极低温度下结果的复杂性或不确定性。

  7. 无法给出“最佳模型”的明确结论:

    • 尽管进行了全面比较，论文最终仍未能给出一个明确的结论，即哪个模型（Hubbard或Emery）在所有情况下都“最优”。这反映了强关联物理的复杂性，以及“最小模型”概念本身的语境依赖性。

总而言之，本研究清晰地揭示了现有模型和方法在描述铜氧化物方面的优点和局限性。它为未来的理论和实验工作指明了方向，鼓励更精细的模型构建、更严格的参数化以及更全面的实验验证。

其他你认为必要的补充

铜氧化物超导研究的重大意义

铜氧化物高温超导体自发现以来一直是凝聚态物理领域的核心研究课题。它们在高于液氮温度下展现超导性，这一特性使其在科学和技术上都极具吸引力。然而，理解其超导机制，特别是其前超导（pseudogap）和奇异金属（strange metal）正常态，仍然是物理学中的一个重大挑战。本研究通过对Hubbard和Emery这两种核心格点模型的深入比较，直接触及了这一领域最根本的问题：如何用最简洁而又准确的模型来描述这些复杂材料的强关联物理？

论文强调，这些模型的成功不仅在于定性地重现实验现象，更在于提供定量预测，这对于揭示超导的微观机制至关重要。通过对光谱、输运、Lifshitz相变和有效质量的综合分析，本工作为我们提供了一个统一的框架，用以评估理论模型与真实材料之间的差距。

对Planckian耗散理论的批判性审视

本研究最引人注目且具有深远影响的发现之一，是对Planckian耗散理论在实验解释中存在的“有效质量混淆”的批判性分析。Planckian耗散假设在奇异金属态中，散射率1/τ与温度T线性相关，且系数a≈1。这一假设被认为是许多非常规超导体的“普适性”特征。然而，论文通过详细的理论推导和计算，明确指出：

  1. 回旋有效质量 (mc) 与输运有效质量 (mtr) 的本质区别: 理论上，mc与费米面拓扑和状态密度有关，而mtr与Drude电导率和动量驰豫有关。两者在非相互作用或弱相互作用系统中可能近似相等，但在强关联系统和Lifshitz相变附近，它们可以相差5倍甚至更多（例如，在δ=26%掺杂下，mc可达20m_e，而mtr仅为2-3m_e）。
  2. 实验解释中的混淆: 论文指出，在一些Planckian耗散的实验研究中，计算散射率时错误地将来自比热（通常与mc相关）的有效质量用于Drude公式（需要mtr）。当作者用正确的m*tr计算时，发现Planckian耗散因子a在LSCO的26%掺杂下约为5，远高于普遍认为的a≈1的“Planckian极限”。

这一发现对Planckian耗散理论的普适性提出了质疑，并可能引发对先前许多实验结果解释的重新审视。它强调了在理论和实验中严格区分物理量定义的重要性，尤其是在处理强关联系统时。对于量子化学研究者来说，这提醒我们，将从一个测量（如光谱、比热）中提取的“有效”参数应用于另一个完全不同的物理过程（如输运）时，必须谨慎核实其物理含义和适用范围。

对未来研究的启示

本研究为未来的理论和实验工作提供了多方面的指导：

实验层面：

  • U值敏感特征的直接探测: 论文识别了费米能级附近的占据态谱重Pwindow的最大值位置（δmax）和准粒子峰宽对U值的强烈依赖性。未来的光电子能谱（ARPES/XPS）实验可以尝试在不同掺杂和温度下精确测量这些特征，以直接约束Hubbard和Emery模型中的有效U值。
  • 精确区分有效质量: 实验物理学家应采用更严谨的方法来测量和区分回旋有效质量和输运有效质量。例如，通过量子振荡（de Haas-van Alphen）提取mc，同时通过光学电导率或霍尔效应测量提取mtr。这将有助于更准确地验证理论模型，并重新评估Planckian耗散等现象。
  • 量化无序效应: 进一步的实验应致力于量化真实样品中的无序强度，并研究其如何影响光谱线宽和输运散射率，以指导理论模型中更精确地纳入无序效应。

理论层面：

  • 改进下折叠方法: 探索更精确的下折叠技术，包括长程跳跃、非局域相互作用或频率依赖的有效耦合，以弥合DFT与格点模型之间的差距。这有助于更准确地定义“最小模型”及其参数。
  • 超越单格点DMFT: 考虑使用簇DMFT（cluster DMFT）、动态簇近似（DCA）或其他包含空间关联的DMFT变体，以捕获铜氧化物中可能存在的短程空间关联和非局域序参数，从而更全面地理解其相图和临界行为。
  • 精确处理无序: 开发和实现更先进的无序处理方法，例如结合相平均或淬火无序平均的DMFT方法，以更真实地模拟样品中的物理环境。
  • 重新评估U值的物理来源: 深入研究Hubbard模型中U值过高的原因。这可能涉及到对Peierls代换的改进、对有效掺杂定义的修正，或者更复杂的相互作用（如非局域或Vinet相互作用）。

对量子化学科研人员的实际启示

对于从事量子化学研究的人员而言，本研究提供了几点重要的实践启示：

  • 模型构建与参数化的精细化: 量子化学在处理分子体系时，也常涉及选择活性空间和定义有效哈密顿量（例如，CASSCF、有效哈密顿量方法）。本研究强调了从第一性原理（DFT）数据“下折叠”到格点模型时，对轨道选择、跳跃积分和有效相互作用参数（如U）的精确定义和验证至关重要。这与量子化学中构建高精度有效哈密顿量的挑战异曲同工。
  • 多体方法与凝聚态物理的交叉: 本文展示了DMFT+NRG作为一种处理强关联凝聚态系统的强大工具。量子化学研究者如果对扩展其研究范围到周期性材料，特别是强关联材料感兴趣，DMFT及其相关方法（如簇DMFT）是一个值得探索的方向。理解这些方法的原理和局限性，有助于在设计和分析材料时做出更明智的选择。
  • 批判性思维与实验-理论对比: 本研究对实验数据解释中的有效质量混淆进行了有力批判。这提醒量子化学家，在将理论计算结果与实验数据进行比较时，必须仔细核查实验测量所基于的理论框架和参数定义。即使是看起来相似的物理量，也可能因为定义方式的细微差别而导致截然不同的结论。
  • 寻找“指纹”特征: 论文识别了费米能级附近的谱重最大值位置（δmax）等对U敏感的“指纹”特征，这为通过实验数据反向推断模型参数提供了新思路。量子化学研究也可以尝试寻找分子或材料中对关键相互作用参数敏感的独特光谱或其他物理特征，以增强理论预测的实验可验证性。

总之，本研究不仅深化了我们对铜氧化物超导体的理解，更为强关联物理研究的方法论和实践提供了宝贵的经验和教训。它鼓励研究人员以更严谨、多维度的方式进行理论与实验的结合，以期最终揭示高温超导这一世纪之谜。