永冈超金属态：粒子掺杂三角哈伯德模型的深度解析

来源论文: https://arxiv.org/abs/2605.13837v1 生成时间: May 14, 2026 07:26

执行摘要

本研究深入探讨了粒子掺杂三角晶格哈伯德模型中一个新颖的量子态——永冈超金属态（Nagaoka Supermetal, NS）。通过结合动态簇近似（DCA）计算和有效的低能模型推导，我们揭示了该态的亚线性直流电阻率、电荷可压缩性和零频谱权重的奇异行为。这些异常特性源于相互作用驱动的高阶范霍夫奇点（HOVHS），并展示了该模型在几何阻挫体系中非费米液体行为的独特视角，为未来的量子模拟实验提供了明确的指导。

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

核心科学问题

在凝聚态物理中，相互作用和几何阻挫的复杂交织为探索奇异量子物态提供了肥沃的土壤。特别是在掺杂的莫特绝缘体中，随着载流子的注入，金属态的性质及其起源仍然是悬而未决的核心科学问题。费米液体理论是描述相互作用费米子的基石，但它在强关联体系中常常失效，导致诸如奇异金属、赝能隙等非费米液体（NFL）行为的出现。最近，一种由高阶范霍夫奇点（HOVHS）驱动的新型NFL态——超金属态——被提出，其特征是态密度（DOS）的幂律发散和亚线性电阻率（ρ ~ T^α，α < 1）。然而，尽管实验研究已在一些固体平台（如 Kagome 金属）中发现了 HOVHS 的存在，但直接观测到超金属态的实验证据仍待证实。

本研究的核心科学问题在于：粒子掺杂的三角晶格哈伯德模型是否能自发地产生这种超金属态？如果是，其微观机制是什么？它是否像传统超金属那样依赖于单粒子能带结构的精细调整，还是由相互作用本身内在驱动？此外，这种超金属态的特性（如电阻率和可压缩性的标度指数）如何随掺杂水平演化？理解这些问题对于揭示几何阻挫体系中强关联电子的丰富物理性质、以及为超冷原子量子模拟实验提供理论指导至关重要。三角晶格因其固有的几何阻挫特性，抑制了传统的反铁磁关联，从而放大了掺杂的作用，使其成为研究巡游铁磁性等新奇物理的理想平台。本文正是旨在通过理论计算，为永冈超金属态的存在提供强有力的证据，并阐明其微观起源。

理论基础

本研究的理论基础主要建立在以下几个关键概念之上：

1. 哈伯德模型 (Hubbard Model): 它是描述强关联电子系统最基本的模型之一，通过考虑电子在晶格上的跳跃（动能项）和同格点上的库仑排斥（相互作用项）来捕捉电子之间的关联效应。在本研究中，具体采用的是三角晶格哈伯德模型，其哈密顿量为 $H = -t \sum_{\langle i,j \rangle, \sigma} (c^\dagger_{i\sigma} c_{j\sigma} + H.c.) + U \sum_i n_{i\uparrow} n_{i\downarrow} - \mu \sum_i (n_{i\uparrow} + n_{i\downarrow})$。三角晶格引入了几何阻挫，这意味着即使在半满填充时，传统的反铁磁长程有序也很难形成，从而为其他新奇的量子态（如巡游铁磁性）的出现提供了条件。

2. 永冈铁磁性 (Nagaoka Ferromagnetism): 这是强关联电子系统中的一个重要现象，最初由永冈洋一提出，认为在近乎半满填充的莫特绝缘体中，单个空穴或粒子可以诱导铁磁性基态。其核心思想是，空穴（或粒子）的引入通过降低电子的动能来稳定铁磁构型，即动能驱动的铁磁性。在三角晶格中，这种机制被认为能促进巡游铁磁关联，其能量尺度由跳跃参数 $t$ 决定，远超超交换作用 $J \sim t^2/U$，因此可以在较高温度下保持稳定。本研究中的“永冈超金属”即暗示了这种铁磁背景下的新奇金属态。

3. 费米液体理论与非费米液体 (Fermi Liquid Theory and Non-Fermi Liquid): 费米液体理论是描述相互作用费米子系统的标准范式，其中激发可以被视为具有良好定义的寿命和有效质量的准粒子。然而，在强关联体系或存在量子临界点附近，准粒子概念可能失效，系统进入非费米液体状态。非费米液体表现出许多异常特性，如违反维德曼-弗朗兹定律、亚线性电阻率、态密度的幂律发散等。本研究关注的超金属态正是一种特定的非费米液体。

4. 高阶范霍夫奇点 (Higher-Order Van Hove Singularity, HOVHS): 范霍夫奇点是能带结构中梯度为零的点，通常导致态密度（DOS）的对数发散。高阶范霍夫奇点则是能带在某个动量点附近不仅梯度为零，而且更高阶的导数也为零，从而导致DOS以更高的幂律发散（例如，DOS ~ |ε|^-γ）。这种强烈的DOS发散可以驱动各种奇异物态，包括超导、电荷密度波和本研究中的超金属态。本研究中，HOVHS不是来自简单的单粒子能带结构，而是通过相互作用驱动的能带重构（具体是通过永冈极化子介导的次近邻跳跃）自发产生的。

5. 有效低能模型 (Effective Low-Energy Model): 为了从微观层面理解永冈超金属态的起源，研究通过二阶微扰理论将哈伯德模型映射到一个简化的有效低能模型。这个模型描述了在极化子自旋背景下，双占据态（或空穴）的动力学。关键在于，相互作用诱导的三位点跳跃过程（即涉及次近邻跳跃的项）被发现是捕捉三角晶格中奇异特性的关键。通过分析这个有效模型的色散关系，可以发现特定的参数条件（如动能重整化因子Z与相互作用U的关系）会导致能带在M点附近发生重构，形成接近四次方的色散关系，从而产生HOVHS。

这些理论基础共同构成了理解本研究中永冈超金属态的关键框架，解释了其在强关联、几何阻挫背景下，通过相互作用驱动的能带重构，如何自发涌现出HOVHS并导致一系列非费米液体特性。

技术难点

