来源论文: https://arxiv.org/abs/2606.20029v1 生成时间: Jun 19, 2026 16:07

克服小体积效应的破局方案：(2+1)维哈密顿量U(1)规范理论中格点阶梯标度连续外推的有限体积方案深度解析

0. 执行摘要

在现代粒子物理与凝聚态物理的交叉领域，格点规范理论（Lattice Gauge Theory, LGT）是研究非微扰强相互作用最强有力的工具。然而，传统的拉格朗日表述格点规范理论严重依赖于马尔可夫链蒙特卡洛（MCMC）抽样，这使它在面对有限化学势的符号问题（Sign Problem）以及实时动力学模拟时举步维艰。此外，在趋近连续极限（Continuum Limit）的过程中，MCMC 还会遭遇致命的“临界减速”（Critical Slowing Down）效应，导致自相关时间呈指数级增长。

作为一种极具前景的替代方案，哈密顿量格点规范理论在变分方法（如张量网络，Tensor Networks）或未来量子计算平台的辅助下，能够彻底免疫符号问题和临界减速。然而，现阶段经典张量网络模拟与量子硬件均受到空间体积的严格限制。在极小的格点体积下，如何精确设定物理标度并平滑地进行连续极限外推，成为了阻碍哈密顿量 LGT 走向定量精确科学的首要技术瓶颈。

近期发表的学术论文《A Finite-Volume Scheme for the Continuum Extrapolation of Lattice Step-Scaling in (2+1)D Hamiltonian U(1) Gauge Theory》针对这一核心痛点，首次提出并成功验证了一种固定比例有限体积外推方案。该方案通过在趋近连续极限的过程中，将物理参考距离与格点纵向尺寸的比例保持为常数，完美地将有限体积效应转化为系统、可控且平滑的修正项，避免了传统外推中因体积不匹配导致的非单调无序振荡。作者以 (2+1) 维紧致 U(1) 规范场论为测试基准，结合对偶哈密顿量表述与矩阵乘积态（MPS）变分模拟，不仅清晰地观测到了静电荷间的 confinement（局域禁闭）通量筒结构，更在极小尺寸的格点上实现了极为精确、稳定的连续极限外推。这一成果为未来在量子模拟器和张量网络上运行非阿贝尔规范场论的阶梯标度（Step-Scaling）计算铺平了道路。

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

1.1 核心科学问题：窗口问题与小体积污染

在强相互作用（如 QCD）的数值模拟中，核心目标之一是确定跑动耦合常数（Running Coupling）并与高能微扰机制进行无缝对接。然而，直接在单一格点模拟中跨越非微扰低能区到微扰高能区会遭遇著名的**“窗口问题”（Window Problem）**。具体而言，为了避免离散化误差，格点步长 $a$ 必须远小于物理探测标度的倒数 $1/\mu$；而为了避免有限体积污染，格点线性尺寸 $L$ 又必须远大于 $1/\mu$，即：

在算力受限的条件下，同时满足这两个不等式极其困难。为了解决这一矛盾，Lüscher、Weisz 和 Wolff 引入了阶梯标度法（Step-Scaling Approach）。该方法通过在一系列较小的有限体积中迭代地改变尺度，每次都在小体积下完成连续极限外推（$a \to 0$），随后通过步长因子 $s$ 将耦合常数过渡到更大的物理尺度，最终将低能参考标度与微扰高能区连接起来。

然而，当我们将阶梯标度法引入哈密顿量 LGT 时，会遇到一个极其棘手的新问题。由于当前的张量网络算法（如 MPS）在处理二维空间时计算量随宽度呈指数增长，我们只能在极小的空间格点（例如宽度 $N_y = 3$ 的狭长条带）上进行计算。在如此微小的体积下，传统的连续外推方法（即直接减小格点步长 $a$ 而不改变格点数）会导致物理体积在每一步模拟中发生剧烈变化。这种多变的有限体积效应会彻底污染格点阶梯标度函数 $\Sigma_s$，导致外推曲线呈现出高度非单调、无规律的振荡（见后文的“反面教材”分析），使得根本无法提取出有物理意义的连续极限值。

