来源论文: https://arxiv.org/abs/2603.24698v1 生成时间: Mar 27, 2026 12:11

0. 执行摘要

核物理和粒子物理的核心问题之一是理解量子色动力学（QCD）中的强相互作用下粒子是如何形成和演化的，特别是被称为“强子化”的过程。弦断裂是强子化模型的核心机制，它描述了色荷（夸克和胶子）如何通过胶子流形成的“弦”在足够长的距离上断裂，产生新的介子或重子。然而，由于“符号问题”（sign problem），基于传统蒙特卡洛的格点规范理论（LGT）计算难以在实时间（闵可夫斯基时间）研究这一复杂现象。

本研究论文提出并应用了一种基于Loop-String-Hadron (LSH) 形式化的张量网络工具包，用于研究具有动力学费米子的 (1+1) 维 SU(2) LGT。LSH 形式化通过将规范不变性内在化到局域自由度中，有效地解决了非阿贝尔规范理论中的高维度和规范约束挑战。该工具包首次被用于在连续统和热力学极限下精确测定静态弦势和弦张力。更重要的是，它深入探讨了在猝灭动力学中动力学弦断裂的实时演化。通过详细分析 LSH 激发态的微观演化，研究揭示了弦膨胀与收缩、端点分裂、粒子簇射、链散射、弦解离与复合以及粒子生成等复杂动力学过程。这些微观洞察与能量传输、纠缠熵生成和关联扩散等宏观特征紧密联系。该工作为未来利用张量网络研究更复杂格点规范理论中的弦断裂和粒子生成开辟了新途径，为最终理解 (3+1) 维 QCD 中的强子化过程奠定了基础。

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

1.1 核心科学问题：弦断裂与强子化

量子色动力学（QCD）是描述夸克和胶子强相互作用的基本理论，也是粒子物理标准模型的重要组成部分。在粒子对撞实验中，高能碰撞产生的夸克和胶子并非孤立存在，而是迅速地形成无色粒子——强子。这一过程被称为“强子化”，是 QCD 动力学中最复杂的现象之一。弦断裂理论是理解强子化机制的关键组成部分。它假设当两个色荷（例如夸克和反夸克）通过胶子场连接形成一个色流管（即“弦”）时，如果弦的长度超过某个临界值，那么真空中的夸克-反夸克对会通过 Schwinger 机制自发产生，从而“断裂”这条弦，形成两个或多个较短的弦或强子。理解弦断裂的静态特性（如弦张力、临界断裂长度）和动态过程（如弦如何随时间演化、断裂产物的生成与传播）对于构建精确的强子化模型至关重要。

然而，从第一性原理出发对 QCD 中的弦断裂动力学进行实时模拟面临巨大挑战。传统的基于蒙特卡洛方法的格点规范理论（LGT）计算，在处理闵可夫斯基时间（实时间）演化时会遭遇严重的“符号问题”。这使得对强子化这一非平衡态、实时动力学过程的直接模拟变得极为困难，甚至无法实现。因此，研究人员转而寻求在简化模型（例如低维度格点规范理论）中利用新的计算工具来深入理解这些基本物理过程。

本研究工作的核心目标就是利用张量网络（Tensor Networks, TNs）这一新兴计算方法，在 (1+1) 维 SU(2) 格点规范理论（一种非阿贝尔规范理论）中，对弦断裂的静态和动态性质进行系统且全面的研究。这一简化模型虽然不是完整的 (3+1) 维 QCD，但它保留了非阿贝尔规范理论的关键特征，并能提供在实时动力学方面的宝贵洞察，为未来向更高维、更复杂模型的拓展奠定基础。

1.2 理论基础：Kogut-Susskind 哈密顿量与 Loop-String-Hadron (LSH) 形式化

Kogut-Susskind 哈密顿量形式化

(1+1) 维 SU(2) 格点规范理论的哈密顿量形式化最初由 Kogut 和 Susskind 提出。在这个形式化中，物理系统被定义在离散的格点上，规范场（胶子）被放置在连接格点的链路上，而物质场（费米子，如夸克）则放置在格点上。SU(2) 规范群的链路上自由度通常用角动量运算符表示，费米子则用费米子算符表示。哈密顿量包含三个主要部分：电场能量项（胶子相互作用）、质量能量项（费米子质量）和相互作用项（费米子与胶子的耦合）。规范不变性通过高斯定律（Gauss’s Law）来保证，要求物理态在高斯定律算符的作用下为零，这使得系统具有严格的局域规范不变性。然而，Kogut-Susskind 形式化中的高斯定律是非阿贝尔的，这意味着它们之间不相互对易，这给构建显式规范不变基态带来了巨大的挑战，尤其是在将这些基态与张量网络方法结合时。

Loop-String-Hadron (LSH) 形式化

为了克服 Kogut-Susskind 形式化中非阿贝尔高斯定律带来的挑战，本研究采用了 Loop-String-Hadron (LSH) 形式化 [114]。LSH 形式化是对 Kogut-Susskind 哈密顿量形式化的一种重构，其核心思想是直接构建内在规范不变的自由度。它通过将 SU(2) 角动量自由度表达为 Schwinger 玻色子双态，并与 SU(2) 费米子双态结合，从而形成局域 SU(2) 单态：

  • 环 (Loop, n_l)： 由两个玻色子组成的局域规范不变单态。它们代表了局域的威尔逊环快照、介子弦的内部区域或背景通量配置。n_l 是一个非负整数，表示环激发的总数。
  • 弦-入/弦-出 (String-in/String-out, n_i, n_o)： 由一个费米子和一个玻色子组成的局域规范不变单态。它们表示费米子在当前格点上，并向左（n_i）或向右（n_o）产生一个单位通量，是扩展弦的端点。n_i, n_o 是二进制量子数，取值 0 或 1。
  • 强子 (Hadron)： 由两个费米子或两个玻色子组成的局域规范不变单态。在 LSH 语境中，n_i = n_o = 0 表示费米子真空态，而 n_i = n_o = 1 表示一个强子激发。