研究强关联费米子系统，特别是具有几何阻挫的系统，面临诸多显著的技术挑战：

1. 强关联多体问题: 哈伯德模型是一个典型的多体问题，电子之间的相互作用导致波函数高度纠缠，无法使用单粒子近似处理。精确求解其基态和激发态需要处理指数级增长的希尔伯特空间，这在计算上是不可行的。尤其是在掺杂情况下，粒子数不再是固定的，这进一步增加了问题的复杂性。

2. 几何阻挫: 三角晶格的几何阻挫特性使得传统的平均场理论或基于序参量的理论难以捕捉其丰富的物理性质。阻挫抑制了长程反铁磁序的形成，导致系统在量子涨落面前更不稳定，可能涌现出多种竞争的量子态。这要求计算方法能够有效地处理非局域关联和量子涨落。

3. 费米子符号问题 (Fermion Sign Problem): 在量子蒙特卡罗（QMC）模拟中，费米子系统的路径积分包含负权重项，导致蒙特卡罗采样的效率极低甚至崩溃，这使得DQMC等精确数值方法在低温度和二维阻挫体系中应用受限。虽然DCA在一定程度上缓解了这一问题，但仍然是一个考量因素，尤其是在大簇和低温度下。

4. 动态簇近似 (DCA) 的计算复杂性: DCA 虽然是解决强关联问题的强大工具，通过将无限晶格映射到有限簇来捕捉短程空间关联，但其计算成本随簇尺寸的增加而呈指数级增长。这意味着能够模拟的簇尺寸有限（本研究中为 Nc=3, 4, 6），可能无法完全捕捉所有尺度的关联效应。选择合适的簇几何形状和大小，以平衡计算成本和物理精度，是一个重要挑战。

5. 解析延拓 (Analytic Continuation): DCA 计算在虚时（Matsubara频率）上进行，但为了获取真实的物理量（如实频谱函数），需要将虚时自能解析延拓到实频轴。这是一个不适定问题（ill-posed problem），即使输入数据存在微小误差，也会导致实频谱函数出现大的波动或伪影。本研究使用了 Padé 近似和最大熵方法进行交叉验证，但其内在的不确定性仍然存在。

6. 有效模型的推导与参数化: 从复杂的哈伯德模型推导出低能有效哈密顿量，并将其参数与原始模型参数关联起来，需要细致的微扰理论分析。确保有效模型能够准确捕捉关键物理（如相互作用驱动的能带重构和HOVHS）是一个精细的建模过程。动能重整化因子Z的确定也依赖于自能的低频极限，这与解析延拓的挑战紧密相关。

方法细节

为了克服上述挑战并深入研究粒子掺杂三角晶格哈伯德模型中的永冈超金属态，本研究采用了动态簇近似（Dynamical Cluster Approximation, DCA）结合连续时间混合杂化展开（Continuous-Time Hybridization Expansion, CTHYB）求解器，并辅以有效的低能模型推导和多项物理量计算。

1. 动态簇近似 (DCA): DCA 是一种将无限晶格映射到有限簇的强关联多体方法，它通过在动量空间中对布里渊区进行粗粒化来恢复自能的动量依赖性。这使得 DCA 能够捕捉单格点DMFT（动态平均场理论）无法描述的短程空间关联，同时在平均场水平上处理长程关联。在本研究中，DCA 的核心是将哈伯德模型问题转化为一个局域的“簇杂质问题”，该簇杂质嵌入一个自洽确定的有效介质中。

   • 簇几何结构: 考虑到计算成本和几何阻挫的特性，本研究采用了 Nc=3、4 和 6 的簇尺寸。特别是，4 站点簇被发现足以捕捉关键物理，因为其动量粗粒化包含 M 点，而 M 点是三角晶格中标准范霍夫奇点和本研究中相互作用诱导能带重构的关键区域。
   • 自洽循环: DCA 迭代求解一个自洽循环：从一个初始的自能猜测开始，计算簇格林函数，然后更新杂化函数，再用杂化函数求解簇杂质问题（得到新的簇自能），最后将簇自能映射回无限晶格自能，重复此过程直至收敛。

2. 连续时间混合杂化展开 (CTHYB) 求解器: 簇杂质问题通过 CTHYB 求解器 [60] 来解决。CTHYB 是一种量子蒙特卡罗方法，它在虚时（Matsubara 频率）上计算簇格林函数和自能。相比于早期的蒙特卡罗方法，CTHYB 能够更有效地处理强关联效应，并且在一定程度上缓解了费米子符号问题，使其适用于研究掺杂的哈伯德模型。本研究利用 TRIQS 库 [61] 实现了 CTHYB 求解器。

3. 解析延拓: 从 CTHYB 求解器得到的虚时自能 $\Sigma(i\omega_n)$ 需要进行解析延拓以获得实频自能 $\Sigma(\omega)$，进而计算实频物理量。本研究主要采用 Padé 近似 [70] 进行解析延拓，并通过 最大熵方法 (MaxEnt) [66, 85, 86] 进行结果的交叉验证，以确保其可靠性。

4. 物理量计算:

   • DC 电阻率 ($\\rho$): 根据 Kubo 公式，DC 电阻率 $\\rho = \lim_{\omega \to 0} [1/\sigma(\omega)]$，其中光学电导率 $\sigma(\omega)$ 由公式 (3) 给出。$\\sigma(\omega) = \sigma_0 \sum_{\mathbf{k}} \int d\epsilon d\nu X(\epsilon) A(\epsilon,\nu) A(\epsilon, \nu+\omega) \frac{n_F(\nu) - n_F(\nu+\omega)}{\omega}$。
   • 近邻自旋关联函数 ($\\langle S_i S_j \rangle_c$): 定义为 $\\langle S_i S_j \rangle - \\langle S_i \rangle \\langle S_j \rangle$，用于表征磁性背景。
   • 电荷可压缩性 ($\\chi_c$): 定义为 $\\partial \\langle n \\rangle / \\partial \\mu$，衡量系统对化学势变化的响应，是探测奇异态的重要热力学量。
   • 零频谱权重 ($A(\\omega=0)$): 通过 Matsubara 格林函数 $G(i\omega_n)$ 在最低频率处外推得到 $A(\\omega=0) = -\frac{1}{\pi} \lim_{\omega_n \to 0^+} \text{Im}G(i\omega_n)$，它与态密度在费米能级附近的强度密切相关。

