拓扑重费米子模型中的关联莫特半金属：哈伯德算符方法与行列式量子蒙特卡洛的深度对标

来源论文: https://arxiv.org/abs/2606.20466v1 生成时间: Jun 20, 2026 00:43

0. 执行摘要

魔角双层石墨烯（Magic-Angle Twisted Bilayer Graphene, MATBG）作为强关联物理与拓扑能带理论交汇的前沿阵地，其核心物理特征在于电中性点附近的半金属导电行为与局部局域磁矩（Localized Moments）的共存。传统的格点哈伯德模型由于拓扑能带的“Wannier阻碍”（Wannier Obstruction），无法直接建立局域化的格点描述。最近提出的**拓扑重费米子模型（Topological Heavy Fermion Model, THFM）**通过将系统解耦为强关联的局域 $f$ 轨道与无相互作用的拓扑瞬态 $c$ 轨道，提供了一个优雅且极具启发性的极简理论框架。

然而，现有的理论解析方法（如 Hubbard-I 近似）在处理局域磁矩与瞬态电子之间的弱杂化耦合时，往往忽略了空间非局域关联，导致在局部磁矩区域给出错误的谱学特征和过高的能隙。动力学平均场理论（DMFT）虽然能够自洽处理动态自能，但在处理低维系统和空间涨落时也存在局限。

本工作发展了哈伯德算符（Hubbard Operator, HO）方法（亦称复合算符方法，Composite Operator Method, COM），并将其从传统的单位置（Single-site）限制推广至**双位置（Dimer）限制。为了对该近似方法的有效性进行无偏（Unbiased）检验，作者在格点正则化的 1D THFM 体系中，利用行列式量子蒙特卡洛（Determinant Quantum Monte Carlo, DQMC）**进行了精确的数值基准测试。结果表明：

1. 传统的本地近似（如 Hubbard-I）无法正确捕捉局域与瞬态自由度之间的非局域关联，高估了单粒子能隙。
2. 引入 Roth 解耦自洽方案的单位置 COM 方法在静态关联函数的描述上与 DQMC 符合极好，但在动力学谱函数上仍有偏差。
3. 基于 Dimer 基础的哈伯德算符方法能够精确地重现 DQMC 的虚时格林函数和实频谱学特征，成功描述了 $\Gamma$ 点能隙闭合、高对称点能隙打开的**莫特半金属（Mott Semi-metal）**行为，证明了空间非局域关联在重费米子物理中的决定性作用。

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

1.1 核心科学问题：Wannier 阻碍与拓扑重费米子模型（THFM）

在强关联电子系统中，电荷的局域化通常通过在实空间构建 Wannier 函数并引入 local Coulomb $U$ 来描述。然而，MATBG 的平带具有非零的微观拓扑性质（如 Chern 数或非平凡的 Berry 曲率分布），这在数学上阻碍了对称性守恒的实空间 Wannier 函数的构建。这种“Wannier 阻碍”使得传统的哈伯德模型不适用于 MATBG。

为了克服这一困难，拓扑重费米子模型（THFM）应运而生。THFM 将 MATBG 的低能物理映射为两部分：

• $f$ 电子（重费米子部分）：局域在实空间 Wigner 晶格的蜂窝格点上，具有完全局域化的特征。它们不直接在格点间跃迁，但具有强烈的局域库仑排斥 $U$。
• $c_1, c_2$ 电子（轻费米子/瞬态部分）：具有狄拉克色散的拓扑导电电子。$c_1$ 电子通过杂化项 $\gamma$ 与 $f$ 电子耦合，而 $c_2$ 电子则不与 $f$ 直接耦合，但通过与 $c_1$ 的跃迁建立联系。

该模型的哈密顿量在格点正则化后可表示为：

其中 $U$ 是 $f$ 轨道的排斥能，$\gamma$ 是 $f$ 与 $c_1$ 的杂化强度，$t_{ij}$ 为轻费米子间的跃迁矩阵元。在电中性点（Half-filling），MATBG 表现为具有 disordered 局部磁矩的特殊半金属。我们的核心任务就是如何精确计算该强关联多体系统的格林函数和谱函数。

1.2 理论基础：哈伯德算符方法（HO）与运动方程（EOM）技术

哈伯德算符方法是一种解析量子集群方法。其核心思想是，对于一个包含强局域关联的哈密顿量，将其分为局域部分 $H_{\text{loc}, i}$ 和非局域部分 $H_{t, ij}$：