因此，每个格点上的局域 Hilbert 空间由 |n_l, n_i, n_o⟩_r 这些量子数来标记。这些局域基态本质上是规范不变的，并且非阿贝尔高斯定律约束已经内在地被满足。这意味着在构建张量网络时，我们只需要强制执行阿贝尔高斯定律（AGL），即电通量的连续性，这大大简化了算法的复杂性。AGL 要求 N_L(r) = N_R(r+1)，其中 N_L(r) 和 N_R(r) 是与格点 r 左右链路相关的总角动量（或电通量）。

LSH 哈密顿量 Ĥ 被重写为包含 n_l, n_i, n_o 算符的电场能量项 (Ĥ_E)、质量能量项 (Ĥ_M) 和相互作用项 (Ĥ_I)，这些项都是由局域算符和最近邻算符构成。该理论还存在两个守恒的全局量子数：总费米子数 Q = Σ_r (n_i(r) + n_o(r)) 和净通量 q = Σ_r (n_o(r) - n_i(r))。

1.3 技术难点与方法细节：张量网络与实时演化

张量网络 (Tensor Networks, TNs)

张量网络，特别是矩阵乘积态（Matrix Product States, MPS）和矩阵乘积算符（Matrix Product Operators, MPO），已成为量子多体物理中模拟一维量子系统基态和实时动力学的强大工具。它们的优势在于能够高效地表示具有局域纠缠结构的量子态，即使物理 Hilbert 空间维度随系统尺寸指数增长，MPS 的计算复杂度也仅为多项式级别。这是因为它们利用了纠缠熵在能量本征态中满足“面积定律”的特性，即纠缠熵不随系统大小线性增长。

LSH-Based MPS 和 MPO 构造

本研究将 LSH 形式化与 MPS/MPO 相结合。一个 MPS 态 |Ψ⟩ = Σ A^(p1...pN) |p1...pN⟩，其中 |pr⟩ 是每个格点 r 上的 LSH 局域基态 |n_l, n_i, n_o⟩_r。MPS 的张量 A 具有虚拟指标，其维度被称为“键维度”（bond dimension, D）。D 限制了系统可以承载的最大纠缠量。MPO 类似地表示算符，例如哈密顿量。

规范不变性与对称性

LSH 形式化的关键优势在于其局域规范不变性。通过 LSH 基态，非阿贝尔高斯定律已经内在地被满足。张量网络实现中只需要确保阿贝尔高斯定律（AGL）的满足。这通过在 MPS/MPO 张量中强制实施局域守恒规则来实现，这些规则确保了虚拟键上的总电荷守恒。此外，系统的全局守恒量子数 Q 和 q（总费米子数和净通量）也被整合到 MPS 结构中，使得张量具有块对角线形式，进一步提高了计算效率和稳定性。

希尔伯特空间截断

为了使计算在实际中可行，需要对局域希尔伯特空间进行截断。主要截断参数是电通量截断 j_max，它限制了 n_l 的最大值（LSH 形式中 N_L/R(r) ≤ 2j_max，对应 n_l(r) 最大值）。文章证明 j_max = 5/2 足以确保动力学计算的收敛性，因为 j_max=1 会导致定性不同的结果。

静态势能计算

  1. 引入静态电荷： 通过修改 AGL N_L(r) - N_R(r+1) = Q(r, r+1) 来引入静态电荷，其中 Q(r1, r1+1) = -1 和 Q(r2-1, r2) = +1 表示在 r1 和 r2 处放置了一对静态色荷。为了强制执行这些修改后的 AGL，在哈密顿量中添加一个惩罚项 Ĥ_p。
  2. 基态计算： 使用密度矩阵重正化群（DMRG）算法计算在存在和不存在静态电荷情况下的基态能量 E(gl) 和 E(0)。DMRG 通过迭代优化 MPS 张量来找到最低能量本征态。
  3. 外推法： 对 DMRG 结果进行外推，以消除有限键维度（D）、有限通量截断（j_max）、有限系统尺寸（N）和有限格点间距（ga）引入的误差，最终获得连续统和热力学极限下的弦张力。

实时动力学计算

  1. 初始态准备： 首先，使用 DMRG 计算出零静态电荷扇区的相互作用真空态 |Ω⟩。然后，通过应用一个 Dirac 介子弦算符 Ŝ_Dirac^(r,∆r)（由 LSH 算符构成）到 |Ω⟩ 来制备初始的介子弦激发态 |Φ(t=0)⟩。这个算符在 r 处产生一个费米子激发，在 r+∆r 处湮灭一个费米子激发，并通过一系列环算符连接。
  2. 时间演化： 初始态 |Φ(0)⟩ 随后通过单点时间依赖变分原理（Single-site Time-Dependent Variational Principle, TDVP）算法进行时间演化，即 |Φ(t)⟩ = e^(-itĤ)|Φ(0)⟩。TDVP 算法结合了 Krylov 基扩展协议，可以在演化过程中动态扩展键维度 D_max，以适应纠缠熵的增长。
  3. 可观测量： 在时间演化过程中，计算各种局域 LSH 激发态的期望值（如总物质激发 ⟨n(r)⟩、弦-入 ⟨n_i(r)⟩、弦-出 ⟨n_o(r)⟩、重子 ⟨n_b(r)⟩、环 ⟨n_l(r)⟩ 和对称通量 ⟨N(r)⟩）、真空扣除的能量密度 h(r;t)、双向纠缠熵 S_VN(r;t) 以及两点关联函数 C_hE(r,s;t) 和 C_n(r,s;t)。
  4. 收敛性分析： 研究了 j_max 和 D_max 对动力学结果的影响，以确保计算的可靠性。对于本文中的轻质量费米子，最大模拟时间 gt=8 处的相对误差为 10%，而重质量费米子则收敛得更好，可以模拟到 gt=16。

