来源论文: https://arxiv.org/abs/2606.12508v1 生成时间: Jun 13, 2026 12:51

混合维度量子蒙特卡洛在M点莫尔材料中的深度应用：突破符号问题的算法创新与多体物理相图解析

0. 执行摘要

在强关联电子学领域，寻找既具有高度实验可调性，又能在理论和数值上进行精确、无偏（unbiased）求解的物理平台，一直是凝聚态物理学家和量子化学家追求的“圣杯”。最近，一类新型的莫尔材料——基于二维三角晶格单层（其低能状态位于布里渊区的三个M点）叠加旋转而成的M点莫尔超晶格系统（如AA堆叠的扭角$SnSe_2$），引起了学术界的广泛关注。

本研究揭示了这类系统在“混合维度”（mixed-dimensional）极限下独特的物理机制与算法突破。由于体系中存在一种非共格的、对称性保护的投影镜面对称性（projective mirror symmetry $\tilde{M}_z$），低能导带电子的单粒子跃迁被严格限制在各谷内部的一维链上，而不同谷的电子则通过二维静电库仑作用发生强烈的电荷和自旋耦合。这种特殊的“混合维度Hubbard模型”在进行随机级数展开（Stochastic Series Expansion, SSE）量子蒙特卡洛（Quantum Monte Carlo, QMC）模拟时，在任意填充数下都完全没有费米子符号问题（sign-problem-free）。

为了克服混合维度极限下因强库仑封锁（Coulomb blockade）导致的链间电荷自洽和热力学平衡极度缓慢的瓶颈，本文深入探讨了作者提出的一种全新的全局“跨链（inter-chain）”更新算法。该更新方案将传统的单链蠕虫算法（worm algorithm）扩展为多链协同的全局更新，并在低温区（$\beta \ge 12$）引入Replica-Exchange（并行回火）机制，成功实现了对AA堆叠扭角$SnSe_2$（AA t-$SnSe_2$）整个关联相图的无偏、高精度数值绘制。研究不仅确立了整数填充下的谷极化Mott绝缘体和自旋二聚体相，还在特定分数填充下发现了打破平移对称性的Wigner-Mott绝缘态，并利用部分平均场（Parton Mean-Field）和强耦合膨胀理论给出了完美的解析阐释。

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

1.1 核心科学问题：M点莫尔体系的独特物理化学特性

传统的莫尔超晶格材料（如扭角双层石墨烯TBG、过渡金属硫族化合物TMDs）其低能有效状态通常位于布里渊区的$K$点或$\Gamma$点。相比之下，M点莫尔材料（M-point moiré materials）展现出截然不同的电子结构与对称性特征：

1. 三谷结构与时间反演对称性：体系在低能下拥有三个独立的谷（对应布里渊区的三个M点），具有谷内时间反演对称性。这阻止了单个谷中单粒子陈数（Chern number）的出现，使得体系的能带结构在拓扑上是平凡的，能够构建定域化极佳的Wannier函数。
2. 投影镜面对称性（Projective Mirror Symmetries）：在AA堆叠的扭角$SnSe_2$（AA t-$SnSe_2$）中，即使考虑晶格弛豫，系统仍然保留了一个投影面内镜面对称性 $\tilde{M}_z$。这一对称性将单粒子哈密顿量在每个谷内约束为严格的准一维（quasi-1D）形式。也就是说，属于 valley-$\eta$ 的电子只能沿着特定的方向 $\mathbf{e}_\eta$ 跃迁，形成相互交叉的、成120度角的三组一维链，而链与链之间的单粒子电荷跃迁被对称性完全禁止。

然而，电子之间的库仑排斥相互作用是完全二维的。这使得系统形成了一个**交叉滑动Luttinger液体（crossed sliding Luttinger liquid）**的二维网络：单粒子运动是严格一维的，而多体关联和集体激发是完全二维的。理解这类强关联混合维度多体系统在弱耦合到强耦合过渡区的行为，是当前凝聚态理论的一大难题。

1.2 理论基础与无符号问题证明

由于系统的一维跃迁特征，我们可以利用独立的Jordan-Wigner（JW）变换将费米子算符转换为硬芯玻色子算符，从而在数学上严格证明符号问题的消失。