5. 有效低能模型推导: 为了揭示永冈超金属态的微观起源，研究采用了二阶微扰理论推导了一个有效低能模型。该模型将哈伯德模型投影到由双占据态构成的低能子空间，其中单占据的自旋背景充当有效真空。核心发现是，相互作用诱导的三位点跳跃过程产生了有效的次近邻（NNN）跳跃项，其强度与背景磁化强度 $m$ 相关。通过分析这个有效模型（如公式 (5) 和 (6) 所示）在动量空间 M 点附近的色散关系，可以发现当动能重整化因子 $Z$ 满足特定条件时，能带会发生重构，导致高阶范霍夫奇点 (HOVHS) 的形成。

这些方法的综合运用使得本研究能够从微观和宏观层面全面刻画永冈超金属态，并揭示其与相互作用驱动的能带重构和HOVHS之间的深刻联系。

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

关键 Benchmark 体系

本研究的理论框架和计算结果通过与现有实验数据及高精度数值模拟（如确定性量子蒙特卡罗，DQMC）进行对比，得到了可靠性验证。这主要包括以下几个方面：

1. 超冷原子量子模拟实验: 论文明确指出，本研究的发现可以直接通过当前超冷原子实验进行验证。特别是，超冷费米原子在光学晶格中为模拟强关联物理提供了原始平台，量子气体显微镜技术可以实现对多体关联的位点分辨测量。此前，这类实验已成功揭示了掺杂载流子的微观结构，包括磁极化子的形成 [26-30] 和赝能隙的出现 [31, 32]。尤其重要的是，在掺杂三角晶格中，永冈铁磁极化子已通过实验直接观测到 [33-39]。本研究的计算参数（如相互作用强度 U/t = 12 和温度 T/t ≈ 1）与当前量子模拟实验的可达范围高度一致 [23, 25, 31-33, 35]，这使得我们的理论预测具有很强的实验可检验性。

2. 确定性量子蒙特卡罗 (DQMC) 模拟: 为了进一步验证动态簇近似（DCA）在掺杂体系中的可靠性，研究将 DCA 结果与在 $L=9 \times 9$ 晶格上进行的 DQMC 模拟结果进行了对比 [62-65]。

   • 双占据率和电荷可压缩性: 如图 S4(a) 和 S4(c) 所示，在较高的温度下，DCA 计算得到的双占据率（$\\langle n_\uparrow n_\downarrow \rangle$）和电荷可压缩性（$\\chi_c$）与 DQMC 模拟结果显示出定量上的一致性。尽管在较低温度下存在细微偏差，但整体趋势高度吻合。
   • 近邻自旋关联函数: 如图 S4(b) 所示，DCA 准确地捕捉了 DQMC 数据中自旋关联函数（$\\langle S_i S_j \rangle_c$）的定性趋势，进一步证实了其在磁性性质描述上的可靠性。

3. 三位点关联实验数据: 研究还将 DCA 结果与来自参考文献 [35] 的实验测量数据（关于双占据率和三位点关联）进行了比较。如图 S3 所示，DCA 框架不仅在双占据率上与实验数据（圆点）取得了定量一致，而且在三位点关联函数（特别是空穴-双占据-自旋-自旋关联函数 $C^{(3)}_{h/d}$ 和空穴-双占据-空穴关联函数 $C^{(3)}_{h/h}$）的定性特征上也表现出良好的一致性。这些关联是捕捉掺杂载流子和自旋背景之间相互作用的关键指标。

通过这些细致的 Benchmark，本研究建立了 DCA 框架在描述三角晶格哈伯德模型中的电荷和磁性自由度方面是可靠的，能够捕捉到其核心的多体物理。

计算所得数据

本研究通过DCA计算获得了大量数据，全面揭示了粒子掺杂三角哈伯德模型中永冈超金属态的特性：

1. 相图 (图 1)：

   • 低掺杂伪能隙 (PG) 和奇异金属 (SM) 态：在粒子掺杂侧（δ > 0），随着掺杂量的增加，系统从莫特绝缘体演变为赝能隙（PG）和奇异金属（SM）态。
   • 永冈超金属 (NS) 态的涌现：在较高粒子掺杂（δ > 0）和低温度下，系统进入永冈超金属（NS）态。该态与永冈铁磁关联并存。
   • 鲁棒性: NS 态表现出显著的鲁棒性，在高达 T/t ≈ 1 的温度下依然存在，这表明其能量尺度由动能 $t$ 而非超交换 $J$ 决定。
   • 空穴掺杂侧: 对比显示，空穴掺杂（δ < 0）侧主要由120° Néel 反铁磁关联主导，并依次经历 PG、SM 和 FL 态，与粒子掺杂侧形成鲜明对比，突出了粒子-空穴不对称性。

2. 直流电阻率 ($\\rho$) (图 2a)：

   • 亚线性温度依赖: 在 NS 态区域，ρ 表现出亚线性温度依赖，即 $\\rho \sim T^\alpha$，其中 $\\alpha < 1$。具体地，对于 δ = 0.5，α ≈ 0.5；对于 δ = 0.4，α ≈ 0.66。这与奇异金属或坏金属的 T-线性标度行为 [3, 4, 74-78] 显著不同。
   • 与磁性关联的交叉: 随着温度升高，系统会经历一个交叉，电阻率恢复到标准的 T-线性依赖。这一交叉与正近邻自旋关联的建立密切相关（图 2b），表明铁磁关联驱动了亚线性标度行为的产生。
   • 低掺杂行为: 对于 δ = 0.1 的低掺杂情况，近邻自旋关联保持反铁磁性，电阻率在整个温度范围内保持线性，没有出现超金属行为。

