来源论文: https://arxiv.org/abs/2503.20566 生成时间: May 28, 2026 12:36

二维量子规范场论的无负符号计算革命：规范规范形式（GCF）与PEPS-VMC方法的深度重构

0. 执行摘要

格点规范场论（Lattice Gauge Theory, LGT）是现代物理学的基石之一。从强相互作用的夸克禁闭机制到凝聚态物理中拓扑序与量子自旋液体的低能有效理论，规范场论提供了不可或缺的数学框架。然而，采用传统的拉格朗日表述下的蒙特卡洛（QMC）方法在面对非零化学势（有限物质密度）或实时动力学（Real-time Dynamics）问题时，会遭遇臭名昭著的“符号问题”（Sign Problem），导致计算复杂度呈指数级增长。

为了绕开符号问题，基于哈密顿量表述的张量网络状态（Tensor Network States, TNS）方法在近年来成为了极具前景的非微扰模拟工具。在(1+1)D体系中，矩阵乘积态（MPS）已取得了巨大成功。然而，当步入(2+1)D及更高维度时，投影纠缠对态（Projected Entangled-Pair State, PEPS）的精确计算却面临双重严峻挑战：一是高维张量收缩的固有高复杂度（#P-hard），二是规范不变量对物理希尔伯特空间强加的局部Gauss定律约束，使得传统的变分优化算法极易陷入局部极小，且计算效率低下。

本研究提出了一种极具开创性的新型张量网络方法。其核心突破在于首次定义了规范不变量张量网络的“规范规范形式”（Gauge Canonical Form, GCF）。通过对格点规范场（Link上的张量B）进行局部规范变换，成功将所有变分参数集中到物质场（Vertex上的张量A）中，使B张量完全退化为无变分参数的静态选择函数。进一步地，研究者将GCF-PEPS与变分蒙特卡洛（Variational Monte Carlo, VMC）完美结合。在VMC采样过程中，由于GCF的引入，每一个物理态构型会自主且唯一地选择物质张量A的单个电荷扇区，使网络收缩的有效键维度从 $D$ 骤降至单个扇区维度 $D_k$，极大地释放了计算算力。

凭借这一创新的方法论，研究团队在多项极具挑战性的(2+1)D阿贝尔格点规范场论体系中实现了超高精度的数值模拟：

1. 纯 $\mathbb{Z}_2$, $\mathbb{Z}_3$ 和 $\mathbb{Z}_4$ 规范场论：在高达 $32 \times 32$ 的大尺寸格点上，基态能量精度达到了与QMC一致的 $10^{-5}$ 级别，并精确确立了相变点。
2. 奇数-$\mathbb{Z}_2$ 规范场论（背景荷 $Q_{\mathbf{x}} = 1$）：首次直接变分得到了非受限相与自发平移对称性破缺的价键固体（VBS）相之间的连续相变行为，临界点 $g_c = 0.64(1)$ 与双重QMC结果完美吻合。
3. 具有动力学物质场的 $\mathbb{Z}_2$ 规范-玻色子耦合模型：首次在二维格点上直接处理了规范场与硬核玻色子（Hard-core Boson）的强耦合，大尺寸有限尺寸外推给出了极其可靠的热力学极限基态能量。
4. 长时间实时动力学模拟：在 $10 \times 10$ 格点（希尔伯特空间维度高达 $2^{81}$）上成功模拟了去禁闭 $\mathbb{Z}_2$ 规范场中单激发（Vison）的无耗散空间传播，且在整个长时间演化轨迹中保持系统能量漂移低于 $0.2\%$。

这项工作不仅彻底奠定了规范不变量PEPS作为模拟二维规范-物质耦合体系的里程碑式工具的地位，其融合VMC与时变变分（t-VMC）的技术路线更将深远地影响量子化学强关联计算、拓扑量子计算设计及非平衡态多体动力学模拟。以下我们将对该工作的科学原理、算法细节、基准评测和实现方案进行全景式的深度解析。

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

1.1 核心科学问题与物理哈密顿量