设 $\hat{d}^{\dagger}_{\mathbf{R},\eta,s}$ 为在胞 $\mathbf{R}$、谷 $\eta$、自旋 $s$ 的费米子产生算符。我们可以引入一维链的全局索引 $\alpha \in [1, 3\mathcal{N}]$（其中 $\mathcal{N}$ 为系统线性尺寸，共 $3\mathcal{N}$ 条独立的单谷跃迁链），以及沿着该链的格点索引 $i \in [1, N]$。重新索引后的费米子算符写为 $\hat{d}^{\dagger}_{\alpha,i,s}$。

由于单粒子跃迁哈密顿量 $H_0$ 不包含跨链跃迁算符，对于给定的自旋 $s$，我们在每条链 $\alpha$ 上独立定义Jordan-Wigner弦算符：

其中 $\hat{N}_{\alpha,s} = \sum_{i=1}^N \hat{n}_{\alpha,i,s}$ 是该链上的总粒子数算符。通过将算符映射为硬芯玻色子 $\hat{b}^{\dagger}_{\alpha,i,s}$：

在代入动能哈密顿量 $H_0$ 时，相邻格点的跃迁项 $\hat{d}^{\dagger}_{\alpha,i,s} \hat{d}_{\alpha,i+1,s}$ 中，JW弦算符在链内的项发生精确抵消。最终，哈密顿量在玻色子表象下写为：

其中，玻色子满足边界条件 $\hat{b}^{\dagger}_{\alpha,i+N,s} = \hat{b}^{\dagger}_{\alpha,i,s}$。为了使所有跃迁矩阵元非正（满足随机蒙特卡洛的Stochastic要求），费米子系统的边界条件必须根据该链上的粒子数奇偶性进行自适应调整：当粒子数为偶数时采用抗周期边界条件（antiperiodic），为奇数时采用周期边界条件（periodic）。这在热力学极限下不改变其物理性质，从而在数学上完美避开了费米子负符号问题。

1.3 技术难点：强静电相互作用下的谷间/链间电荷自洽瓶颈

尽管符号问题在数学上被消除，但在实际的蒙特卡洛采样中存在灾难性的收敛性技术难点：

• 在经典的SSE算法中，所有的非对角更新（对谷内一维链上的玻色跃迁进行蠕虫环更新）都是局域在单条一维链上的。也就是说，一个 Worm 只能在一条确定的链 $\alpha$ 上游走并改变其局域电荷排布和粒子数。
• 然而，不同链之间的电荷耦合是通过二维静电库仑作用 $H_I$ 实现的。在低温极限下，系统由于强烈的各谷互排斥（谷间Hubbard作用 $U'$ 以及最近邻排斥 $V$），会形成极高的能垒。这种现象类似于一维链之间的“库仑封锁”（Coulomb blockade）。
• 在这种状态下，单链更新无法独自完成从一条链向另一条相交链的“电荷转移”。要改变不同链之间的粒子数分布（例如从配置 $s_1 = \{N_1, N_2\}$ 到 $s_2 = \{N_1 - 1, N_2 + 1\}$），在传统的采样中需要经历极长、且统计权重极低、甚至不可逾越的中间高能量状态。这导致不同谷（链）之间的总粒子数被彻底电荷冻结（charge freezing），计算的自相关时间（autocorrelation time）呈指数级发散，数值模拟无法收敛到真正的谷自洽基态。

1.4 方法细节：创新的全局“跨链更新”（Inter-chain Updates）

为了攻克上述核心技术瓶颈，本工作设计了一种极其精妙的**跨链全局联合更新（Inter-chain QMC Update）**方法：