3. 近邻自旋关联函数 ($\\langle S_i S_j \rangle_c$) (图 2b)：

   • 在较高掺杂和低温度下，$\\langle S_i S_j \rangle_c$ 变为正值并显著增强，表明铁磁关联的建立。这种铁磁背景与电阻率的亚线性行为同时出现，证明了 NS 态的永冈驱动起源。

4. 电荷可压缩性 ($\\chi_c$) (图 2c)：

   • 幂律发散: 在 NS 态区域，$\\chi_c$ 表现出幂律发散，$\\chi_c \sim T^{-\gamma}$，其中 $\\gamma > 0$。具体地，对于 δ = 0.5，γ ≈ 0.25；对于 δ = 0.4，γ ≈ 0.16。
   • 标度关系验证: 这些指数与超金属理论预测的 $\\alpha = 1 - 2\\gamma$ 关系高度吻合（例如，δ = 0.5 时，0.5 = 1 - 2 * 0.25），为 NS 态的存在提供了自洽验证。

5. 动能重整化因子 ($Z$) (图 3a)：

   • HOVHS 的判据: 当 Z 满足 0 < Z < 9t/U (对于 U/t = 12，约为 0 < Z < 0.75) 条件时，系统进入 HOVHS 区域。
   • 与电阻率交叉的关联: Z 进入 HOVHS 区域的起始温度与电阻率的交叉温度（图 2a）相匹配，进一步证实了 HOVHS 与 NS 态的关联。

6. 零频谱权重 ($A(\\omega=0)$) (图 3b)：

   • 幂律发散: 在 NS 态区域，A(ω=0) 表现出幂律发散，如 T^-0.25 (δ=0.5) 和 T^-0.16 (δ=0.4)。
   • 与 HOVHS 的关联: 对于 δ = 0.5，A(ω=0) ~ T^-0.25 的标度与非相互作用极限下 HOVHS 的 $|SE|^{-1/4}$ 奇异性精确匹配，有力支持了 HOVHS 是 NS 态微观起源的观点。
   • 与费米液体的对比: 对于 δ = 0.1，A(ω=0) 表现为 $\ln(1/T)$ 标度，这更接近于边际费米液体行为，而非超金属态。

7. 相互作用强度 ($U$) 的影响 (图 4)：

   • 鲁棒性: 亚线性输运在宽泛的 $U$ 范围内保持鲁棒。
   • 非单调起始温度: NS 态的特征起始温度 $T_c$（由电阻率交叉点推断）显示出对 $U$ 的非单调依赖性，在中等耦合强度下达到最大值，这反映了动能驱动的铁磁极化子形成与强关联抑制次近邻跳跃之间的竞争。

这些详尽的计算数据共同描绘了永冈超金属态的独特图景，揭示了其在强关联几何阻挫体系中的涌现机制和关键特性。

性能数据

本研究中的数值计算性能主要体现在所采用方法的计算复杂度、硬件资源利用以及所能达到的模拟规模和精度上。虽然论文没有提供具体的 CPU/GPU 小时数或内存使用量，但可以从方法论和所用计算平台推断出以下性能方面的信息：

1. 动态簇近似 (DCA) 的计算复杂度: DCA 的计算成本随着簇尺寸 ($N_c$) 的增加而呈指数级增长，尤其是在使用连续时间量子蒙特卡罗（CTQMC，如 CTHYB）求解器时。论文指出，由于“求解器计算成本随簇尺寸指数增长”以及“几何阻挫三角晶格上严重的费米子符号问题”，本研究的计算被限制在 $N_c = 3, 4, 6$ 的簇尺寸。这表明即使是相对较小的簇，也需要显著的计算资源。

2. 费米子符号问题: 尽管 CTHYB 求解器在一定程度上缓解了费米子符号问题，但在阻挫体系和低温度下，这一问题仍然是限制计算性能和可达温度范围的关键因素。论文提到“严格的费米子符号问题”，暗示了即使在DCA框架下，低温度区域的模拟仍然具有挑战性，需要更高的计算投入才能获得可靠结果。

3. 解析延拓的计算开销: 将虚时自能解析延拓到实频轴是一个不适定问题，需要大量的计算来获得稳定可靠的结果。Padé 近似相对较快，但对噪声敏感；最大熵方法（MaxEnt）虽然更鲁棒，但通常计算成本更高，并且需要对谱函数进行先验假设。论文通过两种方法的交叉验证来确保结果的可靠性，这无疑增加了整体的计算负担。

4. 计算平台: 论文明确指出，“数值计算是在天津国家超级计算中心的‘天河-1A’集群上进行的”。“天河-1A”是一个大规模并行超级计算机系统，这表明本研究需要利用高性能计算（HPC）资源才能完成。使用这种级别的超算集群，意味着单个计算任务（如一个DCA自洽循环或一次解析延拓）可能需要数百到数千个处理器核心并行运行数小时到数天，才能达到所需的精度和参数范围。

5. 参数空间探索: 为了构建完整的相图、研究不同掺杂和相互作用下的物理量，并进行温度依赖性研究，需要对大量的参数组合进行计算。例如，在 Fig. 2、3 和 4 中，针对多个掺杂点、温度和相互作用强度都进行了详细计算。这种大规模的参数扫描是强关联物理研究的常见需求，也对计算资源提出了极高的要求。

总而言之，本研究的性能数据虽然没有以传统意义上的“运行时间”或“吞吐量”形式给出，但从其采用的计算方法、处理的物理复杂性、克服的计算挑战以及所使用的超算平台来看，这是一个计算密集型项目，依赖于高性能计算资源以在可接受的时间内获得可靠的物理结果。

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

代码实现细节

本研究的数值模拟主要基于动态簇近似（DCA）框架，结合连续时间混合杂化展开（CTHYB）量子蒙特卡罗求解器以及后续的解析延拓技术。其实现细节涉及到多个模块的协同工作：