1.2 理论基础：Kogut-Susskind 哈密顿量及其对偶变换

规范场论在格点上的经典哈密顿量表述为 Kogut-Susskind 形式。对于 $a=1$ 的二维空间格点，其原始哈密顿量由电场项（电能）和磁场项（磁能）两部分组成：

其中，$\hat{L}_{\mathbf{x}, k}$ 是定义在格点链路上的共轭电场动量算符，$\hat{P}_{\mathbf{x}, ij}$ 是由四个链路算符构成的基本面元（Plaquette）算符，代表磁通量。在极弱耦合极限（$g^2 \to 0$）下，磁场项在哈密顿量中占绝对主导地位。由于磁场项是一个四链路的复杂相互作用，它会在系统的基态中产生极其庞大且复杂的量子纠缠。这使得基于纠缠度量控制的张量网络算法（如 DMRG）在趋近连续极限时，所需的键合维度（Bond Dimension）呈指数级暴涨，计算迅速失控。

为了克服这一困难，本研究采用了对偶变换（Duality Transformation）。在施加开边界条件（Open Boundary Conditions）的前提下，通过将格点链路变量重参数化为面元变量，并利用高斯定理（Gauss’ Law）将辅助链路变量完全消去，可以推导出一个完全由面元度规表征的对偶哈密顿量：

在这个对偶表述中，原本复杂的四体相互作用磁场项变成了完全局域的单体面元算符 $\hat{P}_{\mathbf{x}}$，而原本局域的电场项则转变为面元之间的最近邻相互作用。这一奇妙的对偶性使得系统在弱耦合（微扰）区间的纠缠结构被大幅简化，极为适合用张量网络进行高精度变分求解。

1.3 技术难点：静态电荷的引入与高斯定理的修正

为了定义物理耦合常数，我们必须在系统中引入一对静止的外部电荷（正电荷 $+q$ 与负电荷 $-q$），通过测量它们之间的静电势及相互作用力来提取重正化耦合。然而，静态电荷的引入打破了原有高斯定理的齐次形式，高斯定理被修正为：

其中 $Q_{\mathbf{x}} = \pm 1$ 代表在位置 $\mathbf{x}$ 处存在的外部电荷。在对偶变换过程中，这意味着我们无法将所有的辅助链路变量完全消除。为了在对偶哈密顿量中保留这些电荷的影响，必须人为引入一条连接该对电荷的任意非动力学辅助路径 $\gamma$（即狄拉克弦，Dirac String），用以承载它们之间的电通量。最终得到的含静电荷对偶哈密顿量写为：

其中 $n_{\mathbf{x}, i}^{(\gamma)}$ 是一个整数背景场，当且仅当对偶链路被路径 $\gamma$ 穿过时取值为 $\pm 1$，其余地方为 $0$。这一背景场的引入能够精确模拟夸克禁闭中的通量筒（Flux Tube）效应。

1.4 方法细节：固定比例有限体积方案（Fixed-Ratio Scheme）

为了消除连续外推中多变的有限体积效应，本文提出了该研究最核心的方法论创新：固定比例外推方案。该方案的精髓在于，在逐步减小格点步长 $a$（即提高格点分辨率）的过程中，我们不保持格点数不变，而是保持物理参考距离 $r_1$ 与格点空间线性尺寸 $N_x$ 的比例为常数。具体步骤如下：

1. 定义力方案下的重正化耦合常数：首先通过数值模拟计算在不同电荷间距 $r$ 下的基态能量，减去真空背景能量后得到静电势 $V(r)$。利用如下经验公式（完美契合 (2+1) 维 U(1) 理论在长程的线性禁闭和短程的对数库仑行为）进行拟合：

   $$aV(r/a) = \alpha + \sigma (r/a) + \delta \log(r/a)$$

   对其进行解析求导，提取出静电力 $F(r) = dV/dr$，进而定义出跑动耦合常数：

   $$g^2(r) = r^2 F(r)$$