格点规范场论的哈密顿量构建基于 Kogut-Susskind 理论。在 $d$ 维空间晶格上，物质场（如费米子、玻色子）寄生在晶格顶点（Vertex, $\mathbf{x}$）上，规范场（如电场、电磁矢势）则驻留在连接顶点的边（Link, $(\mathbf{x}, \alpha)$）上，其中 $\alpha = 1, \dots, d$ 表示空间晶格的主轴方向。

对于阿贝尔规范群 $\mathbb{Z}_N$（或连续群 $U(1)$），link 上的物理希尔伯特空间由电场算符 $E$ 的本征态 $|n\rangle$ 张成，满足 $E |n\rangle = n |n\rangle$，其中 $n \in \mathbb{Z}_N$。相应的升算符 $U \equiv e^{i\phi}$（满足 $[\phi, E] = i$）负责将电场状态增加一个单位：$U |n\rangle = |n + 1 \pmod N\rangle$。顶点 $\mathbf{x}$ 处的物质场则由产生消灭算符 $c^\dagger_{\mathbf{x}}, c_{\mathbf{x}}$ 描述。

规范不变性（Gauge Invariance）是系统最核心的物理约束，它要求任何物理态 $|\Psi\rangle$ 必须满足局部 Gauss 定律约束。在每个顶点 $\mathbf{x}$ 处，Gauss 定律算符写作：

其中 $\mathbf{e}_\alpha$ 为 $\alpha$ 方向的单位矢量，$Q_{\mathbf{x}}$ 是预先设定的整数，用于决定顶点 $\mathbf{x}$ 上的背景规范电荷扇区（例如奇数-$\mathbb{Z}_2$ 理论对应 $Q_{\mathbf{x}} = 1$，普通理论对应 $Q_{\mathbf{x}} = 0$）。只有满足 $\hat{G}_{\mathbf{x}} |\Psi\rangle = Q_{\mathbf{x}} |\Psi\rangle$ 的态才属于物理希尔伯特空间。

系统总哈密顿量包含物质项 $H_M$、磁场（Plaquette）项 $H_B$ 和电场项 $H_E$：

此处 $m_{\mathbf{x}}$ 代表裸粒子质量（或化学势），$J$ 为规范-物质耦合强度，$h$ 和 $g$ 分别为磁场和电场耦合常数。当系统不包含物质场时，$J=0$，即退化为纯规范场论。

1.2 规范不变张量网络（GI-TNS）的传统瓶颈

为了自然地表示物理希尔伯特空间，张量网络必须显式地满足规范不变性。通常采用的GI-PEPS ansatz 包含两种张量（图1）：

• 物质张量 $A$（位于 Vertex）：拥有 $2d+1$ 个物理/虚拟腿（Leg）。
• 规范张量 $B$（位于 Link）：拥有 3 个腿（2个虚拟腿，1个物理腿）。

规范不变性通过在张量 $A$ 和 $B$ 的虚拟和物理腿上分配特定电荷 $q(j)$，并施加稀疏性对称约束来实现：

这里 $l, r, d, u$ 分别表示左、右、下、上虚拟键的电荷分配。总物理键维度（Bond Dimension）写作各电荷扇区简并度之和 $D = \sum_k D_k$。

传统计算难点： 虽然该方案在数学上无懈可击，但在数值变分优化时会遇到灾难性困难。在传统的收缩（如边界MPS方法）中，由于 $A$ 和 $B$ 张量在所有 link 上相互耦合，收缩复杂度极高。更为致命的是，在进行基于张量重整化群或简单更新（Simple Update）的虚时演化时，由于局部Gauss定律约束的刚性限制，优化极易锁死在高度局域的非物理极小值中，导致算法无法收敛。这使得在 (2+1)D 下利用 PEPS 进行精确的大尺寸计算一直难以实现。

1.3 理论突破一：规范规范形式（GCF）的严格推导

为了消除虚拟键上的冗余规范自由度，作者在连接张量 $A$ 与 $B$ 的 Link 上定义了一个分块对角矩阵 $X$：

其中 $B^{[k]}$ 是一个 $D_k \times D_k$ 的矩阵，它是由在规范张量 $B^n_{lr}$ 中选择物理量子数 $n=k$，并将虚拟指标 $l, r$ 限制在电荷为 $k$ 的扇区内而提取出来的。

