来源论文: https://arxiv.org/abs/2605.15059v1 生成时间: May 15, 2026 18:13

0. 执行摘要

高温超导铜氧化物中，正常态（即非超导态）的物理性质一直是一个核心未解之谜。特别是在空穴掺杂的欠掺杂区域，一种被称为“赝能隙金属”（Pseudogap metal, PG metal）的奇异态普遍存在，其与过掺杂区域的正常Fermi液体（Normal Fermi Liquid, FL）有显著差异。一个长期存在的关键科学问题是：这两种金属态是否由一个量子临界点（Quantum Critical Point, QCP）分隔？本研究通过采用动态簇近似（Dynamical Cluster Approximation, DCA）结合数值重整化群（Numerical Renormalization Group, NRG）作为杂质求解器，对二维Hubbard模型（参数设定为U=7t, t’=-0.3t）的零温正常态相图进行了深入探索。研究明确指出，在临界掺杂 p*=0.14 处，存在一个连续的量子相变。这一转变表现为费米液体能标（从电荷、自旋和d波配对关联函数中提取）的连续塌缩，并且在中间能尺度上呈现出非费米液体（Non-Fermi Liquid, NFL）行为。在相变点p*附近，反节区（antinodal region）的相干低能谱权重急剧消失并被狭窄的金属赝能隙取代，而节区（nodal region）则平滑演化并保持较强的相干性，从而导致了赝能隙金属中的费米弧（Fermi arcs）现象。这项工作不仅为高温超导铜氧化物中的正常态物理提供了重要的理论支持，也展示了DCA+NRG方法在解决强关联系统低能动力学方面的强大能力。

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

核心科学问题：高温超导铜氧化物中的正常态之谜

在过去的三十年里，对高温超导（high-Tc superconductivity）铜氧化物中掺杂莫特绝缘体（Mott insulator）的问题进行了深入研究。其中一个核心未解之谜是空穴掺杂铜氧化物中出现的“正常态”物理。在过掺杂区域，输运、热力学和光谱学测量结果大体与关联Fermi液体（Correlated Normal Fermi Liquid, FL）行为一致。然而，在欠掺杂区域，一种被称为“赝能隙金属”（Pseudogap metal, PG metal）的奇异态出现，其特征是电子谱权重主要在反节区消失，并且费米面重构与Luttinger求和规则不符。一个关键问题是，这种赝能隙金属和过掺杂的费米液体是否由一个量子临界点（Quantum Critical Point, QCP）在某个临界掺杂p*处分隔开[3-6]。

实验证据表明，在多种铜氧化物家族中，赝能隙末端（pseudogap endpoint）的量子临界性正在累积。例如，霍尔数、角分辨磁阻和量子振荡都指向费米面在p处的转变；热力学测量发现比热系数在相同区域显著增强[14-19]。此外，相图中的奇异金属区域（strange-metal region）表现出与量子临界性相关的几个方面，例如 ω/T 标度和非相干准粒子（incoherent quasiparticles）导致的反常金属行为，这些都呈风扇状（fan-like shape）从p处辐射开来。一些实验数据甚至暗示存在一个扩展的量子临界区域[12, 20]，这被认为可能是由无序引起的[21]。

描述这种PG-QCP是一个重大挑战。首先，赝能隙金属本身仍然神秘且未被完全理解；其次，即使存在，这种QCP似乎也不是可由Hertz-Millis-Moriya理论描述的传统有序参数驱动的QCP[22-24]。因此，描述PG金属和PG-QCP的几种方法已被探索，例如具有拓扑序的分数化费米液体（fractionalized Fermi liquid）理论[25-34]；波动传统序（fluctuating conventional order）[35-46]；或动态平均场理论（DMFT）的非局部扩展[47-50]。

理论基础：Hubbard模型与DCA框架

本研究采用二维方格Hubbard模型来描述铜氧化物的物理。Hubbard模型是一个描述强关联电子系统的基本模型，其哈密顿量为：

