来源论文: https://arxiv.org/abs/2607.00134v1 生成时间: Jul 02, 2026 18:21

无对称性破缺的吉布斯态相区分：三维拓扑码有限温度拓扑相变与新型通道不变量的深度解析

0. 执行摘要

拓扑序（Topological Order）在零温（$T=0$）下的稳定性与分类理论已相对完善，但在有限温度（$T>0$）的吉布斯热态（Gibbs States）中，热涨落如何影响拓扑相、以及如何定义和测量热态中的拓扑序，一直是凝聚态物理与量子信息交叉领域的核心难题。传统的朗道对称性破缺范式无法完整描述拓扑相变。更严峻的是，在实际物理系统中，任何微小的外加磁场都会显式地破缺保护拓扑相的广义高阶对称性（Higher-Form Symmetries），这使得有限温度拓扑序的生存空间饱受质疑。

近期，Haruki Watanabe 教授在预印本论文中发表了一项里程碑式的研究，深入探讨了三维 $\mathbb{Z}_2$ 拓扑码（3D Toric Code）在任意磁场下的有限温度拓扑相变。该研究的核心贡献可以概括为以下三个方面：

1. 几何保护机制的发现：研究表明，即使任意高阶对称性被外场显式破缺，3D $\mathbb{Z}_2$ 拓扑码在有限温度下的拓扑相变依然存在。这种稳定性并非由某种具体的哈密顿量对称性保护，而是本质上由晶格本身的几何约束——**毕安基恒等式（Bianchi Identity）**所保护。在去禁闭相（Deconfined Phase）内，拓扑纠缠熵（Topological Entanglement Entropy, TEE） $\gamma$ 完美地保持在 $\ln 2$ 的量子化平台，并在热相变点 $T_c$ 骤降至 0。通过大尺度量子蒙特卡洛（QMC）模拟与有限尺寸标度（FSS）分析，证实该相变属于 3D Ising 普适类。
2. 拓扑纠缠熵的不变性危机：论文通过显式的算符构造证明，随着系统向混合态演化，传统的本体拓扑纠缠熵 $\gamma$ 并不是**拟局部通道（Quasi-Local Channels, QLCs）**下的不变量。一个常数深度的完全正保迹（CPTP）通道可以从一个完全平庸的乘积态中“无中生有”地制造出 $\ln 2$ 的拓扑纠缠熵。这一结论与近期 Yang、Shi 和 Lee 的独立研究遥相呼应，彻底颠覆了“本体熵可以直接作为混合态相不变量”的传统认知。
3. 解码威尔逊圈关联（Decoded Wilson-loop Correlation, $f_W$）的引入：为了克服 $\gamma$ 的非通道不变性缺陷，作者创新性地引入了纠错解码的思想，定义了 $f_W$。该物理量在拓扑去禁闭相中量子化为 1，在平庸相（限制相与希格斯相）中为 0（在热力学极限下）。最关键的是，$f_W$ 具有拟局部通道不变量性。这是物理学界首次提出的、可直接通过大尺度量子蒙特卡洛计算的有限温度通道不变拓扑序参量。

这一研究不仅澄清了热吉布斯态在“三层态分类模型”中的尴尬地位，也为三维自校错量子存储器（Self-Correcting Quantum Memory）的物理机制提供了精确的统计物理图景。

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

1.1 核心科学问题：外场与热涨落双重打击下的拓扑相存活问题

在二维（2D）拓扑码中，任何有限温度 $T>0$ 的扰动都会导致点状的任意子（Anyons）在热涨落驱使下发生扩散与无规行走，从而将系统存储的拓扑信息彻底擦除。这意味着 2D 拓扑码在有限温度下的拓扑纠缠熵会瞬间坍塌，拓扑序仅存在于 $T=0$ 基态。

而在三维（3D）拓扑码中，激发态呈现出不对称的结构：电荷（Electric Charges，对应顶点算符 $\hat{A}_v$ 激发）依然是点状的；而磁通（Magnetic Fluxes，对应面算符 $\hat{B}_p$ 激发）则变成了**线圈状（Loop-like）**的激发。正是因为磁通激发具有“线张力”（Line Tension），在有限温度下它们不会轻易地任意增长并穿过整个系统，只有当温度超过某个临界温度 $T_c$ 时，这些磁通线圈才会发生“汽化”和渗透（Percolation），导致拓扑序崩溃。因此，无磁场时的 3D 拓扑码在有限温度下能够稳定存在拓扑去禁闭相。

然而，当我们加入外磁场（$h^x$ 与 $h^z$）时，问题变得极其棘手：

• $h^x$ 显式破缺了磁 $\mathbb{Z}_2$ 二形式（2-form）对称性，并产生磁通线圈激发；
• $h^z$ 显式破缺了电 $\mathbb{Z}_2$ 一形式（1-form）对称性，并产生电荷激发。

