来源论文: https://arxiv.org/abs/2605.21632v1 生成时间: Jun 08, 2026 18:54

图符蒙特卡洛（DiagMC）在费米子 Rényi 纠缠熵计算中的突破：理论、算法与大规模应用

0. 执行摘要

在现代凝聚态物理与量子化学的前沿研究中，量子纠缠已成为诊断强关联多体系统物态、拓扑序以及量子临界行为的核心工具。然而，由于费米子系统特有的反对称统计规律，计算相互作用费米子的 Rényi 纠缠熵（Rényi Entanglement Entropy, EE）一直面临巨大的理论与计算挑战。传统的辅助场量子蒙特卡洛（AFQMC）等方法在处理二维（2D）强关联费米子系统时，往往受到极其严重的“费米子符号问题”的制约，且在处理空间子区域的纠缠熵时，复杂的时空拓扑边界条件使得计算成本随着子区域尺度的增加而指数级上升。

本研究提出并实现了一种创新的直接图符蒙特卡洛（Diagrammatic Monte Carlo, DiagMC）计算框架，专门用于精确求解相互作用格点费米子系统的 Rényi 纠缠熵。该方案的核心思想在于从费米渐变交换（Fermionic Graded-Swap）算符表象出发，将纠缠熵的求解转化为一个在多叶黎曼面（Replicated Manifold）上、具有混合时间边界条件的非均匀路径积分问题。在此表象下，本方法不仅成功推导出了半移动（Half-shifted）的复制品动量（Replica Momenta），还构建了适用于连接行列式（Connected-Determinant, CDet）求和的行列式展开式。

为了克服高阶微扰展开中的拓扑组合爆炸问题，该框架将高阶 Feynman 图展开与多组态马尔可夫链蒙特卡洛（Many-Configuration Markov Chain Monte Carlo, MCMCMC）采样技术深度融合。在基准测试中，该算法在 $3 \times 3$ 周期性 Hubbard 团簇上展现出了与精确角动量对角化（Exact Diagonalization, ED）完全一致的微扰系数（误差在 $10^{-6}$ 量级）；在实际生产计算中，该算法成功应用于高达 $40 \times 40$ 格点、子区域尺度为 $16 \times 16$ （$|A|=256$）的超大体系。这表明该方法不仅彻底规避了常规配置空间量子蒙特卡洛中的符号问题，而且其计算瓶颈成功从“指数级符号墙”转化为“多项式级内存开销”。本工作为研究高温超导、赝能隙物理、非费米液体以及拓扑费米子系统中的量子纠缠铺平了道路。

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

1.1 核心科学问题与研究背景

在强关联费米子系统中，纠缠熵（特别是第二 Rényi 纠缠熵 $S_2$）不仅满足面积律（Area Law），在存在费米面的非简并系统中还会表现出对面积律的对数违背（Logarithmic Violation）。这种对数违背的系数与费米面的几何拓扑特征（如 Fermi Volume 和 Euler Characteristic）直接相关。然而，如何定量、精确地计算相互作用费米子系统的纠缠熵是多体理论中的著名难题。原因有二：

1. 费米子符号问题：在二维及更高维度下，由于费米子波函数的反对称性，传统的行路径积分 QMC 方法在离开半满（Half-filling）或引入化学势时，权重函数会出现剧烈的正负抵消，导致方差呈指数级发散。
2. 时空拓扑的变化：计算 $S_2[A]$ 需要采用所谓的“复制品方法（Replica Trick）”。在路径积分表象下，这意味着系统在空间子区域 $A$ 内沿着虚时间方向存在跨越不同复制品（Replica Sheets）的循环边界条件，而在其余子区域 $B$ 内保持常规边界条件。这种非均匀、混合的时间边界条件极大地破坏了系统的平移对称性，使得常规格林函数的计算和蒙特卡洛采样变得异常困难。

1.2 理论基础：费米渐变交换算符与路径积分

第二 Rényi 纠缠熵的定义为：

其中 $\rho_A = \text{Tr}_B \rho$ 是子区域 $A$ 的约化密度矩阵，整个体系的密度矩阵为 $\rho = e^{-\beta H}/Z_1$。为了避免直接对约化密度矩阵进行复杂的对角化，引入双复制品空间 $\mathcal{H} \otimes \mathcal{H}$，并将配分函数之比表示为：

