[论文解读] Adaptive gauge cooling for complex Langevin dynamics

Lorenzo Bongiovanni, Gert Aarts+3 more|arXiv (Cornell University)|Nov 5, 2013

High-Energy Particle Collisions Research参考文献 7被引用 9

一句话总结

本文提出自适应规范冷却（adaptive gauge cooling）以稳定非阿贝尔规范理论中具有复数作用量的复数朗之万动力学，例如带有 θ 项的 SU(3) 杨-米尔斯理论。通过动态调整规范冷却参数，该方法能有效最小化复化 SL(N,ℂ) 配置空间中与 SU(N) 的距离，从而实现对真实 θ 值的可靠模拟——通过与 HMC 在虚 θ 情况下的结果一致，以及在 6⁴ 网格上首次直接获得真实 θ 的结果得到验证。

ABSTRACT

In the case of nonabelian gauge theories with a complex weight, a controlled exploration of the complexified configuration space during a complex Langevin process requires the use of SL(N,C) gauge cooling, in order to minimize the distance from SU(N). Here we show that adaptive gauge cooling can lead to an efficient implementation of this idea. First results for SU(3) Yang-Mills theory in the presence of a nonzero theta-term are presented as well.

研究动机与目标

• 为解决具有复数作用量的非阿贝尔规范理论中的符号问题，例如在非零 θ 或化学势下的 QCD。
• 通过控制复化 SL(N,ℂ) 配置空间的探索，改善复数朗之万动力学的收敛性和稳定性。
• 开发并测试一种自适应规范冷却算法，动态调整冷却参数以保持与 SU(N) 的接近性，避免虚假结果。
• 实现 SU(3) 杨-米尔斯理论中真实 θ 项的可靠模拟，该情况此前因符号问题而无法实现。

提出的方法

• 自适应规范冷却在规范变换 Ω_x = exp(−ε α_gf f_ax λ_a) 中使用动态调整的参数 α_gf，以最小化酉性范数 Tr(UU†)/N。
• 在每次朗之万更新后迭代应用冷却步骤，减少与 SU(N) 的距离，同时不破坏规范不变性。
• 通过利用酉性范数在 SL(N,ℂ) 变换下不不变的特性，避免了对规范固定或 Fadeev-Popov 衍生行列式的需求。
• 该算法应用于 SU(2) 和 SU(3) 格点规范理论，通过与 HMC 对虚 θ 的结果对比，验证了平面元和拓扑荷期望值的一致性。
• 冷却参数基于酉性范数的局部梯度进行自适应调整，以优化收敛速度和稳定性。
• 在 6⁴ 网格上对 SU(3) 杨-米尔斯理论进行模拟，作用量为 S = S_W − iθ_L ∑x q_L(x)，使用实数和虚数 θ_L。

实验结果

研究问题

• RQ1 自适应规范冷却是否能稳定具有复数作用量的非阿贝尔规范理论中的复数朗之万动力学？
• RQ2 自适应冷却是否优于固定冷却，能更有效地保持配置空间接近 SU(N)，并确保收敛到正确的物理结果？
• RQ3 在标准重要性采样失效的情况下，复数朗之万动力学结合自适应规范冷却能否对带有真实 θ 项的 SU(3) 杨-米尔斯理论产生可靠结果？
• RQ4 在复数朗之万动力学中，结合自适应冷却时，真实和虚数 θ 下的作用量和拓扑荷分布如何表现？
• RQ5 在虚数 θ 情况下，复数朗之万结果与 HMC 是否存在定量一致，从而验证该方法？

主要发现

• 自适应规范冷却成功稳定了带有 θ 项的 SU(3) 杨-米尔斯理论中的复数朗之万动力学，首次实现了对真实 θ 值的模拟。
• 对于虚数 θ_L，复数朗之万方法在 6⁴ 网格上对 β = 5.9、6.0 和 6.1 的平面元期望值结果与 HMC 结果一致。
• 真实 θ_L 下拓扑荷虚部的运行平均值表现出与 θ_L 的线性依赖关系，与弱耦合区间的预期 O(θ_L) 行为一致。
• 对于真实 θ_L = 2，作用量分布显示在较小 β 时尾部更宽，表明冷却效果减弱，需进一步分析。
• 真实和虚数 θ_L 下的拓扑荷涨落存在显著差异，表明存在尚未理解的非平凡动力学。
• 当配置在规范等价于 SU(N) 时，该方法在单链模型中表现出向 SU(N) 流形的指数收敛，证实了其理论基础。