2. 构建等比例格点序列：选定初始格点体积 $V^{(1)}$，其对应的参考距离为 $r_1^{(1)}$，步进距离为 $r_2^{(1)} = s r_1^{(1)}$。对于后续更精细的格点序列 $V^{(i)}$，我们严格按比例定义其格点单位下的参考距离：

   $$r_1^{(i)} = r_1^{(1)} \frac{V^{(i)}}{V^{(1)}}$$

   从而在整个外推序列中，始终保持比例恒定：

   $$\frac{r_1^{(i)}}{V^{(i)}} = \text{const.}$$

3. ** bare coupling（裸耦合）微调**：在每个格点 $i$ 上，调节其原始哈密顿量中的裸耦合常数 $g^2_{(i)}$，使得在该格点上计算得到的重正化耦合常数在参考距离 $r_1^{(i)}$ 处精确等于一个设定的常数 $u$：

   $$r_1^{(i) 2} F(r_1^{(i)}, g^2_{(i)}) = u$$

4. 外推步进值：在微调好裸耦合后，在相同的格点上测量步进距离 $r_2^{(i)}$ 处的重正化耦合，这就是我们的格点步进函数值：

   $$\Sigma_s^{(i)} = r_2^{(i) 2} F(r_2^{(i)}, g^2_{(i)})$$

   由于物理比例在序列中完全被冻结，$\Sigma_s^{(i)}$ 随格点步长 $a$ 的变化将呈现出极度平滑的二次曲线行为，可以通过标准的齐次多项式 $f(a) = c_1 + c_2 a^2$ 极其干净地外推至连续极限 $a \to 0$。

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

2.1 测试基准体系参数

研究团队选择 (2+1) 维紧致 U(1) 纯规范理论作为数值测试基准。在空间格点设计上，为了适应矩阵乘积态（MPS）的一维拓扑结构，采用了极具代表性的狭长条带几何（Strip Geometry）。

• 格点尺寸：纵向尺寸 $N_x$ 在微调与测试中分别取 $N_x \in \{6, 7, 8, 9\}$ 以及延伸至 $\{10, 11, 12\}$。横向宽度固定为极小的 $N_y = 3$，以此模拟极度受限的小体积环境。
• 物理步进因子：设定为 $s = \sqrt{2}$。
• 基态能求解：使用基于两位置密度矩阵重正化群（2-site DMRG）的 MPS 变分算法。
• 希尔伯特空间截断：由于对偶表述中的面元变量具有无限维的整数谱，在数值模拟中必须对抗局部维度进行截断。定义面元截断参数为 $\Lambda$（即局部态空间维度 $d_O = 2\Lambda + 1$），论文的主要物理结果均在 $\Lambda = 1$（3个局部态）下获得，并系统评估了 $\Lambda = 2$ 到 $\Lambda = 7$ 的系统误差。

2.2 关键物理数据提取

2.2.1 局域禁闭通量筒的观测

在 $8 \times 5$ 的格点、裸耦合 $g^2 = 0.5$ 且电荷间距 $r = 5$ 的系统下，研究团队通过计算相邻面元之间的局域电能分布：

成功绘制了电磁能量密度空间分布图（参见原文 Fig. 4）。结果清晰地显示出：电通量高度局域化地集中在连接正负电荷的单一直线上，并向四周快速衰减。这不仅直观地证实了二维空间下 U(1) 规范理论的弦禁闭（String Confinement）特征，也验证了对偶路径背景场 $n_{\mathbf{x}, i}^{(\gamma)}$ 物理描述的正确性。

2.2.2 静电势拟合与力提取

在同样的 $8 \times 5$ 体系下，对电荷间距 $r/a \in \{2, 3, 4, 5\}$ 的势能数据进行三参数拟合：

提取出的拟合参数为：恒定偏移量 $\alpha = 0.0160$，弦张力 $\sigma = 0.0211$，以及对数库仑系数 $\delta = 0.00594$。计算表明，在所测量的距离范围内，对数项贡献占总势能的 9%~10% 左右，这与 (2+1) 维紧致 QED 的理论预期高度吻合。