其中 $S_A^{\text{gr}}$ 是费米渐变交换算符（Fermionic Graded-Swap Operator）。在同构基底状态下，其作用定义为：

其中 $p(a_i)$ 表示子区域 $A$ 中费米子态的宇称（Parity），即费米子数模 2 的余数。这个特殊的 Koszul 负号 $(-1)^{p(a_1)p(a_2)}$ 是费米子系统区别于玻色子系统的关键。如果错误地将其替换为常规的玻色循环置换，将会在自由费米子极限下直接导致错误的物理相。

通过该算符，可建立相干态路径积分。对于 $n$ 复制品系统，作用量（Action）在虚时间区间 $[0, \beta]$ 内保持局域：

其中 $h_{ij} = -t_{ij} - \mu \delta_{ij}$ 为单粒子哈密顿量矩。所有的空间不均匀性与纠缠区域 $A$ 的信息均被编码到虚时间边界条件中。具体而言，在非纠缠区 $B$ 上，场变量满足标准的温度反周期边界条件：

而在纠缠区 $A$ 上，场变量通过渐变交换算符进行跨复制品循环粘合：

注意，最后一个复制品向第一个复制品过渡时引入的负号，正是费米渐变宇称在相干态边界条件下的数学体现。

1.3 复制品傅里叶变换与格林函数求解

为了解耦上述边界条件，在复制品指数空间中引入离散傅里叶变换。由于非等时反周期的存在，允许的复制品动量（Replica Momenta）被精确地半移动（Half-shifted）：

在此傅里叶扇区内，混合边界条件可写为：

其中 $\Omega_m$ 为定义在单粒子希尔伯特空间上的混合投影算符：

这里 $P_A$ 和 $P_B$ 分别是子区域 $A$ 和 $B$ 的单粒子空间投影算符。不难发现，当 $P_A = 0$（即无纠缠区）时，$\Omega_m = \mathbf{I}$，边界条件退化为常规的反周期条件。

定义自由格林函数满足运动方程：

结合混合边界条件 $\bar{G}^{(m)}(\beta^-, \tau') = -\Omega_m \bar{G}^{(m)}(0^+, \tau')$，我们可以解析求解得到各扇区的精确格林函数：

其中关键算符 $A_m$ 定义为：

由于 $\Omega_m$ 包含了复数相因子 $e^{iq_m}$，格林函数 $\bar{G}^{(m)}$ 在复数域内运算。最终的实体空间-虚时间格林函数通过逆傅里叶变换重构：

这套公式优雅地将复杂的时空不均匀边界条件完全吸收到单粒子算符 $A_m$ 的求逆过程中。这是本方法能够进行大规模计算的数学基石。

1.4 技术难点：微扰级数展开与连接行列式（CDet）

当我们引入相互作用 $U$ 时，配分函数之比 $Z_n[A]/Z_{0,n}[A]$ 可以通过相互作用表象下的 Wick 定理写成行列式展开的形式。在第 $k$ 阶微扰下，自旋 $\sigma$ 的贡献对应一个 $k \times k$ 的格林函数矩阵的行列式：

其中 $X_l = (i_l, \tau_l, \alpha_l)$ 表示第 $l$ 个相互作用顶点的时空和复制品坐标。然而，直接计算 Feynman 图的总和会面临拓扑重叠导致的巨大冗余。更为严重的是，我们需要的是对数量 $\ln Z_n[A]$，这意味着我们必须将展开式限制在**连通图（Connected Diagrams）**上。

本工作采用连接行列式（Connected Determinant, CDet）算法。对于一组给定的时空顶点集合 $\mathbf{X}$，通过子集递归（Subset Recursion）方法，可以在 $O(3^k)$ 的计算复杂度内精确提取出所有连通图的自能与自由能贡献，从而直接对对数配分函数进行微扰级数展开：

这里引入了同伦微扰参数 $\xi$（实际计算中令 $\xi=1$），通过这种方式，原本需要处理的无数 Feynman 拓扑图被巧妙地归纳为高效率的行列式求和。由于不显式绘制 Feynman 图，该方案消除了经典微扰论中的拓扑标记（Labeling）难题。

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

2.1 $3 \times 3$ 周期性 Hubbard 团簇基准测试

为了绝对验证上述“费米渐变交换表象”与“连接行列式蒙特卡洛”结合的正确性，作者首先在一个小尺寸、可进行精确角动量对角化（ED）的系统上进行了严苛的基准测试。

体系参数设置：

