[论文解读] Exploring Complex-Langevin Methods for Finite-Density QCD

D. K. Sinclair, J. B. Kogut|arXiv (Cornell University)|Oct 21, 2015

High-Energy Particle Collisions Research参考文献 12被引用 13

一句话总结

该论文研究了在12⁴格点上使用规范冷却的复数朗之万动力学模拟零温有限密度量子色动力学（QCD）的应用，包含两个轻夸克。尽管存在发散行为的挑战，自适应更新与规范冷却使算法稳定，结果定性正确：当μ < m_N/3时，夸克数密度从零开始上升；在饱和时，手征子凝聚度向零减小；平面子（plaquette）向冻结值增加，但因平衡化问题，接近核物质相变点时未能完全收敛。

ABSTRACT

QCD at non-zero chemical potential ($μ$) for quark number has a complex fermion determinant and thus standard simulation methods for lattice QCD cannot be applied. We therefore simulate this theory using the Complex-Langevin algorithm with Gauge Cooling in addition to adaptive methods, to prevent runaway behaviour. Simulations are performed at zero temperature on a $12^4$ lattice with 2 quarks which are light enough that $m_N/3$ is significantly larger than $m_π/2$. Preliminary results are qualitatively as expected. The quark-number density is close to zero for $μ&lt; m_N/3$, beyond which it increases, eventually reaching its saturation value of $3$ for $μ$ sufficiently large. The chiral condensate decreases as $μ$ is increased approaching zero at saturation, while the plaquette increases towards its quenched value. We have yet to observe the transition to nuclear matter at $μ\approx m_N/3$, presumably because the runs for $μ$ between $m_N/3$ and saturation have yet to equilibrate.

研究动机与目标

• 测试使用规范冷却的复数朗之万动力学是否能为零温有限密度QCD提供可靠结果。
• 研究夸克数密度、手征子凝聚度和平面子随化学势μ的变化行为。
• 确定该方法是否能捕捉到在μ ≈ m_N/3处的核物质相变以及在大μ时的颜色超导态的出现。
• 评估在复数作用量和非紧致规范场存在下，复数朗之万方法的收敛性和正确性。
• 探索将该方法扩展至更轻的夸克质量、更大格点以及有限温区域的可行性。

提出的方法

• 通过将SU(3)规范场复化为SL(3,C)，将复数朗之万动力学应用于有限μ下的QCD，以实现复数费米子行列式模拟。
• 随机演化由带有实高斯噪声的朗之万方程控制，采用部分二阶积分器以提高数值稳定性。
• 通过迭代最小化单位性范数F(U) = (1/4V)ΣTr[U†U + (U†U)⁻¹ - 2] 实现规范冷却，从而抑制规范场的无界涨落。
• 采用自适应时间步长控制漂移力，当力增大时减小dt，以提升稳定性。
• 在每次更新后测量夸克数密度、手征子凝聚度和平面子等可观测量，数据收集覆盖90,000至500,000次更新。
• 模拟初始化时将规范链接置于SU(3)流形上，以探测在μ ≈ m_N/3附近的亚稳态和平衡化行为。

实验结果

研究问题

• RQ1 使用规范冷却的复数朗之万动力学是否能为零温有限密度QCD产生稳定且收敛的结果？
• RQ2 该方法能否正确再现μ ≈ m_N/3以上夸克数密度的增加，表明核物质的出现？
• RQ3 手征子凝聚度是否随μ单调减小，在饱和时趋近于零，如在夸克-胶子等离子体相中所预期的那样？
• RQ4 在μ ≈ m_N/3处是否存在一阶相变的证据？系统是否因平衡化缓慢而被困在亚稳态中？
• RQ5 该方法能否检测到大μ时两味夸克颜色超导体的形成？与重夸克方法或重加权法的预测相比如何？

主要发现

• 当μ < m_N/3时，夸克数密度保持接近零，随后在大μ时上升至饱和值3，表明核物质的出现。
• 手征子凝聚度随μ增加而单调减小，在饱和时趋近于零，与手征对称性恢复一致。
• 平面子随μ增加而上升，趋近于冻结值，表明致密介质中规范场强度增强。
• 当μ > m_N/3时，系统表现出缓慢的平衡化，可能由于一阶相变导致的亚稳态，当前运行中未能实现完全收敛。
• 自适应时间步长和规范冷却稳定了算法，但即使经过冷却，单位性范数仍轻微偏离SU(3)，表明存在残余不稳定性。
• 初步结果定性正确，但由于不完全平衡化，特别是在μ ≈ m_N/3附近，误差条尚未获得。