由于规范变换的自由度，我们可以对张量进行如下的伴随局部变换（Gauge Transformation）：

在这个局部的基底变换下，规范张量 $B$ 得到了极大简化：

这意味着什么？在规范规范形式（GCF）下，所有的规范张量 $B$ 都退化为了完全参数无关（Parameter-free）的 $1/0$ 静态选择器！ 所有系统的物理和变分自由度都被干净地吸纳到了顶点处的物质张量 $A$ 中。现在，变分算法只需集中力量优化物质张量 $A$，这从根本上减少了变分参数的数量，消除了规范冗余，保证了数值梯度的平滑性。这一发现不仅简化了二维PEPS，也直接简化了一维的GI-MPS算法。

   传统的局部张量连接：                   规范规范形式（GCF）简化后：
      |                                      |
   --[A]-- (l) -- [B] -- (r) -- [A]--     --[A]-----------------[A]--
      |            |             |           |         |           |
                 |n>                                 |n>

1.4 理论突破二：GCF 与变分蒙特卡洛（VMC）的超强化学式化学反应

将 GCF 与 VMC 结合，是本工作最精妙、最具开创性的物理设计。在变分蒙特卡洛中，我们计算任意算符 $O$ 的期望值：

其中 $|\mathbf{s}\rangle \equiv |\mathbf{n}, \mathbf{p}\rangle$ 代表格点上全体规范场和物质场的一个特定物理解析构型。由于我们只对严格满足规范不变性（Gauss定律）的物理构型进行采样，因此任何非物理构型的振幅 $\langle \mathbf{s} | \Psi \rangle$ 都精确为零。

此时，GCF 带来的惊人效果显现了：对于任意给定的物理构型 $|\mathbf{s}\rangle = |\mathbf{n}, \mathbf{p}\rangle$，规范场状态 $\mathbf{n}$ 的具体数值已经在每一条 link 上被完全固定。由于 $B$ 张量满足选择函数形式 $B^n_{lr} = \delta_{lr} \delta_{n, q(l)} \delta_{n, q(r)}$，电场值 $n_k$ 会唯一地选择与之相连的物质张量 $A$ 的特定电荷扇区！

具体而言，在网络 $\langle \mathbf{s} | \Psi \rangle$ 的收缩中：

1. 所有规范张量 $B$ 在构型投影后完全消失。
2. 每一个顶点的物质张量 $A$ 被投影为一个单一扇区张量 $\tilde{A}$，其对应的有效虚拟键维度从完整的 $D$ 骤降为单个扇区维度 $D_k$（例如，在 $\mathbb{Z}_2$ 中，$D=2D_k$，相当于键维度直接减半）。

在 PEPS 收缩中，计算复杂度随着键维度呈现极高阶的多项式扩展。将收缩键维度从 $D$ 降低到 $D_k$，使得我们可以轻松处理原本在二维晶格上完全不可企及的高键维度（如 $D=14$），并允许采用更为精细的边界MPS方法（Boundary MPS）来精确收缩单层张量网络。其计算复杂度仅标度为：

其中 $\chi$ 为边界MPS截断维度。本研究中，取 $\chi = 3D_k$ 即可提供极其优异的收敛精度。

1.5 随机重构（SR）与时变变分蒙特卡洛（t-VMC）

对于基态优化和实时动力学演化，本工作采用基于信息几何学的随机重构（Stochastic Reconfiguration, SR）算法。对于动力学演化，即求解含时薛定谔方程，根据时变变分原理（TDVP），参数 $\alpha_m$ 的演化方程满足：

其中 $S_{mn}$ 是由蒙特卡洛采样计算的变分参数重叠矩阵（Fisher信息矩阵）：

梯度算符定义为：

外驱力矢量 $g_m$ 为物理梯度：

为了确保长时间实时演化的数值稳定性与能量守恒，作者摒弃了一阶Euler格式，而创造性地在PEPS中引入了基于Taylor展开的含时变分两步更新方案（Eq. S8）：