2.3 裸耦合微调与普适性验证（Universality Check）

在固定比例方案中，为了验证不同格点尺寸是否确实被匹配到了相同的物理尺度，研究团队进行了极其严格的普适性检验（Universality Check）。

• 微调过程：以 $6 \times 3$ 格点、裸耦合 $g^2 = 0.5$、参考距离 $r_1^{(1)} = 2.0$ 为基准起点。对于 $N_x = 7, 8, 9$ 的格点，其对应的参考距离按比例分别设为 $r_1^{(i)} \in \{2.33, 2.67, 3.00\}$。通过对模拟得到的跑动耦合常数进行二次多项式插值，精确微调出各尺寸格点对应的裸耦合常数 $g^2_{(i)}$（参见原文 Fig. 7）。
• 无量纲势能坍缩：在微调完毕后，绘制无量纲势能 $r_1^{(i)} V(r)$ 随相对距离 $r/r_1^{(i)}$ 的变化曲线。令人振奋的是（参见原文 Fig. 8），来自四个不同尺寸格点的数据点完美无缝地重合在同一条连续曲线上。这有力地证明了该有限体积方案能够完美地消除空间离散化带来的异质影响，确保了多格点系统在物理标度上的一致性。

2.4 连续外推性能数据

在截断参数 $\Lambda = 1$ 下，通过将格点纵向尺寸一路扩展至 $N_x = 12$（对应的最精细格点其裸耦合已被微调至极弱的 $g^2 \sim 10^{-2}$ 级别），获得了连续极限外推的关键曲线（参见原文 Fig. 9）。研究人员对比了三种拟合模型：

1. 纯二次模型 $f_1(a) = c_1 + c_2 a^2$： 在最精细的四个格点数据上表现极佳，外推得到的连续极限值为 $\Sigma_{\sqrt{2}}(0) \approx 0.788$。
2. 对数修正模型 $f_2(a) = c_1 + c_2 a^2 \log a$： 由于引入了对数项，导致外推值对拟合窗口的选择高度敏感，说明在当前距离范围内，高阶微扰对数修正并非主导项。
3. 二次+四次混合模型 $f_3(a) = c_1 + c_2 a^2 + c_3 a^4$： 成功将拟合范围延伸至极粗糙的格点，得出的连续极限值与 $f_1$ 保持高度一致。

这表明，基于固定比例方案的格点阶梯标度函数具有极佳的 $O(a^2)$ 收敛行为。

2.5 局部基截断误差与性能分析

为了评估面元 Hilbert 空间截断 $\Lambda$ 对物理结果的污染，研究团队系统分析了截断误差的收敛速度和耦合常数依赖性：

• 指数收敛特征：在 $4 \times 3$ 格点、裸耦合 $g^2 = 0.5$ 的测试中，随着 $\Lambda$ 从 $1$ 增加到 $7$，重正化耦合常数以及基态能量与极限值（以 $\Lambda=7$ 为参考）的偏差呈现出极其标准的指数衰减行为（参见原文 Fig. 11）。在 $\Lambda = 5$ 时，系统能量的绝对误差已降至 $10^{-14}$ 以下，达到了双精度浮点数的数值极限。
• 最大误差区间：通过扫描不同耦合常数下的截断误差（参见原文 Fig. 12），发现截断误差在中等耦合区间（$g^2 \sim 0.3 - 0.5$）达到峰值。这是因为该区域正处于电磁相互作用的过渡交叠区（Crossover Region），电场项与磁场项竞争最激烈，系统基态需要最多的面元基底来表征。而在强耦合（$g^2 \to 2.0$）和极弱耦合（$g^2 \to 0$）极限下，截断误差均快速收敛至极小值。

2.6 横向有限体积效应（Large-$N_y$ 极限）

