[论文解读] The canonical approach to Finite Density QCD

Slavo Kratochvila, Philippe de Forcrand|ArXiv.org|Sep 26, 2005

High-Energy Particle Collisions Research被引用 9

一句话总结

本文通过固定重子数而非化学势，提出一种有限密度QCD的规范方法，实现了对μ/T ≈ 2以内的QCD相图的精确研究。利用重权重法与费米子行列式傅里叶分解的蒙特卡洛模拟，作者在T–ρ平面中识别出共存区域，并观察到在高化学势下一阶相变边界出现“下弯”现象，与强耦合理论预测一致。

ABSTRACT

We present a canonical approach to study properties of QCD at finite baryon density rho, and apply it to the determination of the phase diagram of four-flavour QCD. For a pion mass of about 350 MeV, the first-order transition between the hadronic and the plasma phase gives rise to a co-existence region in the T-rho plane, which we study in detail. We obtain accurate results for systems containing up to 30 baryons and quark chemical potentials mu up to 2T. Our T-mu phase diagram agrees with the literature when mu/T < 1. At larger chemical potential, we observe a ``bending down'' of the phase boundary. We characterise the two phases with simple models: the hadron resonance gas in the hadronic phase, the free massless quark gas in the plasma phase.

研究动机与目标

• 通过使用固定重子数而非化学势的规范形式，克服有限密度QCD模拟中的符号问题。
• 准确确定四味QCD在T–μ与T–ρ平面中的QCD相图，其中mπ ≈ 350 MeV。
• 通过麦克斯韦构造法实现自由能相等，识别出强子相与夸克-胶子等离子体相之间的共存区域。
• 在标准广义-canonical方法因符号问题失效的高密度区域，检验规范方法的有效性。
• 与现有文献及强耦合理论预测进行比较，特别关注大μ/T下临界线的行为。

提出的方法

• 通过在虚化学势下计算广义-canonical配分函数Z_GC(T, iμ_I T)的傅里叶变换，提取规范配分函数Z_C(T, Q)。
• 在固定β和虚μ_I下进行蒙特卡洛模拟，并通过Ferrenberg-Swendsen重权重法重建重子数投影的配分函数。
• 利用逸度展开与鞍点近似计算重子密度ρ(μ)，并通过麦克斯韦构造法验证一致性。
• 通过条件∫_ρ₁^ρ₂ (f’(ρ) – 3μ) dρ = 0，强制要求强子相与等离子体相的自由能密度相等，从而确定相边界。
• 该方法通过在规范系综中固定重子数，避免了广义-canonical系综中的符号问题。
• 通过与文献数据对比，并检查逸度展开与麦克斯韦构造之间的一致性，验证了该方法的有效性。

实验结果

研究问题

• RQ1 在μ/T > 1的区域，符号问题限制了标准方法的适用性，该规范方法是否能实现QCD相图的可靠计算？
• RQ2 在T–μ平面中，一阶相变边界的形状如何，特别是在μ/T ≈ 1之后？
• RQ3 在高重子密度下，强子相与等离子体相在自由能密度与状态方程方面如何表现？
• RQ4 能否在T–ρ平面中通过麦克斯韦构造法准确确定禁闭相与解禁闭相之间的共存区域？
• RQ5 观察到的临界线“下弯”现象是否与强耦合理论预测及其他格点研究结果一致？

主要发现

• 在T–μ平面中，相边界从μ/T ≈ 1.3开始出现“下弯”行为，偏离线性外推结果，与强耦合理论预测一致。
• 在T–ρ平面中，共存区域的重子密度为常数：对于T ≤ 0.85T_c，ρ_confined = 0.50(5) B/fm³，ρ_QGP = 1.8(3) B/fm³。
• 在T/T_c = 0.92时，临界化学势为μ/T = 1.06(2)，通过麦克斯韦构造法确定，确保两相自由能相等。
• 等离子体相中的状态方程与Stefan-Boltzmann定律高度吻合，ρ/T³的线性项与立方项拟合系数分别为0.82(2)与1.94(6)。
• 有限尺寸效应在逸度展开中掩盖了一阶相变，但麦克斯韦构造法给出的体积修正为指数级小。
• 该规范方法在μ/T ≤ 1时成功复现了文献结果，但在更高μ/T下揭示了新特征，表明标准方法存在局限性。