[论文解读] Lattice quantum chromodynamics at large isospin density: 6144 pions in a box

Ryan Abbott, William Detmold+7 more|arXiv (Cornell University)|Jul 27, 2023

High-Energy Particle Collisions Research被引用 5

一句话总结

本文提出了一种数值稳定的算法，用于在格点QCD中使用单个N×N矩阵分解计算n-单粒子关联函数，实现了最多6,144个带电π介子的模拟。该方法在μ_I ≈ 1.5m_π处揭示了能量密度的峰值，表明发生了玻色-爱斯坦凝聚相变，且在较宽范围的同位旋化学势下，声速超过了共形极限。

ABSTRACT

We present an algorithm to compute correlation functions for systems with the quantum numbers of many identical mesons from lattice quantum chromodynamics (QCD). The algorithm is numerically stable and allows for the computation of $n$-pion correlation functions for $n \in \{ 1, \dots, N\}$ using a single $N imes N$ matrix decomposition, improving on previous algorithms. We apply the algorithm to calculations of correlation functions with up to 6144 $π^+$s using two ensembles of gauge field configurations generated with quark masses corresponding to a pion mass $m_π= 170$ MeV and spacetime volumes of $(4.4^3 imes 8.8)\ { m fm}^4$ and $(5.8^3 imes 11.6)\ { m fm}^4$. We also discuss statistical techniques for the analysis of such systems, in which the correlation functions vary over many orders of magnitude. In particular, we observe that the many-pion correlation functions are well approximated by log-normal distributions, allowing the extraction of the energies of these systems. Using these energies, the large-isospin-density, zero-baryon-density region of the QCD phase diagram is explored. A peak is observed in the energy density at an isospin chemical potential $μ_I\sim 1.5 m_π$, signalling the transition into a Bose-Einstein condensed phase. The isentropic speed of sound in the medium is seen to exceed the ideal-gas (conformal) limit ($c_s^2\leq 1/3$) over a wide range of chemical potential before falling towards the asymptotic expectation at $μ_I\sim 15 m_π$. These, and other thermodynamic observables, indicate that the isospin chemical potential must be large for the system to be well described by an ideal gas or perturbative QCD.

研究动机与目标

• 开发一种用于在格点QCD中计算最多N个相同π介子关联函数的数值稳定算法。
• 实现大同位旋化学势下多π介子系统的第一性原理计算，避免符号问题。
• 利用第一性原理格点QCD模拟探索QCD相图中大同位旋密度、零重子密度区域的性质。
• 使用累积量展开和对数正态近似等统计技术分析跨越多个数量级的关联函数。

提出的方法

• 该算法通过单次N×N矩阵分解计算n-π介子关联函数，显著优于以往需要多次分解的方法。
• 采用关联函数对数的累积量展开，从统计矩估计关联函数，通过截断至N_κ阶累积量实现系统性改进。
• 该方法假设关联函数分布为对数正态分布，其合理性通过二阶与三阶累积量截断结果的一致性得到验证。
• 统计分析采用基于累积量的修正方法，以处理在多个数量级范围内变化的关联函数，降低偏差并有效控制噪声。
• 该技术应用于两个规范场组，其π介子质量为m_π = 170 MeV，体积分别为(4.4³×8.8) fm⁴和(5.8³×11.6) fm⁴。
• 数据可视化采用自定义插值方法，生成平滑的误差带，紧密包裹每个数据点的不确定性椭圆。

实验结果

研究问题

• RQ1 在QCD相图的大同位旋密度、零重子密度区域，能量密度的行为如何？
• RQ2 是否如有效模型所预测的那样，在μ_I ≈ m_π处发生玻色-爱斯坦凝聚相变？
• RQ3 声速如何随同位旋化学势的增加而演化？是否超过共形极限？
• RQ4 多π介子系统的关联函数在多大程度上可由对数正态分布描述？
• RQ5 在微扰QCD预计适用的区域，热力学可观测量（如能量密度和压强）如何演化？

主要发现

• 在μ_I ≈ 1.5m_π处观察到能量密度的峰值，标志玻色-爱斯坦凝聚相的出现。
• 在较宽范围的μ_I内，等熵声速超过共形极限（c_s² ≤ 1/3），表明存在强关联，超出理想气体行为。
• 关联函数可良好地由对数正态分布近似，从而实现从数据中稳健提取能级。
• 累积量展开方法在二阶与三阶截断下结果一致，高阶项虽增加统计噪声，但未显示对数正态性被破坏的迹象。
• 系统在μ_I ≈ 15m_π之后仍保持强相互作用，表明微扰QCD或理想气体模型在中等化学势下不足以描述介质性质。
• 该算法实现了格点QCD中最多6,144个带电π介子系统关联函数的首次计算，使先前研究的规模扩展了一个多数量级。