[论文解读] Gauge cooling for the singular-drift problem in the complex Langevin method --- an application to finite density QCD

Keitaro Nagata, Hideo Matsufuru+2 more|arXiv (Cornell University)|Nov 24, 2016

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

一句话总结

本研究在 $4^3\times8$ 晶格上应用带规范冷却的复数朗之万方法研究有限密度QCD，利用该方法稳定模拟并控制幺正性范数。重子数密度的 onset 出现在 $\mu a \approx 0.55$，与平均场预测一致，且高于相位截断QCD中的结果，其有效性通过漂移项概率分布的指数衰减准则在 $\mu a \lesssim 0.65$ 范围内得到验证。尽管面临奇异漂移问题的挑战，该方法在束缚相QCD模拟中仍展现出前景。

ABSTRACT

We study full QCD at finite density and low temperature with light quark mass using the complex Langevin method. Since the singular drift problem turns out to be mild on a $4^3 imes 8$ lattice we use, the gauge cooling is performed only to control the unitarity norm in this exploratory study. We report on our preliminary data obtained from the complex Langevin simulation up to certain Langevin time. While the data are still noisy due to lack of statistics, the onset of the baryon number density seems to occur at larger $μ$ than half the pion mass, which is the value for the phase quenched QCD. The validity of our simulation is tested by the recently proposed criterion based on the probability distribution of the drift term.

研究动机与目标

• 研究复数朗之万方法在有限重子密度和低温条件下的可行性，特别是在束缚相中的应用。
• 解决有限密度QCD中费米子矩阵出现零特征值所引发的奇异漂移问题。
• 利用近期提出的基于漂移项概率分布指数衰减的准则，检验复数朗之万模拟的有效性。
• 探索即使奇异漂移问题较轻微，规范冷却在控制幺正性范数方面是否仍能有效稳定模拟。

提出的方法

• 将SU(3)规范链接变量复化为SL(3,\mathbb{C})，以通过复数朗之万方程实现随机量化。
• 应用规范冷却以控制幺正性范数，本研究中未直接用于缓解奇异漂移问题。
• 采用平面元作用量和随机动量费米子形式，引入化学势 $\mu$ 以生成复数作用量。
• 通过离散步长 $\epsilon$ 数值模拟朗之万动力学，包含来自复化作用量的高斯噪声和漂移项。
• 通过配分函数的全纯延拓测量可观测量，如重子数密度和手征凝聚。
• 利用漂移项概率分布必须呈指数衰减或更快以确保收敛的准则验证模拟结果。

实验结果

研究问题

• RQ1 在存在奇异漂移项潜在风险的情况下，带规范冷却的复数朗之万方法能否在束缚相中为有限密度QCD提供可靠结果？
• RQ2 在完整QCD中，重子数密度在何种夸克化学势 $\mu$ 值下开始出现？其与相位截断QCD中的结果相比如何？
• RQ3 近期提出的基于漂移项概率分布指数衰减的准则是否适用于有限 $N_s$ 的格点QCD？
• RQ4 随着 $\mu$ 增大，费米子矩阵的特征值分布如何演化？这对奇异漂移问题有何启示？
• RQ5 在奇异漂移问题较轻微的情况下，如本研究中的 $4^3\times8$ 晶格设置，规范冷却在多大程度上能稳定朗之万动力学？

主要发现

• 重子数密度在 $\mu a \approx 0.4$ 到 $0.5$ 之间开始上升，在 $\mu a = 0.5$ 时达到 $\sim 0.01$，对应每晶格体积约一个重子。
• 完整QCD中重子数密度的 onset 出现在 $\mu a \approx 0.55$，高于相位截断QCD中的 $\mu a \approx 0.265$，与平均场预测一致。
• 费米子矩阵的特征值分布显示在 $\mu a = 0.4$ 时存在明显间隙，但在 $\mu a = 0.65$ 时该间隙减弱，表明特征值更接近原点，奇异漂移风险增加。
• 漂移项的概率分布对 $\mu a \lesssim 0.65$ 呈指数衰减，满足近期提出的复数朗之万方法正确收敛准则。
• 自相关时间在 $\mu a \gtrsim 0.4$ 显著增加，表明可能与奇异漂移问题的 onset 相关的动力学问题。
• 在指数衰减准则下，模拟在 $\mu a \lesssim 0.65$ 范围内保持稳定且有效，表明该方法在该参数范围内具有可行性，尽管存在初步噪声和有限统计问题。