[论文解读] Quark Mass Dependence of the QCD Critical End Point in the Strong Coupling Limit

Jangho Kim, Wolfgang Unger|arXiv (Cornell University)|Nov 28, 2016

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

一句话总结

本文通过在格点QCD的对偶表示中使用蒙特卡洛模拟，在强耦合极限下研究了QCD临界端点（CEP）的夸克质量依赖性。通过采用蠕虫算法处理符号问题，绘制了μ-T相图，并确定了有限夸克质量下的CEP位置，表明临界端点随mq变化而移动，且非微扰各向异性ξ(γ)消除了相边界中的非物理回弯现象。

ABSTRACT

Strong coupling lattice QCD in the dual representation allows to study the full $μ$-$T$ phase diagram, due to the mildness of the finite density sign problem. Such simulations have been performed in the chiral limit, both at finite $N_t$ and in the continuous time limit. Here we extend the phase diagram to finite quark masses, with an emphasis on the low temperature first order transition. We present our results on the quark mass dependence of the critical end point and the first order line obtained by Monte Carlo via the worm algorithm.

研究动机与目标

• 将QCD的μ-T相图扩展至手征极限之外的有限夸克质量区域。
• 在强耦合极限下，确定不同夸克质量下的临界端点（CEP）位置。
• 通过应用非微扰各向异性ξ(γ)而非微扰的ξ' = γ²，消除一阶相变边界中出现的非物理回弯现象。
• 与先前的蒙特卡洛模拟和平均场理论结果进行比较，验证该方法的准确性。

提出的方法

• 本研究采用强耦合格点QCD的对偶表示，其中规范场和格拉斯曼变量被积分掉，从而得到介子和重子的世界线表示。
• 配分函数以介子跃迁（kb）、有限夸克质量下的单体（nx）以及具有缠绕数rℓ和符号σ(ℓ)的重子环（ℓ）的形式表达。
• 采用蠕虫算法进行蒙特卡洛采样，以计算诸如手征凝聚〈χ̄χ〉、手征 susceptibility χch、重子密度〈nB〉和重子 susceptibility χB 等可观测量。
• 利用O(2)、Z(2)和过渡（ρ=1）临界指数进行有限尺寸标度，以识别相变并定位CEP。
• 通过非微扰方式推导各向异性参数ξ(γ)，以正确映射时间和空间尺度，替代不合理的ξ' = γ²。
• 通过susceptibility和condensate数据的峰值分析来定位CEP，提取不同mq下的aTE'和aμE'。

实验结果

研究问题

• RQ1 在强耦合极限下，随着夸克质量的增加，μ-T平面上的临界端点（CEP）如何移动？
• RQ2 使用非微扰各向异性ξ(γ)是否能消除使用ξ' = γ²时观察到的一阶相变边界中的非物理回弯现象？
• RQ3 本蒙特卡洛研究中CEP的位置与先前蒙特卡洛结果和平均场理论预测在定量上如何比较？
• RQ4 在有限mq下，单体涨落对手征 susceptibility 的作用是什么，它如何影响CEP确定的精度？
• RQ5 当使用正确的非微扰各向异性时，CEP轨迹是否满足(amq)^{2/5}的三临界标度行为？

主要发现

• 临界端点（CEP）随夸克质量移动：当amq = 0.05时，CEP位于aTE' = 0.73(3) 和 aμE' = 0.7808(5)，在物理单位下为aTE = 0.62(3) 和 aμE = 0.6653(4)。
• 使用非微扰各向异性ξ(γ)消除了使用ξ' = γ²时在低温区域出现的非物理回弯现象，尤其在低温度区域表现显著。
• 有限mq下CEP轨迹与平均场理论存在偏差，包括CEP位置和(amq)^{2/5}斜率均与平均场预测不一致。
• 当amq = 0.05时，手征 susceptibility 峰值出现在aT' = 0.725 至 0.75之间，表明CEP位于此范围内，且由Z(2)临界指数确认。
• amq = 0.01 和 0.02 时的结果与先前蒙特卡洛数据一致，验证了方法的准确性和数值稳定性。
• 在aT' = 0.7处使用ρ = 3进行有限尺寸标度，确认了相变的一阶性质；而在aT' = 0.95处使用ρ = 1进行过渡标度。