其中 $c_{1,2} = (1 \pm i)/2$，该格式具有稳定的阶数 $\mathcal{O}(\Delta t^3)$。配合高精度的共轭梯度（CG）求解器，成功突破了二维PEPS无法进行长时稳定演化的技术壁垒。

2. 关键 Benchmark 体系、计算数据与物理发现

2.1 纯 $\mathbb{Z}_N$ 规范场论（$N=2,3,4$）

研究团队首先在无物质耦合（$J=0$, $Q_{\mathbf{x}} = 0$）的纯规范场论中验证了算法，基键维度取 $D_k=2$（对应总键维度 $D=4, 6, 8$）。

2.1.1 纯 $\mathbb{Z}_2$ 规范场论

在 $16 \times 16$ 的格点系统上，研究团队将 GCF-PEPS 计算的基态能量与经典的无符号问题量子蒙特卡洛（QMC）结果进行了详尽的逐点对比（Table I）。结果表明，在临界点 $g_c \approx 0.3285$ 附近，PEPS的能量精度依然能够保持在 $10^{-5}$ 到 $10^{-6}$ 的极高水准：

这一极高精度的达成，得益于变分蒙特卡洛中的“零能量方差原理”（Zero Energy Variance Principle），当变分态逼近真实本征态时，局部能量方差趋近于零，从而在极小样本量下也能获得超低统计噪声的优异解。

2.1.2 纯 $\mathbb{Z}_3$ 规范场论

对于 $\mathbb{Z}_3$ 场论，研究团队计算了尺寸从 $8\times 8$ 直至 $24\times 24$ 的基态性质。通过对能量的一阶导数 $\partial\langle H\rangle/\partial g$（图2(a)）及二阶差分导数 $\partial^2\langle H\rangle/\partial g^2$（图2(b)）进行尺寸标度分析，清晰地捕获到了相变点的信号。通过有限尺寸外推，精确确定了相变点位于：

该相变展现出明显的一阶相变特征。作为对比，此前无规范约束的无限体系iPEPS（通过简单更新优化）预测的相变点为 $g_c = 0.448(3)$，且能量显著高于本工作的外推值。这有力地证明了，不带显式规范不变量约束且采用局部更新的张量网络，其变分精度在强关联相变区会大打折扣，而本方法的全局VMC变分方案则能彻底逼近真实物理基态。

2.1.3 纯 $\mathbb{Z}_4$ 规范场论

这是张量网络方法在二维晶格上对 $\mathbb{Z}_4$ 规范场论的首次探索（图2(c, d)）。相变点被精确定位在：

有趣的是，这个相变点与 $\mathbb{Z}_2$ 规范场论的相变点（$0.3285$）在误差范围内高度重合。根据三维经典规范场论的对偶理论，量子 $\mathbb{Z}_4$ 规范场论在数学上对偶于 $\mathbb{Z}_2 \times \mathbb{Z}_2$ 场论。本工作通过最直接的非微扰张量计算，强烈暗示并支持了这一等价性猜想在量子晶格规范场论中依然成立。

2.2 奇数-$\mathbb{Z}_2$ 规范场论（Odd-$\mathbb{Z}_2$ Theory）

该系统在每个格点上引入了静态背景电荷 $Q_{\mathbf{x}} = 1$，是研究高温超导、量子自旋液体和二聚体模型（Dimer Models）的核心模型。理论预测，随着耦合常数 $g$ 的增加，该系统会经历一个从去禁闭相（Deconfined Phase）到自发平移对称性破缺的局域禁闭相——**价键固体（VBS, Valence-Bond Solid）**的连续相变。

本研究利用 GI-PEPS 在高达 $32 \times 32$ 的巨大格点上进行了非微扰计算。空间局域不变量——Plaquette 期望值 $\langle P_{\mathbf{x}} \rangle$ 在 $g/h=0.4$ 时表现出完全均匀的空间分布（对应去禁闭相，图3(a)），而在 $g/h=0.8$ 时则呈现出完美的条纹状自发平移对称性破缺图案（对应VBS相，图3(b)）。

为了确定临界点，研究者精确计算了 VBS 序参量 $D_{x/y}$：