[1. 随机选取相交格点 R 及两个关联谷 η, η']
                     │
                     ▼
[2. 分割哈密顿量权重：W(s) = W_d(s) × W_o(s)]
  (W_d 包含支持在 R 点、涉及谷 η 和 η' 的对角相互作用顶点算符)
                     │
                     ▼
[3. 暂时屏蔽 W_d：在无 W_d 的高效权重场 W_o(s) 下进行蠕虫传播]
  (分别在链 η 和 η' 上构造两个连续的 Worm Loop)
                     │
                     ▼
[4. 精确物理跃迁提案：在链 η 上消除一个粒子，在链 η' 上注入一个粒子]
                     │
                     ▼
[5. 延迟 Metropolis 滤波：利用 W_d 的初始与最终权重计算联合接受率]
  (Metropolis 接受度：P_accept = min[1, (p_worm_reverse * W_d(s')) / (p_worm_forward * W_d(s))])

详细步骤解析：

1. 格点与谷对的选择：随机选择一个三链相交的格点 $\mathbf{R}$。在三个谷内随机选择两个参与更新的谷 $\eta$ 和 $\eta'$。
2. 权重分割：将整体配置的统计权重 $W(s)$ 分割为局域对角相互作用权重 $W_d(s)$ 和剩余权重 $W_o(s)$ 乘积的形式。其中 $W_d(s)$ 只包含所有作用在格点 $\mathbf{R}$ 且涉及谷 $\eta$ 和 $\eta'$ 的对角 12-legged 顶点算符。这一步是整个跨链算法的精髓：由于对角相互作用极其阻碍电荷在两条链间的转移，我们必须在生成提案的过程中暂时“屏蔽”它。
3. 无势蠕虫传播：在不考虑 $W_d(s)$ 阻碍的背景下，即在由 $W_o(s)$ 定义的高效场中，运行两个连续的单链 Worm 循环。第一条 Worm 在链 $(\mathbf{R},\eta)$ 上移去一个电荷；第二条 Worm 在链 $(\mathbf{R},\eta')$ 上增加一个电荷。因为移去了最强烈的局域库仑阻碍，这两个 Worm 的闭合概率极高。
4. 延迟接受度计算（Delayed Acceptance）：一旦两条 Worm 顺利闭合，实现了完美的“跨链电荷转移”（其局域格点的物理态完成了从谷 $\eta$ 向 $\eta'$ 的单粒子跳跃），我们重新计算被暂时忽视的 12-legged 对角顶点的贡献。最终的联合接受度由下式给出：

其中 $p^{\text{worm}}_{s \to s'}$ 是生成上述双蠕虫路径的提案概率（其具体的微观状态概率汇总在论文附录的 Table S5 中）。因为扭角 $SnSe_2$ 体系中的局域库仑排斥具有高度的近似 $U(6)$ 对称性（谷内排斥 $U$ 与谷间排斥 $U'$ 非常接近，$U \approx U'$），在电荷完成谷转移后，对角顶点在新态 $s'$ 下的矩阵元之和与旧态 $s$ 几乎完全相等（$W_d(s') \approx W_d(s)$）。这使得这一跨链全局更新的最终接受率极高，彻底粉碎了低温下的库仑封锁，使混合维度QMC模拟得以在极低温度下平稳高效地运行。

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

为了验证算法的正确性、收敛速度以及探究实际 AA t-$SnSe_2$ 材料的强关联物理，作者开展了极其详尽的 Benchmark，并给出了大尺寸体系下的相图计算。系统格点数设置为 $N \times N = 12 \times 12$，共对应 $3 \times 12 \times 12 \times 2 = 864$ 个单粒子量子态。

2.1 18套物理参数体系的设计

实验上，扭角莫尔材料的强关联效应可以通过改变扭角 $\theta$（调节动能带宽 $t$）、介电常数 $\epsilon$（调控库仑相互作用强度）以及栅极屏蔽长度 $\xi$（调控最近邻排斥力 $V$ 与胞内排斥力 $U$ 的比例）来大范围调控。为了与实验直接对接，作者基于第一性原理计算导出的 Wannier 模型，精心设计了 18 套物理参数体系，覆盖了从弱、中等耦合到强耦合的所有典型关联区间（详细参数对应表见下表）：

2.2 与精确对角化（ED）和行列式蒙特卡洛（DQMC）的 Benchmark 验证

为了检验新算法在消除符号问题以及引入跨链更新后的数值精度，作者在小尺寸体系上进行了一系列严格的算法对撞实验。

A. 一维极限与精确对角化（ED）对比（验证动能算符及单链更新）

在严格的一维极限下（设谷间相互作用 $U' = 0$，最近邻相互作用除一维链方向外均置为零），体系退化为 $3\mathcal{N}$ 条互不干扰的一维带排斥力 Hubbard 链。作者在 $\mathcal{N}=6$ 的体系上，对比了精确对角化与 SSE 在胞内排斥力 $U=4$、两种不同一维最近邻相互作用 $V_4=0$ 和 $V_4=1.5$ 时的单链能级。正如论文图 S36 所示：

• 在整个化学势 $\mu \in [-1, 3.5]$ 的扫描区间内，计算得到的每链平均电荷数 $\langle\hat{N}\rangle$ 和每链平均能量 $\langle H\rangle$ 的 SSE 曲线与 ED 离散格点数据实现了完美的无缝重合。
• 这直接证明了基于奇偶性修正的自适应一维边界条件在消除符号问题上的数学严密性。

B. 二维强关联极限与 ED 对比（验证跨链更新）

在胞内排斥 $U=U'=5$，所有最近邻相互作用均设为 $V_i = 0.5$ 的二维强耦合极限下（系统尺寸 $\mathcal{N}=2$），作者计算了总能量 $\langle H\rangle$、总粒子数 $\langle\hat{N}\rangle$ 以及其热力学涨落 $\langle\hat{N}^2\rangle_c$ 随温度（等效倒温度 $\beta \in [1, 10]$）的变化趋势（论文图 S37）：

• 在包含跨链更新的完整 SSE 算法下，能量、电荷以及敏感的二阶涨落 $\langle\hat{N}^2\rangle_c$ 的模拟值与 ED 结果在所有温度区间完全一致，误差线远小于线宽。
• 这无可辩驳地证实了跨链更新方案满足细致平衡条件（detailed balance），没有引入任何物理系统偏差。

C. 与行列式蒙特卡洛（DQMC）在半填充下的对比（性能及收敛速度验证）

半填充 $\nu=3$ 处，该体系具有隐藏的二部图（bipartite）晶格对称性，此时 DQMC 同样没有符号问题（见论文图 1 示意）。作者在 $\mathcal{N}=4$ 尺寸、胞内排斥 $U=12, V_i=0$ 的超强关联极限下，将不含/包含跨链更新的两种 SSE 算法与 DQMC 数据进行了对撞（论文图 6(d) 和图 S39）：

• 当不包含跨链更新时（图 6(d) 橙色曲线），由于低能下发生电荷冻结，谷不平衡涨落 $\chi_{\eta\eta'}$ 随 $\beta$ 升高而剧烈偏离 DQMC 真实解（甚至在 $\beta \ge 6$ 后发生数值崩塌，无法正确收敛）。
• 当引入跨链更新后（图 6(d) 紫色曲线），SSE 结果与 DQMC 在所有温度下实现了完美的无缝契合，自相关时间降低了2 到 3 个数量级，成功使低温区域的计算成为可能。

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

由于该算法涉及混合维度一维链与二维库仑作用的特殊拓扑结构，以及独特的跨链全局非对角更新，现有的标准开源QMC框架（如 ALPS 库、standard Stochastic Series Expansion 代码）无法直接运行。必须在基础 QMC 代码上进行高度定制化编写。以下是基于 C++ 或 Julia 复现该核心算法的技术架构与实操指南。

3.1 核心数据结构设计

为了实现高效的对角与非对角蠕虫更新，核心状态必须包含：

1. 格点状态数组（State Array）：大小为 $3\mathcal{N}^2 \times 2$，其中 $3\mathcal{N}^2$ 代表总 Wannier 轨道数（每个莫尔原胞包含 3 个谷），$2$ 代表自旋。每个元素的值为 $\{0, 1\}$，分别代表空态和占用态。
2. 算符串数组（Operator String）：大小为固定截断长度 $L$ 的一维数组（在热力学平衡期间自适应调整）。数组的每个元素为一个二元组 [type, bond]：
   • [0, 0]：恒等算符 $H_{0,0}$，用于对齐级数长度。
   • [1, b]：对角相互作用算符 $H_{1, b}$（12-legged 顶点，作用于两个胞的 6 个轨道上）。
   • [2, b]：自旋为 $\uparrow$ 的非对角跃迁算符 $H_{2, b}$。
   • [3, b]：自旋为 $\downarrow$ 的非对角跃迁算符 $H_{3, b}$。
   • [4, b]：对角常量算符 $H_{4, b}$（4-legged 顶点，用于非对角更新的转换）。
3. 顶点索引网络（Vertex Link List）：由于 12-legged 和 4-legged 顶点的存在，必须在非对角更新前构建一条格点在虚时（传播指数 $p \in [1, L]$）方向的链表，以便蠕虫算符（Worm）在遇到顶点时能快速决定其出射腿（exit leg）。

3.2 跨链更新算法的 C++ 核心代码实现框架

以下是跨链更新的典型伪代码复现实现：

#include <vector>
#include <random>
#include <cmath>

struct Operator {
    int type; // 0: identity, 1: diagonal_12leg, 2: hop_up, 3: hop_down, 4: diagonal_4leg
    int bond; 
};

// 模拟体系核心类
class MixedDimQMC {
public:
    int N_size; // 晶格线性尺寸 N
    int L_cutoff; // SSE截断长度 L
    std::vector<int> state; // 局域电荷状态：[site][valley][spin]
    std::vector<Operator> op_string; // 算符级数串
    
    // 核心：跨链全局更新提案
    void inter_chain_update(double beta, std::mt19937& rng) {
        std::uniform_int_distribution<int> select_site(0, N_size * N_size - 1);
        std::uniform_int_distribution<int> select_valley(0, 2);
        
        int focus_site = select_site(rng);
        int v1 = select_valley(rng);
        int v2 = (v1 + 1 + (select_valley(rng) % 2)) % 3; // 确保选择两个不同的谷
        
        // 1. 临时移除对角相互作用项，计算初始权重贡献 W_d
        double log_W_d_initial = calculate_local_diagonal_weight(focus_site, v1, v2);
        
        // 2. 构造两个 Worm，分别在谷 v1 和 v2 对应的一维链上游走
        bool worm_v1_success = run_worm_on_chain(focus_site, v1, -1, rng); // 移除一个电荷
        if (!worm_v1_success) { 
            rollback_changes(); 
            return; 
        }
        
        bool worm_v2_success = run_worm_on_chain(focus_site, v2, +1, rng); // 注入一个电荷
        if (!worm_v2_success) { 
            rollback_changes(); 
            return; 
        }
        
        // 3. 计算完成转移后的新对角相互作用项权重 W_d_final
        double log_W_d_final = calculate_local_diagonal_weight(focus_site, v1, v2);
        
        // 4. 获取转移提案概率纠正因子
        double transition_factor = get_worm_proposal_probability_ratio(focus_site, v1, v2);
        
        // 5. 进行延迟 Metropolis 滤波接受测试
        double accept_prob = std::exp(log_W_d_final - log_W_d_initial) * transition_factor;
        std::uniform_real_distribution<double> dist_zero_one(0.0, 1.0);
        
        if (dist_zero_one(rng) < accept_prob) {
            // 成功接受更新，状态持久化
            commit_changes();
        } else {
            // 拒绝更新，状态回滚
            rollback_changes();
        }
    }

private:
    double calculate_local_diagonal_weight(int site, int valley1, int valley2) {
        double log_weight = 0.0;
        // 扫过整个 op_string，只对支持在 site、且属于 valley1 或 valley2 上的对角 12-legged 算符累加能量贡献
        for (const auto& op : op_string) {
            if (op.type == 1 && bond_contains_site_and_valleys(op.bond, site, valley1, valley2)) {
                log_weight += std::log(get_matrix_element(op));
            }
        }
        return log_weight;
    }
    
    bool run_worm_on_chain(int site, int valley, int charge_change, std::mt19937& rng) {
        // 标准一维特定谷链上的 Worm 传播算法
        // 在传播时不考虑与另一条链的 12-legged 耦合强度，遇到 12-legged 对角算符直接穿透而过
        return true; 
    }
    
    double get_worm_proposal_probability_ratio(int site, int v1, int v2) {
        // 根据 Table S5 矩阵元素，提取 p_worm(s'->s) / p_worm(s->s') 的精确解析比值
        return 1.0;
    }
    void rollback_changes() {} 
    void commit_changes() {}
    bool bond_contains_site_and_valleys(int b, int s, int v1, int v2) { return true; }
    double get_matrix_element(const Operator& op) { return 1.0; }
};

3.3 并行回火（Parallel Tempering）配置复现指南

为了使在强 Mott 区（例如参数 15、18 对应的超低温极限 $\beta \approx 40$）下的量子相图得到精确绘制，系统需要配置高效的并行回火温度阶梯。

• 极值设定：极高温端设定为 $\beta_{\text{min}} = 0.1$（此时体系处于完全各向同性的局域高温电荷涨落相，没有收敛阻碍），极低温端设为目标观察温度 $\beta_{\text{max}} = 12 \sim 40$。
• 阶梯分布策略：不要采用传统的线性温度梯度，这会导致交换概率在低温端雪崩式下跌。必须采用自适应的平方根温度梯度分布（论文附录公式 S3.39）：

• 交换接受率：两个相邻温度 $\beta$ 和 $\beta'$ 的交换尝试接受率由下式给出（论文公式 S3.32）：

其中 $n$ 和 $n'$ 分别是相邻 Replica 此时的非恒等算符级数级数截断数。平方根阶梯分布可以确保在整个温度区间内相邻 Replica 的电荷激发分布谱有稳定的重叠，从而在全温度段实现稳定大于 $20\%$ 的高效 Replica 交换率。

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

4.1 关键经典文献背书

本工作建立在凝聚态量子蒙特卡洛和强关联多体理论的诸多基石研究之上：

1. Stochastic Series Expansion (SSE) 方法奠基：Sandvik 在其经典文献 [Phys. Rev. B 59, R14157 (1999)] 和 [Stochastic Series Expansion Methods (2019), arXiv:1909.10591] 中，系统建立了将量子哈密顿量进行高氏虚时幂级数展开并进行无偏采样的完备数学架构，这是本算法的底座。
2. Rotor Parton 理论：Florens & Georges 在文献 [Phys. Rev. B 66, 165111 (2002)] 和 [Phys. Rev. B 70, 035114 (2004)] 中提出了强关联多轨道系统的量子转子（Rotor） Parton 拆分方案，这一解析理论为本研究的相图物理规律（特别是 $\nu=1$ 处 Mott 绝缘体最稳定、临界作用力最小的性质）提供了强有力的解析微观机理解释。
3. M点莫尔物理的提出：Călugăru 等人于 [Nature 643, 376 (2025)] 中首次揭示了 M 点莫尔超晶格系统中投影对称性引发的准一维跃迁行为，本工作正是针对该基础材料物理的首个精确量子数值模拟。

4.2 局限性深度剖析与批判性讨论

尽管本工作在数值物理学和量子算法层面上取得了非凡的突破，但作为一项具有开拓性的研究，它依然包含若干不可忽视的局限性：

1. 链间单粒子跃迁 $t_\perp$ 的直接丢弃对相图相边界的潜在扰动

为了满足极其干净的“混合维度极限”，从而在数学上通过独立链 Jordan-Wigner 映射彻底干掉符号问题，作者在模拟哈密顿量中直接抹去了胞间次近邻跨链跃迁 $t_\perp$（以及面内弛豫引起的 $\tilde{M}_z$ 对称性微破缺项 $t'$）。根据第一性原理 Wannier 投影的实际物理参数（见论文附录表 S4），对于真实的 AA t-$SnSe_2$ 晶格：

• 在扭角为 $\theta = 7.34^\circ$ 时，胞间跨链单粒子跃迁强度 $t_\perp / t \approx 0.254$，占了主导动能带宽的四分之一，这绝非是一个微扰小量。
• 实验上，如此强度的跨链跃迁会导致一维物理链之间的相干杂化，使“交叉滑动 Luttinger 液体”向完全各向同性的二维 Fermi 液体发生转化。虽然在强相互作用下由于局域化效应 $t_\perp$ 的有效带宽会被极大地重整化缩减（Renormalized bandwidth reduction），但它会直接破坏 1D 电荷守恒，导致系统在除了半填充 $\nu=3$ 外的区域重新产生费米子符号问题。因此，本工作基于 $t_\perp = 0$ 绘制的 $\nu=1$ Mott 相边界和 Wigner-Mott 相相边界，在实验对照中可能会存在由于单粒子杂化带来的温区及相互作用阈值的系统性偏移。

2. 对自旋轨道耦合（SOC）与时间反演对称性破缺相互作用的简化

在真实的过渡金属及硒化物超晶格中，强自旋轨道耦合（Spin-Orbit Coupling, SOC）和非共面局域磁结构导致的非局域相互作用往往非常活跃：

• 论文为了保持 QMC 顶点的正定性，忽略了这些会带来复数矩阵元的非局域自旋跃迁项和非共振相互作用项。
• 这些被忽略的项可能会诱导出拓扑手性超导态（Chiral Superconductivity）或陈绝缘体（Chern Insulators）等非凡相，这些相在本算法框架下因为符号问题重新出现而无法被无偏模拟。

3. 大尺寸扩展下的计算开销瓶颈

由于跨链全局更新算法包含了延迟 Metropolis 滤波，每一次接受测试都需要扫过整个算符串来计算 $W_d(s')$ 强度的指数矩阵元变化：

• 当系统线性尺寸 $\mathcal{N} \to 24$ 甚至更大时，算符串的长度 $L$ 会随着 $\beta N^2$ 线性发散（通常在极低温下达到 $L \ge 10^5$），这会导致 $W_d(s)$ 的搜索和求和过程变得极其昂贵，算法的单步复杂度可能会从标准 SSE 的线性格点标度退化为超线性（甚至是二次标度），限制了其向极低能热力学极限的大尺寸外推。

5. 其他必要的补充

5.1 解析互补：Parton 平均场理论对 Mott 相稳定性等级的论证

为了让物理图像超越单纯的蒙特卡洛数值拟合，作者引入了强大的量子转子 Parton 理论（Rotor Parton Theory）对相图行为进行了完美的物理透视。

将电子算符拆分为携带电荷的量子转子 $e^{i\hat{\varphi}_\mathbf{R}}$ 和不带电荷的自旋子（spinon） $\hat{\gamma}_{\mathbf{R},\eta,s}$：

利用自洽约束 $\sum_{\eta,s} \hat{\gamma}^{\dagger}_{\mathbf{R},\eta,s}\hat{\gamma}_{\mathbf{R},\eta,s} - 3 \equiv \hat{Q}_{\text{ch}}$，其中 $\hat{Q}_{\text{ch}}$ 代表局域总电荷算符。在自洽平均场下，量子转子的相干重整化因子（Quasiparticle spectral weight $Z = Q^2 = \langle \cos \hat{\varphi} \rangle^2$）决定了系统是金属还是绝缘体。当相互作用足够强，重整化因子发生凝聚并骤降至 $Z=0$ 时，转子发生相位去相干，系统发生 Mott 金属-绝缘体转变（MIT）。

通过计算不同整数填充下的基态局域简并度，可以完美预测这一转变的敏感度。根据 U(2)$^{\otimes 3}$ 对称性哈密顿量，不同填充下的格点微观电子排布简并度 $d_\nu$ 表现出巨大的差异：

• 对于 $\nu = 1$：简并度 $d_1 = \binom{6}{1} = 6$。系统可以通过最少的配型方式在单格点定位一个电子，电荷涨落最小，因而转子相干性最容易被破坏。
• 对于 $\nu = 2$：因为胞内排斥满足 $U > U'$，两个电子必须占据两个不同的谷。其配型退化简并度为 $d_2 = \binom{6}{2} - 3 = 12$。
• 对于 $\nu = 3$：三个电子占据三个不同的谷，简并度为 $d_3 = 2^3 = 8$。

根据重整化激发熵理论，较小的简并度对应着较大的关联稳定区间。因此，自洽 Parton 计算精确地预测出 Mott 绝缘体形成的稳定层级为：

这一结论与无偏量子蒙特卡洛模拟得到的逆压缩率峰值高度及电荷能隙大小（论文图 7 第四、六列）具有惊人的一致性，形成了“数值模拟-微观物理机制-自洽解析图景”的闭环验证。

5.2 分数填充下的物理：打破莫尔平移对称性的 Wigner-Mott 绝缘体

当最近邻排斥力 $V$ 较大时（对应参数集 1, 4, 7 等低屏蔽长度 $\xi=2.5\text{ nm}$ 体系），由于相邻原胞间的强静电排斥，电子倾向于在空间中形成电荷密度波（CDW），以实现库仑能的最小化。

在混合维度 QMC 模拟中，作者在分数填充 $\nu = 1/3, 5/3, 2$ 处观测到了极其强烈的动量解析电荷易受性 $\chi(\mathbf{k})$ 峰值，且该峰值精确锁死在莫尔布里渊区的布拉格点 $K_M$ 上：

这表明体系的实空间电荷密度打破了原有的莫尔晶格平移对称性，自发凝聚为等效的 $\sqrt{3} \times \sqrt{3}$ 空间超胞（Supercell）。

实空间电子定域化构型解析（对应论文图 8(g)-(j)）：

• 在 $\nu = 1/3$ 处：每三个原胞包含一个电子。电子精准定位在 $\sqrt{3} \times \sqrt{3}$ 超胞的顶点上，胞内其余格点为空，形成了完美的经典三角形 Wigner 晶体。
• 在 $\nu = 2$ 处：通常人们认为 $\nu = 2$ 在高 $V$ 下应该是一个均匀的 Mott 绝缘体（每个原胞均匀分布 2 个电子）。但在这里，强烈的最近邻相互作用驱动系统进入了一个全新的 Wigner-Mott 绝缘态：超胞中的三个格点发生了电荷自发歧化，分别展现出 $\{2, 2, 2\}$（三个格点各占 2 个电子，简并度极高）向极低温下的空间不均匀共振排布转化。这表明二维长程静电相互作用能够精细地操控混合维度系统中的电荷空间拓扑织构。

5.3 强耦合极限下的低温自旋相行为

当温度进一步降低，电荷激发被彻底冻结在整数填充的绝缘体内部时，系统的剩余物理完全由自旋自由度主导。通过对模拟得到的虚时自旋-自旋关联函数 $S(i)$ 展开高阶级数拟合与熵 $S(T)$ 的积分提取（论文图 9）：

1. 对于半填充 $\nu=3$：每个胞的三个谷各占一个自旋 $1/2$ 的电子。三个相互平行的 1D 自旋链方向通过 $120^\circ$ 交叉贯穿整个二维平面。在极低温下，由于谷内单粒子跃迁的虚时过程，产生了谷内 Heisenberg 反铁磁超交换作用（exchange coupling $J = 4t^2/U$）。模拟结果显示，在 $\beta \approx 15$ 以下，电荷涨落熵骤降至零，而自旋 Stiffness 则保持有限，体系物理退化为三组相互重叠、自旋非共格的、经典的一维反铁磁 Heisenberg 自旋链，激发谱为无能隙的反铁磁自旋子（spinon）。
2. 对于 $\nu=1$ 和 $\nu=2$：基态系统为了使交换能最小化，自发选择了局域自旋二聚化（Spin dimerization）。实空间分析表明，体系被一维自旋单态二聚体（Spin-singlet dimer tiles）以经典的无序方式（对于 $\nu=1$）或打破平移对称性的价键固体（Valence-Bond Solid, VBS）方式（对于 $\nu=2$）铺满，在温区 $\beta \ge 35$ 展现出明显的自旋能隙。

本研究所展现的“混合维度量子蒙特卡洛”算法创新以及对 M 点莫尔强关联相图的无偏绘制，不仅为 AA t-$SnSe_2$ 等材料在实验上的输运和扫描隧道显微镜（STM）测量提供了极其精确、可直接对撞的理论基准，也为探索轨道选择性、强各向异性强关联量子化学体系开辟了一条崭新且极具潜力的数值研究路径。