[论文解读] Canonical approach to the finite density QCD with winding number expansion

Atsushi Nakamura, Shotaro Oka+1 more|arXiv (Cornell University)|Apr 16, 2015

High-Energy Particle Collisions Research被引用 3

一句话总结

本文提出了一种绕数展开方法，用于在格点QCD中进行化学势展开，实现了无符号问题的有限密度可观测量计算。通过在跳跃参数中展开狄拉克行列式，并按时间方向的绕数重新分组各项，作者利用零化学势的规范配置计算了规范系综的正则系综配分函数，成功再现了手征对称性守恒和夸克数密度，且在高化学势和低温下表现出类似相变的行为。

ABSTRACT

The canonical partition function is related to the grand canonical one through the fugacity expansion and is known to have no sign problem. In this paper we perform the fugacity expansion by a method of the hopping parameter expansion in temporal direction for the lattice QCD: winding number expansion. The canonical partition function is constructed for Nf=2 QCD starting from gauge configurations at zero chemical potential. After derivation of the canonical partition function we calculate hadronic observables like chiral condensate and quark number density and the pressure at the real chemical potential.

研究动机与目标

• 通过利用无符号问题的规范系综，解决有限密度格点QCD中的符号问题。
• 通过时间方向的绕数计数，发展一种基于跳跃参数展开的系统性化学势展开方法。
• 利用零化学势模拟数据，在真实化学势下计算强子可观测量（如手征对称性守恒和夸克数密度）的广义系综。
• 通过观察在高化学势和低温下手征对称性守恒的类似相变行为，验证该方法的有效性。

提出的方法

• 广义系综配分函数通过化学势展开与规范系综关联：$ Z_G = \sum_n Z_C(n) \xi^n $，其中 $ \xi = e^{\mu/T} $。
• 狄拉克行列式在跳跃参数 $ \kappa $ 中展开，得到 $ \text{Tr} \log D_W = -\sum_{n=1}^\infty \frac{\kappa^n}{n} \text{Tr}(Q^n) $。
• 将展开式中的每一项按其时间绕数 $ n $ 重新分组，得到 $ \text{Tr} \log D_W = \sum_n w_n \xi^n $，其中 $ \xi^n = e^{n\mu/T} $。
• 行列式重写为 $ \text{Det} D_W = \sum_n z_n(U) \xi^n $，其中系数 $ z_n(U) $ 通过绕数分辨迹的傅里叶变换获得。
• 通过 $ Z_G(e^{i\theta}, T, V) $ 的傅里叶变换提取规范系综配分函数 $ Z_C(n) $，避免了数值符号问题。
• 强子可观测量通过 $ \langle O \rangle_C(n) = O_C(n)/Z_C(n) $ 计算，其中 $ O_C(n) $ 通过类似化学势展开方法获得算符夸克收缩项。

实验结果

研究问题

• RQ1 能否在时间方向通过跳跃参数展开，系统地执行广义系综配分函数的化学势展开？
• RQ2 能否通过绕数计数，从零化学势的规范配置中准确重构规范系综配分函数？
• RQ3 该方法是否能在无符号问题下，准确再现真实化学势下的物理可观测量（如手征对称性守恒和夸克数密度）？
• RQ4 在高化学势和低温下，手征对称性守恒是否表现出类似相变的行为？

主要发现

• 规范系综中的手征对称性守恒在大重子数时迅速下降，表明存在类似相变的行为。
• 在广义系综中，手征对称性守恒随化学势增加而急剧减小，尤其在低温下，表明手征对称性可能被恢复。
• 夸克数密度随化学势增加而上升，在高 $ \mu/T $ 时表现出明显上升趋势，与相变行为一致。
• 该方法成功利用仅来自零-μ 规范配置的数据，在真实化学势下计算了可观测量，避免了符号问题。
• 不同 $ \beta $ 值下的结果一致，且在较低温度下相变区域向更高 $ \mu/T $ 移动。
• 该方法在 $ 8^3 \times 4 $ 格点上得到验证，表明其在 $ n_B \sim 30 $ 量级的重子数下具有可行性。