1. DCA 框架: DCA 的核心是一个自洽循环，它将无限晶格问题映射到有限簇上的杂质模型。

   • 布里渊区粗粒化: 首先，将布里渊区划分为 $N_c$ 个动量块（patches），每个块由一个簇动量 $Q$ 代表。簇尺寸 $N_c$ 的选择（本研究中使用 3、4、6 站点簇）决定了动量空间的分辨率。对于三角晶格，确保簇几何能够覆盖关键的高对称点（如 M 点和 K 点）至关重要。
   • 有效介质构建: 根据格林函数和自能的迭代，构建一个有效的杂化函数 $\Delta_Q(i\omega_n)$，它描述了簇与周围无限晶格的耦合。
   • 簇自能计算: 在每个自洽迭代中，通过 CTHYB 求解器解决 $N_c$ 站点簇的杂质问题，得到簇的格林函数 $G_c(Q, i\omega_n)$ 和自能 $\Sigma_c(Q, i\omega_n)$。
   • 自洽收敛: 迭代过程直到簇自能或平均粒子数达到预设的收敛标准。

2. CTHYB 求解器: 作为 DCA 的核心组件，CTHYB 求解器在 Matsubara 频率轴上对簇的格林函数和自能进行采样。

   • 高阶展开: CTHYB 通过对杂化函数在相互作用自由哈密顿量中的展开来实现，其主要优点是可以在任意相互作用强度下工作，并且对高阶关联有较好的描述。
   • 测量: 在蒙特卡罗模拟过程中，测量格林函数和自能，并通过统计平均得到它们在虚时轴上的值。
   • 温度与噪声: CTHYB 的计算在虚时进行，并且其结果会带有统计噪声，尤其是在低温区域。这些噪声对后续的解析延拓构成了挑战。

3. 解析延拓: 从虚时自能到实频自能的转换是关键一步，其准确性直接影响到物理量的可靠性。

   • Padé 近似: 本研究主要使用 Padé 近似进行解析延拓。这是一种数学技术，通过将函数表示为两个多项式的比值来近似函数。Padé 近似对噪声敏感，但由于其较高的低温分辨率而被选用。
   • 最大熵方法 (MaxEnt): 作为交叉验证，MaxEnt 是一种更鲁棒的正则化方法，可以处理噪声数据，但通常需要先验信息且计算成本较高。

4. 物理量计算: 根据解析延拓得到的实频自能，计算各种物理量：

   • 光学电导率: 使用 Kubo 公式，涉及到动量空间积分和谱函数的计算。
   • 电荷可压缩性: 通过粒子数对化学势的导数计算。
   • 零频谱权重: 从 Matsubara 格林函数的外推得到。
   • 自旋关联: 直接从 CTHYB 测量得到。

5. 有效低能模型推导: 这是一个理论推导过程，通过二阶微扰理论从哈伯德模型得到一个简化的有效哈密顿量。这涉及到对原始哈密顿量的投影，以及对中间态的求和。

整个代码实现是多学科的交叉，结合了量子场论、统计物理、数值线性代数和蒙特卡罗方法。

复现指南

要复现本研究中的主要结果，需要熟悉强关联数值方法、量子蒙特卡罗技术以及高性能计算环境。以下是复现的关键步骤和考虑因素：

1. 环境设置:

   • 操作系统: 通常为 Linux 或 HPC 环境。
   • 编程语言与库: 熟悉 C++ 和 Python，因为 TRIQS 库主要基于这些语言。需要安装 TRIQS 及其相关组件（如 triqs_cthyb 求解器）。
   • 高性能计算: 如果模拟较大簇或极低温度，可能需要访问超算集群（如论文中提到的天河-1A），并了解 MPI 等并行计算技术。

2. 模型与参数配置:

   • 哈伯德模型: 定义在三角晶格上，近邻跳跃 $t$，同格点库仑排斥 $U$，化学势 $\mu$。
   • 交互参数:
     • 相互作用强度 $U/t$: 论文主要使用 $U/t = 12$。
     • 温度 $T/t$: 需要在一定范围内扫描，例如从 $T/t = 0.1$ 到 $1$。
     • 掺杂水平 $\delta$: 粒子掺杂侧 (δ > 0)，例如 δ = 0.1, 0.4, 0.5。
     • 化学势 $\mu$: 通过自洽迭代调整，以达到目标掺杂水平。

3. DCA 模拟:

   • 簇几何: 选择三角晶格的簇几何。论文指出 $N_c=4$ 站点簇足以捕捉核心物理，且包含 M 点。
   • 自洽循环:
     • 初始化: 从一个初始自能开始，通常为自由电子自能或上一轮计算的收敛自能。
     • CTHYB 求解: 在每个DCA迭代中，使用 triqs_cthyb 求解器计算簇杂质问题。需要配置蒙特卡罗模拟的步数、热身步数、测量步数和马尔可夫链的长度。确保蒙特卡罗统计量的误差足够小。
     • 杂化函数更新: 根据新的自能更新杂化函数，并执行DCA的自洽条件。
     • 收敛标准: 设定收敛标准，例如自能或粒子数的相对变化小于 $10^{-3}$ 到 $10^{-5}$。
   • 数据输出: 保存虚时格林函数和自能，以及各种物理量（如粒子数、双占据率、自旋关联）。

4. 解析延拓:

   • 数据准备: 提取 CTHYB 结果中的虚时自能 $\Sigma(i\omega_n)$。
   • Padé 近似: 使用 Padé 近似程序（可自行实现或使用现有库）将 $\Sigma(i\omega_n)$ 延拓到实频轴。需要选择合适的 Padé 阶数。
   • 最大熵方法: 作为对比或验证，也可以使用 MaxEnt 方法。MaxEnt 通常需要一个谱函数的先验，并涉及一个优化过程。

5. 物理量后处理:

   • 电阻率: 根据实频自能计算光学电导率 $\sigma(\omega)$，然后取 $\omega \to 0$ 极限得到 $\rho$。
   • 电荷可压缩性: 根据 $\mu$ 的变化计算 $\delta \\langle n \\rangle / \\delta \\mu$。
   • 零频谱权重: 从虚时格林函数外推。
   • 自旋关联: 直接从 CTHYB 测量得到。

复现过程中可能会遇到费米子符号问题导致的噪声，尤其是在低温下。需要仔细调整蒙特卡罗模拟参数和自洽循环的容差，以确保结果的稳定性和可靠性。