由于模拟采用的是固定 $N_y=3$ 的条带几何，研究团队在 $6 \times N_y$ 的格点上（横向宽度 $N_y$ 从 3 扫描至 10）深入评估了横向有限体积效应（参见原文 Fig. 13 和 Fig. 14）：

• 随着横向尺寸 $N_y$ 的增加，重正化耦合常数展现出清晰的指数级收敛行为，可以用如下公式进行拟合： $$f(N_y) = g^2(\infty) + A e^{-\kappa(N_y - N_0)}$$ 外推得到的无限大宽度极限值为 $g^2(\infty) = 0.2905$。而在 $N_y=3$ 处测得的值约为 $0.395$，说明极窄的条带几何带来了约 26% 的常数系统偏差。这一系统偏差虽然是一个显著的加性修正，但由于其在连续外推中是结构性固定的，因此完全不会破坏连续极限外推的平滑度与稳定性。

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

3.1 技术栈与核心软件包

本研究所采用的数值模拟计算链条高度规范，核心库均采用国际物理学界开源的成熟项目：

1. ITensor 库（C++版本 / Julia版本）：用于构建哈密顿量的张量算符表示，并执行高精度的矩阵乘积态（MPS）二位置 DMRG 变分基态求解。其优异的奇异值分解（SVD）控制能够支持将截断阈值压低至 $10^{-15}$，这对于捕捉极其微弱的格点离散化效应至关重要。
2. SciPy (scipy.integrate.quad)：用于对连续 U(1) 群的 Haar 测度进行数值积分，从而高精度地计算哈密顿量中各对偶算符（如 $\hat{L}, \hat{L}^2, \hat{P}$）在面元表征下的矩阵元。

3.2 2D 拓扑向一维 MPS 链的映射（蛇形映射）

由于 MPS 在本质上是一维链状张量网络，为了处理二维空间格点，必须采用**“蛇形路径映射”（Snake Path Mapping）**。对于大小为 $N_x \times N_y$ 的物理面元格点，我们为其分配一维 MPS 站点索引 $i_{\text{MPS}}$ 如下：

这种映射方式使得算符的局部相互作用在纵向（$N_x$ 方向）保持温和的局域性，而横向（$N_y$ 方向）的相互作用则会跨越较长的 MPS 距离。这解释了为什么该算法的计算复杂度随 $N_x$ 呈线性增长，而随横向宽度 $N_y$ 呈指数暴涨，也从算法层面论证了为什么在当前经典算力下必须限制 $N_y=3$。

3.3 DMRG 多阶段收敛控制策略

由于在弱耦合区系统极易陷入亚稳态（局部能量极小值），研究团队设计了一套三阶段渐进 DMRG 优化方案：

1. 预调理阶段（Preconditioning Phase）：
   • 执行至少 10 次完整变分扫射（Sweeps）。
   • 设置较为宽松的奇异值截断阈值（$10^{-5}$）。
   • 人为注入一定强度的数值噪声（Noise），强制系统跳出局域势阱，确保探索到全局真正的基态空间。
2. 高精度过渡阶段（High-Precision Phase）：
   • 将 SVD 截断阈值收紧至 $10^{-14}$。
   • 逐渐减小并最终关闭噪声注入，锁定物理基态的大致方向。
3. 极限收敛阶段（Convergence Phase）：
   • 将 SVD 阈值推向极限的 $10^{-15}$。
   • 限制最大键合维度上限为 $\chi_{\max} \sim 500 N_x N_y (2\Lambda + 1)$（在实际模拟中，由于对偶表述的优越性，该上限并未饱和，系统在极低维度下即可完全收敛）。

3.4 核心复现指南与伪代码示例

复现本研究的关键在于如何构建含狄拉克弦背景场的对偶哈密顿量算符（MPO）。以下提供一份基于 Julia 语言和 ITensor 库风格的伪代码指南，用以说明如何构建该对偶哈密顿量：

# 伪代码：构建 (2+1)D U(1) 静态电荷对偶哈密顿量 (以 Lambda=1 为例)
using ITensors

