[论文解读] Complex Langevin for Lattice QCD at $T=0$ and $μ\ge 0$

D. K. Sinclair, J. B. Kogut|arXiv (Cornell University)|Nov 7, 2016

Quantum Chromodynamics and Particle Interactions参考文献 11被引用 6

一句话总结

该论文将复 Langevin 方程（CLE）结合规范冷却方法应用于零温、有限化学势 μ 的两味格点 QCD 模拟，测试其再现正确相结构与物理量的能力。在 β=5.6 的 12⁴ 格点上，CLE 在 μ≈mN/3 附近显示出相变，但对有效场强（如手征子凝聚与平面子）存在系统性偏差；然而，在耦合更弱的 β=5.7 的更大 16⁴ 格点上，与已知结果的符合度显著提高，表明 CLE 在弱耦合、连续极限下可能给出正确结果。

ABSTRACT

QCD at finite quark-/baryon-number density, which describes nuclear matter, has a sign problem which prevents direct application of standard simulation methods based on importance sampling. When such finite density is implemented by the introduction of a quark-number chemical potential $μ$, this manifests itself as a complex fermion determinant. We apply simulations using the Complex Langevin Equation (CLE) which can be applied in such cases. However, this is not guaranteed to give correct results, so that extensive tests are required. In addition, gauge cooling is required to prevent runaway behaviour. We test these methods on 2-flavour lattice QCD at zero temperature on a small ($12^4$) lattice at an intermediate coupling $β=6/g^2=5.6$ and relatively small quark mass $m=0.025$, over a range of $μ$ values from $0$ to saturation. While this appears to show the correct phase structure with a phase transition at $μ\approx m_N/3$ and a saturation density of $3$ at large $μ$, the observables show departures from known values at small $μ$. We are now running on a larger lattice ($16^4$) at weaker coupling $β=5.7$. At $μ=0$ this significantly improves agreement between measured observables and known values, and there is some indication that this continues to small $μ$s. This leads one to hope that the CLE might produce correct results in the weak-coupling -- continuum -- limit.

研究动机与目标

• 测试带规范冷却的复 Langevin 方程（CLE）在有限重子密度下模拟格点 QCD 的有效性，标准方法因符号问题而失效。
• 确定 CLE 是否能正确再现 QCD 的相结构，包括 μ≈mN/3 处的强子物质到核物质的相变及大 μ 下的饱和行为。
• 评估小格点上 CLE 物理量（尤其是手征子凝聚与平面子）的系统性误差，并评估在更大格点与更弱耦合下的改进效果。
• 研究尽管费米子行列式存在极点导致漂移项的全纯性被破坏，CLE 是否能在弱耦合与连续极限下收敛至正确分布。
• 探索 CLE 在研究高 μ 相（如颜色超导性）及有限温度下的临界点方面的潜力。

提出的方法

• 通过将 SU(3) 规范场推广至 SL(3,C)，将标准 Langevin 方程推广至复数场，从而在有限 μ 下实现具有复费米子行列式的 QCD 模拟。
• 实施规范冷却，通过在演化过程中将配置投影回 SU(3) 流形，抑制发散行为，防止场强无界增长。
• 使用复作用 S(U,μ) = β∑□{1−1/6Tr[UUUU+(UUUU)⁻¹]} − Nf/4 Tr{ln[M(U,μ)]} 描述系统，其中 M(U,μ) 为重叠 Dirac 算符。
• 通过带复噪声项 ηₗ 的随机微分方程演化规范场，确保在 Langevin 时间中满足正确归一化与高斯白噪声条件。
• 在 μ=0 与 μ>0 条件下对 12⁴ 与 16⁴ 格点进行模拟，将物理量（如平面子、手征子凝聚与夸克数密度）与 RHMC 的参考结果进行比较。
• 系统性地改变 β 与格点尺寸，评估有限尺寸与有限时间步长效应，并趋近弱耦合与连续极限。

实验结果

研究问题

• RQ1 在 T=0 的两味 QCD 中，带规范冷却的 CLE 是否能正确再现 μ≈mN/3 处从强子物质到核物质的相变？
• RQ2 在 β=5.6 的小格点（12⁴）上，CLE 模拟在多大程度上表现出与已知物理量（如平面子与手征子凝聚）的系统性偏差？
• RQ3 在 β=5.7 条件下，格点尺寸增大至 16⁴ 时，CLE 结果与 μ=0 条件下的参考值之间的符合度如何改善？
• RQ4 尽管漂移项因费米子行列式零点而存在极点，CLE 是否仍能在弱耦合与连续极限下收敛至正确分布？
• RQ5 对规范冷却或漂移项进行何种修改可进一步减少 CLE 失败现象，并改善向正确结果的收敛？

主要发现

• 在 β=5.6、m=0.025 的 12⁴ 格点上，CLE 模拟在 μ≈0.33（mN/3）附近显示出相变，与预期一致，但物理量存在系统性偏差。
• 手征子凝聚随 μ 增加单调减小，而非在相变前保持恒定，表明与预期行为存在显著偏差。
• 在小 μ 下，平面子表现出微小但非零的偏差，而夸克数密度在大 μ 下正确饱和于 3。
• 在 β=5.7 的 16⁴ 格点上，CLE 的平面子（0.42374(4)）与手征子凝聚（0.1738(11)）与 RHMC 结果（分别为 0.42305(1) 与 0.1754(2)）表现出显著改善的符合度。
• 在 μ=0 时，平面子的系统性误差减少约两倍，手征子凝聚的误差减少近一个数量级，表明 CLE 可能在弱耦合与连续极限下给出正确结果。
• 在 β=5.7 的 16⁴ 格点上，μ>0 的初步结果表明与已知值的符合度持续改善，支持 CLE 在连续极限下可能有效的希望。