• 几何结构：$3 \times 3$ 二维正方格点，具有周期性边界条件（PBC）。
• 相互作用强度：$U = 4$，处于典型的中等偏强关联区域。
• 物理参数：化学势 $\mu = 1$，反温度 $\beta = 6$（中低温区）。
• 自旋对称性相关系数：$\alpha_\uparrow = \alpha_\downarrow = 1.7$（用于移频同伦路径中）。
• 纠缠子区域 $A$：选取第一行三个格点（即子区域大小为 $3 \times 1$，占总格点数的 $1/3$）。
• 切比雪夫逼近：使用 $25 \times 25$ 阶切比雪夫多项式来逼近和插值格林函数的虚时间依赖性，以保证高精度。

MCMC采样设置：

• 预热步数（Warmup Steps）：$2 \times 10^7$
• 生产计算步数（Production Steps）：$4 \times 10^8$ steps，分布于 128 条独立马尔可夫链。

数据对比：

通过精确对角化（ED），对配分函数进行有限差分（Finite-difference）求导，提取出前三阶微扰系数。表1展示了 ED 结果与本工作图符蒙特卡洛（CDet）估算值的惊人一致性：

由表1可见，CDet 方法在极高精度上复现了精确对角化的微扰级数系数。第一阶微扰系数由于其非随机部分的确定性，实现了完美对齐；第二阶和第三阶的统计误差分别控制在 $2.1 \times 10^{-5}$ 和 $6.9 \times 10^{-6}$ 以内。这一结果强有力地证明了该算法在数学构架和代码实现上的双重无误性。

2.2 $40 \times 40$ 超大格点生产计算

在基准测试通过后，作者展示了该算法在常规方法（如 ED 或传统 AFQMC）根本无法企及的宏观尺度下的计算能力。

体系参数设置：

• 几何结构：$40 \times 40$ 二维周期性正方格点（总格点数 $N = 1600$）。
• 相互作用与温度：$U = 3$（弱到中等关联），$\mu = 1$，$\beta = 8$（极低温区，对于 $40 \times 40$ 体系具有很高的物理挑战性）。
• 移频系数：$\alpha_\uparrow = \alpha_\downarrow = 1.3$。
• 纠缠子区域 $A$：中心正方形区域，半径 $r = 8$，即边长为 16 的正方形，包含格点数 $|A| = 256$。在这个尺度下，传统的边界粘合配分函数法会因为巨大的子区域周长而彻底失效。
• 计算资源：生产跑使用 $2 \times 10^8$ 采样步，同样在 128 条独立链上并行分发。

阶数展开系数数据分析：

在图 1（论文第 4 页）的右侧面板中，展示了直至第 8 阶的微扰系数 $c_k / N$：

• 第 0 阶（无相互作用项，自由费米子贡献）：$f_A^{(0)} = -0.04089329$。该值是通过对自由格林函数矩阵直接进行行列式求解得到的确定性常数。
• 奇数阶项与偶数阶项的交替演化：微扰系数随阶数 $k$ 的增加呈现出阻尼振荡特征。第二阶系数 $c_2/N \approx -4.7 \times 10^{-3}$，处于明显的主导地位，后续高阶项（$k=3$ 到 $k=8$）快速衰减，表明微扰展开级数在 $U=3$ 时具有优异的收敛半行为，可以通过 Padé 逼近或简单的直接求和获得高精度纠缠熵。

2.3 性能数据与瓶颈转移

传统的量子蒙特卡洛方法（如无偏的辅助场 BSS-QMC）在计算 Rényi 纠缠熵时，计算复杂度随子区域周长 $L_{\partial A}$ 以及 $U$ 和 $\beta$ 呈指数级增长。最核心的障碍是符号问题（Sign Problem）。而本算法的性能特征呈现出完全不同的格局：

