非阿贝尔规范场与轴子场耦合的动力学研究
1. 非阿贝尔规范场与轴子场耦合的物理背景
在粒子物理和宇宙学研究中，轴子场作为一种赝标量场，与非阿贝尔规范场的耦合动力学一直是一个引人入胜的课题。这种耦合不仅存在于QCD轴子模型中，也出现在许多宇宙学场景中，特别是所谓的"自然膨胀"模型中。这类模型利用轴子场的平移对称性来保证膨胀势的平坦性，从而为解决标准暴胀模型中的精细调节问题提供了可能途径。
轴子场与规范场的耦合通常通过拓扑荷密度项实现，形式为φFF̃，其中φ是轴子场，F是规范场强张量。在阿贝尔规范场（如电磁场）情况下，这种耦合会导致 tachyonic不稳定性，产生一系列可观测的宇宙学信号。然而，非阿贝尔规范场的情况更为复杂，因为规范协变性会引入动量依赖的三点自相互作用，这些相互作用会驱动系统快速趋向平衡。
2. 研究模型与理论框架
2.1 基本拉格朗日量
我们考虑的系统由一个赝标量场φ∈R和一个非阿贝尔规范场组成。假设φ对于规范群是中性场。系统的拉格朗日量密度可写为：
L = 1/2 ∂μφ∂μφ - V(φ) - 1/4 FμνFμν - φχ/fa
其中χ是拓扑荷密度，定义为： χ ≡ εμνρσ g² FμνFρσ / (64π²)
这里我们采用(+---)度规签名；Fμν是Yang-Mills场强，a=1,...,N²c-1标记生成元；ε0123=1；g²=4πα是规范耦合常数。为减少计算成本，本研究主要关注Nc=2的情况，尽管对于QCD轴子的暖暴胀研究，Nc=3才是相关情形。
2.2 运动方程与约束条件
从上述拉格朗日量导出的运动方程在部分规范固定Aa0=0后，空间规范势Aai作为正则变量，Aa0的欧拉-拉格朗日方程给出需要在初始时刻满足的高斯约束：
D·E = κ(∇φ)·B
其中E和B既是规范代数空间也是物理空间的矢量，D是伴随表示中的协变导数。E和B由规范势定义导致两个Jacobi恒等式：
D·B = 0 Ḃ = D×E
最后，φ和Aai的欧拉-拉格朗日方程与Eai的定义一起给出了动力学变量的演化方程：
φ̈ = ∇²φ - V'(φ) - κE·B Ȧ = E Ė = -D×B + κ(φ̇B - ∇φ×E)
3. 数值实现方法
3.1 空间离散化
为了数值研究动力学，我们用格点间距a>0对空间方向进行离散化。在此过程中，保持规范协变性至关重要。我们采用标准格点规范理论的策略，将空间规范势Aai替换为幺正链接矩阵Ui∈SU(2)作为动力学变量。
3.2 时间离散化与有限体积效应
在模拟中，方程需要通过两种方式进一步近似：离散化时间坐标和将系统置于有限体积中。对于时间坐标，我们选择步长δt̃≪ã，并用三阶Runge-Kutta算法求解方程。这个算法不是辛算法，意味着能量不能精确守恒。然而，通过限制δt̃≤0.15ã，我们发现对于t̃≤50且κ̃≈2.0的情况，能量守恒的相对精度可达10^-4，这对我们的目的已经足够。
对于空间体积，我们采用立方盒子，所有方向采用周期性边界条件。测试已在N=256,512,1024下进行，验证了体积依赖性的不显著性。我们的SU(2)生产运行使用固定体积N³=756³。
3.3 初始条件与功率谱
我们保持凝聚态的初始条件固定，确保它携带大部分初始能量密度，并测试了多种为凝聚态添加噪声的方法。在默认设置中，φ̃的噪声在连续情况下由方程给出，而电场根据方程初始化；磁场设为零以满足高斯约束。在替代方案中，我们只在磁场中放入噪声（链接已从规范势的指数映射获得），设置φ̃和Ẽai为零。我们也测试了同时初始化电场和磁场的情况，尽管这样不能满足高斯约束。
一旦演化开始，我们就根据方程测量等时功率谱作为时间的函数。在格点单位中，功率谱由方程给出。
4. 平衡化早期阶段的动力学
4.1 规范场能量密度的指数增长
线性扰动理论表明，如果κ̃φ̇̃≠0，那么波数k̃<κ̃|φ̇̃|max的红外模式可能经历指数增长。我们的数值模拟结果验证了这一预期。图2展示了几个κ̃值下，电场和磁场的能量密度随时间t̃的演化。灰色带表示指数增长区域。在灰色带之间，增长会暂停一段时间；这对应于|φ̇̃|变小，导致不稳定带变窄。经过几个不稳定带后，规范场能量密度达到饱和；对于小κ̃，这需要更长时间，因为不稳定带更窄。
在灰色带中，规范场能量密度的相对增长率φ̇gauge/egauge近似为常数，我们将行为拟合为：
egauge(t̃) ≈ egauge,0 eΓ̃Δt̃, Δt̃≡t̃-t̃0
4.2 轴子振荡的阻尼
对于轴子场，我们从测量数据中通过空间平均提取轴子凝聚态：
φ̃(t̃) ≡ 1/N³ Σx φ̃(t̃,x)
它的演化如图3所示，揭示了对于κ̃≪1.0的长周期振荡和对于κ̃≫1.0的过阻尼状态。回顾图2，阻尼是由于能量从凝聚态转移到规范自由度和凝聚态碎裂造成的。
仔细观察发现，阻尼模式是非平凡的，并且只在一段时间后才开始。在SU(2)情况下，我们发现数据可以很好地表示为：
φ̃(t̃) ≈ { A0 sin(ω̃1t̃), t̃<t̃0 { sin(ω̃2t̃-φ2){exp -Υ̃(t̃-t̃0) +A}, t̃≥t̃0
4.3 功率谱与能量均分
类似于方程，我们定义电场和轴子功率谱为：
P̃E ≡ dA k̃³/(2π²) (|Ẽk,+(t̃)|² + |Ẽk,-(t̃)|²) P̃φ ≡ k̃³/(2π²) |φ̃k(t̃)|²
其中dA=3(1)表示规范代数的生成元数量，Ẽk,±和φ̃k是模式函数。图6展示了κ̃=1.0时从初始时间t̃=0到最终时间t̃=50的时间演化。深蓝色曲线对应早期时间，深红色曲线对应晚期时间。初始谱的尖锐截断是由于方程中施加的动量截断Λ。
功率谱在不稳定带k̃<κ̃|φ̇̃|max内显示初始增长，用虚线垂直线表示。随后能量被转移到更高动量的模式。这个过程在SU(2)中比在U(1)中明显更平滑，后者显示出类似于图1的振荡特征。我们将SU(2)功率谱的平滑性归因于该情况下存在的规范自相互作用。
5. 结论与展望
本文的目的是开始数值研究轴子型暴胀子场(φ)和非阿贝尔规范系综如何相互作用。为此，我们改变了表征相互作用强度的无量纲耦合κ̃=αm/(2πfa)，并在756³的格点上模拟了系统。大部分初始能量密度由空间均匀的φ凝聚态携带。添加少量噪声会开启到规范系综的衰变通道，凝聚态的能量密度迅速转移到其他自由度（规范涨落和φ的空间变化）。
我们观察到的SU(2)动力学显示出平滑的参数依赖性，这是线性响应区域的特征。电场和磁场的能量密度随时间呈指数增长，增长率近似线性正比于κ̃。轴子振荡受到阻尼，在SU(2)情况下，阻尼率近似正比于κ̃²。然而，轴子阻尼只有在规范系综获得足够能量密度后才开始，这个延迟因子与κ̃成反比。最终，规范区域的总能量密度大约是轴子场的2(N²c-1)倍，表现为能量均分。
从物理上讲，将我们观察到的能量均分与平衡化（或热化）联系起来是很诱人的。然而，严格来说，后者不能用经典格点规范理论研究，因为热动量尺度k∼T已被非物理的格点截断尺度k∼π/a所取代。然而，考虑到总能量守恒且哈勃率为零，我们设置的物理类比对应于瞬时再加热，其中所有凝聚态能量密度最终都转移到均分化的UV自由度。
转向与U(1)情况的比较，我们注意到尽管它与SU(2)有许多相似之处，但我们也展示了非常显著的差异。在U(1)中，对κ的依赖并不平滑，但即使很小的非零值也会导致快速均分。我们怀疑这是由于在没有规范场自相互作用的情况下，tachyonic不稳定性没有太大衰减。
在本研究中，我们保持非阿贝尔规范耦合固定，g²=1.0。原则上，在保持κ固定的同时改变规范自相互作用的强度是很有趣的。由于经典理论的尺度不变性，这并不完全简单。然而，探索如何以物理上有意义的方式实现这种变化将是很有趣的。
第二个重要的下一步是在模拟中包含宇宙的膨胀。膨胀选项已内置在CosmoLattice平台中，因此我们希望在不久的将来回到这个话题。
最后，包含费米子代表了另一个有趣的方向。即使基本费米子不允许经典描述，费米子系综的平均确实可以通过时间依赖的化学势来包含，这可以偏向sphaleron过程。我们将这个扩展留待未来的研究。