综上所述，本研究通过将 LSH 形式化与先进的张量网络方法相结合，提供了一个强大且系统化的框架，用于在非阿贝尔格点规范理论中研究复杂的强子化物理，尤其是在实时动力学领域，这是传统方法难以触及的。

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

本研究在 (1+1)D SU(2) 格点规范理论中，针对静态弦势和实时弦断裂动力学进行了详细的计算与分析。以下将介绍关键的 Benchmark 体系、计算获得的数据及其物理意义，以及对性能和收敛性的评估。

2.1 静态弦势与弦张力

Benchmark 体系设置

为了研究静态弦势，研究引入了一对静态色荷放置在格点 r1 和 r2 处。这些静态电荷通过修改阿贝尔高斯定律 (AGL) 来实现，即 N_L(r) - N_R(r+1) = Q(r, r+1)，其中 Q(r1, r1+1) = -1 且 Q(r2-1, r2) = +1，其他位置为零。为了强制系统保持在这一特定的静态电荷扇区，在哈密顿量中添加了一个大的惩罚项 Ĥ_p。计算在固定费米子质量 m/g = 0.5 的情况下进行，通过 DMRG 算法获得存在和不存在静态电荷时的基态能量 E(gl) 和 E(0)。静态势能定义为 V(gl) = E(gl) - E(0)，其中 gl = ga(r2 - r1) 是物理弦长。

计算所得数据与物理现象

  • 弦势曲线 (图 3a)： 静态势能 V(gl)/g 作为物理弦长 gl 的函数绘制。曲线呈现出经典的禁闭行为：在短距离 gl 处，势能大致呈线性增长，这表明静态色荷之间存在一个连续的通量管（弦）。随着 gl 的增加，当达到某个临界分离距离时，势能曲线开始平坦化，形成一个平台。这标志着弦断裂的发生：未断裂弦的能量与断裂弦配置（由动态色荷形成较小的通量管并屏蔽中心区域的电场）的能量变得可比。对于较小的系统体积 gL=14，平台在 gl > 7 后略微向下倾斜，这是有限体积效应。

  • 通量截断 j_max 的影响： 结果显示，对于最低截断 j_max = 1，静态势能的大小和线性部分的斜率对格点间距 ga 有显著依赖性。然而，对于 j_max = 3/2 和 j_max = 2，结果几乎重合，且对 ga 的依赖性很小。这表明能量在 j_max = 3/2 或 j_max = 2 时已经接近收敛，对于更大的 j_max 值，结果不会有显著变化。临界分离距离似乎对 ga 不敏感，但对于 j_max = 1，临界分离距离较短。

  • 弦张力 κ/g^2： 弦张力定义为静态势能线性部分的斜率，它是理论的一个自然质量尺度，也是禁闭的诊断。通过在 gl ≈ 3 到 gl ≈ 7 范围内的线性拟合，并在热力学和连续统极限下进行外推（如图 3b 所示），获得了费米子质量 m/g = 0.5 时的连续统弦张力估计值 κ/g^2 = 0.330(3)。这个值与 Ref. [109] 中使用高斯变分方法在较大 ag 值下估计的 0.3-0.4 的范围一致。

  • 微观激发态 (图 4)： 通过计算真空扣除的总物质激发 ⟨n(r)⟩_s、弦-入 ⟨n_i(r)⟩_s、弦-出 ⟨n_o(r)⟩_s、重子 ⟨n_b(r)⟩_s、环 ⟨n_l(r)⟩_s 和对称通量 ⟨N(r)⟩_s 的期望值，揭示了未断裂弦和断裂弦基态的微观差异。对于较短的弦（gl=3），在弦两端附近观察到弦-入和弦-出激发交替模式，这表明产生了最小长度的介子弦对。在弦内部，弦-入激发主导。对于较长的弦（gl≈8），激发数幅度更大，且通量在中间区域几乎耗尽，这对应了断裂弦的概念。

2.2 实时弦断裂动力学

Benchmark 体系设置

为了研究实时弦断裂动力学，首先计算了参数集 N=128, x=16, j_max=5/2 下的相互作用真空态 |Ω⟩。然后，通过将 Dirac 介子弦算符 Ŝ_Dirac^(r,∆r) 作用于 |Ω⟩ 来制备初始态 |Φ(t=0)⟩，其中 r=55, ∆r=16。这个初始态并非哈密顿量的本征态，因此会随时间非平凡地演化。研究重点关注两种不同的裸费米子质量 m/g = 0.2（轻质量）和 m/g = 2.0（重质量），在 gt=8（轻质量）和 gt=16（重质量）的总演化时间内进行模拟。

