[论文解读] Probability distribution functions in the finite density lattice QCD

Shinji Ejiri, Yoshiyuki Nakagawa+6 more|arXiv (Cornell University)|Dec 4, 2012

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

一句话总结

该论文提出一种基于直方图的多点重加权方法，用于有限密度格点QCD，研究高温高化学势区域的QCD相结构，重点分析了Polyakov环分布。研究确定了重夸克区域中一阶相变与交叉相变的临界面，表明在临界点附近，夸克行列式复相位的影响可忽略不计。

ABSTRACT

We study the phase structure of QCD at high temperature and density by lattice QCD simulations adopting a histogram method. We try to solve the problems which arise in the numerical study of the finite density QCD, focusing on the probability distribution function (histogram). As a first step, we investigate the quark mass dependence and the chemical potential dependence of the probability distribution function as a function of the Polyakov loop when all quark masses are sufficiently large, and study the properties of the distribution function. The effect from the complex phase of the quark determinant is estimated explicitly. The shape of the distribution function changes with the quark mass and the chemical potential. Through the shape of the distribution, the critical surface which separates the first order transition and crossover regions in the heavy quark region is determined for the 2+1-flavor case.

研究动机与目标

• 使用格点QCD模拟研究有限温度与化学势下的QCD相结构。
• 通过分析Polyakov环的概率密度函数（直方图）解决有限密度QCD中的符号问题。
• 确定在重夸克条件下，2+1味QCD中一阶与交叉相变的临界面。
• 估计夸克行列式复相位对临界点位置的影响。
• 验证多点重加权技术在有限密度QCD中用于宽范围直方图重建的适用性。

提出的方法

• 使用直方图方法，计算Polyakov环概率密度函数随夸克质量和化学势的变化。
• 应用多点重加权技术，结合不同β值的模拟，实现对序参量的宽范围采样。
• 采用夸克行列式的领先阶跳跃展开，将其表示为plaquette和Polyakov环算符的函数。
• 通过相位因子的累积量展开定义有效势，其贡献来自相位角的二阶、四阶和六阶累积量。
• 以相位压制的有效势为基准，通过累积量展开计算复相位因子的修正。
• 通过有效势出现双阱结构的条件，结合夸克质量与化学势的关系，确定临界点。

实验结果

研究问题

• RQ1 在重夸克极限下，Polyakov环概率密度函数的形状如何随夸克质量和化学势变化？
• RQ2 多点重加权技术能否可靠地重建宽范围耦合参数下的完整分布函数？
• RQ3 在重夸克区域，夸克行列式复相位对临界点位置有何影响？
• RQ4 由直方图导出的有效势是否能可靠地作为识别一阶相变的信号？
• RQ5 在(κ_ud, κ_s, μ)参数空间中，临界面如何随化学势增加而演化？

主要发现

• Polyakov环分布函数的形状随夸克质量和化学势显著变化，在一阶相变点表现出双模态行为。
• 通过直方图方法确定了在重夸克区域中2+1味QCD一阶与交叉相变的临界面。
• 在重夸克区域，对于24³×4格点，夸克行列式复相位对临界点位置的影响可忽略不计，二阶累积量主导相位修正。
• 相位因子对有效势的修正在Polyakov环值较小时较大，但对临界点位置无显著影响。
• 通过将N_f=2的结果外推，确定了N_f=2+1的临界点，临界条件为 2κ_ud^Nt cosh(μ_ud/T) + κ_s^Nt cosh(μ_s/T) = 2[κ_cp^N_f=2(0)]^Nt。
• 结果表明，在重夸克极限下，相位压制近似对临界点定位是可靠的，支持在探索性研究中结合相位压制模拟与重加权方法。