在物理学传统观点中，对称性一旦显式破缺，由这些对称性定义和保护的拓扑相就应该不复存在。那么：

1. 在没有高阶对称性保护的情况下，3D $\mathbb{Z}_2$ 拓扑码的有限温度拓扑相变是否依然保持为一个不可逾越的、锐利的物理相边界（即非终止边界，Non-terminating boundary）？
2. 在热混合态分类（CPTP 通道等价类）的框架下，我们应该用什么数学和物理量来精确刻画这个热态拓扑相，使其既能对磁场/电场扰动保持鲁棒，又不能被局部局域物理过程伪造？

1.2 理论基础与三层态分类架构

为了给上述问题提供严谨的回答，必须首先理清有限温度吉布斯态在量子物理状态分类体系中的位置（见论文 Fig. 1）：

1. 第一层：基态（纯态，Ground States）：等价关系为有限深度局部酉变换（Finite-Depth Local Unitary, FDLU）。在此框架下，拓扑序由基态简并度、非阿贝尔规范结构、以及基态拓扑纠缠熵（TEE） $\gamma$ 完全刻画。
2. 第二层：混合态（Mixed States）：等价关系为双向拟局部通道（Two-way Quasi-Local Channels, QLCs）。只要两个混合态可以通过有限深度的局部量子通道（CPTP）相互转化，它们就属于同一个混合态相。最近的研究表明，这种分类过于粗糙（例如它能把去相干的拓扑态与平庸乘积态等同起来），因而需要引入带有局部可恢复性（Local Recoverability）约束的拓扑通道连通性。
3. 第三层：吉布斯态（Gibbs States, $\rho_\beta = e^{-\beta \hat{H}}/Z$）：介于前两者之间。目前的共识是，如果在连续改变哈密顿量参数和温度的过程中，系统不跨越任何热力学奇点（即关联长度或马尔可夫长度不发散），则这两个吉布斯态属于同一个热态相。由于温度的变化不能由 FDLU 算符直接生成（因为它要重新分配能级权重，属于非酉的虚时演化过程），吉布斯态的相分类难度极大。

1.3 技术难点：毕安基恒等式的几何本征约束与 TEE 的“伪造”危机

毕安基恒等式的约束力

3D 拓扑码哈密顿量中，每一个面（Plaquette）关联着一个规范算符 $\hat{B}_p = \prod_{l \in \partial p} \hat{\sigma}^z_l$。对于晶格中的任意一个基本立方体（Elementary Cube） $c$，由于每个棱边（Link） $l$ 恰好被相邻的四个面所共享，在同一个立方体的六个面求积时，每个 $\hat{\sigma}^z_l$ 均会出现两次。因此我们有如下的**毕安基恒等式（Bianchi Identity）**算符恒等关系：

这是一个不依赖于任何物理哈密顿量、任何磁场强度的严格晶格几何代数恒等式。它保证了磁通线圈在三维空间中绝不可能暴露出点状的“单极子”（Monopoles）。磁通激发必须以闭合圈的形式存在。这就是拓扑相能够在有限温度下抵御任意磁场扰动、保持锐利相边界的根源：它是几何保护的，而非对称性保护的。

拓扑纠缠熵的伪造机制

尽管毕安基恒等式确保了本体物理相边界的存在，但在定量刻画上，著名的 Levin-Wen / Kitaev-Preskill 拓扑纠缠熵 $\gamma$ 在有限温度下遇到了巨大的理论挑战。作者通过显式的量子通道构造，证明了 $\gamma = \ln 2$ 可以被非拓扑的微观操作完全“伪造”出来。

具体地，考虑三维立方晶格上的平庸初始乘积态（所有自旋处于 $z$-方向本征态）：

显然该状态的 TEE $\gamma = 0$。现在定义一个常数深度的量子通道（QLC），其操作如下：