计算所得数据与物理现象

  • 相互作用真空态 (图 5)： 相互作用真空态在两端存在边界效应，但在中间区域保持平移不变。轻质量真空的总激发数 ⟨n⟩ 远大于重质量真空，这表明轻质量系统更容易产生激发。对于轻质量和重质量，真空中的弦-入和弦-出激发数 ⟨n_i⟩ 和 ⟨n_o⟩ 相等，符合平移不变性。

  • 初始介子弦态 (图 6, gt=0 面板)： 初始介子弦态通过 Dirac 介子弦算符的施加而形成。它在弦的两端（r≈55, 71）产生局域的物质激发，这些激发在相互作用真空中是真空扣除的。弦内部区域是一个非平衡区域，具有背景通量。与静态弦不同，动态弦的环激发在动态边缘附近耗尽，但在内部未耗尽，因为动态弦的端点会扩散。静态弦的环激发则在静态边缘附近达到峰值，在内部耗尽。

  • 微观动力学：LSH 激发态演化 (图 6 & 8)

    • 弦膨胀与收缩 (图 7a)： 初始动力学主要由弦-出和弦-入态在介子弦两端的内向运动（收缩）和外向运动（膨胀）驱动。弦的收缩导致弦内部通量从边缘开始耗尽。外向运动将通量分布到弦的外侧区域。
    • 粒子簇射 (Particle Showers)： 费米子在整个演化过程中持续地左右跳跃，导致每个粒子轨迹反复分裂，产生粒子“簇射”。后续的“簇射”携带较少的平均物质激发，但速度大致相同。
    • 重子运动与生成： 弦-入/弦-出态的跳跃也伴随着重子的输运。对于重质量，向外运动的弦-入/弦-出和重子波前继续无干扰地向边界移动，内部碰撞后净弦-入/弦-出差值 (n_i-n_o) 几乎消失，但重子仍然持续追踪端点运动。对于轻质量，当弦-入和弦-出波前碰撞时 (gt≈2)，(n_i-n_o) 差值消失，同时伴随着更多重子的产生，这表明弦解离过程 (图 7c) 在发挥作用。
    • 弦解离与复合 (图 7c)： 对于轻质量，内部波前碰撞后，中心格点的通量继续快速耗尽，这是弦解离的标志。随后的介子生成过程在有任意数量环存在的情况下，产生净弦-入/弦-出激发，这些激发也会移动和分裂。这导致了激发簇射的增殖及其在体内的后续散射。
    • 质量依赖性差异： 轻质量和重质量动力学存在显著差异。轻质量情况下，电通量快速耗尽，弦解离普遍，重子大量增殖，发生非弹性链散射事件，能量呈弹道式输运。重质量情况下，电通量耗尽较慢，弦解离受到抑制，粒子生成减少，电通量部分屏蔽，能量传输呈非弹道和非扩散性。

  • 能量、熵与关联演化 (图 9, 10, 11, 12)

    • 能量扩散 (图 9a,e & 10a,i)： 初始时，介子弦在两端具有能量峰值，内部区域能量密度较小且均匀。随着时间推移，这些能量包向弦的内外区域扩散。轻质量的峰值能量较低，扩散速度较快，电场能量密度 h_E(r;t) 在弦内部耗尽并扩散到外部，质量能量密度 h_M(r;t) 在内部累积，总电能接近零 (gt=8)。重质量的峰值能量较高，扩散速度较慢，电场能量密度初期耗尽，但随后扩散缓慢，总电能和质量能都趋于 ≈25 (gt=16)。
    • 纠缠熵生成 (图 9d,h & 10d,l,h,p)： 真空扣除的双向纠缠熵 S_VN(r;t) 在初始时在弦两端迅速生成。对于轻质量，缠结熵在中间区域迅速增加 (gt≈2 处主散射事件)，最终趋于饱和，总纠缠熵约为重质量的两倍。对于重质量，纠缠熵生成相对较慢，在中心区域的生成也较少。轻质量的缠结熵边界扩张速度 ≈ 3.58 格点/gt，与能量前沿速度接近。
    • 对称性分辨 Schmidt 权重 w(q) (图 11)： w(q) 表示半系统出射净通量为 q 的概率。初始时，q=1 扇区占主导。对于轻质量，q=0 扇区在 gt≈2 时期（主散射事件时）最终超过 q=1 扇区，表明弦断裂（形成分离的介子对）。对于重质量，q=1 扇区在整个模拟时间内仍然占主导，与部分电场屏蔽相符。
    • 两点关联函数 (图 12)： 电场能量关联函数 C_hE(r,s;t) 和物质激发关联函数 C_n(r,s;t) 提供了关联传播的视图。初始时，关联集中在弦的对角线附近和两端。对于轻质量，关联结构迅速致密并扩散到整个内部区域，在 gt≈2 和 gt≈4 处产生显著关联。对于重质量，关联传播更有序，在主碰撞点附近没有产生显著关联，反映了连接通量管的持续存在。

2.3 性能数据与收敛性分析

本研究详细讨论了计算结果在不同参数下的收敛性，这对于确保结果的可靠性至关重要。

  • 通量截断 j_max 的收敛性 (图 17, 18)： 动态模拟结果对通量截断 j_max 有显著依赖。当 j_max 从 1/2 增加到 3/2 时，动力学行为（如能量密度、纠缠熵）发生剧烈变化，甚至呈现出截然不同的定性行为（例如，j_max=1/2 显示出振荡行为和受限传输，而 j_max=5/2 显示出连续扩散和弦断裂）。然而，当 j_max 从 3/2 增加到 5/2 时，结果差异很小，表明 j_max=5/2 足以实现收敛。这强调了在模拟非阿贝尔规范理论时选择足够大的 j_max 的重要性。

  • 键维度 D_max 的收敛性 (图 19, 20)： 张量网络方法中的键维度 D_max 限制了 MPS 可以表示的纠缠量。在实时演化中，纠缠熵通常会线性增长，因此 D_max 必须足够大才能精确捕捉动力学。研究使用了三个不同的 D_max 值（100, 200, 600）进行了计算。

    • 重质量 m/g=2： 对于重质量情况，两个最大的 D_max 值（200 和 600）在整个模拟时间内（gt=16）几乎重合，总电能的相对误差在 gt=16 时仅为亚百分比。这意味着对于重质量，D_max=200 已经足以获得收敛结果。
    • 轻质量 m/g=0.2： 对于轻质量情况，相对误差更大，在 gt=8 时达到百分之几。这表明轻质量动力学中纠缠熵的增长更快，需要更大的 D_max 才能在更长时间内保持高精度。因此，轻质量的模拟时间被限制在 gt=8。
    • 结论： 尽管轻质量结果在 gt=8 处积累了一些误差，但总的来说，张量网络方法在所选 D_max 范围内成功捕捉了 LSH 激发态的演化，并提供了对弦断裂动力学的深入洞察。

总体而言，本研究在 (1+1)D SU(2) LGT 中，不仅首次精确测量了弦张力，而且通过 LSH 形式化和张量网络工具，首次在微观层面详细揭示了实时弦断裂动力学的复杂过程，并对关键参数的收敛性进行了严格的评估，为未来更复杂、更精确的 QCD 模拟奠定了坚实基础。

3. 代码实现细节，复现指南与所用软件包