1. 无符号问题困扰：由于 DiagMC 在虚时间-空间连续表象中对已求和的连通 Feynman 图进行采样，蒙特卡洛采样的权函数不依赖于多体状态的相位，因而从根本上规避了由于热力学极限下波函数反对称所导致的指数符号墙。
2. 瓶颈转化为内存限制（Memory Bottleneck）：随着晶格尺寸 $N_1 \times N_2$ 增大以及子区域 $A$ 形状变得复杂，计算并存储格林函数张量 $\bar{G}^{(n)}_{\alpha\beta}(i, \tau; j, \tau')$ 需要消耗巨量内存。对于 $40 \times 40$ 的格点，虚时间网格精细度（由于 Chebyshev 阶数高达 $25 \times 25$）会导致所需存储的中间体数据达到数十吉字节（GB）。为了防止高频调用的内存溢出，本工作设计了极其巧妙的LRU（Least Recently Used）缓存机制与分块表（Chunked Table）构建方案。这使得物理学家能够根据可用内存动态调节计算开销，从而将不可逾越的“计算时间符号墙”转化为可通过硬件升级或软件优化解决的“内存带宽多项式问题”。

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

3.1 算法架构设计

一个完整的连接行列式图符蒙特卡洛（CDet-DiagMC）代码由以下四个核心模块构成：

1. 自由格林函数解算器（Free G-Solver）：

   • 输入：格点尺寸 $N_1 \times N_2$， hopping $t$，化学势 $\mu$，温度 $\beta$，子区域 $A$ 的几何定义。
   • 步骤：
     1. 构造单粒子 Hamilton 矩阵 $h$。
     2. 对 $h$ 进行精确对角化（以利于后续矩阵指数 $e^{-\tau h}$ 的高效计算）。
     3. 针对每一个复制品动量 $q_m = (2m+1)\pi / n$，构造边界条件矩阵 $\Omega_m = P_B - e^{iq_m} P_A$。
     4. 求解关键的中间矩阵 $A_m = [\mathbf{I} + \Omega_m^{-1} e^{-\beta h}]^{-1}$。
     5. 计算各 Fourier 扇区格林函数 $\bar{G}^{(m)}(\tau, \tau')$ 并通过傅里叶反变换合并，将结果存入高速缓存中。

2. Chebyshev 时间插值模块（Chebyshev Interpolator）：

   • 虚时间 $\tau$ 跨越 $0$ 到 $\beta$，为了避免在蒙特卡洛步骤中实时进行矩阵指数乘法，需预先在 Chebyshev 节点上计算好格林函数并构建多维插值表。

3. CDet 递归引擎（Connected Determinant Recursion Engine）：

   • 采用高效的子集递推。设当前顶点集为 $V$。对于任意子集 $S \subseteq V$，根据下述递归式提取连通部分阻尼权重 $C(S)$：
   $$C(S) = D(S) - \sum_{S' \subset S, v_0 \in S'} C(S') D(S \setminus S')$$

   其中 $D(S)$ 是该子集对应的满行列式，通过标准高斯消元或 LU 分解在 $O(|S|^3)$ 时间内完成。固定某一特征顶点 $v_0$ 可以有效避免重复计数。

4. MCMCMC 采样器（Many-Configuration Sampler）：

   • 采用信封配置（Envelope Configuration） $V = (X_1, \dots, X_N)$，使多个微扰阶数在同一次马尔可夫链中被同时采样。包含增加/减少顶点、修改顶点时空坐标等更新步骤。

3.2 复现指南与伪代码流程

以下是实现核心格林函数计算及 CDet 部分的关键伪代码流程：

# 核心算法伪代码：求解复制品格林函数
import numpy as np
from scipy.linalg import expm, inv

def compute_replica_green_function(N_sites, t, mu, beta, subregion_A, n_replica, tau_1, tau_2, alpha, beta_index):
    # 1. 构造单粒子哈密顿量
    h = construct_tight_binding_hamiltonian(N_sites, t, mu)
    
    # 2. 构造投影算符
    P_A = np.zeros((N_sites, N_sites))
    for idx in subregion_A:
        P_A[idx, idx] = 1.0
    P_B = np.eye(N_sites) - P_A
    
    # 3. 循环计算扇区格林函数
    G_m = []
    for m in range(n_replica):
        qm = (2 * m + 1) * np.pi / n_replica
        Omega_m = P_B - np.exp(1j * qm) * P_A
        
        # 计算 A_m = [I + Omega_m^-1 * exp(-beta * h)]^-1
        exp_beta_h = expm(-beta * h)
        inv_Omega = inv(Omega_m)
        A_m = inv(np.eye(N_sites) + inv_Omega @ exp_beta_h)
        
        # 根据 tau_1 与 tau_2 的大小计算特定时间的格林函数值
        if tau_1 > tau_2:
            G_temp = expm(-tau_1 * h) @ A_m @ expm(tau_2 * h)
        else:
            G_temp = expm(-tau_1 * h) @ (A_m - np.eye(N_sites)) @ expm(tau_2 * h)
        G_m.append(G_temp)
        
    # 4. 执行反傅里叶变换重建实体空间格林函数
    G_real = np.zeros((N_sites, N_sites), dtype=complex)
    for m in range(n_replica):
        qm = (2 * m + 1) * np.pi / n_replica
        phase = np.exp(1j * qm * (beta_index)) # beta_index 为复制品指数差 
        G_real += G_m[m] * phase
    return G_real / n_replica

3.3 推荐开源软件包及 Repo 链接

虽然本论文作者暂未公开其针对此纠缠熵计算的专用私有生产代码，但由于该方法深度基于常规的连接行列式（CDet）蒙特卡洛框架，读者可以通过参考并修改以下著名开源 DiagMC 与多体计算库来实现复现：

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

   • 简介：由法国原子能署（CEA）和弗拉特龙研究所（Flatiron Institute）共同维护的高性能多体量子计算框架，内置极佳的 C++ 矩阵库与虚时间格林函数插值算法，并提供 Python 接口。
   • GitHub 链接：https://github.com/TRIQS/triqs

2. CDet-inspired Repositories by Riccardo Rossi

   • 简介：连接行列式算法的原创提出者 Riccardo Rossi 公开过关于高阶图符蒙特卡洛和 CDet 的核心算法例程，这是理解微扰级数 subset recursion 逻辑的最佳范本。
   • 参考学术检索：检索 GitHub 上关于 CDet 或 Diagrammatic Monte Carlo 的开源实现（例如：github.com/riccardorossi 相关多体微扰公开仓库）。

3. ALPSCore (Algorithms and Libraries for Physics Simulations)

   • 简介：专为经典与量子蒙特卡洛模拟设计的高性能核心计算库，提供优秀的马尔可夫链更新、自动误差分析（Autocorrelation analysis）及数据序列化功能。
   • GitHub 链接：https://github.com/ALPSCore/ALPSCore

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

4.1 关键引用文献分析

本研究的成功确立在以下几项里程碑式研究的肩膀上：

1. 连接行列式算法的提出：
   • R. Rossi, Phys. Rev. Lett. 119, 045701 (2017) [文献 [25]]。这篇文献首次提出了 Connected-Determinant (CDet) Diagrammatic Monte Carlo，彻底改变了 DiagMC 依赖于漫长 Feynman 图标记搜索的现状，使高阶展开在计算复杂度上变得高度可行。
2. 多组态 MCMC 采样法：
   • F. Šimkovic and R. Rossi, arXiv:2102.05613 [文献 [26]]。该工作提供了在一次蒙特卡洛运行中同时提取多个不同微扰阶数系数的“信封配置”采样技术（MCMCMC），极大地加速了本算法中高阶修正系数 $c_k$ 的收敛性能。
3. 费米子纠缠熵的边界条件构建：
   • T. Grover, Phys. Rev. Lett. 111, 130402 (2013) [文献 [8]]。Grover 探讨了在辅助场蒙特卡洛（DQMC）中如何利用格林函数矩阵构造 Rényi 纠缠熵的方法。本研究可以看作是将这种空间边界设计推广到了微扰图符展开的无限连续虚时间空间中。

4.2 局限性与前沿评论

尽管本工作在规避符号问题、实现大体系纠缠熵计算方面取得了里程碑式的进展，但在面对更复杂的实际强关联物理场景时，仍暴露出以下不容忽视的局限性：

1. 极强关联区域的收敛发散难题（The Convergence Limit in Strong-Coupling Regime）：

   • 本质上，本方法是一个基于温和同伦路径下的微扰级数展开法。尽管采用了 CDet 和很多加速求和技术，但对于极强关联区域（如典型的 $U/t > 8$ 且处于掺杂区域的 2D Hubbard 模型），其微扰级数可能会变成非渐近收敛的。在这类区域，Dyson 级数的复数奇点可能极度靠近实轴，导致即便求到第 10 阶甚至更高阶，级数求和也会发生剧烈震荡，必须依赖极为精细的解析延拓技术（如博雷尔求和 Borel Summation 或 conformal maps），这会带来难以控制的系统误差。

2. 内存消耗随温度与尺度的爆炸性增长（The Memory Wall at Low Temperatures）：

   • 为了在低温（大 $\beta$）下保证切比雪夫插值精度，切比雪夫多项式的阶数必须相应提高（如从 $25 \times 25$ 提升到 $100 \times 100$）。格林函数张量的数据尺寸正比于 $N_{\text{replica}}^2 \times N_{\text{sites}}^2 \times N_{\text{Cheby}}^2$。在超大系统下，即使有 LRU 缓存，频繁的磁盘读写与缓存失效重构也会导致计算时间显著变长。内存带宽成为了限制该方法走向 $100 \times 100$ 格点极低温物理研究的最大绊脚石。

3. 对 Rényi 指数 $n$ 的泛化困难：

   • 本文主要探讨了第二 Rényi 纠缠熵 $S_2$（$n=2$）。虽然理论上可以自然推广到更高的整数 $n$（如 $S_3, S_4$ 等），但随着 $n$ 的增大，复制品动量 $q_m$ 的数量增多，边界条件矩阵 $\Omega_m$ 的复相因子求逆变得更加频繁，计算耗时将线性增加。更重要的是，物理学家最关心的往往是冯·诺依曼纠缠熵（Von Neumann Entropy, 即 $n \to 1$ 极限）。如何在图符蒙特卡洛框架下平滑地进行 $n \to 1$ 的解析延拓，目前在理论上仍是一个悬而未决的开放问题。

5. 其他必要补充：物理洞察与未来展望

5.1 费米子渐变交换宇称的深层拓扑学诠释

本研究最令人瞩目的理论亮点在于其边界条件设计中对 Koszul 符号 完美的代数处理。在统计物理中，由于多粒子波函数的波色/费米统计性质不同，其多叶 Riemann 面的互连方式在代数结构上存在本质差异。在玻色子系统里，跨页过渡是平凡的。而在费米子系统中，非局域交换算符带有天然的 Jordan-Wigner 弦（String）。

本研究通过引入半移动动量 $q_m = (2m+1)\pi/n$，在数学上相当于在复制品空间中注入了一个大小为 $\pi$ 的有效规范磁通（Effective Gauge Flux）。正是这个规范磁通在虚时间演化边界上自动补偿了费米子交换所产生的负号。这个物理图像极其优美：量子统计带来的波函数相位差，可以在多叶复制空间中被等效转化为一个几何规范场带来的 Aharonov-Bohm 效应。这为我们理解多体费米子纠缠的代数本源提供了崭新的物理视角。

5.2 探测强关联多体物理新效应的有力武器

借助本方法能够处理大尺度、大子区域纠缠熵的独特优势，未来可以在以下几个极具争议的凝聚态物理前沿方向大显身手：

1. 赝能隙（Pseudogap）本质与费米面重构（Fermi Surface Reconstruction）的诊断：

   • 在掺杂的 2D Hubbard 模型中，随着空穴掺杂浓度的变化，系统会经历从大费米面到小费米口袋（Fermi Pockets）的重构。由于费米子纠缠熵的对数违背系数直接正比于费米面的投影面积（$S_2 \sim L_A \ln L_A$），通过精确测定不同温度和掺杂下的纠缠熵对数违背系数值，可以判定系统在特定温区是否发生了非费米液体行为，并定量给出费米面重构的临界掺杂点。

2. 量子自旋液体（Quantum Spin Liquids）与拓扑序的判定：

   • 虽然本文以 repulsive Hubbard 模型为例，但该方案可以无缝推广到具有强自旋-轨道耦合的 Kitaev 相互作用系统。在存在拓扑序的二维费米子体系中，纠缠熵含有不依赖于边界周长的常数修正项——拓扑纠缠熵（Topological Entanglement Entropy, $\gamma$）。本算法超强的无符号大尺寸模拟能力，使得在足够大的几何结构下准确外推并分离出这个极微小的拓扑常数成为可能，从而提供量子自旋液体最直接、最无可争议的数值证据。

5.3 结语

Boyuan Shi 的这项工作是一次完美的理论物理形式美学与高性能数值算法的联姻。它不仅证明了我们能够超越传统的配置空间蒙特卡洛的视野，转而在高阶微扰的“虚空图符空间”中自由驰骋，还指明了通过精细设计代数结构，费米子符号问题完全可以被局部化、多项式化。在量子计算和冷原子气体显微镜技术（Quantum Gas Microscopy）突飞猛进的今天，这种高精度的数值基准计算将为我们在实验室中操纵、理解和利用量子多体纠缠提供最坚实、最可靠的理论引路灯塔。