通过二次多项式有限尺寸外推（图3(c)），得到热力学极限（$L \rightarrow \infty$）下的转变点为：

这与双重空间QMC计算得到的对偶模型（完全挫败的横场伊辛模型）临界值 $g_c \simeq 0.634$ 达成了高度一致。

此外，研究团队还计算了 $32 \times 32$ 格点上的非局域算符——**威尔逊圈（Wilson Loop, $\langle W \rangle$）**的周长-面积标度行为。通过提取弦张力 $\sigma$（图3(d)），发现在 $g \lesssim 0.6$ 时 $\sigma$ 几乎为零（代表周长定则，Deconfined），而在 $g > 0.6$ 后 $\sigma$ 陡然线性上升（代表面积定则，Confined），从非局域拓扑可观测量的角度完美互补印证了相变的连续发生。

   Wilson Loop 标度律示意：
   ln<W> 
     |       / (面积定则 Area Law, 禁闭相 Confined)
     |      /
     |     /
     |----/--- (周长定则 Perimeter Law, 去禁闭相 Deconfined)
     |__________________ Area S

2.3 $\mathbb{Z}_2$ 规范场论耦合硬核玻色子

这是数值计算难度极高、此前在 (2+1)D 下处于完全空白状态的规范-物质耦合模型。首先，研究者在小尺寸 $3 \times 3$ 系统（包含2个硬核玻色子，耦合参数 $h=1, g=0.33, J=0.5$）上与精确对角化（ED）进行了严苛的比对。在 $D=6$ 的 PEPS 计算下，能量达到了极其惊人的收敛：

• ED 精确能量：$-0.4707135061$
• GI-PEPS 能量：$-0.470713502(3)$ （在 $10^5$ 采样下精度即达 $10^{-8}$ 级别！）

随后，研究团队将系统拓展至 $16 \times 16$ 半满玻色子的大尺寸系统。能量随键维度 $D$ 呈现出单调向下收敛（图4(a)），由于引入了活跃的费米/玻色子物质场，纠缠熵显著增加，收敛所需的键维度从纯规范场的 $D=4$ 提升至 $D=12$。利用不同核心 bulk 尺寸 $L_b \times L_b$ 的外推（图4(b)），在热力学极限下给出了极其可信的基态能量 $-1.3322(4)$。

物理发现：在引入动态物质场后，即使在原本纯规范场属于“禁闭”的强电场区域（如 $g = 0.4$），由于硬核玻色子的存在，电荷屏蔽机制（Charge Screening）开始起主导作用。这导致原本表现为“面积定则”的威尔逊圈，在任意非零物质耦合下都完全退化回了“周长定则”标度（图4(d)），从数值上直观呈现了规范禁闭与物质屏蔽之间的物理纠缠。

2.4 实时动力学：Vison激发的空间传播

在去禁闭的 $\mathbb{Z}_2$ 基态（$g=0.1$）上，通过在左下角连接处施加局部 Pauli 算符 $\sigma_z$，构造一个受激的准粒子——Vison。实时含时演化采用 $\Delta t = 0.005$ 步长。

在 $6 \times 6$ 的小晶格上，GI-PEPS 与 ED 计算的各时变 Plaquette 演化曲线 $\langle P_{\mathbf{x}}(t) \rangle$ 完全重合，展现了变分动力学的超高保真度（图5(a)）。

而在 ED 彻底失效、希尔伯特空间高达 $2^{81}$ 的 $10 \times 10$ 庞大系统上，GI-PEPS 成功观察到了 Vison 激发波包极其清晰地朝着右上方向进行空间无耗散传播的全过程（图5(b-d)）。在长达 $T=18$（对应高达 7200 步变分时空演化）的超长演化轨迹中，系统的总能量漂移依然被牢牢锁死在 $0.2\%$ 之内，展现了时变变分蒙特卡洛（t-VMC）算法配合Taylor双步更新方案在解决高维张量演化中“能量漂移/雪崩失效”难题上的绝对成功。

3. 代码实现细节、复现指南与开源生态

3.1 核心算法实现流程与伪代码

