三角格子M点莫尔材料中的隐抗铁磁性、持久谷涨落与U(6)交叉行为：行列式量子蒙特卡洛深度解析

来源论文: https://arxiv.org/abs/2606.12530v1 生成时间: Jun 13, 2026 10:32

0. 执行摘要

莫尔（Moiré）超晶格材料为探索强关联电子态提供了一个前所未有的高度可调平台。绝大多数已知的莫尔系统（如转角双层石墨烯、过渡金属二硫族化合物TMD等）的低能物理主要由布里渊区（BZ）的 $\Gamma$ 点或 $K/K'$ 点控制。然而，最近理论提出了一类全新的莫尔材料平台——通过将低能态位于三个 $M$ 点的二维单层材料进行转角堆叠，可以实现一类独特的三谷哈伯德模型（Three-Valley Hubbard Model）。

这类材料（以转角 AA 堆叠的 $\text{SnSe}_2$ 为代表，即 AA $t-\text{SnSe}_2$）具有以下极其优异且反常的物理特性：

1. 谷选择性的准一维跳跃（Quasi-1D Hopping）：受体系中涌现的动量空间非共形（Nonsymmorphic）对称性 $\tilde{M}_z$ 的约束，电子在特定的谷内仅能沿着特定的方向进行强跳跃，形成了空间上相互交织但动量上解耦的准一维链状结构。
2. 高度对称的交互作用：在无需精细调节的情况下，由于三个谷的 Wannier 轨道在实空间高度重合，且谷间洪特耦合（Hund’s coupling）在小转角下受到极大抑制，其库仑排斥相互作用表现出近乎完美的 $U(6)$ 对称性（结合了 2 个自旋自由度与 3 个谷自由度）。
3. 天然无符号问题（Sign-Problem-Free）：尽管体系处于几何受挫的三角格子上，但在半填充（每个莫尔胞 3 个电子，即 $\nu=3$）时，由于动量空间跳跃的特殊二分（Bipartite）结构，该系统天然地允许**无符号问题的行列式量子蒙特卡洛（DQMC）**数值模拟。

本文将对该系统的核心理论框架、DQMC 算法实现、基准测试（Benchmark）、关键关联物性（如“隐”奈尔反铁磁、持久谷涨落、自旋-电荷-谷相图的 $U(6)$ 交叉行为）以及 Parton 平均场理论进行全方位的技术性深度解析。

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

1.1 体系哈密顿量与对称性分析

考虑一个定义在三角格子上的三轨道（代表三个谷 $\eta \in \{0, 1, 2\}$）紧束缚哈密顿量。电子的自旋记为 $s \in \{\uparrow, \downarrow\}$。动能项 $H_t$ 描述了电子在同谷内的各向异性跳跃：

其中 $d^{\dagger}_{\mathbf{R},\eta,s}$ 在位置 $\mathbf{R}$ 处创建了一个自旋为 $s$、属于谷 $\eta$ 的电子。其独特的跳跃结构（如图 1 所示）由空间群对称性决定。在谷 $\eta$ 中，最强的跳跃项 $t$ 沿着主轴方向 $C_{3z}^{\eta} \mathbf{a}_{M_2}$，而垂直于该方向的次近邻跳跃项 $t_{\perp}$ 较弱。最关键的物理在于，最邻近链间跳跃项（通常记作 $t'$）由于非共形对称性 $\tilde{M}_z$ 的约束被严格禁止（即 $t' = 0$）。因此，整个动能项在实空间上分裂成解耦的、相互穿插的矩形网格（Rectangular Lattices）。