1. 构建局域算符基底 $\Psi_i$： 由于局域希尔伯特空间 $H_i$ 是有限维的，我们可以选择一组有限的算符集合 $\Psi_i = (\psi_{1,i}, \psi_{2,i}, \dots)^T$，使其在与局域哈密顿量 $H_{\text{loc},i}$ 的对易子运算下闭合（即运动方程）：

   $$\mathbf{J}_{\text{loc}, i} = [\Psi_i, H_{\text{loc},i}] = E_{\text{loc},i}\Psi_i$$

   其中 $E_{\text{loc},i}$ 是描述局域激发能的矩阵。

2. 投影非局域项： 非局域项 $H_t$ 的引入会使运动方程不再闭合，产生格点外的算符 $\delta \mathbf{J}_i$：

   $$[\Psi_i, H_t] = \delta \mathbf{J}_i$$

   HO 方法的关键步骤是将这些高阶算符 $\delta \mathbf{J}_i$ 线性投影回我们选择的局域算符基底上：

   $$\delta \mathbf{J}_i \approx \sum_j E_{t, ij} \Psi_j$$

   结合局域与非局域贡献，完整的运动方程近似为：

   $$\mathbf{J}_i \approx \sum_j (\delta_{ij}E_{\text{loc},i} + E_{t, ij})\Psi_j$$

3. 自洽求解能谱与谱权重： 为了求解系统的格林函数 $G_{ij}(\omega) = \langle \langle \Psi_i; \Psi_j^{\dagger} \rangle \rangle_{\omega}$，我们需要定义两个核心矩阵：

   • 能量矩阵 $E$：定义为 $E = M I^{-1}$
   • 谱权重矩阵 $I_{ij}$：$I_{ij} = \langle \{\Psi_i, \Psi_j^{\dagger}\} \rangle$
   • 动力学激发矩阵 $M_{ij}$：$M_{ij} = \langle \{\mathbf{J}_i, \Psi_j^{\dagger}\} \rangle$

   其中 $\{\dots\}$ 代表费米子对易子，矩阵元中的期望值在等时（Equal-time）自洽评估。最终，频域格林函数可写为：

   $$G_{ij}(\omega) = \left[ (\omega + i\varepsilon - E)^{-1} I \right]_{ij}$$

   这里的 $\varepsilon$ 是唯象的准粒子展宽，而非规范的 $I$ 矩阵则是该系统偏离无相互作用自由费米子行为的关键指标。

1.3 技术难点一：Roth 解耦自洽方案