要在数值上复现这一方案，其核心在于正确构建对称分块物质张量 $A$，利用 GCF 对物理构型进行采样，并采用边界MPS技术收缩无 $B$ 参与的单层对称网络。以下提供核心演化与收缩逻辑的伪代码实现框架：

import numpy as np
# 假定使用基于张量块对称性的科学计算库（如 TenPy 或自定义的分块稀疏张量库）

class GCF_PEPS_VMC:
    def __init__(self, L_x, L_y, gauge_group="Z2", D_k=2, chi=6):
        self.L_x, self.L_y = L_x, L_y
        self.gauge_group = gauge_group
        self.D_k = D_k # 单电荷扇区的键维度
        self.chi = chi # 边界MPS收缩截断键维度
        
        # 初始化物质张量 A（仅在顶点），其结构严格满足 Gauss 定律约束的多分块形式
        self.A_tensors = self.initialize_block_sparse_A()
        
    def initialize_block_sparse_A(self):
        """
        构建满足 A^p_{lrdu} 对称稀疏性的变分张量块。
        对于 Z2 规范群，变分参数只存在于 p + q(l) + q(d) - q(r) - q(u) = Q_x (mod 2) 的非零块中。
        """
        A = {}
        # 遍历所有顶点，初始化每个扇区的参数
        for x in range(self.L_x):
            for y in range(self.L_y):
                # 分配随机数作为初始变分实/复数参数
                A[(x, y)] = create_symmetric_tensor(phys_dim=2, bond_dims=[self.D_k]*4)
        return A

    def get_amplitude(self, config_s):
        """
        突破点：投影物理构型 s = |n, p>
        此时 Link 上的 B 消失，A 投影为 A~ (有效键维从 D 降为 D_k)
        """
        A_projected = {}
        for vertex in self.A_tensors:
            # 根据 config_s 中该顶点周围的 link 规范场值(n_l, n_r, n_d, n_u)
            # 以及顶点的物理物质态 p，提取对应的唯一分块，使其键维只有 D_k
            n_l, n_r, n_d, n_u, p = config_s.get_local_fields(vertex)
            A_projected[vertex] = self.A_tensors[vertex].get_block(n_l, n_r, n_d, n_u, p)
            
        # 使用传统的边界MPS压缩算法收缩这个仅有 D_k 的纯物质二维张量网络
        amplitude = contract_2D_peps_boundary_mps(A_projected, chi=self.chi)
        return amplitude

    def metropolis_step(self, current_config):
        """
        规范不变性保持的 Metropolis 物理构型更新：
        1. Plaquette flip: 随机选择一个面，在其四周 link 施加规范升降算符 P_x^t，完全不破坏 Gauss 定律。
        2. 局部 Link flip (若存在物质场耦合): 施加 c*_x U_{x,a} c_{x+a} 协同翻转。
        """
        proposed_config = current_config.propose_plaquette_flip()
        
        # 计算振幅之比（VMC 核心项）
        p_ratio = np.abs(self.get_amplitude(proposed_config) / self.get_amplitude(current_config))**2
        
        if np.random.rand() < p_ratio:
            return proposed_config
        return current_config

    def update_parameters_SR(self, learning_rate=0.01):
        """
        随机重构 (Stochastic Reconfiguration) 变分参数更新
        """
        # 1. 蒙特卡洛采样收集 M 个物理构型
        configs = self.run_monte_carlo_sampling(samples=10000)
        
        # 2. 计算局部能量 Eloc(s) 以及对数梯度 O_m(s)
        E_loc, O = self.compute_observables(configs)
        
        # 3. 组装重叠矩阵 S_mn 和梯度矢量 g_m
        S = self.assemble_fisher_matrix(O)
        g = self.assemble_gradient(O, E_loc)
        
        # 4. 使用共轭梯度法（CG）求解线性系统: (S + epsilon * I) * delta_alpha = g
        delta_alpha = conjugate_gradient_solver(S, g, regularizer=1e-3)
        
        # 5. 更新 A 张量的参数
        self.apply_parameter_update(delta_alpha, learning_rate)

3.2 开源依赖与复现资源