本研究的核心在于将 Loop-String-Hadron (LSH) 形式化与张量网络 (Tensor Networks, TNs) 方法相结合，以实现对 (1+1)D SU(2) 格点规范理论的静态和动态模拟。以下将详细阐述其代码实现细节、复现指南以及所用的软件包。

3.1 代码实现细节

3.1.1 LSH 哈密顿量与基态表示

  • 局域 Hilbert 空间： LSH 形式化的关键在于其局域规范不变的自由度。每个格点 r 的局域 Hilbert 空间由三个量子数 (n_l, n_i, n_o) 定义，分别代表环、弦-入和弦-出激发数。n_l 是非负整数（由 j_max 截断），n_i, n_o 为 0 或 1。这些量子数构成了每个格点上有限维度的基。
  • 哈密顿量算符： Kogut-Susskind 哈密顿量被重写为 LSH 算符的形式，即 Ĥ = Ĥ_E + Ĥ_M + Ĥ_I。这些项由局域的 n_l, n_i, n_o 数算符和阶梯算符 (Â^±, Ŝ^±, χ^±) 构成。这些算符的具体作用方式（如对 n_l 加减一，对 n_i, n_o 进行翻转等）在附录 2 中详细定义。这些算符构成了作用在 MPS 上的 MPO 的基本构建块。
  • 规范不变性与对称性： LSH 基态本质上是规范不变的，因此无需显式处理非阿贝尔高斯定律。代码实现中，只需要强制执行阿贝尔高斯定律（AGL），即 N_L(r) = N_R(r+1)。这通过在构建 MPS/MPO 张量时引入局域守恒规则来实现，确保虚拟键上的总电荷守恒。此外，全局守恒量子数 Q (总费米子数) 和 q (净通量) 也被整合到张量网络结构中，使得张量具有块对角线形式，这不仅提高了计算效率，也确保了物理态的对称性。

3.1.2 张量网络结构：MPS 和 MPO

  • 矩阵乘积态 (MPS)： 量子态 |Ψ⟩ 被表示为 MPS 的形式，即一系列张量 A_r 的乘积。每个 A_r 张量有一个物理腿（对应局域 LSH 基态 |p_r⟩）和两个虚拟腿（连接相邻格点）。这些张量 A_r 内部的元素在存储时会组织成块对角线形式，以体现全局守恒量子数 Q 和 q 的守恒。
  • 矩阵乘积算符 (MPO)： 哈密顿量 Ĥ 以及其他可观测算符（如介子弦算符 Ŝ_Dirac）被表示为 MPO 的形式。MPO 同样由一系列张量 W_r 构成，每个 W_r 有两个物理腿（对应输入和输出的局域基态）和两个虚拟腿。MPO 的张量同样利用块对角线结构来编码规范不变性和守恒量。

3.1.3 算法选择

  • 基态计算： 采用密度矩阵重正化群 (DMRG) 算法来寻找哈密顿量的基态。DMRG 通过迭代优化 MPS 张量来收敛到最低能量本征态。在静态弦势计算中，为了强制静态电荷，哈密顿量中添加了惩罚项，并使用 DMRG 寻找这个“修正”哈密顿量的基态。
  • 实时演化： 采用单点时间依赖变分原理 (Single-site Time-Dependent Variational Principle, TDVP) 算法进行实时演化。TDVP 是一种高效的算法，用于在 MPS 流形上演化量子态。它结合了 Krylov 基扩展协议，允许在演化过程中动态增加 MPS 的键维度 D，以适应纠缠熵的增长，从而在保持计算效率的同时提高精度。论文指出，对于轻质量情况，DMRG 初始基态的截断误差设置为 ε = 10^-8，TDVP 最大键维度 D_max = 600。

3.1.4 希尔伯特空间与键维度管理

  • 通量截断 j_max： LSH 形式中，玻色子环量子数 n_l 理论上可以取任意非负整数。为了进行数值计算，必须施加一个上限 j_max。论文中 j_max 被设置为 5/2，这对应于 LSH 形式中最大环量子数 n_l,max = 5，导致每个玻色子自由度的局域 Hilbert 空间维度为 6。对于 n_i, n_o (0或1)，最终局域 Hilbert 空间维度为 d = 8j_max + 2，在 j_max = 5/2 时为 22。文章强调了 j_max 的选择对动力学结果的定性影响。
  • 键维度 D_max： MPS 的键维度 D 是控制计算精度和资源消耗的关键参数。在 DMRG 中，通过设置截断误差 ε 来控制 D。在 TDVP 实时演化中，D 可以动态增加，最大不超过 D_max。论文中对 D_max 进行了收敛性研究，发现 D_max = 600 对于轻质量在 gt=8 内，对于重质量在 gt=16 内都能提供足够好的收敛性。

3.2 复现指南