所用的软件包及开源 Repo Link

本研究的数值模拟和数据分析高度依赖于以下先进的软件工具和库：

1. TRIQS (Toolbox for Research on Interacting Quantum Systems):

   • 功能: TRIQS 是一个功能强大的开源 C++/Python 库，旨在为研究强关联电子系统提供一套统一的框架和工具。它包含了各种量子杂质模型求解器（如 CTHYB、CT-INT），以及处理格林函数、自能和自洽循环的模块。
   • 链接: TRIQS 的官方网站是其主要的信息来源和下载地址。可以在其 GitHub 页面找到源代码。
     • GitHub Repo: https://github.com/TRIQS (这是一个组织，里面有多个子项目，例如 triqs_cthyb 在 https://github.com/TRIQS/triqs_cthyb)
   • 在本研究中的应用: TRIQS 库被用来实现 DCA 框架，并作为 CTHYB 求解器的后端。它提供了处理 Matsubara 频率格林函数和自能数据结构的基础设施，极大地简化了强关联计算的实现。

2. triqs_cthyb 求解器:

   • 功能: 这是 TRIQS 生态系统中的一个核心量子杂质模型求解器，专门用于解决连续时间混合杂化展开（CTHYB）算法。它基于量子蒙特卡罗方法，能够精确求解在 DCA 或 DMFT 框架下得到的簇杂质问题。
   • 链接: 通常通过 TRIQS 安装包一起安装，也可以在其独立的 GitHub 仓库找到。
     • GitHub Repo: https://github.com/TRIQS/triqs_cthyb
   • 在本研究中的应用: 用作 DCA 框架中的杂质求解器，计算虚时格林函数和自能。

3. Padé 分析延拓方法:

   • 功能: Padé 近似是一种将虚时数据解析延拓到实频轴的数学技术。它不依赖于统计采样，计算速度快，但对原始数据的噪声敏感。
   • 链接: 论文引用了 [70] S. P. Vidberg and J. W. Serene, “Solving the Eliashberg equations by means of n-point Padé approximants,” Journal of Low Temperature Physics 29, 179 (1977)。此方法通常需要自行实现或从现有数值库中调用。没有一个标准化的开源库，但相关实现可以在 NumPy, SciPy 等科学计算库的用户贡献代码中找到。
   • 在本研究中的应用: 主要用于将虚时自能转换为实频自能，以计算实频物理量如光学电导率。

4. 最大熵方法 (MaxEnt):

   • 功能: 最大熵方法是另一种常用的解析延拓技术，通过最大化谱函数的熵并结合贝叶斯推断来处理虚时数据中的噪声，获得最可能的实频谱函数。它比 Padé 近似更鲁棒，但通常计算成本更高且需要先验信息。
   • 链接: 论文引用了 [66, 85, 86]。MaxEnt 通常也需要自行实现或使用专门的库。一些研究组会发布自己的 MaxEnt 实现，但没有一个广泛认可的官方开源库。
   • 在本研究中的应用: 用作 Padé 近似结果的交叉验证，以增强对实频谱函数可靠性的信心。

需要强调的是，尽管 TRIQS 及其求解器是开源的，但本研究中的具体 DCA 脚本、数据分析流程以及有效低能模型推导的代码（特别是与永冈超金属态相关的独特分析部分）并未在论文中直接提供开源链接。研究者需要根据论文中详细描述的方法论和公式，结合 TRIQS 库的功能，自行编写DCA的驱动代码和后处理脚本。

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

关键引用文献

本研究建立在大量前人工作的基础上，并引用了诸多关键文献，这些文献涵盖了强关联物理、非费米液体、永冈铁磁性、范霍夫奇点以及数值计算方法等多个领域。以下列举并简要分析其中一些核心引用：

1. 超金属态与高阶范霍夫奇点 (HOVHS) [8-13]:

   • [8] L. Classen and J. J. Betouras, Annual Review of Condensed Matter Physics 16, 229 (2025): 这篇综述为 HOVHS 及其与平带的联系提供了广阔的背景。
   • [9] H. Isobe and L. Fu, Phys. Rev. Res. 1, 033206 (2019): 首次提出了超金属的概念，并将其与 HOVHS 驱动的奇异金属态联系起来。
   • [10] N. F. Yuan, H. Isobe, and L. Fu, Nature Communications 10, 5769 (2019): 进一步探讨了 HOVHS 的“魔力”及其物理效应。
   • 这些文献为本研究中永冈超金属态的理论框架和特性（如幂律发散的态密度和亚线性电阻率）提供了直接的理论支持，并区分了由单粒子能带精细调整产生的HOVHS与本研究中相互作用驱动的HOVHS。

2. 永冈铁磁性与磁极化子 [33-39]:

   • [33] M. Xu et al., Nature 620, 971 (2023): 这篇文献（共同作者包括本研究作者 M. Xu 和 M. Greiner）揭示了在费米-哈伯德模拟器中挫折和掺杂诱导的磁性。
   • [34] M. Lebrat et al., Nature 629, 317 (2024): 观测到费米-哈伯德模拟器中的永冈极化子，提供了实验基准。
   • [35] M. L. Prichard et al., Nature 629, 323 (2024): 直接成像了动力学阻挫哈伯德系统中的自旋极化子。
   • [36] J. P. Dehollain et al., Nature 579, 528 (2020): 在量子点阵列中观测到永冈铁磁性。
   • 这些实验进展构成了本研究的实验动机和关键基准，证明了在三角晶格中永冈铁磁极化子的存在，并为本研究中永冈驱动的NS态提供了背景。

3. 动态簇近似 (DCA) 与量子杂质求解器 [58-61]:

   • [58] T. Maier et al., Rev. Mod. Phys. 77, 1027 (2005): DCA 综述，详细介绍了其理论框架和应用。
   • [59] M. H. Hettler et al., Phys. Rev. B 61, 12739 (2000): DCA 的早期工作。
   • [60] E. Gull et al., Rev. Mod. Phys. 83, 349 (2011): CTHYB 求解器的综述，介绍了其原理和应用。
   • [61] O. Parcollet et al., Computer Physics Communications 196, 398 (2015): TRIQS 库的介绍，本研究使用的核心软件平台。
   • 这些文献是本研究方法论的基础，为 DCA 框架的选择、CTHYB 求解器的使用以及 TRIQS 库的运用提供了理论和实践指导。

