[论文解读] Numerical methods for the sign problem in Lattice Field Theory

Lorenzo Bongiovanni|arXiv (Cornell University)|Mar 10, 2016

Theoretical and Computational Physics参考文献 16被引用 3

一句话总结

本博士论文研究了三种主流方法以克服格点场论中的符号问题：带规范冷却的复 Langevin 动力学、Lefschetz 柱面以及基于 LLR 算法的态密度方法。论文展示了带规范冷却的复 Langevin 动力学在 SU(3) Yang-Mills 理论中含 θ 项时的可行性，并在模型上将其与柱面方法进行比较；同时将态密度方法应用于有限化学势下的相对论性玻色气体，结果与解析解高度一致。

ABSTRACT

The great majority of algorithms employed in the study of lattice field theory are based on Monte Carlo's importance sampling method, i.e. on probability interpretation of the Boltzmann weight. Unfortunately in many theories of interest one cannot associated a real and positive weight to every configuration, that is because their action is explicitly complex or because the weight is multiplied by some non positive term. In this cases one says that the theory on the lattice is affected by the sign problem. An outstanding example of sign problem preventing a quantum field theory to be studied, is QCD at finite chemical potential. Whenever the sign problem is present, standard Monte Carlo methods are problematic to apply and, in general, new approaches are needed to explore the phase diagram of the complex theory. Here we will review three of the main candidate methods to deal with the sign problem, namely complex Langevin dynamics, Lefschetz thimbles and density of states method. We will first study complex Langevin dynamics, combined with the gauge cooling method, on the one-dimensional Polyakov line model, and then we will apply it to pure gauge Yang-Mills theory with a topological theta-term. It follows a comparison between complex Langevin dynamics and the Lefschetz thimbles method on three toy models, which are the quartic model, the U(1) one-link model with a mu dependent determinant, and the SU(2) non abelian one-link model with complex beta parameter. Lastly, we introduce the density of state method, based on the LLR algorithm, and we will employ it in the study of the relativistic Bose gas at finite chemical potential.

研究动机与目标

• 解决格点场论中的符号问题，该问题导致标准蒙特卡罗方法在作用量为复数或权重非正定情况下无法应用。
• 评估并比较三种有前景的数值方法——复 Langevin 动力学、Lefschetz 柱面及态密度方法——在处理具有符号问题的理论时的表现。
• 测试带规范冷却的复 Langevin 动力学在非阿贝尔规范理论中的鲁棒性与准确性，特别是 SU(3) Yang-Mills 理论中含 θ 项的情形。
• 将基于 LLR 算法的态密度方法应用于有限化学势下的相对论性玻色气体，该系统具有严重的符号问题。
• 对复 Langevin 动力学与 Lefschetz 柱面方法在基准模型上的系统性比较，包括四次模型、U(1) 单链接模型及 SU(2) 单链接模型。

提出的方法

• 通过将场变量复化，并根据 Fokker-Planck 方程随机演化，实现复 Langevin 动力学，其中在非阿贝尔理论中使用规范冷却以稳定动力学。
• 通过在复化流形中寻找作用量的临界点来构建 Lefschetz 柱面，并沿虚部作用量恒定的流形积分，确保测度为实数且正定。
• 态密度方法利用 LLR 算法从一系列受限模拟中计算态密度，从而实现从高度振荡积分中重构配分函数。
• 对于相对论性玻色气体，通过拟合配分函数对控制参数的对数导数，利用 LLR 算法计算自由能差。
• 在复 Langevin 模拟中，通过梯度流监测拓扑荷，以提取拓扑易失性并验证收敛性。
• 通过分析场构型分布及各柱面对路径积分的贡献，实现复 Langevin 与柱面方法的比较。

实验结果

研究问题

• RQ1 带规范冷却的复 Langevin 动力学能否在 SU(3) Yang-Mills 理论中含 θ 项（符号问题严重）时产生可靠结果？
• RQ2 复 Langevin 动力学的结果与 Lefschetz 柱面方法在四次模型、μ 依赖性行列式存在的 U(1) 单链接模型及 SU(2) 单链接模型上的结果相比如何？
• RQ3 基于 LLR 算法的态密度方法在有限化学势下的相对论性玻色气体中，对配分函数与自由能的计算精度如何？
• RQ4 规范冷却在非阿贝尔规范理论中稳定复 Langevin 动力学方面起什么作用？它如何影响拓扑荷分布的收敛性？
• RQ5 不同 Lefschetz 柱面的贡献如何影响路径积分？在复作用量存在时，哪些固定点主导测度？

主要发现

• 带规范冷却的复 Langevin 动力学在 SU(3) Yang-Mills 理论中含 θ 项时，成功再现了预期的拓扑荷分布，且在虚 θ 情况下与 HMC 结果一致。
• 实 θ 与虚 θ 值下的拓扑荷分布表现出相反的曲率，证实了预期的 CP 奇性行为，并验证了该方法的一致性。
• 在四次模型中，复 Langevin 动力学正确捕捉了场构型的分布，且稳定柱面与模拟中观察到的主导贡献相匹配。
• 在 U(1) 单链接模型中，该方法识别出正确的贡献柱面，且复 Langevin 结果与基于柱面的预测一致，尤其在存在 μ 依赖性行列式时表现良好。
• 在复 β 的 SU(2) 单链接模型中，带规范冷却的复 Langevin 模拟收敛至正确分布，柱面分析证实单一稳定柱面占主导。
• 基于 LLR 算法的态密度方法在相对论性玻色气体中准确计算了自由能差，结果在统计误差范围内与解析的平均场理论预测一致。