要复现本研究结果，需要遵循以下步骤，并熟悉张量网络的基本概念和编程。

  1. 环境设置：

    • 安装 Julia 编程语言（ITensors.jl 是一个 Julia 包）。
    • 安装 ITensors.jl 包：在 Julia 环境中运行 Pkg.add("ITensors")。
    • 理解 LSH 形式化：仔细阅读论文附录 2，理解 LSH 基态、算符及其作用规则，以及如何将 Kogut-Susskind 哈密顿量映射到 LSH 形式。

  2. LSH 哈密顿量构建：

    • 定义局域 Hilbert 空间： 根据 j_max 设置每个格点 r 的局域 Hilbert 空间 sites[r]，该空间需要能够表示 (n_l, n_i, n_o)。ITensors.jl 允许用户自定义 SiteType 来实现这一点。
    • 构建哈密顿量 MPO： 根据论文 Eq. (4) 和 Eq. (5) 中定义的 LSH 哈密顿量，构建对应的 MPO。这需要实现所有基本 LSH 算符（如 ĥ_l, ĥ_i, ĥ_o, Ŝ^++ 等）及其在局域基态上的作用规则。注意处理 1/√(N_L(r)+1) 等因子，这些因子会导致算符的 MPO 结构稍微复杂。
    • 整合对称性： 在定义 SiteType 和构建 MPO 时，确保将全局量子数 Q 和 q 的守恒规则整合进去，使张量自动具有块对角线结构。

  3. 静态弦势计算 (参考 Sec. III)：

    • 基态计算（无静态电荷）：
      • 使用 DMRG 算法计算 LSH 哈密顿量在 Q=N, q=0 扇区（相互作用真空态）的基态 |Ω⟩。设置合适的 D_max 和收敛精度 ε。
      • 记录基态能量 E(0)。
    • 基态计算（有静态电荷）：
      • 修改 AGL： 定义一个修改后的哈密顿量，其中包含惩罚项 Ĥ_p (Eq. (25))，以及在 r1 和 r2 处通过 AGL 强制引入静态电荷。这通常通过在 sites 定义中为 r1 和 r2 指定特殊的 Q(r, r+1) 值来实现。
      • 使用 DMRG 算法计算这个修改后的哈密顿量的基态 |Φ_static⟩。
      • 记录基态能量 E(gl)。
    • 计算静态势能： V(gl) = E(gl) - E(0)。
    • 外推法： 对不同 N, ga, j_max, D_max 的结果进行外推，以获得连续统和热力学极限下的弦张力。

  4. 实时动力学计算 (参考 Sec. IV)：

    • 初始态制备：
      • 构建 Dirac 介子弦算符 Ŝ_Dirac^(r,∆r) 的 MPO (Eq. (30))。这个算符作用在真空态 |Ω⟩ 上，产生初始态 |Φ(t=0)⟩ = (1/N_0) Ŝ_Dirac^(r,∆r) |Ω⟩。
      • N_0 是归一化因子。可以通过 inner(Ŝ_Dirac^(r,∆r) |Ω⟩, Ŝ_Dirac^(r,∆r) |Ω⟩) 来计算。
    • 时间演化：
      • 使用 TDVP 算法对 |Φ(t=0)⟩ 进行时间演化。选择合适的 dt 和总演化时间 gt。
      • 在 TDVP 循环中，允许键维度动态增长，但要设置 D_max 的上限（例如 600）。
    • 可观测量计算：
      • 在每个时间步，计算各种局域 LSH 激发态的期望值 ⟨n(r)⟩ 等，以及能量密度 h(r;t)、纠缠熵 S_VN(r;t) 和两点关联函数 C_hE, C_n。这通常通过构建这些算符的 MPO，然后计算 inner(Φ(t), MPO_O, Φ(t)) 来完成。

本研究主要依赖于 ITensors.jl 软件包 [124]。

  • ITensors.jl： 这是一个用 Julia 语言编写的开源张量网络库。它提供了强大的工具来构建和操作 MPS、MPO 以及其他张量网络，并实现了多种算法，包括 DMRG 用于寻找基态，以及 TDVP 用于实时演化。ITensors.jl 的设计允许用户方便地定义自定义的 SiteType 和物理算符，这对于实现 LSH 形式化中的复杂局域自由度和算符至关重要。其对对称性的内置支持（通过 qn (量子数) 系统）也使得将 LSH 的全局守恒量整合到张量网络中变得相对简单。

开源 Repo Link：

论文中并未提供直接指向其具体代码实现的开源仓库链接。然而，ITensors.jl 本身是一个活跃开发的开源项目，其官方文档和 GitHub 仓库是理解如何使用该库实现类似计算的最佳资源。

  • ITensors.jl GitHub 仓库： https://github.com/ITensors/ITensors.jl
  • ITensors.jl 官方文档： https://itensors.github.io/ITensors.jl/stable/

用户可以在 ITensors.jl 的文档中找到关于自定义 SiteType、构建 MPO、使用 DMRG 和 TDVP 的详细教程和示例。通过这些资源，结合论文中对 LSH 形式化、哈密顿量算符和算法细节的描述，理论上可以复现本文中的计算。虽然没有直接的代码，但论文提供了足够的细节来指导熟悉 ITensors.jl 的研究人员进行实现。

4. 关键引用文献与工作局限性评论

4.1 关键引用文献综述

本研究工作在 (1+1)D SU(2) 格点规范理论中对弦断裂的静态与动态性质进行了深入探讨，其成果建立在多个关键领域的前沿进展之上，并引用了大量重要文献。

4.1.1 LSH 形式化

Loop-String-Hadron (LSH) 形式化是本研究的理论基石，它将非阿贝尔规范理论的规范不变性内在地编码到局域自由度中，极大地简化了张量网络的构建。核心引用文献包括：

  • [114] Indrakshi Raychowdhury and Jesse R. Stryker, “Loop, String, and Hadron Dynamics in SU(2) Hamiltonian Lattice Gauge Theories,” Phys. Rev. D 101, 114502 (2020)。 这篇论文首次详细介绍了 SU(2) LSH 形式化，为本研究提供了理论框架。
  • [117] Saurabh V. Kadam, Indrakshi Raychowdhury, and Jesse R. Stryker, “Loop-string-hadron formulation of an SU(3) gauge theory with dynamical quarks,” Phys. Rev. D 107, 094513 (2023)。 LSH 形式化被扩展到 SU(3) 规范群，展示了其普适性。
  • [118, 119] 进一步阐述了 LSH 形式化在更高维度（如 (2+1)D）中的应用和算符表示。

4.1.2 张量网络在格点规范理论中的应用

张量网络方法（特别是 MPS 和 MPO）是实现本研究计算的工具。早期和相关工作奠定了基础：

  • [23, 27] Steven R. White, “Density matrix formulation for quantum renormalization groups,” Phys. Rev. Lett. 69, 2863-2866 (1992) 和 Ulrich Schollwoeck, “The density-matrix renormaliza- tion group in the age of matrix product states,” An- nals Phys. 326, 96-192 (2011)。 DMRG 作为 MPS 方法的基础，是计算基态的核心算法。
  • [31] T. Byrnes, P. Sriganesh, R. J. Bursill, and C. J. Hamer, “Density matrix renormalization group approach to the massive Schwinger model,” Phys. Rev. D 66, 013002 (2002)。 这标志着首次将 DMRG 应用于格点规范理论（Schwinger 模型）。
  • [32-38] 一系列关于 TNs 在 LGTs 中应用的进展，包括 Abelian 和 non-Abelian 规范群，以及高维度拓展。