4. 量子蒙特卡罗 (QMC) 模拟 [62-65]:

   • [62] R. Blankenbecler et al., Phys. Rev. D 24, 2278 (1981): DQMC 的开创性工作。
   • [63] J. E. Hirsch, Phys. Rev. B 28, 4059 (1983): 费米子哈伯德模型的离散哈伯德-斯特拉托诺维奇变换。
   • 这些 DQMC 工作的引用用于对 DCA 结果进行基准测试，确保了在特定参数下 DCA 结果的可靠性。

5. 解析延拓方法 [70, 85, 86]:

   • [70] H. J. Vidberg and J. W. Serene, Journal of Low Temperature Physics 29, 179 (1977): Padé 近似方法。
   • [85] G. J. Kraberger et al., Phys. Rev. B 96, 155128 (2017): 最大熵方法的综述。
   • 这些文献为本研究中将虚时自能转换为实频物理量的技术提供了理论和应用基础。

6. 三角晶格哈伯德模型相关研究 [51-57]:

   • [51] D. Pereira and E. J. Mueller, Phys. Rev. B 112, 245120 (2025): 动能磁性在方格和三角格哈伯德模型交叉。
   • [53] J. Fournier et al., SciPost Phys. 17, 072 (2024): 三角格哈伯德模型中两个 T-线性散射率区间。
   • [79] D. Yudin et al., Phys. Rev. Lett. 112, 070403 (2014): 在三角晶格哈伯德模型中范霍夫奇点附近的费米凝聚。
   • 这些文献展示了三角晶格哈伯德模型丰富的物理性质，包括赝能隙、奇异金属和动能磁性，为本研究提供了重要的上下文和对比。

对这项工作局限性的评论

尽管本研究在揭示永冈超金属态方面取得了重要进展，但任何理论和数值工作都存在一定的局限性。以下是对这项工作的几点评论和潜在局限性：

1. 簇尺寸限制:

   • 计算成本: 动态簇近似（DCA）的计算成本随簇尺寸 $N_c$ 的增加呈指数级增长。本研究的计算限制在 $N_c = 3, 4, 6$ 的小簇。尽管论文指出 $N_c=4$ 在捕捉 M 点物理方面是足够的，但更大的簇可能能够捕捉到更长程的关联效应或更精细的相变特征。在某些情况下，小簇可能无法完全重现连续晶格的物理性质，尤其是在接近量子临界点时。
   • 费米子符号问题: 几何阻挫和低温度会加剧费米子符号问题，即使在 DCA 框架下也限制了可达到的簇尺寸和温度范围。这意味着在更低的温度或更大的簇尺寸下，结果的可靠性可能需要更强的计算能力来保证。

2. 解析延拓的不确定性:

   • 不适定问题: 将虚时自能解析延拓到实频自能是一个不适定问题，固有的噪声和方法选择（Padé 近似或最大熵方法）会引入一定的不确定性。尽管本研究使用了两种方法进行交叉验证，但其结果仍然是数值近似，而非精确解。对于更复杂的谱函数结构，解析延拓的精度可能成为限制。

3. 忽略顶点修正:

   • 泡泡近似: 在光学电导率的计算中，本研究使用了泡泡（bubble）近似，即忽略了顶点修正。论文指出，对于三角晶格，较高的连通性和磁阻挫促进了自能的局域性，从而压制了顶点修正的影响。然而，在某些强关联或临界区域，顶点修正可能仍然是重要的，其忽略可能导致定量上的一些偏差。

4. 有效模型的近似性:

   • 微扰理论: 论文推导的有效低能模型是基于二阶微扰理论。这种近似在相互作用强度不是极端强的情况下是有效的。然而，对于极强关联区域，更高阶的微扰项可能变得重要，从而影响有效模型的准确性。

5. 缺乏实验直接证实:

   • 理论预测: 尽管论文强调其发现“可直接通过当前超冷原子实验进行验证”，但截至目前，永冈超金属态的直接实验观测尚未实现。理论预测的验证需要未来的实验工作。

6. 粒子-空穴对称性与超金属的特异性:

   • 本研究明确指出，永冈超金属态仅在粒子掺杂侧出现，而空穴掺杂侧则表现出更传统的费米液体或奇异金属行为。这突出了粒子-空穴不对称性在三角晶格哈伯德模型中的重要性。但这也意味着该态的涌现具有一定的特异性。

这些局限性并不减损本研究的重要性，而是指出了未来研究可以改进和拓展的方向，例如利用更先进的数值方法（如张量网络态或全相互作用蒙特卡罗）处理更大簇和更低温度，或开发更精确的解析延拓技术。

其他你认为必要的补充

永冈超金属态的深远意义与前沿影响

本研究关于永冈超金属态在粒子掺杂三角哈伯德模型中的发现，具有深远的理论意义和潜在的前沿影响，为强关联物理领域注入了新的活力：

1. 非费米液体物理学的新范式:

   • 超越费米液体: 传统费米液体理论在描述许多强关联系统中的异常金属行为时遭遇困境。永冈超金属态的发现，提供了一个清晰的、由相互作用驱动的非费米液体示例，其特征是亚线性电阻率和幂律发散的谱权重及电荷可压缩性。这拓展了我们对非费米液体相的理解。
   • 与奇异金属的区分: 奇异金属（Strange Metal）是凝聚态物理中的一个热点，通常表现出 T-线性电阻率。永冈超金属态的亚线性电阻率 $\rho \sim T^\alpha (\alpha < 1)$ 明确地区分了这两种非费米液体行为，为分类和理解不同的异常金属相提供了新的维度。
   • 高阶范霍夫奇点的内生机制: 许多关于HOVHS的研究往往依赖于精细调整能带结构来使其出现在费米能级附近。本研究的关键突破在于，它展示了HOVHS可以在强关联相互作用驱动下内生地涌现，而非仅仅是单粒子能带结构的巧合。这种“相互作用驱动的能带重构”机制，为理解其他复杂材料中的HOVHS和相关物理现象提供了新的思路。

