1. Z2格点规范理论中的强子动力学研究概述在量子多体物理研究中格点规范理论作为连接高能物理与凝聚态物理的重要桥梁为我们理解强相互作用系统提供了独特视角。Z2格点规范理论因其离散对称性和相对简单的数学结构成为研究规范场论非微扰特性的理想模型。这项研究聚焦于该理论框架下的强子动力学行为特别是3介子态和四夸克态之间的相互作用与转化机制。传统观点认为强子间的束缚主要源于吸引力作用但我们的研究表明在某些参数区间排斥相互作用同样能导致稳定的多粒子束缚态形成。这一发现挑战了常规认知为理解复杂量子系统中的多体效应开辟了新思路。通过精确的数值模拟我们系统考察了不同耦合强度下强子态的演化特性揭示了排斥束缚态的动力学稳定性及其与系统参数的依赖关系。研究采用了基于矩阵乘积态(MPS)的时间演化块解耦(TEBD)算法这是一种针对一维量子多体系统的高效数值方法。我们特别关注了初始态分别为3介子态和四夸克态两种情况下的时间演化行为通过分析粒子数密度分布、束缚态占据数等物理量建立了完整的强子动力学图像。数值结果得到了有效模型分析和精确对角化计算的验证形成了从微观机制到宏观现象的多层次理解。2. 理论基础与模型构建2.1 Z2格点规范理论的基本框架Z2格点规范理论是在离散时空格点上构建的规范场论模型其规范群为最简单的离散群Z2。该模型虽然数学结构相对简单却包含了规范理论的核心特征如规范不变性、禁闭现象等。在本文研究的自旋链实现中物质场由位于格点上的自旋1/2粒子表示规范场则体现在链接变量上。系统的完整哈密顿量包含三个主要项H -J∑i,jσ^z_iσ^z_j - K∑iσ^x_i - h∑iσ^z_i - m∑ini其中第一项表示最近邻自旋相互作用(J项)第二项是粒子对产生湮灭项(K项)后两项分别代表规范场(h项)和物质场(m项)的能量贡献。这个哈密顿量在Z2规范变换下保持不变体现了系统的规范对称性。特别值得注意的是K项的作用它允许局域粒子对的产生和湮灭是连接不同粒子数态的关键。当hm时系统处于共振条件此时3介子态和四夸克态能量简并K项将在它们之间产生有效耦合。这一特性为研究多粒子束缚态的形成提供了理想平台。2.2 强子态的表示与分类在我们的研究中强子态被定义为系统中局域的能量激发。根据包含的粒子数不同主要分为以下几类1介子态单个粒子-反粒子对对应自旋链中的↑↓或↓↑构型3介子态三个相邻的粒子-反粒子对构型为↑↓↑↓↑↓四夸克态四个粒子组成的束缚态构型为↑↑↓↓↑↑这些态在自旋链中的空间分布如图1所示。值得注意的是四夸克态并非简单的两个独立介子而是通过规范场相互作用形成的关联态。在参数适当时这种关联可以导致稳定的多粒子束缚态即使存在净排斥相互作用。关键点强子态的识别不能仅基于粒子数还需考虑空间关联特性。实际分析中我们通过计算两点关联函数来确认束缚态的形成。3. 数值方法与计算细节3.1 矩阵乘积态与TEBD算法矩阵乘积态(MPS)是一类针对一维量子多体系统的高效表示方法其核心思想是将多体波函数分解为一系列局部张量的乘积。对于系统尺寸为L的链MPS表示为|ψ⟩ ∑{s} A^{s1}A^{s2}...A^{sL}|s1s2...sL⟩其中A^{si}是位置i处的秩-3张量si表示该位置的物理自由度。MPS的表达能力由键维数χ决定χ越大表示能承载的纠缠越多但计算成本也越高。时间演化块解耦(TEBD)算法是模拟MPS时间演化的有效方法。其基本步骤包括将哈密顿量分解为奇数位和偶数位两部分对每个小时间步长δt应用二阶Suzuki-Trotter分解对每个两体门作用后的MPS进行奇异值分解并截断我们采用的参数为Jδt0.025最大键维数χ32。如图S1所示这个参数设置已经能保证结果的收敛性。特别对于研究的强子动力学由于纠缠增长相对缓慢较小的χ即可准确捕捉系统行为。3.2 初始态制备与背景扣除研究选取了两种典型初始态中心局域的3介子态|↑↓↑↓↑↓⟩位于链中心中心局域的四夸克态|↑↑↓↓⟩位于链中心这些初态都是乘积态可用键维数χ1的MPS精确表示。实际操作中我们发现真空背景涨落会引入高频噪声。这是由于真实的Z2规范理论真空态并非简单的无物质态而是包含由K项引起的粒子数涨落。为消除这一影响我们采用了背景扣除技术单独演化真空初态(全↓态)记录其粒子数密度分布从主计算结果中减去背景演化结果图S3展示了背景扣除前后的对比可见该方法有效消除了均匀分布的背景涨落使束缚态的动力学特征更加清晰。3.3 收敛性验证与误差控制为确保数值结果的可靠性我们系统考察了各参数的收敛性键维数χ的影响如图S1所示当χ≥32时四夸克数的时间演化曲线已基本收敛。这表明系统的纠缠熵增长有限中等大小的χ已足够。Trotter步长δt的影响图S2显示当Jδt≤0.025时时间离散化误差可忽略。我们因此固定Jδt0.025进行所有计算。系统尺寸效应通过比较不同链长(L50-100)的结果确认边界效应不影响中心区域的动力学。此外我们还计算了累积平均的收敛性(图S4)确保观测到的束缚态特征不是暂态现象。所有检查共同保证了研究结论的可靠性。4. 强子动力学的主要结果4.1 3介子初态的演化特性当系统初始制备为中心3介子态时时间演化展现出丰富的动力学行为。根据K/J的不同可观察到三种典型区域大K/J区域(K/J0.5)清晰的3介子-四夸克振荡(图2)长时平均四夸克数趋于稳定值(图3)形成分离的吸引和排斥束缚态中等K/J区域(0.2K/J0.5)振荡幅度减小但仍可辨识部分态与连续谱杂化排斥束缚态仍占主导小K/J区域(K/J0.2)快速弛豫到准稳态仍保持有限的四夸克数纯排斥束缚机制图S7中的密度分布清晰展示了这些特征。特别值得注意的是即使在小K/J区域系统仍保持局域的四夸克成分这明确表明了排斥束缚态的存在。4.2 四夸克初态的互补研究为验证结果的普适性我们也考察了四夸克初态的演化行为(图S7)。主要发现包括大K/J区域等幅的3介子-四夸克振荡长时平均各占50%(图S8)与3介子初态对称小K/J区域部分解离为1介子连续谱但仍保持有限束缚概率杂化程度高于3介子初态这种不对称性源于排斥束缚态中3介子和四夸克成分的不等权重。当K→0时四夸克态与连续谱的耦合强度(J²/2h)与带宽(4J²/h)相当导致不完全解离表现出准局域化行为。4.3 共振条件的影响研究还考察了偏离共振条件(h≠m)的情况(图S9)。当|h-m|增大时3介子-四夸克共振耦合减弱束缚态特征逐渐消失3介子表现出禁闭动力学四夸克态解离为1介子图S10定量展示了长时平均四夸克数随|h-m|的变化在共振点(hm)出现明显峰值。这一结果证实了共振条件对形成四粒子束缚态的关键作用。5. 有效模型与机制分析5.1 低能有效哈密顿量为深入理解数值结果我们通过二阶微扰理论导出了有效模型。在hm的共振条件下主要考虑三个子空间所有可能的3介子态{|m₃^c⟩}所有可能的四夸克态{|q₄^c⟩}分离的1介子态{|m₁^c,r⟩}有效哈密顿量矩阵元包括一阶项K连接|m₃⟩和|q₄⟩二阶跃迁-J²/2h连接|q₄⟩和|m₁^{r2}⟩能移项来自虚跃迁过程完整的有效模型如图1(b)所示其能谱与精确对角化结果高度一致(图S12 vs 图S13)。5.2 束缚态的形成机制基于有效模型可以明确区分两种束缚机制大K/J下的常规束缚态由K项主导形成远离连续谱的能级态为|m₃⟩±|q₄⟩叠加全K范围的排斥束缚态源于能移与杂化的平衡位于连续谱上方以|m₃⟩成分为主图S12的能谱分析清晰展示了这些特征。随着K/J减小吸引束缚态逐渐消失但排斥束缚态始终保持。5.3 精确对角化验证为验证有效模型的可靠性我们对L12系统进行了精确对角化(图S11)。主要发现大K/J时分离的吸引和排斥束缚带与连续谱有明显能隙小K/J时吸引束缚态融入连续谱排斥束缚态保持分离这些结果与TEBD模拟完全一致共同支持了排斥束缚态的稳定性。有效模型仅在高激发态附近出现微小偏差但对研究的物理现象捕捉准确。6. 讨论与展望本研究系统考察了Z2格点规范理论中多粒子束缚态的动力学行为揭示了排斥相互作用导致束缚的新机制。这一发现不仅丰富了我们对强相互作用系统的认识也为量子模拟实验提供了明确的理论指导。在实际材料体系中类似的排斥束缚机制可能存在于以下场景强关联电子系统中的电荷密度波超冷原子系统中的束缚分子态量子磁性材料中的自旋激发未来的研究方向包括推广到更高维格点系统考虑非阿贝尔规范群情况研究有限温度效应探索与拓扑序的联系特别值得关注的是近期超冷原子实验技术的进步使得Z2规范理论的量子模拟成为可能。我们的结果为这类实验提供了明确的理论预测和可观测信号将有力推动格点规范理论的实验研究。