4.1.3 (1+1)D SU(2) LGT 的 TN 模拟

本研究直接拓展了先前在 (1+1)D SU(2) LGT 上的张量网络工作：

  • [106] Stefan Kühn, Erez Zohar, J. Ignacio Cirac, and Mari Carmen Bañuls, “Non-Abelian string breaking phenomena with Matrix Product States,” JHEP 07, 130 (2015)。 这是使用 MPS 研究 SU(2) 弦断裂的早期工作，但可能在系统尺寸和希尔伯特空间截断方面有所限制。
  • [108] Mari Carmen Bañuls, Krzysztof Cichy, J. Ignacio Cirac, Karl Jansen, and Stefan Kühn, “Efficient basis formu- lation for 1+1 dimensional SU(2) lattice gauge theory: Spectral calculations with matrix product states,” Phys. Rev. X 7, 041046 (2017)。 探讨了 SU(2) LGT 的高效基组表示。
  • [109] P. Sala, T. Shi, S. Kühn, M. C. Bañuls, E. Demler, and J. I. Cirac, “Variational study of U(1) and SU(2) lattice gauge theories with Gaussian states in 1+1 dimensions,” Phys. Rev. D 98, 034505 (2018)。 使用高斯态研究了 SU(2) LGT 的静态势能。

4.1.4 实时动力学与弦断裂

对实时动力学和弦断裂的研究一直是 TNs 的重要应用方向：

  • [39, 61, 81] 突出显示了 TNs 在实时模拟弦断裂动力学中的重要性。
  • [19, 20] Julian S. Schwinger, “On gauge invariance and vacuum polarization,” Phys. Rev. 82, 664-679 (1951)。 Schwinger 机制是粒子对产生的基础，与弦断裂紧密相关。
  • [126, 127, 128] 描述了 TDVP 及其 Krylov 基扩展协议，这是实时演化的核心算法。

4.1.5 ITensors.jl 库

  • [124] Matthew Fishman, Steven R. White, and E. Miles Stoudenmire, “The ITensor software library for tensor network calculations,” (2020)。 本研究所使用的具体张量网络库。

4.2 本工作局限性评论

尽管本研究在利用 LSH 形式化和张量网络模拟 (1+1)D SU(2) 格点规范理论的弦断裂方面取得了显著进展，但仍存在一些固有的局限性：

  1. 低维度限制： 本研究集中在 (1+1) 维系统。虽然 (1+1)D LGT 能捕捉非阿贝尔规范理论的核心特征，并提供 QCD 的宝贵物理洞察，但它与真实的 (3+1) 维 QCD 存在显著差异。许多重要的 QCD 现象，如三维空间中的色通量管几何结构、非平凡的拓扑效应等，无法在 (1+1) 维中完全体现。将这些方法扩展到更高维度（如 (2+1)D 或 (3+1)D）是未来的研究方向，但会引入更大的计算挑战（如 PEPS 等高维 TNs 的复杂性）。

  2. 动力学计算的极限外推缺失： 论文虽然在静态弦势的计算中实现了键维度、系统尺寸和格点间距的极限外推，以获得连续统和热力学极限下的弦张力，但在实时动力学方面并未进行这样的系统性外推。作者明确指出，动态可观测量的热力学和连续统极限外推在计算上非常耗时，并留作未来的工作。这意味着目前的动态结果是在固定格点间距和系统尺寸下获得的，可能仍受到有限尺寸效应和离散化误差的影响。

  3. 有限键维度对动力学的影响： 尽管 TDVP 算法允许动态扩展键维度 D_max，但在长时间演化中，纠缠熵会线性增长，最终可能超出 D_max 的承受能力。研究表明，对于轻质量费米子，计算结果在 gt=8 处开始积累约 10% 的相对误差，这限制了模拟的可靠时间尺度。虽然重质量结果收敛性更好，但有限 D_max 仍然是动态模拟的根本限制。

  4. 希尔伯特空间通量截断 j_max 的影响： 论文强调了 j_max 的选择对动力学结果的定性影响。当 j_max 取值较小时 (例如 j_max=1/2)，动力学行为与较大 j_max 值有显著差异。虽然论文声称 j_max=5/2 足以确保收敛，但本质上这仍是一个截断，可能会遗漏某些高能激发态对动力学的贡献。

  5. 渐进强子内容的识别： 本研究详细分析了 LSH 激发态的演化，包括粒子生成和输运。然而，它并未直接计算时间演化态与渐进强子态（如介子和重子）之间的重叠。因此，虽然微观过程被详细描述，但无法直接量化强子化率和最终产物的强子谱。这对于理解强子化模型的预测能力至关重要，但超出本工作的范围。

  6. LSH 形式中扩展弦的模糊性： 论文在附录 6 中指出，在 LSH 形式中，跨多个格点的扩展规范不变介子弦的解释并非唯一。虽然局域 LSH 激发态（环、弦、强子）是明确的，但如何将它们组合成宏观的“介子弦”可能存在多种方式，这给介子内容的识别带来了模糊性。

  7. 计算资源需求： 尽管张量网络方法相较于蒙特卡洛具有优势，但 TDVP 算法（特别是 2-site 版本）在计算上仍然非常昂贵。特别是在需要较大键维度和长时间演化的情况下，计算资源（内存和 CPU 时间）仍然是一个重要的制约因素。

  8. 缺乏与其他方法的直接比较： 虽然论文与Ref. [109] 中的静态弦张力进行了比较，但对于实时弦断裂动力学，本文的设置是独特的，因此缺乏与其他现有方法的直接定量比较，无法全面评估其在特定动力学方面的优势和局限性。