1. 在每个晶格顶点 $v$ 上附加一个处于最大混合态的辅助比特（Ancilla Qubit）：$\rho^{\text{anc}}_v = \frac{1}{2}\mathbb{I}$，其计算基状态可视为一个经典随机硬币 $\tau_v = \pm 1$；
2. 对于每一个棱边 $l = \langle vv' \rangle$，应用受控非门（CNOT），将两端顶点 $v, v'$ 上的辅助比特作为控制位，棱边 $l$ 上的系统比特作为目标位。在 $\sigma^z$ 基下，该操作将棱边的自旋状态映射为两端硬币的乘积：$z_l \to \tau_v \tau_{v'}$；
3. 抛弃（Trace out）所有的辅助比特。

由于该量子通道的所有操作都是完全局部的、且各门算符对易，这显然是一个深度为 $O(1)$ 的完全正保迹拟局部通道。经过该通道作用后，系统的状态变为：

这恰好是无通量经典随机闭合圈气体（Flux-free classical loop gas）的均匀混合态！对于该状态，由于每个面上的 $\hat{B}_p$ 均满足 $\prod_{l \in \partial p} z_l = \prod_{l \in \partial p} \tau_v \tau_{v'} = +1$，这完美符合 3D 拓扑码去禁闭相的局域特征。利用 Levin-Wen 的四区域纠缠熵组合（Fig. 5 中的环形 donut 方案）对其计算纠缠熵，由于环形区域包含非平凡同调圈的闭合条件，计算结果为：

这一构造惊人地表明：一个平庸的乘积态，仅需经过一轮常数深度的局部量子通道操作，就能在宏观上产生与拓扑相一模一样的、量子化的 $\ln 2$ 拓扑纠缠熵！ 这意味着，本体拓扑纠缠熵 $\gamma$ 在混合态的分类框架中无法单独作为区分子系统拓扑性质的“不变量”。它不是拟局部通道不变量。

1.4 方法细节：解码威尔逊圈关联 $f_W$ 的构建与物理实质

为解决上述危机，必须引入全局相干的性质。作者指出，虽然局部通道可以伪造局域面算符的零通量约束，但它绝对无法在环面上建立起宏观的、具有同调非平凡性质的全局规范环关联。

因此，作者设计了“解码威尔逊圈关联”（Decoded Wilson-loop Correlation） $f_W$：

其中，$\bar{W}_1, \bar{W}_2$ 是在三维环面空间中，两个沿着同调方向 $a$、相互平行的非平凡大圈 $C_1, C_2$ 上的 bare 威尔逊圈算符：

直接测量裸威尔逊圈算符 $\bar{W}_i$ 会因为热激发的电荷和磁通涨落而导致其期望值以面积律衰减。为了提取真正的拓扑关联，必须首先对测量得到的自旋快照进行通量清洁（Flux Cleaning，即解码纠错）。

通量清洁通道 $\mathcal{C}$ 的精确定义

解码过程可以形式化为一个量子通道 $\mathcal{C}$：

其中，$\eta = \{\eta_p = \pm 1\}$ 是系统中的磁通缺陷结构（即被违反的面算符，满足 $\hat{B}_p = -1$）。投影算符 $\hat{P}_\eta = \prod_p \frac{\mathbb{I} + \eta_p \hat{B}_p}{2}$ 提取通量缺陷分布。根据毕安基恒等式，这些缺陷点在对偶晶格上必然构成闭合的环线。解码纠错算符 $\hat{X}_{R(\eta)}$ 则是应用最小权重完美匹配（Minimum-Weight Perfect Matching, MWPM）算法，在对偶晶格上寻找一个以缺陷环为边界的最小权重膜表面 $R(\eta)$，并在该表面上施加 $\hat{\sigma}^x_l$ 操作，从而将缺陷对消，使系统重回零缺陷的去禁闭流形：

解码威尔逊圈算符 $W_i$ 即为裸算符在清洁通道下的海森堡绘景对偶作用：$W_i = \mathcal{C}^\dagger (\bar{W}_i)$。 由于 $\bar{W}_i$ 与所有的投影算符对易，解码后的两点关联在每个快照上被简化为一个经典的布尔二值判定问题：

其中 $P_{\text{fail}}$ 指的是纠错码解码器错误地将测量得到的线圈划分到了错误的同调类（Homology Class）中的概率。

同时，由于对称性破缺（特别是外加 longitudinal 场 $h^z$ 会打破系统的膜对称性），我们必须减去单点期望值的平方。单点期望值 $m_a = \langle W_i \rangle$ 对应于系统的同调磁化强度（Holonomy Magnetization）：

$f_W$ 对三个物理相的完美响应图景

1. 拓扑去禁闭相：热激发涨落极其稀疏，解码器工作极度精准，$P_{\text{fail}} \to 0$；由于电荷未凝聚，同调磁化未建立， $m_a = 0$。因此 $f_W \to 1$。
2. 磁平庸限制相：磁通发生凝聚，缺陷线圈无处不在。解码器彻底失效，判定结果如同抛硬币，导致 $P_{\text{fail}} \to 1/2$。由于无电荷凝聚，$m_a = 0$。因此 $f_W \to 0$。
3. 电平庸希格斯相： longitudinal 场 $h^z$ 导致电荷发生凝聚。虽然由于缺陷稀疏，解码器的保真度很高（$P_{\text{fail}} \to 0$），但电荷凝聚强行将系统的同调状态锁定在了特定方向，使得同调磁化饱和：$m_a \to \pm 1$。因此 $f_W \to (1 - 0) - 1^2 = 0$。

由此可见，$f_W$ 在三种不同的物理相界中表现出了惊人的灵敏度与区分能力。更为重要的是，由于 $f_W$ 在定义上严格执行了连接关联函数的差值，任何通过局部通道制造的平庸伪造态（如前述 $\rho_{\text{out}}$）由于没有长程相干，其单点磁化必然与两点关联完全抵消（即 $\langle W_1 W_2 \rangle = m_a^2 = 1$），从而产生 $f_W = 0$。这铁证了 $f_W$ 的通道不变量鲁棒性。

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

为了验证上述理论和 $f_W$ 的实际测量能力，作者基于大尺度量子蒙特卡洛（QMC）模拟，针对三维简立方晶格上的拓扑码体系开展了系统性的计算。计算选取的体系尺寸参数为 $L \in \{6, 8, 10, 12, 14, 16, 18\}$，周期性边界条件下的自旋数（棱边数）为 $N = 3 L^3$。

2.1 经典基准体系：无场极限点 $(h^x, h^z) = (0,0)$

在无磁场极限下，3D 拓扑码可以通过 Wegner 对偶（Wegner Duality）严格映射到经典的 3D Ising 模型上。这为量子蒙特卡洛算法的精度提供了一个无参数可调的、极高精度的基准系统。在该极限点，理论预测的具体热容 $C_v$ 表现出电、磁两个完全解耦的物理扇区（见附录 A 2）：

1. 电扇区 $C_v^A$（顶点算符）：表现为平滑的经典肖特基异常（Schottky Anomaly），无任何热力学相变，峰值位于 $T \approx 0.833$ 处： $$\frac{C_v^A}{N} \approx \frac{1}{3} \beta^2 \text{sech}^2 \beta$$
2. 磁扇区 $C_v^B$（面算符）：精确继承 3D Ising 模型的临界相变特征。根据对偶变换，其精确相变温度为： $$T_c^{(0,0)} = [\text{arctanh}(e^{-2 J_c})]^{-1} \approx 1.31335$$ 其中 $J_c \approx 0.2216544$ 为 3D Ising 模型的临界耦合常数。

量子蒙特卡洛对基准点的复现精度（见论文 Fig. 10）

作者在 QMC 算法中引入了两个至关重要的修正：算符数纠正（Operator-count correction）和有限长链自相关去偏（Finite-chain autocorrelation de-bias）。

• 峰值位置精度：对于 $L \ge 10$ 体系，QMC 提取的具体热容峰值位置 $T_p(L)$ 与经典精确极限的偏差小于万分之七：$|\Delta T_p| \le 7 \times 10^{-4}$。
• 峰值高度精度：在所有模拟尺寸下，QMC 测得的峰值高度与经典解析曲线的偏差控制在 $1.5\%$ 以内（等价于 1 倍标准差以内 $1 \sigma$）。
• 临界指数拟合：在无场点，自由三参数拟合得到的临界温度和指数为： $$T_c(\infty) = 1.3168 \pm 0.0008, \quad \nu = 0.567 \pm 0.025$$ 这与标准 3D Ising 普适类临界指数 $\nu_{\text{Ising}} \approx 0.6301$ 具有极佳的相容性，验证了模拟方法的可靠性。

2.2 扰动物理体系：代表性测试点 $(h^x, h^z) = (0.5, 0.1)$

该参数点处于强物理扰动区（外加磁场已经显式破缺了所有的系统哈密顿量对称性），是检验几何保护拓扑序相变是否依然存在的终极试金石。通过具体热容 $C_v$ 的有限尺寸标度标定（论文 Fig. 10(b)），在引入磁场扰动后，相变点发生了红移：

拟合得到的关联长度指数为 $\nu = 0.573 \pm 0.040$，依然完全符合 3D Ising 普适类（无任何一阶相变不连续性特征），表明磁场只改变了相边界的位置，并未瓦解相变的物理性质。

2.3 核心不变量性能比对：$\gamma$ 平台 versus $f_W$ 分辨力

为了直观展示这两个拓扑参量在热力学极限下的性能，下表汇总了论文在三个正交参数扫描切片（T-scan, $h^x$-scan, $h^z$-scan）中得到的关键数值结果（基于 $L=8$ 生产数据集与 $L=10, 12, 14$ 有限尺寸验证系列）：

特别性能数据：高场下的有限尺寸几何限差（Aperture Limit）

在强场扰动测试点 $(h^x, h^z) = (0.75, 0.15)$ 下，纠缠熵标度曲线（论文 Fig. 11）展示了一个极具启发性的物理限制： 在 $T=0.5$ 下，$\gamma(L)$ 分别输出为 $\gamma(8) = 0.625(9)$，$\gamma(10) = 0.612(15)$，$\gamma(12) = 0.621(19)$，$\gamma(14) = 0.615(28)$。这些值高度稳定地收敛于一个略低于 $\ln 2$ 的非平凡平台（约合 $0.89 \ln 2$）。

这并非意味着物理拓扑序发生了消退，其本质是由 Levin-Wen 探测几何的“固定孔径限制”（Fixed Aperture Limit） 造成的：由于蒙特卡洛采样中探测环的几何尺寸在物理上是保持不变的（其内径 $R_{\text{in}}=1.25$、外径 $R_{\text{out}}=2.75$ 固定），当外场增加导致系统真实的物理关联长度 $\xi$ 变大并逐渐逼近探测环的几何线宽时，由于边界处局域涨落泄漏导致区域间的消去不完全，从而系统性地下调了外推的本体拓扑值。这证实了：在大场下利用纠缠熵提取拓扑特征时，尺寸的不匹配会导致系统性的偏差。相比之下，$f_W$ 则是通过物理纠错（解调）在同调层面上进行判定，完全免受这种局域几何泄漏的影响。

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

3.1 核心算法：连续时间世界线量子蒙特卡洛（CT-WLQMC）

为了处理含有横场与纵场的量子自旋晶格系统，作者采用了无符号问题（Sign-problem-free）的连续时间世界线蠕虫算法（Continuous-Time Worldline Worm Algorithm）。在世界线表示中，配分函数 $Z = \text{Tr} e^{-\beta \hat{H}}$ 被展开为关于局部哈密顿量非对易项插入的世界线路径积分形式（见论文附录 C 1）：

其中，系统的状态由虚时演化 $\tau \in [0, \beta)$ 上的自旋配置文件 $s(\tau)$ 完全确定。空间中的顶点算符 $\hat{A}_v$ 对应于六条系统棱边上的自旋协同翻转事件，而横场 $h^x$ 则对应于单自旋在虚时轴上的翻转。由于哈密顿量中所有非对角项的矩阵元均为非负：

这保证了世界线权重 $W(\omega) > 0$ 处处为正，不存在任何量子相位带来的负符号问题，从而可利用马尔可夫链蒙特卡洛进行极高精度的概率采样。

3.2 开源软件包资源

本研究大尺度计算的核心程序基于以下优秀的开源软件生态进行开发和扩展：

1. QMC 引擎：ParaToric
   这是本研究所使用的、高度优化的三维拓扑码世界线蒙特卡洛采样专用软件包。由 S. M. Linsel 和 L. Pollet 共同开发。
   • 开源仓库链接：https://github.com/sm-linsel/ParaToric (或者参考 SciPost Codebases 官方发布)
   • 主要特点：采用现代 C++20 标准编写，提供了高效的蠕虫更新（Worm Updates）以及专门针对拓扑码多体顶点对易子设计的局部块更新（Block Updates）。
2. 纠错解码核心：PyMatching
   用于计算 $P_{ ext{fail}}$ 判定中的最小权重完美匹配（MWPM）算法。由 Oscar Higgott 编写，底层的图匹配算法经过了高度的 C++ 优化。
   • 开源仓库链接：https://github.com/quantumlib/PyMatching
   • 作用：读取 QMC 输出的虚时 $\tau=0$ 切片上的通量缺陷三维坐标，将其投影到二维切片中构建综合纠错图，计算非平凡纠错回路的拓扑同调。

3.3 数值复现指南与核心算子伪代码

为了复现论文中的核心结论（即在吉布斯态中计算 $f_W$ 以及拓扑纠缠熵 $\gamma$），研究人员可以按照以下步骤搭建复现管道（Pipeline）：

第一阶段：QMC 热平衡采样与快照生成

1. 初始化三维周期性立方晶格（尺寸 $L=8$，棱边上放置比特，共 $3 \times 8^3 = 1536$ 个物理自旋）。
2. 设定哈密顿量参数 $(J_e = J_m = 1, h^x = 0.5, h^z = 0.1)$，以及目标倒温度 $\beta = 1/T = 2.0$。
3. 运行 ParaToric 采样器，执行不少于 $10^7$ 步热平衡扫掠（Thermalization sweeps）。
4. 在平衡后，每隔 $N_{bs} = 200$ 步测量扫掠，记录 $\tau=0$ 虚时切片上所有物理自旋在 $\sigma^z$ 计算基下的本征值快照 $s \in \{\pm 1\}^{1536}$。至少收集 $N_s = 5000$ 个独立的自旋位形。同时记录所有面算符的本征值状态 $\{b_p \in \{\pm 1\}\}$。

第二阶段：最小权重完美匹配解码计算 $P_{ ext{fail}}$

对于获得的每个自旋快照，其上被破坏的面算符（$\hat{B}_p = -1$）的位置即为三维空间中的磁缺陷分布。由于计算 $f_W$ 需要沿着二维剖面进行解码分析，我们将缺陷分布投影至与其垂直的 $xy$ 剖面上，得到一个二维平面上的零维缺陷集合（即点状缺陷）。

以下是利用 PyMatching 进行解码计算的 Python 核心逻辑：

import numpy as np
import pymatching

def calculate_decoded_wilson_loop_correlation(snapshots, L):
    """
    复现 3D 拓扑码解码威尔逊圈关联 f_W 的核心算法
    snapshots: 包含自旋本征值和面本征值快照的数据集
    L: 晶格线性尺寸
    """
    # 1. 初始化纠错图 (Decoding Graph)
    # 节点对应于二维对偶晶格点，边对应于对偶晶格的边（跨越系统物理棱边）
    proj_graph = pymatching.Matching()
    
    # 构建二维周期性对偶格点，定义其边权重与纠错纠缠概率相关
    for x in range(L):
        for y in range(L):
            # 物理上对应的两个水平/垂直棱边
            edge_id_h = get_edge_index_h(x, y, L)
            edge_id_v = get_edge_index_v(x, y, L)
            
            proj_graph.add_edge(node1=get_node_id(x, y, L), 
                                 node2=get_node_id((x+1)%L, y, L), 
                                 qubit_id=edge_id_h, 
                                 weight=1.0)
            proj_graph.add_edge(node1=get_node_id(x, y, L), 
                                 node2=get_node_id(x, (y+1)%L, L), 
                                 qubit_id=edge_id_v, 
                                 weight=1.0)

    logical_outcomes = []
    
    for snap in snapshots:
        # 2. 提取当前快照下的缺陷分布（Syndrome）
        syndrome = snap["plaquettes"] == -1
        # 降维：将三维通量闭合线投影为二维平面点状缺陷
        projected_syndrome = project_syndrome_to_2D(syndrome, L)
        
        # 3. 使用 MWPM 寻找最优校正链
        prediction = proj_graph.decode(projected_syndrome)
        
        # 4. 判断校正后在大圈 C_1、C_2 上的非平凡同调绕数（Winding Number）
        # 若预测校正算符的逻辑类与实际测量自旋态不符，则判定纠错失败
        actual_winding = compute_actual_winding(snap["spins"], L)
        decoded_winding = actual_winding ^ np.sum(prediction[logical_qubits_along_loop]) % 2
        
        logical_outcomes.append(1 if decoded_winding == 0 else -1)
        
    # 计算 1 - 2 * P_fail
    two_point_fW = np.mean(logical_outcomes)
    return two_point_fW

第三阶段：通过热力学积分计算同调磁化强度 $m_a$

由于在拓扑去禁闭相中，系统的对称性内存势垒高达 $e^{c \beta L}$，任何局域更新算法都会被锁死在特定的同调分支中，无法直接采样 $m_a$。必须通过**热力学积分（Thermodynamic Integration）**计算自由能之差 $\Delta F$：

1. 在哈密顿量中，针对大膜 $\Sigma$ 上的电场项施加连续参数 $\lambda$： $$\hat{H}_{\lambda} = \hat{H} + 2 \lambda h^z \sum_{l \in \Sigma} \hat{\sigma}^z_l$$
2. 设计 6 点梯形积分网格 $\lambda \in \{0.0, 0.2, 0.4, 0.6, 0.8, 1.0\}$。
3. 在每个网格点运行 QMC，测量膜自旋期望值 $\langle M^z_\Sigma \rangle_{+, \lambda} = \langle \sum_{l \in \Sigma} \hat{\sigma}^z_l \rangle$。
4. 执行梯形积分： $$\Delta F = 2 h^z \int_{0}^{1} d\lambda \langle M^z_\Sigma \rangle_{+, \lambda}$$
5. 同调磁化强度最终输出为： $$m_a = \tanh \left( \frac{\beta \Delta F}{2} \right)$$

将上述两个步骤结果结合，即可输出精准复现的 $f_W = \langle W_1 W_2 \rangle - m_a^2$ 标度曲线。

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

4.1 关键引用文献

本文建立在以下数个关键性理论基石之上，这些文献对于深入理解该工作至关重要：

1. 有限温度 TEE 的奠基性研究：
   • C. Castelnovo and C. Chamon, “Entanglement and topological entropy of the toric code at finite temperature”, Phys. Rev. B 77, 184442 (2007).
     首次在数学上推导了热吉布斯态下的拓扑纠缠熵，并发现二维拓扑熵在任何有限温度下都会坍塌。
2. 经典 Z2 规范-希格斯模型相图：
   • E. Fradkin and S. H. Shenker, “Phase diagrams of lattice gauge theories with Higgs fields”, Phys. Rev. D 19, 3682 (1979).
     奠定了希格斯相与局域限制相在没有对称性破缺时属于同一个物理相的宏观图像（即 Fradkin-Shenker 连续性定理）。
3. 混合态拓扑相的通道分类新理论：
   • T.-H. Yang, B. Shi, and J. Y. Lee, “Topological connectivity of mixed states”, arXiv:2506.04221 (2025).
     与本论文并行发展的顶尖工作，独立阐明了为什么经典的 Kitaev-Preskill TEE 在混合态拓扑分类中会失效，并提出了通道纠缠分类学。
4. 三维拓扑码经典相变极限：
   • D. A. Reiss and K. P. Schmidt, “Quantum phase transitions in the three-dimensional toric code in a transverse field”, SciPost Phys. 6, 078 (2019).
     提供了高精度基态相变参数基准。

4.2 局限性与独家学术评论

尽管本工作在定义有限温度拓扑不变量以及证实“几何保护相变”方面取得了突破性进展，但作为一个面向量子信息和凝聚态理论最前沿的研究，以下几个关键局限性与争议点依然值得深入剖析：

1. 解码器依赖性与“假阴性”问题

$f_W$ 的不变量属性在很大程度上依赖于纠错通道 $\mathcal{C}$ 中解码算法的选择。论文在蒙特卡洛计算中采用了全局最优的 最小权重完美匹配（MWPM） 算法。然而，正如论文在第四节中所坦承，MWPM 解码器的固有纠错阈值温度 $T_{\text{MWPM}}$ 一般而言是严格小于系统真正的去禁闭物理相变温度 $T_c$ 的（即 $T_{\text{MWPM}} \le T_c$）。

这产生了一个极其尴尬的**“假阴性”区间 $T_{\text{MWPM}} < T < T_c$：在此温区内，系统在物理上依然处于真拓扑相（磁通线圈尚未发生全局渗透，具有非零线张力），但由于 MWPM 解码器能力不足而过早发生了算法失效（Algorithmic failure），导致 $f_W$ 提前降为 0。这意味着，虽然 $f_W \to 1$ 确实是拓扑相的充分条件**，但它并不能作为充要条件。如何设计一个在有限温度下纠错能力逼近物理相边界极限的最优通道（比如基于最大似然对数张量网络解码器），依然是一个未决的技术瓶颈。

2. “几何保护”概念的普适性与边界划分

作者提出利用毕安基恒等式构建“几何保护”相变，并声称这超越了传统的朗道广义规范对称性。然而，从现代高阶 generalized 规范场论的视角来看，毕安基恒等式本身就是非阿贝尔/阿贝尔一形式规范对称性的本征约束。任何物理晶格上的规范理论，其物理态流形必然只能限制在满足毕安基恒等式的规范子空间中。

因此，将“几何保护机制”与“涌现对称性保护机制”割裂开来，在物理学概念上略显牵强。它更像是从一形式自发对称性破缺语言到晶格动力学约束语言的等价翻译，而非一个独立于高阶群论之外的新范式。

3. 局限在阿贝尔（Abelian）格点模型，向非阿贝尔与非厄米系统的推广困难

本工作的所有推导和蒙特卡洛实现都深度依赖于 $\mathbb{Z}_2$ 的二值阿贝尔代数结构。在 $\mathbb{Z}_2$ 中，由于两个逆操作的叠加依然回归到平凡状态，纠错边界判定被简化为简单的同调类划分。然而：

• 对于更具有应用价值的、可实现非平凡拓扑量子计算的非阿贝尔（Non-Abelian）拓扑模型（如三维双子模型 Double Semion），由于其激发态具有非平凡的熔合规则（Fusion Rules），毕安基恒等式的形式变得极其复杂且不再具有加性对消性。
• 解码算法也会从多项式复杂度的最小匹配算法退化为极其高昂的计算难题。这限制了 $f_W$ 向通用非阿贝尔有限温度拓扑体系推广的可行性。

5. 其他必要补充：物理深度探讨与量子化学交叉启示

5.1 有限温度下的经典记忆 versus 量子记忆（长程纠缠的消亡）

对于三维拓扑码，一个常被混淆的物理事实是：它在有限温度下是否依然是一个真正的量子存储器？

答案是否定的。3D 拓扑码在 $T>0$ 时只能作为经典存储器，而无法存储量子相干。

这一结论可以通过论文在第五节中提到的拓扑纠缠负度（Topological Entanglement Negativity, $\gamma_N$） 得到最深刻的阐明：

• 长程量子纠缠的判定：对于一个纯态拓扑相，拓扑纠缠熵 $\gamma$ 精确刻画了系统长程纠缠（Long-Range Entanglement, LRE）的规模。但在有限温度吉布斯态（混合态）中，由于系统与环境发生了能量交换，总纠缠熵中混入了大量的经典热力学熵（Gibbs 经典涨落）。Levin-Wen 的四区域差值虽然消去了局域面积律和经典热熵，但无法单独区分剩下的是量子相干（纠缠）还是经典的宏观统计自旋锁定。这一判定必须求助于不随经典热涨落增加的物理参量——纠缠负度 $\gamma_N$。
• 3D 拓扑码的热相边界分析：在三维拓扑码中，由于电荷激发（星算符 $\hat{A}_v$ 激发）是没有线张力约束的零维点状缺陷，任何有限温度 $T>0$ 都会引发电荷的无成本激增与任意扩散。这导致在任何 $T > 0$ 下，系统的 $\gamma_N \equiv 0$，即所有的长程量子纠缠在哪怕 infinitesimal 极小的热涨落下一瞬间便蒸发殆尽，仅剩经典的统计长程规范关联。因此，3D 拓扑码在有限温度下仅能用来存储一个经典的比特（即磁通线圈不扩散对应的经典二进制状态），这对应于经典的自校错自旋记忆存储器（如三维经典 Ising 界面）。
• 真正的三维量子记忆—— fermionic 拓扑码： 为了在有限温度下保有真量子相干（$\gamma_N > 0$），必须引入费米子拓扑码（Fermionic Toric Code，见论文 Table I 最后一列）。在 fermionic 码中，激发态表现为带有费米统计的超对称结构。即使引入了外场扰动，由于费米弦（Fermionic string）代数结构的特殊反常，任何局域拟局部通道都无法在不付出巨大自由能代价的情况下将这种长程费米子关联破坏。这在 3D 空间中实现稳定的真量子有限温度拓扑序开辟了崭新的道路。

5.2 规范场论视角：与三维经典 $\mathbb{Z}_2$ 规范-希格斯模型的对偶性

本研究所展现的三维拓扑码在磁场和电场下的有限温度拓扑相变，其背后的物理根源在于一门经典的格点场论—— 3D $\mathbb{Z}_2$ 规范-希格斯模型（格点规范场）。

通过对量子哈密顿量进行 铃木-特罗特分解（Suzuki-Trotter Decomposition）：

在 Trotter 极限下，其配分函数路径积分可以映射到一个四维经典的 $\mathbb{Z}_2$ 规范-希格斯模型。当我们在有限温度下进行一维的尺寸缩减（Dimensional Reduction，由于虚时方向的紧致化，虚时轴的周长 $\beta$ 成为一个固定的短轴），该物理系统在长波红外极限下会退化为一个三维经典规范-希格斯物理体系。其中：

• 量子横磁场 $h^x$ 的作用对应于经典规范场中的规范相互作用耦合常数（$\beta_g$）；
• 量子纵电场 $h^z$ 的作用则对应于经典规范-希格斯耦合（$\kappa$）。

在格点规范场论的语言中，$f_W$ 在去禁闭相中的量子化行为，本质上是格点规范理论处于非局域限制相（Deconfined / Free Charge phase）的数学表征。而 $f_W \to 0$ 则是规范电荷发生凝聚导致系统退化为限制相或希格斯相的标志。根据著名的 Fradkin-Shenker 定理，限制相与希格斯相可以通过绕过临界点（即绕过去禁闭相这个“舌状”突起，见 Fig. 9）进行平滑的、没有任何相变边界的交叉（Crossover）过渡。$f_W$ 的不变量边界刚好完美合围并定义了这个去禁闭“舌状相区”，这极大地丰富了格点规范场相分类的数学工具。

5.3 对量子化学与强关联电子理论的交叉启示

对于量子化学与凝聚态物理交叉领域的科研人员而言，本研究蕴含着两个深刻的启示：

1. 混合态电子结构计算的参考不变量选择：
   在强关联体系（如包含过渡金属的多核活性中心、双自由基等）的量子化学计算中，我们经常使用多组态自洽场方法（MCSCF）或密度矩阵重整化群（DMRG）来处理强静态关联。在处理高压、高温或超快激光泵浦下的非平衡态、吉布斯热态电子结构时，传统的分子轨道或自由能面很难给出干净的电子结构相边界。本论文提出的“解码威尔逊圈关联 $f_W$”的计算框架（即：大尺度随机采样 + 经典自适应纠错解码）实际上提供了一种全新思路。例如，在处理非厄米开系统或带有强环境耗散的电子热态时，我们可以构造类似的、带有“活性空间错误自恢复”性质的纠缠不变量，用以比传统的经典关联函数更干净地勾勒强关联电子相图。
2. 电子分馏（Fractionalization）与高温超导机制（RVB 态）：
   热态拓扑序在 3D $\mathbb{Z}_2$ 拓扑码中的稳定性直接关联着高温超导中的共振价键理论（Resonating Valence Bond, RVB）。在 RVB 理论中，高温超导的正常态可以被视为一个电荷去禁闭、并自发发生自旋-电荷分离（Spin-charge separation）的拓扑自旋液体相。本研究以极高的数值精度证明，这种自旋液体流形在有限温度和较强物理扰动磁场下，依然能够通过毕安基恒等式（对应于电子局域双占据约束——古茨维勒投影 Gutzwiller Projection 的几何必然性）保持其热力学边界不终止。这为 RVB 热态拓扑超导稳定性的微观唯象理论提供了坚实的统计物理学支撑。