function build_dual_hamiltonian(Nx, Ny, g2, Lambda, charge_path)
    # 1. 初始化面元物理指标，局部维度为 2*Lambda + 1
    d_local = 2 * Lambda + 1
    sites = siteinds("Qdit", Nx * (Ny - 1); dim = d_local)
    
    # 2. 初始化哈密顿量算符 MPO
    ampo = OpSum()
    
    # 3. 构造磁场项 (磁能): 单体局域算符
    # H_M = (1 / 2*g^2) * (2 - P - P^dagger)
    inv_g2_half = 1.0 / (2.0 * g2)
    for i in 1:length(sites)
        ampo += 2.0 * inv_g2_half, "Id", i
        ampo += -1.0 * inv_g2_half, "P", i      # P 算符需提前定义其矩阵元
        ampo += -1.0 * inv_g2_half, "Pdag", i   # P^dagger 算符
    end
    
    # 4. 构造电场项 (电能): 最近邻面元相互作用，并引入背景场 n_gamma
    # H_E = (g^2 / 2) * (L_i - L_j + n_ij)^2
    # 展开后为: (g^2/2) * (L_i^2 + L_j^2 + n_ij^2 - 2*L_i*L_j + 2*n_ij*L_i - 2*n_ij*L_j)
    g2_half = g2 / 2.0
    for x in 0:(Nx-1)
        for y in 1:(Ny-1)
            # 获取当前面元的 MPS 索引以及相邻面元索引
            idx_curr = x * (Ny - 1) + y
            
            # 横向相邻 (y 方向最近邻)
            if y < (Ny - 1)
                idx_next_y = idx_curr + 1
                n_gamma = get_background_field(charge_path, idx_curr, idx_next_y)
                
                ampo += g2_half, "L2", idx_curr
                ampo += g2_half, "L2", idx_next_y
                ampo += g2_half * (n_gamma^2), "Id", idx_curr
                ampo += -2.0 * g2_half, "L", idx_curr, "L", idx_next_y
                ampo += 2.0 * g2_half * n_gamma, "L", idx_curr
                ampo += -2.0 * g2_half * n_gamma, "L", idx_next_y
            end
            
            # 纵向相邻 (x 方向最近邻)
            if x < (Nx - 1)
                idx_next_x = idx_curr + (Ny - 1)
                n_gamma = get_background_field(charge_path, idx_curr, idx_next_x)
                
                ampo += g2_half, "L2", idx_curr
                ampo += g2_half, "L2", idx_next_x
                ampo += g2_half * (n_gamma^2), "Id", idx_curr
                ampo += -2.0 * g2_half, "L", idx_curr, "L", idx_next_x
                ampo += 2.0 * g2_half * n_gamma, "L", idx_curr
                ampo += -2.0 * g2_half * n_gamma, "L", idx_next_x
            end
        end
    end
    
    # 5. 转化为 MPO 格式
    H = MPO(ampo, sites)
    return H, sites
end

3.5 相关开源仓库与软件链接

• ITensor 官方网站及源码仓库：
  • 官方文档：itensor.org
  • GitHub 仓库：github.com/ITensor/ITensors.jl (Julia 版本推荐)
• SciPy 数值计算库：
  • scipy.org

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

4.1 关键引用文献

本项研究工作植根于格点规范理论与重正化群流的经典理论大厦，其最关键的奠基性文献包括：

1. Kogut & Susskind (1975)：《Hamiltonian Formulation of Wilson’s Lattice Gauge Theories》(Phys. Rev. D 11, 395)。 奠定了格点规范理论哈密顿量表述的根基，定义了电磁项的基本演化形式。
2. Lüscher, Weisz, & Wolff (1991)：《A numerical method to compute the running coupling in asymptotically free theories》(Nucl. Phys. B 359, 221)。 首次正式提出并阐明了阶梯标度法（Step-Scaling），成功解决了非阿贝尔规范场论模拟中的“窗口问题”。
3. Sommer (1994)：《A New way to set the energy scale in lattice gauge theories and its applications to the static force…》(Nucl. Phys. B 411, 839)。 提出了著名的 Sommer Scale 设定方案，基于静态电荷作用力定义重正化耦合常数，是本研究“力方案”的理论源头。
4. Unmuth-Yockey (2019)：《Gauge-invariant rotor Hamiltonian from dual variables of 3D U(1) gauge theory》(Phys. Rev. D 99, 074502)。 为三维 U(1) 理论在对偶表述下的张量网络构建提供了直接的方法学启示。