尽管存在这些局限性，本研究通过引入 LSH 形式化，为非阿贝尔 LGT 的张量网络模拟提供了强大的新范式，并对 (1+1)D SU(2) LGT 中的弦断裂动力学提供了前所未有的微观洞察。这些局限性也为未来的研究指明了方向。

5. 其他必要的补充

5.1 对 QCD 强子化模型的深远影响与未来方向

本研究工作不仅在格点规范理论的量子模拟领域取得了重要的技术突破，更重要的是，它为理解粒子物理中 QCD 的强子化过程提供了新的视角和工具，并为未来的研究开辟了广阔前景：

  • 理解强子化模型： 弦断裂和强子化是描述夸克-胶子等离子体冷却后如何形成强子的核心物理过程。现有现象学模型（如 Lund 模型、Cluster 模型）在预测强子产额方面取得了成功，但其参数并非直接来源于 QCD 第一性原理。本研究通过在简化 LGT 模型中对弦断裂动力学的实时模拟，为这些模型提供了微观层面的量子力学基础。未来，通过对弦断裂动力学中粒子生成、散射和输运过程的深入量化分析，有望为强子化模型的构建提供更坚实的第一性原理输入，甚至可以用于测试现有模型中的假设和局限性。

  • 拓展至高维 LGTs： LSH 形式化的一大优势在于其固有的局域规范不变性，这使得它易于推广到更高维度的格点规范理论。本研究为在 (2+1)D 甚至 (3+1)D 系统中构建基于 LSH 的张量网络（如 Projected Entangled-Pair States, PEPS 或 Tree Tensor Networks, TTNs）奠定了基础。这将是未来模拟更接近真实 QCD 系统的关键一步，有望揭示高维系统中弦的几何形状、拓扑效应和更复杂的强子化模式。

  • 量子模拟器与量子计算的基准测试： 本研究中发展的张量网络工具和方法可作为未来量子模拟实验和量子计算实现的基准。通过在超冷原子、离子阱等量子平台上实现 LGT 哈密顿量并模拟弦断裂过程，可以验证量子模拟器的能力，并推动量子计算在解决粒子物理难题上的应用。

  • 更精细的动力学可观测量： 除了本研究中使用的 LSH 激发态期望值、能量密度、纠缠熵和两点关联函数外，未来可以计算更多精细的动力学可观测量，如电流-电流关联函数，以揭示介质的输运性质。通过计算时间演化态与渐进强子态的重叠，可以更直接地量化强子化率和粒子产额。这些量化的研究将进一步加深对弦断裂物理的理解。

  • 系统化散射模拟： 本研究揭示了弦断裂动力学中复杂的链散射事件，包括弹性与非弹性散射。未来可以进行受控的散射模拟，以更好地表征散射结果及其概率，并研究其与初始条件的关系。这对于理解高能碰撞中粒子相互作用的最终态至关重要。

5.2 LGT 方法在量子化学领域的潜在启示

尽管格点规范理论主要应用于粒子物理，但本研究中使用的张量网络方法（特别是 DMRG 和 TDVP）最初来源于量子化学和凝聚态物理领域。LGT 领域的进展可以为量子化学提供一些潜在的启示和新的研究思路：

  • 强关联电子系统： 量子化学在处理强关联电子系统（例如过渡金属化合物、重元素体系）时面临巨大挑战。LGT 强调了局域规范不变性，这在某种程度上与量子化学中局域轨道的占据和相互作用有关。LSH 形式化将规范不变性内化到局域自由度中的思想，或许可以启发量子化学中对强关联电子系统构建更有效的局域基组或守恒量驱动的计算方法。

  • 精确激发态和动力学： 本研究对实时弦断裂动力学的精确模拟，特别是对能量传输、纠缠熵生成和关联扩散的分析，与量子化学中对分子激发态、非绝热动力学、能量转移和光化学反应的研究高度相关。LSH 方法揭示了复杂的微观过程如何导致宏观动力学，这可以为量子化学中理解复杂化学反应机理提供新的建模框架和诊断工具。

  • 张量网络方法的新发展： LGT 领域对张量网络，尤其是高维张量网络（如 PEPS、TTNs）和动态张量网络算法的不断探索，可能会反过来推动量子化学中张量网络方法的进步。例如，处理 LGT 中快速增长的纠缠熵的策略，或将对称性更高效地整合到张量网络中的技术，都可能被借鉴到量子化学中。

  • 规范对称性与分子对称性： LGT 中对规范对称性的严格处理与量子化学中对分子对称性（如点群对称性、自旋对称性）的利用有异曲同工之处。LSH 形式化通过重构自由度来满足规范定律，这种思想或许可以启发量子化学中构建对称性适应的量子态或算符，从而简化计算并提高精度。

  • 量子计算的交叉应用： LGT 正在成为量子计算的重要应用场景。量子化学同样是量子计算的潜在杀手级应用。LGT 中开发出的量子算法（如用于基态搜索、动态模拟）以及对量子硬件误差的鲁棒性研究，可以直接或间接为量子化学领域的量子算法开发提供经验和工具。

  • 场论与分子物理的对话： 粒子物理中的场论方法，特别是格点场论，提供了描述基本相互作用的通用框架。量子化学可以从这种场论视角中获得启发，例如在描述键断裂/形成或分子间相互作用时，可以尝试引入一些场论中的概念，如“弦”或“通量”，来提供更深层次的物理理解。当然，这需要细致的类比和映射，因为 LGT 关注的是基本粒子和基本力，而量子化学关注的是原子核和电子的复杂多体系统。

总而言之，本研究不仅是对 (1+1)D SU(2) LGT 弦断裂动力学的全面研究，更是通过其创新的 LSH 形式化和先进的张量网络方法，为粒子物理和量子化学这两个看似不同却又共享底层量子物理工具的领域，开辟了相互借鉴和共同发展的广阔空间。