自洽计算矩阵 $M$ 和 $I$ 时，不可避免地会引入高阶静态关联函数。例如单位置的复合算符方法（COM）中，会产生如下非局域关联参数 $p_{\alpha\alpha'}$：

该参数包含电荷密度-电荷密度、自旋-自旋以及库珀对-库珀对的非局域关联。传统的 Hubbard-I 近似简单地将其设为高温极限下的解耦值（在半满时 $p = 0.25$），从而完全忽略了空间关联。COM 通过引入 Roth 解耦方案，将这些高阶关联函数通过格林函数的极点结构进行自洽求解，使其成为温度 $T$ 和相互作用 $U$ 的自洽函数。这是超越本地近似的第一步。

1.4 技术难点二：双位置算符基（Dimer Basis）的构建与维数灾难

尽管 COM 引入了 Roth 解耦，但其算符基依然限制在单格点内，无法精确描述强局域 $f$ 电子与近邻轻电子的短程量子干涉（类 Kondo 杂化）。为了解决这一问题，本工作发展了 Dimer 算符基。我们将相邻的两个格点 (1, 2) 视为一个整体 Dimer，其局域哈密顿量精确包含了它们之间的跳跃：

在 Dimer 基础下，为了使运动方程在原子极限下严格闭合，其算符基底 $\Psi_d$ 的维数急剧膨胀，如公式 (18) 所示，其包含了多种复杂的复合激发态（例如带自旋翻转的辅助算符 $\xi_{1\downarrow}S^{-}_2$ 等共 24 个算符）。这导致 $I$ 和 $M$ 矩阵的大小变为 $24 \times 24$。手工推导这些高阶矩阵元的 Roth 解耦极度困难。为此，作者采用了一种创新的混合方案：利用精确的有限温 DQMC 数值模拟直接计算等时关联函数，进而构建 $I$ 和 $M$ 矩阵，避免了繁琐的解析解耦，同时保证了多体关联的无偏性。

2. 关键 Benchmark 体系、计算所得数据与谱学分析

为了检验不同解析近似（Hubbard-I, COM, Dimer HO）的精确度，本工作采用一维格点正则化拓扑重费米子模型进行系统性的 DQMC 对标。

2.1 静态关联函数 $p$ 的基准测试

首先考察刻画 $f$ 轨道与传导 $c_1$ 轨道自旋及电荷耦合的静态关联函数。在局部磁矩区（$W_2 \gg U, \gamma$），我们将计算得到的自旋关联 $\langle S^{+}_{fi} S^{-}_{c1i} \rangle$ 和电荷关联 $\langle n_{fi\sigma} n_{c1i\sigma} \rangle$ 随跳跃常数 $t$ 的变化绘制于下图（对应论文 Fig. 2）：

物理分析：

• Hubbard-I 的完全失效：Hubbard-I 假设 $p = 0.25$，认为 $f$ 轨道磁矩与 $c_1$ 传导电子完全解耦（红线）。在 $t \le 3.0$ 的物理感兴趣区域，这与 DQMC 的实际物理行为存在巨大偏差。
• COM 的优异表现：得益于 Roth 解耦方案，单位置 COM（虚线）在静态关联上与 DQMC 实现了高达 $10\%$ 以内的定量符合，成功捕捉到了杂化诱导的局部反铁磁关联和电荷屏蔽效应。

2.2 虚时格林函数 $\mathcal{G}(\tau, k = \pi/2)$ 的动态分析

虚时格林函数在 $\tau = \beta/2$ 的行为提供了谱能隙的直接测度：$\mathcal{G}(\beta/2) \sim e^{-\beta \Delta / 2}$。在极高对称点 $k = \pi/2$（谱能隙最明显的区域），作者给出了详细的虚时演化（对应论文 Fig. 3）：

• 参数组 A ($U=1.5, t=3.0$)：
  • 在 $\tau = 10.0$ 处：
    • $\mathcal{G}_{\text{DQMC}} \approx 1.2 \times 10^{-2}$
    • $\mathcal{G}_{\text{Dimer}} \approx 1.0 \times 10^{-2}$（高度契合）
    • $\mathcal{G}_{\text{COM}} \approx 1.5 \times 10^{-3}$（低估了一个数量级）
• 参数组 B ($U=3.0, t=5.0$)：
  • 在 $\tau = 10.0$ 处：
    • $\mathcal{G}_{\text{DQMC}} \approx 8.0 \times 10^{-3}$
    • $\mathcal{G}_{\text{Dimer}} \approx 6.5 \times 10^{-3}$（极其接近）
    • $\mathcal{G}_{\text{COM}} \approx 2.0 \times 10^{-4}$（严重偏离）

物理机制解释： COM 近似由于算符基限制在单格点内，无法容纳双位置上的协同相干散射。它高估了局域磁矩的自囚禁能力，从而给出了一个过大的、近乎原子极限的莫特能隙（表现为虚时格林函数衰减过快）。而 Dimer HO 精确地引入了短程跃迁杂化，释放了局域电子的动能，极大地压低了能隙，从而与 DQMC 的精确趋势完美重合。

2.3 实频谱函数（Spectral Function）与莫特半金属相的演线

论文 Fig. 5 展示了利用随机解析延拓后的 DQMC 谱函数与三种近似方法的对比：

1. Hubbard-I 近似：给出完全对称、刚性的分裂能带，能隙大小恒为 $U$。在 $\Gamma$ 点（$k=0$），能带表现为平凡的线性色散交叉，忽略了相互作用对狄拉克锥的重整化。
2. COM (单位置) 方法：虽然引入了自旋涨落重整化，但在偏离 $\Gamma$ 点时仍然表现出过宽的莫特能隙，且能带极其尖锐，缺乏寿命衰减特征。
3. Dimer HO 与 DQMC 结果对比：
   • $\Gamma$ 点 ($k=0$)：能隙闭合，展现出由轻费米子 $c_2$ 轨道主导的无能隙结构。
   • $k = \pi/2$ 点：由于 $f$ 轨道强关联，展现出显著重整化后的关联能隙。DQMC 显示其真实 gap 约为 $U/3$（由于 $c_2$ 轨道的间接屏蔽和短程反铁磁关联），Dimer HO 准确地捕捉到了这一显著的能隙收缩（能隙从 $1.5$ 收缩至 $\approx 0.5$）。
   • 这完美证实了该系统处于 莫特半金属（Mott Semi-metal） 相。

3. 代码实现细节、复现指南与开源软件包

本工作中的所有基准数据均基于著名的开源量子蒙特卡洛软件库 ALF (Algorithms for Lattice Fermions) 进行计算。以下是针对关联物质科研人员的复现指南及核心实现逻辑：

3.1 ALF 软件包及开源 Repo 链接

• ALF 官方仓库：https://github.com/ALF-Project/ALF
• ALF 官方文档：https://alf.physik.uni-wuerzburg.de/
• 本文所用的具体正则化拓扑重费米子模型（THFM）可通过在 ALF 中自定义哈密顿量模组（Hamiltonian Module）实现。

3.2 DQMC 复现关键参数配置

为了在电中性点实现无费米子符号问题的精确模拟，需在自定义输入文件 parameters 中进行如下设置：

! ALF 关键参数配置示例
&VAR_LATTICE
L1 = 100         ! 沿x方向的格点数 (热力学极限近似)
L2 = 1           ! 1D 链状拓扑体系
Lattice_Type = "Chain"
/
&VAR_MODEL
Model_Name = "THFM_1D"
N_Fermions = 2    ! 考虑 SU(2) 自旋对称，Nf = 2
U_Coupling = 1.5  ! 局域库仑排斥 U
Gamma_Hybrid = 1.0! f-c1 杂化强度
t_hopping = 3.0   ! c1-c2 跳跃常数
/
&VAR_QMC
Beta = 100.0      ! 极低温极限，保证基态物理，beta = 1/T
dtau = 0.005      ! Trotter 虚时切片大小，控制系统误差
Nsweep = 10000    ! 蒙特卡洛扫射步数
Warmup = 2000     ! 预热步数
/

3.3 核心计算算法逻辑：从虚时到实频

1. 虚时格林函数测量： ALF 通过辅助场量子蒙特卡洛（AFQMC）在虚时轴 $\tau_k = k \Delta \tau$ 上测量 $f$ 电子和 $c$ 电子的混合格林函数：

   $$\mathcal{G}_{\alpha\beta}(\tau, k) = -\langle T_{\tau} c_{\alpha}(k, \tau) c^{\dagger}_{\beta}(k, 0) \rangle$$

2. 随机解析延拓 (Stochastic Analytic Continuation, SAC)： 由于虚时向实频的反演极其不稳定，本工作采用 ALF 自带的 SAC 算法（基于 Sandvik 法或最大熵方法 Maximum Entropy Method, MEM）求解积分方程：

   $$\mathcal{G}(\tau) = -\int d\omega \frac{e^{-\tau\omega}}{1 + e^{-\beta\omega}} A(\omega)$$

   通过设定合理的默认模型（Default Model）并对虚时协方差矩阵进行精确加权，最终得到高质量的实频谱函数 $A(k, \omega)$。

4. 关键引用文献与局限性评论

4.1 关键引用文献

1. THFM 模型的提出：
   • Song, Z.-D. and Bernevig, B. A. Magic-angle twisted bilayer graphene as a topological heavy fermion problem. Phys. Rev. Lett. 129, 047601 (2022).
     (奠定了本文的模型基础，首次完成了 MATBG 的拓扑-重费米子映射。)
2. 复合算符方法（COM）的理论框架：
   • Mancini, F. and Avella, A. The hubbard model within the equations of motion approach. Advances in Physics 53, 537 (2004).
     (提供了算符运动方程和投影方法的技术核心。)
3. Roth 解耦方案：
   • Roth, L. M. Electron correlation in narrow energy bands. I. The two-pole approximation in a narrow s band. Phys. Rev. 184, 451 (1969).
     (提供了自洽计算静态关联函数、超越简陋局域近似的关键解耦思想。)

4.2 对本工作局限性的客观评论

虽然本工作成功证明了 Dimer HO 方法在描述拓扑重费米子模型的非局域关联和莫特半金属相方面的巨大优势，但从严苛的量子化学与关联物理视角来看，它仍存在以下无法忽视的局限性：

1. 维数灾难导致多轨道扩展困难（The Scalability Bottleneck）： 本文的基准测试仅限于最简单的 $N_f = 2$（SU(2) 自旋物理）的一维模型。然而，真实的 MATBG 体系具有 Spin-Valley 简并，展现出 SU(8) 高度对称性（公式 19）。若要在 Dimer 算符基下构建 SU(4) 或 SU(8) 的算符闭合集，其基底维数将呈指数级爆炸，矩阵大小将超越解析推导和数值自洽的处理极限。这意味着，目前该方法无法直接、完整地定量预测真实 MATBG 2D 样品的物理参数。
2. 一维基准测试的局限性（Dimensionality Gap）： 为了获取无符号问题的精确 DQMC 参考数据，作者妥协选择了一维链状拓扑正则化体系。一维体系由于强烈的量子涨落，根据 Mermin-Wagner 定理不具备真正的长程磁有序，这使得 paramagnetic 状态下的 HO 方法异常成功。但在真实的二维 MATBG 中，空间几何阻挫、非平凡的 Berry 磁通分布以及可能存在的自发对称性破缺（如铁磁相或 Chern 绝缘体相）会使非局域关联函数极度复杂，1D 基准测试的成功并不能完全保证 2D 复杂相图下的同等精度。
3. 准粒子寿命描述的唯象性（Phenomenological Broadening）： 在有限的算符截断阶数下，HO 方法本质上描述的是具有无限寿命的相干准粒子（Green’s 极点无虚部）。为了得到平滑的谱函数，作者在公式 (10) 中人工引入了唯象的拓宽因子 $\varepsilon = 0.025$。这与自洽求解自能虚部 $\text{Im}\Sigma(\omega)$ 的 DMFT 相比，在处理非相干散射、热力学阻尼和电阻率等输运性质的精细计算时，物理严谨性稍显逊色。

5. 补充物理分析：莫特半金属相的演变与几何屏蔽效应

为了更直观地理解非局域关联如何雕琢“莫特半金属”这一奇特物态，我们需要深入解构 $f$ 轨道与 $c_1, c_2$ 导电通道之间的协同屏蔽机制。

5.1 杂化诱导的“几何屏蔽”（Geometric Screening）

在传统的单点安德森晶格模型（Anderson Lattice Model）中，局域磁矩与传导电子直接发生近邻杂化。当温度低于 Kondo 温度 $T_K$ 时，局域磁矩被完全屏蔽，系统在费米面处打开能隙，转变为强关联的 Kondo 绝缘体。

但在 THFM 中，空间几何结构带来了全然不同的物理：

• $f$ 电子只直接与 $c_1$ 杂化，而 $c_2$ 则在格点间自由穿梭。
• 当 $U$ 较大且非局域关联被正确引入时（如 Dimer HO 和 DQMC 所示）：$f$ 电子的电荷涨落被极大地抑制。通过 Dimer 单元内部的相干跃迁，$f$ 电子与近邻 $c_1$ 电子建立起局域的反铁磁自旋关联（$\langle S^{+}_{fi} S^{-}_{c1i} \rangle < 0$）。
• 然而，这种局部反铁磁纠缠并不是孤立的。由于 $c_1$ 与具有狄拉克锥色散的 $c_2$ 强力混合，狄拉克电子的量子几何（Quantum Geometry）和浓缩在 $\Gamma$ 点附近的 Berry 曲率在动量空间中对这种自旋极化进行了选择性屏蔽。
• 结果是，在离开 $\Gamma$ 点的区域（如 $k \approx \pi/2$），电子展现出强烈的局域化特征，打开了由于 Coulomb $U$ 和自旋涨落共同驱动的莫特关联能隙。而在 $\Gamma$ 点附近，由于拓扑保护和 $c_2$ 轨道的非局域散射，关联自能迅速衰减，能隙闭合，形成了独特的莫特半金属状态。

5.2 准粒子权重 $Z_{QP}$ 的约束分析

辅助材料（Appendix C.2）中提出了一项关于准粒子峰强度的精妙物理判定。通过虚时格林函数在 $\beta/2$ 的精确值，我们可以严格界定准粒子相干权重 $Z_{QP}$ 的上限：

在论文 Fig. 12 中，作者展示了不同跃迁常数 $t$ 下的 $\mathcal{G}_{\text{DQMC}}(\beta/2)$ 数值：

• 对于 $U=3.0$（强关联区）：随着 $t$ 从 $1.0$ 增加到 $5.0$，$\mathcal{G}(\beta/2)$ 迅速从 $\approx 3 \times 10^{-2}$ 坠落至 $\approx 1 \times 10^{-3}$ 以下。
• 这意味着准粒子相干权重 $Z_{QP}$ 被压制得极低，表明系统在强关联和强跃迁混杂下，传统的相干重费米子液体准粒子图像已经趋于瓦解，非相干的莫特散射占据了主导地位。

5.3 展望：超越 Dimer 走向 Quadrimer 基底

为了在未来精确解决真实 MATBG 的 SU(4)/SU(8) 莫特绝缘相和超导配对机制，HO 方法必须寻求新的数学妥协。一个可行的方向是放弃完全实空间的算符展开，转而采用混合动量-实空间算符基，或者利用张量网络（Tensor Networks）对 Dimer 乃至 Quadrimer（四位置）的算符稀疏矩阵进行自适应截断。这将是一项极具挑战性但也极具应用前景的理论计算化学/物理工程。