4.2 局限性深度评论

尽管该研究在方法论上取得了极其漂亮的突破，但作为一项先驱性的概念验证（Proof-of-Concept）工作，它在实际定量精度和拓扑拓展上仍存在不可忽视的局限性：

1. 二维条带几何（Strip Geometry）的局限与真实的 2D 连续极限差距

由于一维 MPS 张量网络在表达二维空间纠缠时的物理瓶颈，本研究在进行连续极限外推时，横向宽度被牢牢限制在 $N_y = 3$（以格点步长 $a$ 为单位）。这意味着，在趋近连续极限 $a \to 0$ 的过程中，由于没有相应按比例增加 $N_y$，系统的物理横向宽度 $a N_y$ 实际上是在不断萎缩的。因此，该方案最终外推得到的连续极限，本质上是一个无限长但无限窄的“一维极细条带”上的连续极限，而非真正的二维体相（Bulk）连续极限。虽然研究证明了这一几何局限不会破坏外推的数学平滑度，但若要提取真正无偏的 2D 物理跑动耦合，横向体积必须满足 $a N_y = \text{const.}$，而这在当前的一维 MPS 框架下是完全无法实现的计算灾难。

2. U(1) 阿贝尔群的对偶优越性难以平移至 SU(2) 或 SU(3) 非阿贝尔规范群

本研究所采用的对偶哈密顿量表述之所以能够彻底免除磁场项的量子纠缠瓶颈，得益于 U(1) 规范群的阿贝尔性质——其一维和二维变量可以完美对偶。然而，对于描述真实世界强相互作用（如 QCD）的 SU(2) 和 SU(3) 非阿贝尔规范群，由于其非对易性，无法进行如此干净、无辅助变量的面元对偶。这意味着，在处理非阿贝尔规范场论时，弱耦合极限下的磁场纠缠屏障依然存在，本研究所依赖的低纠缠变分优势将大打折扣。

3. 数值积分精度对极弱耦合区的硬封顶

在微调过程中，为了深入更精细的格点（如 $N_x = 12$），裸耦合常数必须下调至 $g^2 \sim 10^{-2}$ 以下。然而，在此区间内，哈密顿量矩阵元的数值表现极其微弱，已逼近 SciPy 数值积分器的绝对容差极限（epsabs = 1e-8）。若要继续探索更精细的连续区域，积分器的数值不稳定性将成为硬性瓶颈，需要发展全新的解析极限展开或高精度专门积分算法。

5. 其他必要的补充：双重对偶哈密顿量推导与未来量子计算展望

5.1 深入探秘：含静电荷对偶哈密顿量的配分函数导出逻辑

为了让物理和化学工作者更直观地理解对偶哈密顿量中背景场 $n_{\mathbf{x}, i}^{(\gamma)}$ 的数学来源，我们在此将 Appendix A 的物理脉络进行概念级梳理。

在欧氏时空下，引入一对静态电荷 $q$ 与 $-q$ 的规范场配分函数写为：

其中，由于静态电荷不随时间运动，其外部电流只有时间分量非零：$J_0(\vec{x}, t) = q(\delta_{\vec{x}, \vec{x}_1} - \delta_{\vec{x}, \vec{x}_2})$。在时空中，这对电荷分别勾勒出两条平行的 Polyakov 圈（参见原文 Fig. 15）。

利用贝塞尔函数展开式将指数上的余弦作用量转化为离散的面元整数流 $n_{x, \mu\nu}$，并对规范场 $A_{x, \mu}$ 进行分部求和（Summation by Parts），可以得到高斯约束条件：