1. 分块稀疏对称张量底层： 该研究的核心代码和算法构建在 Python 开源张量网络生态之上。作者明确提及，算法的测试和部分对称性分块逻辑依赖了著名的 TenPy (Tensor Network Python) 代码库。

   • TenPy 官方仓库：https://github.com/tenpy/tenpy
   • 复现建议：可以利用 TenPy 中的 SymmetricTensor 模块来声明满足阿贝尔规范群不变量约束的 block-sparse 张量结构。

2. 变分蒙特卡洛（VMC）收缩框架： 由于单层网络投影后的收缩属于经典的无物理指标PEPS收缩，复现者可参考并引入高级PEPS-VMC收缩算法：

   • PEPS-VMC 算法实现参考：https://github.com/vanderstraeten/peps-vmc （提供了有关边界MPS在VMC中的高效压缩参考流程）。

3. (1+1)D Schwinger 模型复现（End Matter部分）： End Matter 中提及的一维 Schwinger 模型利用基于 GCF 的传统 MPS 时空演化（TEBD），可以直接通过标准的 $U(1)$ 对称性 MPS 代码包实现，复现时无需显式写出规范场的 $B$ 张量，只需在 TEBD 每次局域演化分裂双位张量（SVD）时，手动修正对称荷的分配即可。这是一种极其优雅的一维LGT实现方案。

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

4.1 关键里程碑文献

本工作建立在格点规范场与张量网络数十年的交叉发展历史之上，其最关键的立论支撑文献包括：

1. Kogut & Susskind (1975) [[39]]：奠定了晶格规范场论的哈密顿量构建（KS哈密顿量），是本工作所有物理模型的出发点。
2. Verstraete & Cirac (2004) [[29]]：提出了投影纠缠对态（PEPS）理论，为二维强关联物理提供了张量网络的基础语言。
3. Buyens, Haegeman, et al. (2014) [[14]]：在 (1+1)D 下首次完整阐述了矩阵乘积态（MPS）中的规范不变对称性分配。本工作的 GCF 形式正是对该文一维对称性机制在二维高阶情况下的重大推广与提炼。
4. Robaina, Banuls, Cirac (2021) [[38]]：使用非规范不变量约束下的无限 PEPS（iPEPS）加简单更新（Simple Update）演化来研究纯 $\mathbb{Z}_3$ 规范场论。该文的能量数据在强关联和临界相变区由于变分受限而偏高，本工作通过 GCF-PEPS 全局 VMC 优化对其基准数据进行了决定性的超越和修正。
5. Sandvik & Vidal (2007) [[40]] / Carleo et al. (2012/2017) [[57, 70]]：开创了张量网络变分蒙特卡洛（VMC）与含时变分蒙特卡洛（t-VMC）先河，为高维波函数计算提供了极佳的非定域优化范式。

4.2 局限性与前沿评论

尽管本工作在二维阿贝尔规范场论的模拟上取得了里程碑式的飞跃，但作为一个新兴的技术方向，该方案仍存在一些不容忽视的物理与技术局限，这也是未来该领域的前沿主攻方向：

1. 非阿贝尔规范群（Non-Abelian Gauge Groups）的拓展壁垒： 本工作目前的“规范规范形式（GCF）”高度依赖于阿贝尔群 $\mathbb{Z}_N$ 或 $U(1)$ 的对角化与一维群表示。对于高能物理中更为关键的非阿贝尔规范场论（如含有自相互作用的 $SU(2)$ 或 $SU(3)$ 规范群，即真实的 QCD 模拟），其 link 上的规范张量 $B$ 具有非平庸的非阿贝尔不可约表示耦合（Clebsch-Gordan 结构），其代数结构极度复杂。如何定义非阿贝尔规范群的 GCF，并在物理采样投影时保持有效键维度的有效缩减，是本方法向强相互作用物理进行终极跨越的巨大理论难点。