图 1 所示的实空间和动量空间几何结构表明：每个谷内的有效布里渊区并非三角形，而是矩形。这一特征为规避三角形格子的几何受挫、进而产生长程磁有序奠定了几何基础。

    [实空间跳跃示意图 (t' = 0)]              [矩形 Brillouin 区 (各谷解耦)]
          ● (η=0)                                     M2 ─── M3
         / \                                          │     │
     t  /   \ t_⊥                                     │  Γ  │
       /     \                                        │     │
      ● ───-─ ●                                       M1 ─── K
    (η=1)    (η=2) 

相互作用项由局部各向异性的哈伯德排斥以及最邻近排斥 $V$ 组成：

其中 $\hat{n}_{\mathbf{R},\eta} = \sum_{s} d^{\dagger}_{\mathbf{R},\eta,s} d_{\mathbf{R},\eta,s}$ 且 $\hat{n}_{\mathbf{R}} = \sum_{\eta} \hat{n}_{\mathbf{R},\eta}$。参数 $\alpha$ 刻画了谷间的各向异性：

• 当 $\alpha = 1$ 时，相互作用表现出完美的 $U(6)$ 对称性。
• 当 $\alpha < 1$ 时，对称性降低为 $[U(2)]^{\otimes 3}$，此时物理上更倾向于每个谷占据一个电子的态。

1.2 费米子符号问题的消除机制

在量子蒙特卡洛（QMC）模拟中，符号问题是阻碍强关联费米子计算的核心障碍。对于一般的三角格子哈伯德模型，在非半填充或存在几何受挫跳跃时，蒙特卡洛权重会出现负值或复数值，导致方差呈指数级发散。

本模型的妙处在于，当满足以下条件时，系统天然不存在符号问题：

1. 体系处于半填充状态（$\nu = 3$）。
2. 引入一个非平庸的反共轭粒子-空穴对称性（Particle-Hole Symmetry, PHS）算符 $\hat{\mathcal{P}}$，其在基矢上的作用定义为：

其中，符号函数 $\text{sgn}(\mathbf{R}, \eta) = (-1)^m$（对于格点 $\mathbf{R} = n \mathbf{a}_{M1} + m \mathbf{a}_{M2}$）在各个矩形链方向交替变换。由于 $t' = 0$，动能项在每个谷内都是严格二分的（Bipartite），从而保证了动能项在该 PHS 变换下保持不变。

在 DQMC 的配分函数展开中，自旋向上（$\uparrow$）和自旋向下（$\downarrow$）的费米子行列式权重 $W_{\uparrow}$ 和 $W_{\downarrow}$ 通过该反共轭 PHS 建立了如下共轭关系：

因此，总权重 $W(\{\phi\}) = W_{\uparrow} W_{\downarrow} = |W_{\uparrow}(\{\phi\})|^2 \ge 0$ 恒大于等于零。这在数学上严格证明了该模型在半填充下的 DQMC 模拟完全免于符号问题。

1.3 相互作用项的 Hubbard-Stratonovich (HS) 分解

为了实施 DQMC，必须将非二次型的费米子相互作用算符转化为与辅助场耦合的二次型。由于模型中存在局部和非局部的电荷密度排斥，有两种分解方案：键分解（Bond Decomposition）与三角形分解（Triangle Decomposition）。

1.3.1 谷空间基矢变换

首先引入谷空间的对角化算符：

利用这些算符，相互作用能写成如下平方和的形式（省去常数项）：

1.3.2 三角形分解的优势

对于最邻近排斥 $V$，由于三角格子的几何特性，每个键属于两个共享的三角形。采用三角形分解（如图 2 所示）可以将最邻近相互作用写为：

其中 $g = V/2$。这种分解在数学上比直接的键分解提供了更宽的无符号问题参数区间：

相比之下，键分解的要求为 $V/U \le (2\alpha + 1)/18$。因此，后续所有的物理计算均采用三角形 HS 分解方案。

2. 关键 Benchmark 体系与数值物性计算

为了确保所开发的无符号 DQMC 算法的绝对正确性，论文在小尺寸体系上进行了详尽的严格对角化（ED）和随机级数展开（SSE）基准测试。

2.1 严格对角化 (ED) 基准测试

由于包含 3 个谷和 2 个自旋通道，一个单一单元胞的局部希尔伯特空间维度高达 $2^6 = 64$。即使是极小的 $2 \times 2$ 晶格，其总空间维度也达到了 $2^{24} \approx 1.67 \times 10^7$。利用体系的粒子数守恒、自旋守恒以及空间群对称性（$C_{3z}, C_{2x}$ 以及平移对称性），可以将总矩阵块划分为独立的等价类（Equivalence Classes）。

图 2 展示了在 $L=2$ 晶格上，ED 与 DQMC（在不同的虚时步长 $\Delta\tau$ 下）关于能量 $E/t$ 和电荷涨落 $\text{var}(\hat{N})$ 的对比：

    [ED 与 DQMC 能量 Benchmark 曲线对比示意]
    E/t
    12 ┼───────────────────* (ED)
    10 ┼                  *  * (DQMC, Δτ=0.1)
     8 ┼                 *    * (DQMC, Δτ=0.05)
     6 ┼                *      * (DQMC, Extrapolated Δτ→0)
       └─────────────────────────────
       10^{-1}         10^0          10^1  βt

当 $\Delta\tau \to 0$ 时，DQMC 的外推结果与 ED 在所有温度区间（整个 $\beta t$ 轴）均取得了在数值误差（$10^{-5}$ 数量级）范围内的完美符合。

2.2 隐反铁磁性与长程序外推

体系的一个核心物理发现是隐奈尔反铁磁性（Hidden Néel Antiferromagnetism）。尽管实空间格子是三角形的（这通常会由于受挫而导致非共线 120 度磁结构），但由于 $t' = 0$，自旋关联在每个谷内具有矩形二分性。其磁结构因子 $S^{(\eta)}(\mathbf{q})$ 在矩形布里渊区的边缘交点 $\mathbf{Q}^{(\eta)}$ 处展现出极强的峰值：

为了确定在零温极限下系统是否存在长程反铁磁序（LRO），需进行热力学极限（$L \to \infty$）外推。使用自旋波理论给出的有限尺寸标度公式：

其中 $m_{\infty}$ 为自旋序参量。外推结果（如图 3 所示）显示：

• 在各向异性显著（$\alpha \le 0.9$）且耦合较强时，$m_{\infty}^2 > 0$，表明系统在基态具有真正的长程奈尔反铁磁序。
• 在接近 $U(6)$ 对称点（$\alpha \to 1$）且处于中等强度耦合（$U/W \approx 2.4$）时，$m_{\infty}^2 \le 0$（在物理上对应无长程序），此时长程磁有序被强烈的谷涨落和低能电荷激发完全抑制，形成了宽广的非均匀金属/交叉区（Crossover Regime）。

    [磁性长程序标度外推图]
    S(Q)/L^2
    0.4 ┼                                  ● α = 0.5 (强 AFM)
    0.3 ┼                            ▲     ▲ α = 0.9 (弱 AFM)
    0.2 ┼                      ■     ■     ■ α = 1.0 (无 AFM, m_∞=0)
    0.1 ┼                ◆     ◆     
    0.0 ┼─────────●─────┴─────┴─────┴─────────────────
        0 (L→∞)       1/12        1/8        1/4      1/L

3. 代码实现细节、复现指南与开源工具

3.1 模拟所用开源软件包

本工作中的行列式量子蒙特卡洛计算完全基于开源的量子多体计算框架 ALF (Algorithms for Lattice Fermions)。这是一个高度模块化的、基于 Fortran 和 Python 的通用辅助场 QMC 软件包。

• ALF 官方主页/源码仓库：https://alf.physik.uni-wuerzburg.de
• 解析延拓工具：包内自带的基于经典最大熵方法（Maximum Entropy Method）的 MaxEnt 模块，用于处理虚时格林函数 $G(\mathbf{k}, \tau)$ 并提取实频谱函数 $A(\mathbf{k}, \omega)$。

3.2 核心 Hamiltonian 模块的代码构建逻辑

要在 ALF 中实现本模型，必须在自定哈密顿量子程序（通常位于 Hamiltonian_main.F90）中定义以下物理结构：

! ====================================================================
!  ALF Fortran90 核心代码片段：定义 M 点莫尔三角形格子的动能跳跃和相互作用顶点
! ====================================================================
Subroutine Ham_Set_Hop()
    Implicit None
    Integer :: I, J, n, s, nv
    Real (Kind=Kind(0.d0)) :: t, t_perp
    
    ! 初始化跳跃矩阵 
    ! 3个谷 (nv=3)，每个格点有2个自旋通道
    Call Hop_Matrix_Init()
    
    t = 1.d0 ! 主跳跃方向强跳跃
    t_perp = 0.25d0 ! 链间次近邻跳跃
    
    Do I = 1, N_coord
        Do nv = 1, 3
            ! 依据主轴方向添加跳跃项，避免触发 t' 产生符号问题
            Call Op_make(Op_T(I, nv), N_sites)
            Op_T(I, nv)%O(I, I_next) = -t
            Op_T(I, nv)%O(I, I_perp) = -t_perp
            Op_T(I, nv)%g            = 1.d0
            ! 确保其在 PHS 算符下的自共轭性
        End Do
    End Do
End Subroutine Ham_Set_Hop

Subroutine Ham_Set_Interaction_Triangle()
    Implicit None
    Integer :: I, J, n_tri, nv
    Real (Kind=Kind(0.d0)) :: U, alpha, g_tri
    
    ! 定义电荷与谷电荷算符在实空间三角形上的分解
    ! 对于每个三角形 lacktriangle 构造一个 HS 相互作用顶点
    Do n_tri = 1, N_triangles
        Op_V(n_tri)%g = -SQRT(g_tri * dtau) ! 虚时演化系数
        Op_V(n_tri)%type = 2 ! 离散 HS 场类型 (ISPIN=2)
        ! 填充格点 R1, R2, R3 的电荷密度项 
        Do nv = 1, 3
            Op_V(n_tri)%O(R1, R1) = 1.d0
            Op_V(n_tri)%O(R2, R2) = 1.d0
            Op_V(n_tri)%O(R3, R3) = 1.d0
        End Do
    End Do
End Subroutine Ham_Set_Interaction_Triangle

3.3 物理复现参数设置指南

为了完全复现论文中第三节与第四节的关键结果（例如图 2 中的磁化强度与图 3 中的谱函数），建议在 parameters 输入文件中配置以下核心控制参数：

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

4.1 核心参考文献

1. 莫尔多谷系统与 $U(6)$ 对称性起源：
   • Călugăru, D., Vasiliou, K., Hu, H., Bernevig, B. A., Krauth, W., & Parameswaran, S. A. Stochastic Series Expansion for Mixed-Dimensional Moiré Systems, Companion Paper, [arXiv:2606.XXXXX]. (建立了链状极限下全填充范围的 SSE 理论。)
2. SnSe2 莫尔能带与模型构建：
   • Pi, H., Kwan, Y. H., Hu, H., Jiang, Y., Călugăru, D., Shan, J., Mak, K. F., Ugeda, M. M., Vergniory, M. G., & Bernevig, B. A. Engineering topological flat bands in $\Gamma$-valley moiré systems, arXiv:2605.13984 (2026). (提供了精确的第一性原理及 Wannier 轨道能带参数。)
3. 多轨道哈伯德理论与 Kugel-Khomskii 模型：
   • Kugel, K. I., & Khomskii, D. I. The Jahn-Teller effect and magnetic properties of transition metal compounds, Phys. Uspekhi 25, 231 (1982). (强关联多轨道磁性关联的理论基石。)

4.2 局限性深入剖析

尽管本文利用无符号 DQMC 取得了漂亮且极具启发性的物理结论，但对于真实的三谷莫尔系统，该方案依然存在明显的局限：

1. 掺杂（Doping）下的符号问题灾难

本工作的无符号证明极其依赖于半填充（$\nu = 3$）处的粒子-空穴对称性（PHS）。一旦体系引入电荷掺杂（即偏离 $\nu = 3$），PHS 被立即破坏，费米子行列式将重新出现极具毁灭性的符号问题。这意味着，最令超导领域感兴趣的掺杂莫尔哈伯德超导电性无法在当前 DQMC 框架下进行无偏求解。

2. 对非共形对称性破缺敏感度不足

模型中严格假设了链间最邻近跳跃 $t' = 0$。但在真实的 AA 堆叠 $t-\text{SnSe}_2$ 材料中，晶格松弛（Lattice Relaxation）和微小的扭曲变化会诱导出微小但非零的 $t'$。虽然论文通过一阶扰动理论论证了 $t'$ 的引入相当于弱的反铁磁链间耦合 $J' \sim (t')^2/U$，但在数值计算中，任何微小的 $t' \neq 0$ 都会在数学上破坏 PHS 的无符号基础。数值稳定性对模型细节的这种极度敏感降低了算法在处理更贴近实际的非理想材料时的鲁棒性。

3. 谷间洪特耦合 $J_H$ 的忽略

为了维持 $U(2)^{\otimes 3}$ 对称性以保证 DQMC 的运行，本模型忽略了真实的谷间自旋交换项（即谷间洪特耦合 $J_H$）。$J_H$ 项（形式为 $J_H \sum_{\eta > \eta'} \mathbf{S}^{(\eta)} \cdot \mathbf{S}^{(\eta')}$）在物理上具有极为关键的作用：它会强行锁死不同谷内的反铁磁自旋取向（自旋锁死，Spin-Locking），从而产生复杂的非共线磁结构。忽略该项使得 DQMC 无法对强耦合下真实的六自由度自旋-轨道耦合相图进行完全客观的刻画。

4. 最大熵（MaxEnt）解析延拓的内在多值性

利用虚时 $G(\tau)$ 反演实频谱 $A(\omega)$ 是一个数学上极端不适定的逆问题。尽管作者使用了 ALF 框架中最为先进的 MaxEnt 方法，并引入了空间群对称性平均来减小噪声，但在极低能激发的精细结构解析上（例如确认 Mott 间隙内是否存在微弱的准粒子残留结构），由于核函数的指数衰减特性，反演结果仍具有一定的人为拟合色彩。

5. 补充物理图像深化：强耦合有效自旋模型、Order-by-Disorder 与 Parton 理论

为了更好地辅助量子化学与强关联领域科研人员理解该系统，本节补充推导三个关键的辅助解析理论，它们与 DQMC 数值结果互为补充。

5.1 强耦合极限下的有效自旋模型推导

当相互作用极强（$U \gg t, t_{\perp}$）时，电荷自由度在低能处被冻结，每个晶格格点被严格限制为占据 3 个电子（一谷一电子）。利用 Rayleigh-Schrödinger 二阶微扰论，可以推导出有效的低能自旋哈密顿量：

其中 $P$ 为低能子空间的投影算符，$Q = I - P$。激发能为 $E_{\text{exc}} = U$。带入跳跃算符的显式：

利用 Fierz 自旋恒等式 $\sum_{s,s'} d^{\dagger}_{\alpha,s} d_{\beta,s} d^{\dagger}_{\beta,s'} d_{\alpha,s'} = 2 \mathbf{S}_{\alpha} \cdot \mathbf{S}_{\beta} + \frac{1}{2} n_{\alpha} n_{\beta}$，最终可以优雅地写为各谷完全解耦的反铁磁海森堡链哈密顿量：

其中超交换相互作用强度为 $J^{\eta}_{\Delta\mathbf{R}} = \frac{4 (t^{\eta}_{\Delta\mathbf{R}})^2}{U}$。由于 $t^{\eta}_{\Delta\mathbf{R}}$ 的各向异性，有效相互作用是沿链方向的强反铁磁 $J = 4t^2/U$ 与链间弱反铁磁 $J_{\perp} = 4t_{\perp}^2/U$ 的组合。

5.2 零温 Order-by-Disorder 现象与自旋波理论

对于一谷情形（例如 $\eta=0$），有效模型由两套高度独立的矩形子晶格交织而成。在经典力学极限下，两套子晶格的奈尔矢量 $\mathbf{n}_A$ 和 $\mathbf{n}_B$ 之间的相对夹角 $\phi$ 是任意的，整个系统在经典层面上存在一个连续对易的简并流形 $U(1)$。

然而，当考虑量子涨落带来的零点能修正时，简并会被强行解除。引入 Holstein-Primakoff 玻色化方案：

在动量空间中对哈密顿量进行二次型对角化，并执行 Bogoliubov 变换：

系统的量子零点能为：

其中 $A_{\mathbf{k}}(\phi)$ 与 $B_{\mathbf{k}}(\phi)$ 为自旋波动力学矩阵元。通过对简并参数 $\phi$ 的导数分析（计算结果见图 4 所示），量子零点能在 $\phi = 0$ 和 $\phi = \pi$ 处取得极小值。这属于非常典型的由无序诱导有序（Order-by-Disorder）的物理过程。量子涨落通过选择共线的条纹状（Stripe）磁结构，消除了系统的经典对易简并度。

    [Order-by-Disorder 零点能随角度变化示意]
    E_ZP(φ)
    2.1355 ┼───╮             ╭───╮             
           │   │             │   │             
    2.1345 ┼   ╰───╮     ╭───╯   ╰───╮
    2.1340 ┼         ╰─────╯             ╰─────
           └─────────┴─────────┴─────────┴─────────
           0         π/2        π        3π/2       2π  φ

5.3 Parton 平均场理论与 Quantum Rotor（量子转子）表征

为了透彻理解为什么在 $U(6)$ 对称点（$\alpha = 1$）附近会存在一个极其反常的金属相（金属相被一直拓宽到很强的临界作用力 $U_c/W \approx 6$），作者采用了一种将物理电子分解为带电转子（Charge Rotor）和中性自旋子（Neutral Spinon）的 Parton 方法：

其中 $f^{\dagger}$ 为费米型自旋子（Spinon），满足约束条件；$\hat{\varphi}_{\mathbf{R},\eta}$ 为共轭的相位转子，其角动量算符 $\hat{\ell}_{\mathbf{R},\eta} = -i\partial/\partial\varphi_{\mathbf{R},\eta}$ 表示局域电荷激发。

利用该分解，哈密顿量在平均场近似下写为：

其中自旋子的准粒子权重 $Z = Q^2$ 刻画了系统的金属相相干度。其自洽方程解（如图 5 所示）给出了如下物理图景：

• 在一般情况下，电荷涨落的刚度（Stiffness）很大（$\sim U$），迫使系统在较小的 $U$ 下绝缘化。
• 但当 $\alpha \to 1$ 时，相互作用能矩阵的本征值发生简并，谷涨落模式的激发能量降为 0（即谷激发变软）。通过转子与自旋子的动能耦合作用，原本非常坚硬的电荷刚度被软化的谷涨落完全拖垮，使得局域电荷难以被有效局域化。这也从解析层面上优雅地阐明了为什么高度对称的谷涨落能极大地稳定并拓宽金属相，导致了 DQMC 数值中所观测到的宽广的 $U(6)$ 关联交叉行为。