$$H = \sum_{rr'\sigma} -t_{rr'} c_{r\sigma}^\dagger c_{r'\sigma} - \mu \sum_{r\sigma} n_{r\sigma} + U \sum_r n_{r\uparrow} n_{r\downarrow}$$

或在动量空间表示为：

$$H = \sum_{k\sigma} (\epsilon_k - \mu) n_{k\sigma} + U \sum_r n_{r\uparrow} n_{r\downarrow}$$

其中，$c_{r\sigma}^\dagger$ 和 $c_{r\sigma}$ 是在 $r$ 位置产生和湮灭自旋 $\sigma$ 电子的算符，$n_{r\sigma} = c_{r\sigma}^\dagger c_{r\sigma}$ 是数算符。$t_{rr'}$ 表示跳跃矩阵元，包括最近邻（$t$）和次最近邻（$t'$）跳跃。$\mu$ 是化学势，$U$ 是在位库仑相互作用。本文将 $t=1$ 作为能量单位，并设定 $t'=-0.3t$ 和 $U=7t$，这些参数适用于描述铜氧化物物理。通过调节化学势 $\mu$ 来调整空穴掺杂 $p$（填充 $n=1-p$）。

为了解决Hubbard模型，本研究使用了动态簇近似（DCA）[48-50, 63, 96]。DCA是一种流行的簇DMFT方法，它通过将完整晶格问题映射到一个自洽求解的有限尺寸量子杂质簇来处理非局部关联。在DCA中，自能 $\Sigma_k(\omega)$ 被近似为在布里渊区（Brillouin zone, BZ）的几个“补丁”（patch）内是常数，从而将动量空间分解为有限数量的代表性区域。这种方法有效地捕捉了短程空间关联。

技术难点：低能分辨率与零温动力学

描述量子临界性需要极高的低能和低温分辨率，这是大多数计算方法难以企及的。传统的DMFT方法或基于量子蒙特卡洛（QMC）的簇方法通常受限于有限温度和无法直接获得实频动力学。具体来说：

  1. 低能尺度： QCP的标志是特征能标趋近于零。这要求求解器能够解析指数级低的能量尺度，才能精确捕捉这些行为。
  2. 零温动力学： 量子相变发生在零温。虽然可以通过在有限低温下外推来近似零温，但这可能引入误差，并且难以分辨真正的零温效应。
  3. 实频谱函数： 实验上观察到的许多现象，例如费米弧和能隙，都体现在实频谱函数中。直接计算实频谱函数对于理解物理至关重要，但通常比 Matsubara 频率上的计算更具挑战性。

方法细节：DCA+NRG的结合

为了克服上述挑战，本研究将DCA与数值重整化群（Numerical Renormalization Group, NRG）作为杂质求解器相结合。NRG是一种强大的技术，特别擅长：

  • 高低能分辨率： NRG能够解析从高能到指数级低的能量尺度上的动力学，这对于研究QCP至关重要。
  • 零温物理： NRG最初设计用于解决Kondo问题，能够直接在零温下获得精确结果。
  • 实频动力学： NRG可以直接计算实频谱函数，避免了Matsubara频率到实频的解析延拓问题。

然而，NRG的局限性在于其计算成本随希尔伯特空间大小呈指数增长，限制了它能处理的杂质模型复杂性（例如，单个杂质或小簇）。

本研究的具体实现细节如下：

  1. 四补丁DCA方案： 使用“星形补丁几何”（star patching geometry）[63]，将布里渊区划分为四个补丁：K1=(0,0)（包含节区），K2=(0,π)，K3=(π,0)（反节区），K4=(π,π)。这种划分允许捕获费米面的动量依赖性。
  2. MuNRG求解器： 采用MuNRG实现[80, 81]来求解DCA映射得到的量子杂质簇。MuNRG基于QSpace张量库[87-91]，并利用U(1)电荷和SU(2)自旋对称性来降低计算成本。
  3. 参数设置： 离散化参数设定为 $\Lambda=8$，保留多重态数量 $N_{keep}=3 \times 10^4$，截断能量 $E_{trunc} \geq 5$。
  4. Wilson链几何： 采用交错Wilson链几何[92, 93]，浴模式按K1=(0,0)，K2=(π,0)，K3=(0,π)，K4=(π,π)顺序排列。尽管这种排列形式上打破了(π,0)和(0,π)补丁之间的对称性，但计算结果显示自能差异小于10^-3，因此在DMFT循环中强制执行对称性Σ(0,π)(ω) ← Σ(π,0)(ω)。
  5. 数据展宽： NRG得到的离散谱函数通过log-Gaussian方案[80]进行展宽，展宽参数 $\sigma=0.7$。
  6. 自能计算： 使用Kugler的对称估计器方法[86]计算自能 $\Sigma_K(\omega)$。
  7. 动态关联函数： 计算本地或短程算符的实频动态关联函数 $\chi[O](\omega)$。能标通过 $\chi[O](\omega)$ 在log-log尺度上的拐点（inflection points）来提取，这些拐点对应于谱函数曲率的变化。
  8. 无对称性破缺： 本研究故意排除了对称性破缺解，所有结果都对应于正常态，以聚焦于量子临界点的内在物理。这与以前的一些簇DMFT研究形成对比，后者可能允许或发现超导或反铁磁序的出现。

通过DCA与NRG的结合，本研究能够以前所未有的分辨率，在极低的能量和温度尺度上直接访问Hubbard模型的实频动力学，从而精确捕捉量子临界行为及其在相图中的具体体现。

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

本研究的核心在于通过DCA+NRG方法，对二维Hubbard模型在特定参数下的零温正常态相图进行精确刻画。所选参数为 $U=7t$ 和 $t'=-0.3t$，这些参数被认为与高温超导铜氧化物相关。以下是计算所得的关键数据及其物理意义：

2.1 相图和能标演化 (图1)

(a) 交叉能标 (Crossover Energy Scales)

图1(a)展示了掺杂驱动的零温正常态相图。我们通过动态关联函数 $\chi[O](\omega)$ 的频率依赖性提取了电荷（红色）、自旋（蓝色）和 $d_{x^2-y^2}$ 波配对（绿色）的交叉能标。这些能标的几何平均值被标记为 $T_{FL}$ 和 $T_{NFL}$。主要发现如下：

  • 量子临界点 p = 0.14：* 这是一个关键的临界掺杂点，标志着赝能隙金属（p < p*）与正常Fermi液体（p > p*）之间的连续量子相变。在此点，$T_{FL}$ 能标被抑制至趋近于零，预示着长程关联的消失和非费米液体行为的出现。
  • FL和NFL区域： 在高于 $T_{FL}$ 的能量下表现出非费米液体（NFL）行为，而在低于 $T_{FL}$ 的能量下表现出费米液体（FL）行为。
  • 赝能隙金属的内部区域： PG金属区域进一步细分为PG_ 和 PG+ 两个子区域，它们不是由相变分隔，而是由交叉点 $P_{LS}$ 和 $p_s$ 标记。这些交叉点对应于反节区谱权重和自能极点位置的变化，而非FL能标的塌缩。

(b) 静态关联函数 (Static Susceptibilities)

图1(b)展示了选择的静态关联函数 $\chi'[O](\omega=0)$ 随掺杂 $p$ 的变化。电荷、自旋和 $d_{x^2-y^2}$ 波配对关联函数在 $p^*$ 附近都表现出显著的峰值，这暗示着在QCP处可能存在发散，与能标的塌缩相符。这表示系统在QCP附近对这些序的波动非常敏感。

2.2 动态d波配对谱函数 (图2)

图2展示了局部 $d_{x^2-y^2}$ 波单重态配对谱 $\chi''[P_{x^2-y^2}](\omega)$ 随掺杂的演化，采用log-log坐标轴。这一谱函数是识别QCP的关键证据：

  • FL行为： 在低于 $T_{FL}$ 的频率下，$\chi''(\omega)$ 表现出与 $\omega$ 成正比的线性行为，这是费米液体的典型特征。
  • NFL平台： 在 $p^*$ 附近，$\chi''(\omega)$ 在中频区域出现一个平台，即 $\chi''(\omega) \approx \text{常数}$，这与Kondo模型等非费米液体中的玻色子谱相似。
  • $T_{FL}$ 塌缩： 随着掺杂 $p$ 从两侧趋近 $p^*$，FL交叉能标 $T_{FL}$ 被推向更低的频率，平台区域向下延伸。在QCP处，数据显示平台延伸到 $\omega \to 0$，这意味着静态关联函数 $\chi'(\omega=0)$ 可能对数发散。

2.3 单粒子谱 (图3)

单粒子谱也清晰地揭示了这一相变：

(a1, a2) 簇谱函数 $A_K(\omega)$

  • 反节区谱 ($K=(0,\pi)$ 或 $K=(\pi,0)$)： 在 $p > p^*$ 的正常FL区，反节区谱在 $\omega=0$ 处有一个狭窄的共振峰。随着 $p$ 减小趋近 $p^*$，这个峰变得更尖锐。然而，在 $p < p^*$ 的赝能隙金属区，这个共振峰消失，取而代之的是一个狭窄的赝能隙。这源于反节区自能在 $p < p^*$ 时出现的极点。
  • 节区谱 ($K=(0,0)$)： 节区谱随着掺杂的变化相对平滑，始终保持相干。
  • $(K=(\pi,\pi))$ 补丁： 在所有考虑的掺杂下都保持带隙，因此未显示。

(b) 准粒子权重 $Z_K$

准粒子权重 $Z_K = [1 - \partial_\omega \text{Re}\Sigma_K(\omega)|_{\omega=0}]^{-1}$ 进一步确认了QCP：

  • 反节区 $Z_{(0,\pi)}$： 随着掺杂从两侧趋近 $p^*$，反节区 $Z_{(0,\pi)}$ 被抑制至趋近于零。这与反节区谱函数结构的改变以及簇关联函数中FL能标的塌缩一致。
  • 节区 $Z_{(0,0)}$： 节区 $Z_{(0,0)}$ 在整个相变过程中保持有限值，这可能暗示在QCP处存在动量选择性的NFL行为。然而，论文指出，这可能需要更高精度的计算来确认是否存在对数抑制。

2.4 动量分辨谱函数与Fermi弧 (图4)

图4展示了动量分辨谱函数 $A_k(\omega)$ 随掺杂穿过QCP的变化：

  • 反节区切片 (图4a)： 在FL侧（$p > p^*$），低能谱权重很高，但在PG侧（$p < p^*$）急剧耗尽，与反节区赝能隙和 $Z_{(0,\pi)}$ 的抑制一致。
  • 节区切片 (图4b)： 节区的色散低能特征随掺杂变化很小，与 $A_{(0,0)}(\omega)$ 的平滑演化一致。
  • Fermi面重构 (图4c1, c2)：
    • 在正常FL区域（$p = 0.155 > p^*$），谱权重沿费米面均匀分布，费米面体积符合Luttinger求和规则，形状类似于自由电子系统。
    • 在PG+区域（$p = 0.125 < p^*$），谱权重强烈依赖于费米面上的位置：节区具有显著谱权重，而反节区谱权重极小，从而形成费米弧，似乎终止于反铁磁区边界。

2.5 零频谱权重角分布 (图5)

图5量化了PG+区域中零频谱权重沿费米面的角分布。它对应于ARPES实验中测量的准粒子权重：

  • 结果显示在 $\theta = \pi/4$（节区）附近达到最大值，并随着 $|\theta - \pi/4| \geq \pi/8$（反节区）急剧衰减至零。这与实验观察到的节区-反节区二分法（nodal-antinodal dichotomy）一致。

2.6 性能数据（隐含）

虽然论文没有明确提供具体的CPU小时数或计算耗时，但所采用的DCA+NRG方法本身就暗示了其高计算成本和资源密集性：

  • NRG的计算复杂度： NRG虽然能处理指数级低的能标，但其计算复杂度随希尔伯特空间大小指数增长。虽然使用了对称性（U(1)电荷和SU(2)自旋）来降低复杂度，并设置了 $N_{keep}=3 \times 10^4$ 的截断，但在探索多个掺杂点和多种算符的动力学时，这仍然需要大量的计算资源。
  • DCA迭代的收敛： DCA方法需要自洽迭代，每次迭代都涉及NRG求解杂质模型。这种迭代过程的收敛可能需要数百甚至数千次迭代，每一轮迭代都需要重新计算杂质格林函数和自能。
  • 参数空间探索： 为了构建相图，研究了多个掺杂点 $p$。此外，还探索了次最近邻跳跃 $t'$ 的多个值（虽然主文本集中于 $t'=-0.3t$，但SM中提到了其他值），这进一步增加了计算量。
  • 高精度需求： 为了捕获QCP的连续性并区分微小效应（例如准粒子权重 $Z_K$ 的对数抑制），需要极高的计算精度，这通常意味着更长的计算时间和更大的 $N_{keep}$。

论文中提到了在GCS超级计算机SUPERMUC-NG和Arnold Sommerfeld Center for theoretical physics的计算资源支持，这进一步证实了这项研究对高性能计算资源的依赖性。这些计算通常会在数万到数十万CPU小时的量级，甚至更高，具体取决于迭代次数和收敛标准。

3. 代码实现细节，复现指南，所用的软件包及开源 repo link

本研究的核心在于通过DCA结合NRG求解器来处理二维Hubbard模型，因此其代码实现涉及两个主要部分：DCA框架和NRG杂质求解器。虽然论文并未直接提供GitHub仓库链接（这在学术论文中很常见），但详细描述了所用方法和软件包，从而可以推断其实现方式和潜在的开源工具。

3.1 DCA框架实现细节

DCA框架主要负责处理晶格问题到簇问题（杂质模型）的映射，以及随后的自洽迭代。核心步骤包括：

  1. 晶格哈密顿量定义： 实现了二维Hubbard模型的哈密顿量，包括最近邻 $t$ 和次最近邻 $t'$ 跳跃项，化学势 $\mu$ 和在位库仑相互作用 $U$。动量空间色散 $\epsilon_k = -2t(\cos k_x + \cos k_y) - 4t' \cos k_x \cos k_y$ 是基础。
  2. 补丁（Patch）划分： 采用了“星形补丁几何”[63]，将布里渊区划分为四个动量补丁 $K_1=(0,0)$、$K_2=(0,\pi)$、$K_3=(\pi,0)$ 和 $K_4=(\pi,\pi)$。每个补丁定义了其积分区域，用于计算补丁分辨的态密度 $\rho_K(\epsilon)$。
  3. 簇哈密顿量构建： DCA将晶格问题映射到一个有效簇杂质模型。簇哈密顿量 $H_{imp}$ 形式上与晶格哈密顿量相似，但只包含簇内相互作用和与浴的耦合。这个浴参数（混合函数 $\Delta_K(z)$）是通过自洽条件从晶格格林函数中提取的。
  4. DCA自洽循环：
    • 步骤一： 从一个猜测的自能 $\Sigma_K(z)$ 开始（通常是零或来自前一个掺杂点的收敛解）。
    • 步骤二： 利用 $\Sigma_K(z)$ 和晶格哈密顿量计算平均的晶格格林函数 $G_K(z)$。公式为 $G_K(z) = \int d\epsilon \frac{\rho_K(\epsilon)}{z + \mu - \epsilon - \Sigma_K(z)}$。
    • 步骤三： 从 $G_K(z)$ 中提取有效的混合函数 $\Delta_K(z)$，作为杂质模型浴的参数。本文使用了非对称估计器（Sup. S-IA）来稳定计算 $\Delta_K(z)$，避免了Dyson方程形式中数值误差对极点附近的放大。
    • 步骤四： 将 $\Delta_K(z)$ 传递给NRG求解器，求解簇杂质模型，得到新的簇自能 $\Sigma_K^{new}(z)$。
    • 步骤五： 将新的自能 $\Sigma_K^{new}(z)$ 与旧的自能 $\Sigma_K(z)$ 进行混合（例如线性混合），以更新 $\Sigma_K(z)$。
    • 步骤六： 重复步骤二到五，直到自能收敛到预设精度。
  5. 对称性处理： 在DMFT循环中，强制执行了 $\Sigma_{(0,\pi)}(\omega) \leftarrow \Sigma_{(\pi,0)}(\omega)$ 的对称性，以保证物理结果的合理性。

3.2 NRG杂质求解器 MuNRG 实现细节

NRG求解器负责在给定混合函数 $\Delta_K(z)$ 的情况下，计算簇杂质模型的格林函数和自能。本研究采用了MuNRG实现[80, 81]：

  1. QSpace张量库： MuNRG是基于QSpace张量库[87-91]实现的。QSpace是一个用于处理具有非阿贝尔对称性的张量网络的开源库，特别适合NRG这种依赖于对称性加速的方法。
  2. 对称性利用： MuNRG利用了U(1)电荷和SU(2)自旋对称性。这些对称性在NRG的对角化过程中大幅减少了需要处理的希尔伯特空间大小，从而使得计算在指数级低能尺度上变得可行。
  3. Wilson链构建： 杂质模型被映射到Wilson链，这是一个半无限链状结构，其耦合强度呈对数衰减。本研究采用了交错（interleaved）Wilson链几何[92, 93]，其中浴模式按照 K1=(0,0), K2=(π,0), K3=(0,π), K4=(π,π) 的顺序进行排列。
  4. 对角化过程： NRG的核心是对Wilson链进行迭代对角化，从链的末端开始，逐步包含更多浴态。每次对角化后，保留一定数量 ($N_{keep}$) 的最低能量多重态。
  5. 谱函数计算： NRG能够直接计算实频谱函数。MuNRG使用Kugler的改进估计器方法[86]来计算自能，以提高数值稳定性。
  6. 数据展宽： NRG输出是离散谱线。为了得到连续的谱函数，采用log-Gaussian方案[80]进行展宽，参数 $\sigma=0.7$。
  7. 无对称性破缺： 计算中严格禁止了对称性破缺，如超导或磁序，以专注于正常态的物理。

3.3 复现指南（高层次概述）

复现这项工作需要对强关联物理、DCA和NRG方法有深入理解，并具备高性能计算环境。以下是高层次的复现步骤：

  1. 环境设置：
    • 安装C++/Python开发环境。
    • 获取并安装QSpace张量库及其依赖。
    • 获取MuNRG代码库。如果不是公开的，可能需要自行实现NRG求解器或联系作者。
    • 实现DCA框架，包括晶格格林函数计算、补丁积分、混合函数更新等模块。
  2. 代码实现：
    • Hubbard模型： 实现二维Hubbard模型的参数设置和晶格格林函数计算。
    • DCA迭代： 编写DCA主循环，集成NRG求解器。
    • NRG求解器： 将MuNRG（或自行实现的NRG）集成到DCA框架中，使其能接收混合函数并返回簇自能。
    • 数据分析： 实现谱函数、关联函数、能标提取等数据后处理和分析工具。
  3. 参数配置：
    • 设定Hubbard模型参数：$U=7t$, $t'=-0.3t$。
    • 设定DCA参数：四补丁几何、掺杂点范围。
    • 设定NRG参数：$\Lambda=8$, $N_{keep}=3 \times 10^4$, $\sigma=0.7$。
    • 设置温度为极低值 ($T=10^{-10}t$)。
  4. 运行计算：
    • 对一系列掺杂 $p$ 值运行DCA自洽迭代。
    • 确保自能收敛。
    • 计算收敛后的单粒子谱、电荷/自旋/配对关联函数等。
  5. 数据分析和可视化：
    • 提取 $T_{FL}$、$T_{NFL}$ 等能标。
    • 绘制相图、谱函数、准粒子权重等。
    • 比较结果与论文中的图表。

虽然论文没有给出直接的GitHub链接，但明确提到了使用的关键软件包：

  • QSpace张量库 (QSpace tensor library)： 这是一个开源的C++库，用于处理具有U(1)和SU(2)等非阿贝尔对称性的张量网络。它是MuNRG的基础。可以在以下引用中找到其论文，通常这些论文会包含项目主页或代码仓库链接。

    • [87] A. Weichselbaum, Non-Abelian symmetries in tensor networks: A quantum symmetry space approach, Annals of Physics 327, 2972 (2012).
    • [90] A. Weichselbaum, QSpace - An open-source tensor library for Abelian and non-Abelian symmetries, SciPost Phys. Codebases, 40 (2024).
    • [91] A. Weichselbaum, Codebase release 4.0 for QSpace, SciPost Phys. Codebases, 40 (2024).
    • 根据这些引用，可以推断QSpace的最新版本可以在SciPost Codebases上找到，并可能通过Google Scholar或作者主页找到其开发仓库。

  • MuNRG： 文中提到了 “MuNRG implementation [80, 81] of NRG [82-86]"。MuNRG很可能是一个基于QSpace的内部实现或对现有NRG代码的修改和优化版本。虽然其代码本身可能未直接公开，但其理论和数值细节在[80, 81]中有所描述。对于复现者来说，可能需要基于QSpace和NRG的原理自行实现或使用其他公开的NRG代码（例如NRG-MPS、open-NRG等）进行修改。

  • DCA框架： DCA框架本身可能是一个作者团队内部实现的代码库，因为它将特定的物理模型、补丁几何和NRG求解器集成在一起。一些通用的DCA框架（例如TRIQS/DCA）存在，但针对本研究的特定参数和NRG求解器集成可能需要定制化开发。

总结： 最核心的开源组件是QSpace张量库。MuNRG和DCA框架很可能是研究团队的内部代码。对于复现者而言，理解论文中描述的算法细节，并利用现有的开源NRG和DCA工具进行定制化开发，是实现复现的关键路径。

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

4.1 关键引用文献及其贡献

这项工作建立在强关联物理和高温超导研究的丰富文献基础之上，以下是一些关键引用及其在本研究中的作用：

  1. 高温超导铜氧化物和赝能隙：

    • [1] P. A. Lee et al., Rev. Mod. Phys. 78, 17 (2006): 对掺杂莫特绝缘体和高温超导物理的综述，奠定了研究背景。
    • [2] B. Keimer et al., Nature 518, 179 (2015): 综述了从量子物质到高温超导的进展，强调了铜氧化物中的核心问题。
    • [3-6] C. Castellani et al., Phys. Rev. Lett. 75, 4650 (1995); C. M. Varma, Phys. Rev. B 55, 14554 (1997); C. M. Varma, Phys. Rev. Lett. 83, 3538 (1999); J. L. Tallon et al., Physica Status Solidi (b) 215, 531 (1999): 这些早期工作提出了赝能隙和量子临界点的概念，为本研究的核心问题提供了基础。
    • [14-19] S. Badoux et al., Nature 531, 210 (2016) 等： 提供了实验上关于铜氧化物中QCP的证据，包括载流子密度变化、霍尔数、比热增强等，证实了理论研究的重要性。

  2. DCA和簇DMFT方法学：

    • [47] A. Georges et al., Rev. Mod. Phys. 68, 13 (1996): DMFT的开创性综述，奠定了强关联系统理论研究的基础。
    • [48-50] M. H. Hettler et al., Phys. Rev. B 58, R7475 (1998); M. H. Hettler et al., Phys. Rev. B 61, 12739 (2000); T. Maier et al., Rev. Mod. Phys. 77, 1027 (2005): DCA方法的早期发展和综述，解释了DCA如何将非局部关联纳入考虑。
    • [63] E. Gull et al., Phys. Rev. B 82, 155101 (2010): 详细讨论了DCA中的动量空间各向异性和赝能隙，包括“星形补丁几何”等细节，为本研究的DCA实现提供了指导。

  3. NRG作为杂质求解器：

    • [80] S.-S. B. Lee and A. Weichselbaum, Phys. Rev. B 94, 235127 (2016): 描述了自适应展宽技术在NRG中提高谱分辨率的方法，对本研究处理实频动力学至关重要。
    • [82-86] K. G. Wilson, Rev. Mod. Phys. 47, 773 (1975); H. R. Krishna-murthy et al., Phys. Rev. B 21, 1003 (1980); F. B. Anders and A. Schiller, Phys. Rev. Lett. 95, 196801 (2005); R. Bulla et al., Rev. Mod. Phys. 80, 395 (2008): NRG方法的开创性工作和综述，提供了NRG理论和计算的基础。
    • [92, 93] A. K. Mitchell et al., Phys. Rev. B 89, 121105 (2014); K. M. Stadler et al., Phys. Rev. B 93, 235101 (2016): 讨论了广义Wilson链和交错NRG，这些是本研究高级NRG实现的关键。

  4. 与Sordi et al.工作的对比：

    • [68] G. Sordi et al., Scientific Reports 2, 547 (2012): 提出在某些参数下（$t'/t=0, U/t \leq 6$），临界性与一阶相变的有限温度端点相关联。这与本研究发现的连续相变形成了对比，是后续研究需要弥合的差异。

  5. 边缘费米液体 (MFL)：

    • [74] C. M. Varma et al., Phys. Rev. Lett. 63, 1996 (1989): 描述了铜氧化物正常态的MFL现象学，为本研究在QCP处节区准粒子权重Z(0,0)的可能对数抑制提供了理论背景。

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

尽管这项工作取得了显著进展，但仍存在一些局限性，值得在未来的研究中加以改进和探索：

  1. 准粒子权重 $Z_K$ 的精确性： 论文指出，节区准粒子权重 $Z_{(0,0)}$ 在QCP附近似乎保持有限值，这与FL能标的完全塌缩形成对比，可能暗示动量选择性NFL。然而，论文也承认，要精确解析这种对数抑制（如边缘FL理论预测的），需要极高精度的数值计算，并且当前数据可能不足以给出明确的零温结论。未来的工作需要进一步提升计算精度，以确定 $Z_{(0,0)}$ 是否真的保持有限或仅是对数抑制。

  2. 有限簇尺寸效应： DCA是一种簇近似方法，本研究使用了2x2（四补丁）簇。虽然论文指出“定性结果不强烈依赖于簇方案或周期化”，但更大尺寸的簇可以提供更精确的动量分辨信息，尤其是在费米面重构和Luttinger面拓扑结构等问题上。例如，论文中提到，“改善k分辨精度，使用更大簇可能导致Luttinger面表现出不同行为”。这暗示了当前簇尺寸可能对一些精细的拓扑特征存在限制。

  3. 零温近似的普适性： 尽管NRG在零温下表现出色，但实际的高温超导铜氧化物是在有限温度下运行的。虽然本研究在零温下确立了QCP的存在，但理解有限温度下NFL区域的量子临界标度行为（如论文中提及的后续工作[75]）对于与实验进行更直接的比较至关重要。

  4. 与一阶相变关联的差异： 论文明确指出，与Sordi等人[68]在不同参数区域（$t'/t=0, U/t \leq 6$）发现的与一阶相变端点相关的临界性不同，本研究在 $U=7t, t'=-0.3t$ 参数下发现的是连续量子相变。虽然论文认为这两种情况可能代表了相图的不同区域，但需要更系统的工作来桥接这些差异，并理解Hubbard模型相图中各种临界行为的完整图景，包括是否存在隐藏的一阶相变或极低温度下的临界终点。

  5. 无对称性破缺： 本研究故意排除了对称性破缺解，所有分析都限制在正常态。然而，在真实铜氧化物中，超导序、电荷序或磁序都可能与赝能隙和QCP相互作用。未来的研究需要将这些有序态纳入考虑，以提供更全面的相图，并理解有序态与量子临界性之间的竞争与协同。

  6. 计算成本与参数空间探索： 尽管DCA+NRG是强大的方法，但其计算成本仍然很高。这限制了在更大参数空间中（例如更宽范围的U和t’，或探索更复杂的哈密顿量）进行系统探索的可能性。更高精度的计算（更大的 $N_{keep}$）也会进一步增加成本。

  7. 无序效应： 论文提及一些实验数据暗示存在扩展的量子临界区域可能由无序引起[21]。本研究排除了无序效应，因此其结果主要描述了理想晶格下的内在量子临界行为。将无序纳入理论模型将是一个重要的扩展方向。

这些局限性并非对本研究质量的否定，而是科学探索中常见的里程碑，为未来的工作指明了方向。通过解决这些问题，我们可以对高温超导体的复杂物理有更全面、更深入的理解。

5. 其他你认为必要的补充

5.1 该研究在高温超导领域的里程碑式意义

这项工作代表了高温超导（high-Tc superconductivity）铜氧化物研究领域的一个重要里程碑，主要体现在以下几个方面：

  1. 解决了长期争论的核心问题： 论文明确指出，在二维Hubbard模型中，赝能隙金属（PG metal）与正常Fermi液体（FL）之间存在一个连续的量子临界点（QCP）。这一发现为长达三十年的争论提供了明确的理论答案，即这两种截然不同的正常态物理是相互关联并通过一个QCP转变的，而不是简单地互不相干或通过一阶相变分隔。
  2. 提供了PG-QCP的微观机制： 研究详细揭示了QCP的特征，包括FL能标的连续塌缩，动态关联函数中非费米液体（NFL）行为的出现，以及单粒子谱在反节区处的急剧变化。这些微观层面的洞察对于理解QCP的本质至关重要。
  3. 理论上复现了实验观察： 论文成功地在理论模型中复现了铜氧化物实验中观察到的关键现象，如费米弧（Fermi arcs）、节区-反节区二分法（nodal-antinodal dichotomy）、以及反节区谱权重在欠掺杂区域的耗尽。这些理论与实验的契合极大地增强了Hubbard模型作为描述铜氧化物基石的可靠性。
  4. 深化了对正常态相图的理解： 通过详细绘制掺杂-温度相图（尽管本文侧重于零温），本研究为理解高温超导相图中的“正常态”区域（包括赝能隙金属、奇异金属和Fermi液体）提供了坚实的理论框架，有助于统一不同区域的物理。

5.2 DCA+NRG方法学的突破

本研究的方法学本身就是一项重要的成就，展示了DCA+NRG在强关联计算中的独特优势：

  1. 克服了传统方法的局限性： 传统的量子蒙特卡洛（QMC）方法受限于符号问题和有限温度，难以获得精确的零温实频动力学。大多数DMFT方法虽然能处理无穷维极限，但在捕捉非局部关联方面存在不足。DCA+NRG的结合完美地弥补了这些不足，使得研究能够在极低的能量和温度尺度上直接访问实频动力学。
  2. 高低能分辨率的兼顾： NRG以其在解析指数级低的能量尺度上的高分辨率而闻名，这对于识别QCP的能标塌缩至关重要。DCA则负责处理短程非局部关联。两者的结合使得系统能够同时捕捉高能和低能的物理特征，这是理解强关联系统多尺度行为的关键。
  3. 实频谱函数的直接获取： 避免了从Matsubara频率到实频的解析延拓问题，提供了更可靠的谱学信息，这对于与角分辨光电子能谱（ARPES）等实验数据进行直接比较至关重要。
  4. 为未来研究提供了强大工具： DCA+NRG的成功应用，为未来研究其他强关联系统中的量子相变、非费米液体行为和复杂的正常态物理提供了一个强大的数值工具。

5.3 对非费米液体行为的深入探讨

本研究在QCP附近揭示的非费米液体（NFL）行为是其另一大亮点：

  1. NFL区域的明确界定： 动态关联函数中出现的中频平台（$\chi''(\omega) \propto \text{const.}$）清晰地界定了NFL区域，其范围在FL能标 $T_{FL}$ 和更高能标 $T_{NFL}$ 之间。随着接近QCP，$T_{FL}$ 的塌缩使得NFL区域扩展到更低的能量尺度，最终在QCP处可能延伸到零频。
  2. 与边缘FL理论的关联： 论文提及了节区准粒子权重 $Z_{(0,0)}$ 的可能对数抑制，这与边缘费米液体（Marginal Fermi Liquid, MFL）现象学[74]相吻合。MFL是描述铜氧化物正常态中奇异金属行为的重要理论框架，其特征是准粒子寿命呈现线性衰减，导致谱函数在低能处平坦化。本研究为这种理论在微观Hubbard模型中的实现提供了证据。
  3. 量子临界标度行为： 尽管本文侧重于零温相图，但其发现的NFL区域为后续进行量子临界标度分析奠定了基础。论文提及的后续工作[75]将更深入地探讨有限温度下的NFL区域的标度行为，这将进一步揭示QCP的普遍性类别。

5.4 费米弧和节区-反节区二分法的解释

本研究通过理论模型成功解释了铜氧化物中的关键实验现象：

  1. 费米弧的形成机制： 在赝能隙金属区域（p < p*），计算结果显示节区具有显著谱权重，而反节区谱权重急剧耗尽，从而形成了费米弧。这与ARPES实验观察到的费米弧现象[8-12, 76]高度一致，并且为理解其微观起源（即反节区赝能隙的出现）提供了直接证据。
  2. 节区-反节区二分法： 零频谱权重沿费米面的角分布（图5）清晰地展现了节区（$\theta = \pi/4$ 附近）谱权重最大，而反节区谱权重急剧衰减的现象。这种节区-反节区二分法是铜氧化物正常态物理的标志性特征，本研究提供了其在Hubbard模型中的理论解释。

5.5 未来研究方向

这项工作不仅解决了重要问题，也为未来的深入研究指明了方向：

  1. 更高精度计算 $Z_K$： 为了确定节区准粒子权重 $Z_{(0,0)}$ 在QCP处的精确行为，需要更高精度的数值计算，这可能涉及更大的 $N_{keep}$ 或更优化的NRG算法，以分辨对数抑制与有限值之间的细微差别。
  2. 更大簇尺寸的DCA： 探索更大尺寸的簇（例如4x4簇）可以提供更精细的动量分辨信息，尤其是在费米面拓扑结构和Luttinger面演化方面，可能揭示当前四补丁近似未捕获的物理细节。
  3. 有限温度下的量子临界标度分析： 论文提及后续工作[75]将进行有限温度下的量子临界标度分析。这将是理解QCP在非零温下的表现以及与实验数据更紧密联系的关键一步。
  4. 引入对称性破缺： 本研究故意排除了对称性破缺。未来的工作可以考虑在DCA框架中允许超导、电荷或磁序的出现，从而研究这些有序态与赝能隙以及量子临界性之间的竞争和协同作用，构建更完整的Hubbard模型相图。
  5. 无序效应的研究： 考虑到实验中无序可能导致扩展的量子临界区域[21]，将无序效应纳入DCA+NRG框架将是一个重要的研究方向，以理解无序对量子临界行为的影响。
  6. 应用于其他强关联材料： 将DCA+NRG方法应用于其他具有类似挑战性正常态物理的强关联材料，例如镍氧化物（nickelates）等，可以拓展其普适性，并为这些新材料的研究提供理论指导。

通过这些未来研究，我们可以进一步完善对Hubbard模型的理解，并为最终揭示高温超导的奥秘奠定坚实基础。