2. 真正的二维费米子物质耦合挑战： 在本工作的 (2+1)D 物质场耦合模拟中，仅处理了硬核玻色子。而在高能物理或真实的凝聚态关联体系中，物质场通常是费米子（如 Schwinger 费米子，狄拉克费米子）。在二维 PEPS 中直接处理费米子需要引入复杂的 Swap Gate（交换门）或者通过 Grassmann 积分重构。费米子算符带来的“负正负号”如何完美地融合到 VMC 的概率采样机制中，而不引发采样统计噪声的雪崩，仍是一个亟待解决的工程与理论痛点。

3. 临界区边界 MPS 截断维度 $\chi$ 的收敛稳定性： 在相变临界区，系统的关联长度发散，边界 MPS 收缩所需的截断维度 $\chi$ 往往需要极速膨胀以保证收缩精度。虽然 VMC 已经通过投影将单层有效键维降到了 $D_k$，但在强电荷屏蔽或去禁闭的临界点，由于量子纠缠极强，$\chi$ 的增加依然会使计算成本显著上升。如何在高度纠缠区更稳定、自适应地控制 $\chi$ 的系统误差，需要引入更精密的重整化群方法。

5. 跨界碰撞：格点规范场 PEPS-VMC 对量子化学强关联计算的启示

作为面向现代计算化学与多体物理的技术评论，我们需要深刻认识到：格点规范场论的 GCF-PEPS-VMC 方案，与量子化学中解决过渡金属配合物或多活性中心强关联电子结构问题的核心思想，有着异曲同工、甚至可深度互补的精妙逻辑。

5.1 Gauss定律约束与量子化学自旋/空间对称性的对偶性

在量子化学中，精确构建电子波函数最困难之处在于活性空间（Active Space）的维度随着轨道数和电子数呈指数爆炸（即 Full CI 极限）。通常，量子化学家通过引入自旋自适应轨道波函数或利用分子点群空间对称性来强行将哈密顿矩阵限制在特定的自旋对称区（如单重态 $S^2=0$、三重态 $S^2=2$）。

本工作解决格点规范物理的核心，就是处理“Gauss定律约束”。Gauss定律实际上是一种强烈的、局域的、顶点级的动力学电荷守恒约束。本研究所引入的规范规范形式（GCF），本质上是一种完美的“物理轨道/活性空间伴随对角化”技术。量子化学张量网络算法（如 qc-DMRG 或二维 qc-PEPS）在处理分子体系时，由于库仑相互作用的非定域性，经常遇到参数冗余的问题。如果能够借鉴 GCF 的思想，在化学键连接（Link）上引入局部变换，将轨道间的纠缠冗余和化学不变量（如原子轨道局域电荷、自旋多重度配对）通过静态选择张量 $B$ 予以固定，使得所有的优化仅集中在定域分子轨道上，将极大精简张量网络的变分参数库，使超大过渡金属簇波函数的 PEPS 精确变分成为可能。

5.2 变分蒙特卡洛（VMC）与张量网络的融合是解决高维大活性空间计算的终极出路

传统的量子化学二维 PEPS 计算面临的最大瓶颈，在于收缩复杂度。由于真实分子体系通常具有多体非定域库仑力，其 PEPS 的有效键维度 $D$ 往往需要非常大（如 $D \ge 16$），导致直接数值收缩几乎是不可完成的任务。

本工作给出了最完美的示范：用 VMC 绕过高维张量网络的收缩死劫。在量子化学中，我们同样可以对电子在定域分子轨道上的占据数构型（即行列式构型 $|\mathbf{n}\rangle$）进行马尔可夫链采样。由于电子守恒约束，特定的占据数构型会直接对分子 PEPS 张量进行强制性的电荷/自旋扇区裁剪，使有效收缩键维大幅下降。将此处的 GCF-PEPS-VMC 框架移植至多活性中心催化剂、铁硫簇（Fe-S Clusters）的电子能级与非平衡激发态动力学模拟，不仅能在多项式时间内完美保持物理自旋对称性，其时变变分动力学（t-VMC）更有望直接精确捕获多发色团、人工光合作用系统中的电荷、激子超快传输动力学过程。

综上所述，这一项源自高能-凝聚态物理最前沿的张量计算革命，不仅为格点规范场论的无负符号计算推开了全新的大门，也为未来多维强关联化学物理体系的高保真模拟提供了极其强大的跨界技术范式。