2. 理解强关联磁性与电荷耦合:

   • 永冈机制的广延性: 永冈铁磁性长期以来是一个重要的理论概念。本研究将其与一种新颖的金属态——超金属态——联系起来，揭示了动能驱动的铁磁关联不仅可以稳定磁性序，还可以显著影响电荷输运和热力学性质。这深化了我们对强关联体系中自旋和电荷自由度之间复杂耦合的理解。
   • 几何阻挫的重要性: 三角晶格的几何阻挫在此过程中扮演了关键角色，通过抑制传统反铁磁序，为永冈铁磁极化子的形成和超金属态的涌现创造了条件。这强调了晶格几何对强关联物理相图的决定性影响。

3. 量子模拟的直接实验可及性:

   • 超冷原子平台的价值: 论文明确指出，本研究的理论预测“可直接通过当前超冷原子实验进行验证”。超冷费米原子在光学晶格中提供的可控性和纯净环境，使其成为检验此类复杂多体理论的理想平台。这不仅为超冷原子实验指明了新的研究方向（如探测亚线性电阻率、电荷可压缩性和谱权重的幂律发散），也为凝聚态物理的“量子模拟”提供了又一个成功范例。
   • 潜在的实验探测方法: 实验可以通过测量扩散系数（与电阻率相关）、电荷可压缩性（通过粒子数涨落）和光谱函数（通过射频谱或角分辨光电子能谱的量子模拟对应）来直接探测超金属态的特征。

4. 对高温超导及其他奇异物态的启示:

   • 关联与临界现象: 许多奇异金属和非费米液体相与高温超导的正常态紧密相关。永冈超金属态的发现，特别是其HOVHS的起源和多体性质，可能为理解高温超导等更复杂的强关联体系中的奇异金属行为和临界现象提供新的线索。
   • 新的材料探索方向: 理论上预测的这种相互作用驱动的HOVHS机制，也可能指导我们去探索具有相似晶格几何和强关联特性的新材料，以期发现类似的超金属或其他新奇量子物态。

综上所述，永冈超金属态的发现不仅是一个有趣的理论现象，更是一个可能改变我们对非费米液体、强关联磁性和电荷输运理解的关键里程碑。它连接了理论建模、数值计算和实验量子模拟的前沿，为未来的研究开辟了广阔的道路。

对未来研究方向的展望

本研究不仅揭示了永冈超金属态的丰富物理，也为未来的研究提出了多个引人入胜的方向：

1. 更全面的相图探索:

   • 超越单层三角晶格: 本研究关注的是二维单层三角晶格。未来可以扩展到多层三角晶格或更复杂的几何结构（如 Kagome 晶格），探究几何阻挫与相互作用如何共同塑造更多样化的相图。
   • 更宽泛的参数空间: 虽然本文已探索了不同掺杂和相互作用强度，但系统性地研究更广泛的 $U/t$、温度 $T/t$ 范围以及其他次近邻跳跃 $t'$ 等参数的影响，可以揭示相边界的精细结构和可能的其他新奇相。
   • 非平衡态动力学: 研究永冈超金属态在非平衡条件下的动力学行为，如响应外部电场或磁场的弛豫过程，可以提供对其动力学特征更深入的理解。

2. 精细化微观起源与理论建模:

   • 高阶微扰理论: 论文中的有效低能模型是基于二阶微扰理论。在极强关联区域，高阶微扰项可能变得重要。发展更高阶的微扰理论或更精确的有效模型推导，可以提高对HOVHS起源的定量描述精度。
   • 顶点修正的作用: 尽管本文对三角晶格中忽略顶点修正的合理性进行了讨论，但在某些参数区域或对更精确的定量结果要求下，纳入顶点修正的计算是必要的。这可能需要更复杂的 DCA 扩展或利用其他多体方法。
   • 自旋极化态的深入研究: 论文提到未来工作将扩展到自旋极化区域 [30, 38]。这将有助于澄清集体自旋激发如何影响电荷载流子并重塑输运图景。在外部磁场下，系统的铁磁关联和HOVHS如何演化是一个值得深入探讨的问题。

3. 实验验证与量子模拟的进一步推动:

   • 直接实验观测: 最紧迫的方向是实现永冈超金属态的直接实验观测。超冷原子实验应专注于探测亚线性电阻率、电荷可压缩性的幂律发散以及零频谱权重的奇异行为。这可能需要开发新的实验技术，例如精确测量光学晶格中费米子的扩散系数或热力学响应。
   • 多维度探测: 除了输运和热力学测量，还需考虑其他探测手段，如非弹性光散射、中子散射或射频谱，以获取自旋、电荷和轨道自由度的信息，从而全面刻画超金属态的微观特征。
   • 与真实材料的联系: 尽管本文主要针对模型系统，但将其理论预测与实际具有三角晶格或类似几何阻挫特征的材料（如某些过渡金属氧化物或层状材料）联系起来，将是未来研究的重要方向。

4. 超金属态与其他量子物态的关联:

   • 与超导性的竞争与合作: HOVHS 通常与超导性密切相关。在超金属态中，HOVHS可能与超导涨落共存甚至竞争。研究掺杂范围内是否存在超导序，以及其与永冈超金属态之间的相互作用，是值得深入探索的课题。
   • 量子临界性: 超金属态可能位于量子相变附近，其非费米液体行为可能与量子临界点相关。研究超金属态附近的量子临界现象及其临界指数，可以揭示更普遍的强关联物理规律。
   • 量子自旋液体 (Quantum Spin Liquid, QSL): 论文提及未来研究将探索量子自旋液体相中的输运特征 [84]。这可能涉及更强阻挫或更高阶相互作用的模型，将超金属态的理解扩展到更广阔的量子物态领域。

这些未来研究方向的探索，将进一步丰富我们对强关联电子系统、非费米液体和量子相变的理解，并有望指导新材料的设计和量子技术的进步。