[論文レビュー] 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

ひとこと要約

本稿では、1つのN×N行列分解を用いて、格子QCDにおける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相図を探索する。
• 相関関数が多数のオーダーにわたり変動する状況を、積率展開や対数正規近似などの統計的手法を用いて分析する。

提案手法

• アルゴリズムは、1つのN×N行列分解を用いてn粒子相関関数を計算し、従来の複数回の分解を要する手法よりも著しく改善されている。
• 相関関数の対数の積率展開を用いて、統計的モーメントから相関関数を推定し、N_κ個の積率で切り捨てることで系統的改善が可能である。
• 相関関数の分布が対数正規分布に従うと仮定し、2次および3次積率切り捨ての結果が一貫していることから妥当性が検証された。
• 統計的解析では、積率に基づく補正を用いて、多数のオーダーにわたる変動を示す相関関数を扱い、バイアスを低減しノイズを管理する。
• 手法は、m_π = 170 MeVで、体積が(4.4³×8.8) fm⁴および(5.8³×11.6) fm⁴の2つのゲージ場アンサンブルに適用された。
• データ可視化には、各データポイントの不確実性楕円をきっちり包む滑らかな誤差帯を生成するためのカスタム補間法が用いられた。

実験結果

リサーチクエスチョン

• RQ1 QCD相図の、大アイソスピン密度・零バリオン密度領域におけるエネルギー密度の振る舞いは何か？
• RQ2 有効モデルの予測通り、μ_I ≈ m_πでボーズ・アインシュタイン凝縮転移が発生するか？
• RQ3 音速はアイソスピン化学ポテンシャルの増加に伴いどのように変化するか？また、共形限界を超えるか？
• RQ4 多数粒子系の相関関数は、どの程度まで対数正規分布で記述可能か？
• RQ5 微小なQCDが適用可能と予想される領域において、エネルギー密度や圧力などの熱力学的観測量はどのように変化するか？

主な発見

• μ_I ≈ 1.5m_πでエネルギー密度にピークが観測され、ボーズ・アインシュタイン凝縮相の始まりを示唆する。
• 広いμ_Iの範囲で断熱的音速が共形限界（c_s² ≤ 1/3）を超えており、理想気体の振る舞いを越えた強い相関が示唆される。
• 相関関数は対数正規分布でよく近似され、データからエネルギー準位を安定して抽出可能である。
• 積率展開法では、2次および3次切り捨ての結果が一貫しており、高次項は統計的ノイズを増加させるが、対数正規性の破綻を示唆しない。
• μ_I ≈ 15m_πをはるかに超えて、系は強い相互作用的性質を示し続け、微小なQCDや理想気体モデルでは中間の化学ポテンシャル領域を記述できないことが示された。
• このアルゴリズムにより、6,144個の電荷を帯びたπ粒子を含む系の相関関数の最初の格子QCD計算が可能となり、従来の研究を1桁以上超えた。