为了求解这个含源的散度方程，我们将面元流分解为两部分：

1. 第一项 $\epsilon_{\mu\nu\rho} \Delta_\rho m(x)$ 是无源的格点旋度项（其中 $m(x)$ 是定义在对偶格点上的标量场），其散度恒等于 $0$。
2. 第二项 $n_{\mu\nu}^{(q)}(x)$ 是一个特定的含源特解。为了满足电荷源的约束，我们必须选择一个以两条 Polyakov 圈为边界的二维时空片（即狄拉克片，Dirac Sheet，参见原文 Fig. 16）。在该狄拉克片内部，$n_{\mu\nu}^{(q)}(x)$ 取值为 $\pm q$，其余区域为 $0$。

当我们进行 Hodge 对偶（Hodge Duality）将二维时空的面元流映射为对偶格点上的链路流时，这一时空狄拉克片在纯空间格点上就投影为一条连接两静态电荷的非动力学路径 $\gamma$（狄拉克弦）。在连续时间极限（$\beta_t \to \infty$）下，这一路径直接转化为了哈密顿量中耦合在电场项上的整数背景场 $n_{\mathbf{x}, i}^{(\gamma)}$（参见原文 Eq. 16），从而在哈密顿量形式下完美再现了配分函数的全部物理物理效应。

5.2 对比分析：“正面方案”与“反面教材”的戏剧性反差

为了凸显固定比例外推方案的决定性优势，原作者在附录 B 中专门展示了如果不进行比例锁定、采用常规方案外推所带来的毁灭性后果：

• 常规方案行为（Fig. 18）：在不改变参考距离格点数 $r/a$ 的情况下，直接在 $3 \times 3$ 到 $7 \times 7$ 的格点上测量步进函数。此时，随着格点精细化，物理体积急剧改变。计算得到的 $\Sigma_s(a)$ 曲线呈现出完全杂乱无章、非单调的折线分布。在极小体积效应的严重干扰下，连续极限趋势被彻底掩盖，外推尝试完全失效。
• 固定比例方案行为（Fig. 9）：相同格点尺寸下，一旦开启固定比例锁定，有限体积效应被规整为单一、单调的高阶离散项。外推数据点以极高精度收敛于一条平滑的二次抛物线，展现出无可挑剔的数学美感。

这一剧烈对比确立了该方案在未来小体积物理模拟中的统治性地位。

5.3 跨界展望：从张量网络走向 NISQ 与 FTQC 量子计算时代

格点规范理论的哈密顿量变分模拟，不仅是经典张量网络算法竞技的舞台，更是未来通用量子计算机（Quantum Computers）最具杀手级应用前景的领域。在当前的含噪中等规模量子（NISQ）时代，量子比特的数量和门保真度极其有限，我们同样只能在极小的空间格点上运行规范场论。

本研究所提出的固定比例有限体积外推方案，具有以下极其适合量子硬件的优良属性：

1. 极低的空间需求：方案被证实即使在纵向尺寸 $N_x = 6$、宽度 $N_y = 3$ 的超微型格点上（仅需数十个量子比特），依然能获得质量极高的连续外推趋势。这与当前超导、离子阱等主流昂贵的硬件资源规模高度契合。
2. 完全免疫符号问题：直接在量子态上进行能量观测，能够直接外推包含静态物质流乃至未来动态费米子物质流的物理耦合，这是经典蒙特卡洛算法无法企及的领域。
3. 更先进的张量网络架构支撑：为了解决 $N_y=3$ 窄条带局限性，未来在经典端可以使用更先进的 genuinely 二维张量网络架构，如投射纠缠配对态（PEPS）、等距张量网络（isoTNs）或树状张量网络（TTNs）。这些新架构与本文提出的固定比例方案相结合，将彻底打通（2+1）维乃至真实（3+1）维非阿贝尔规范场论在极小体积下的高精度重正化研究道路，让强相互作用的非微扰定量计算迎来真正的黄金时代。