Critical net-proton number fluctuations with hydrodynamics

技術要約：流体力学を用いた正の陽子数揺らぎの決定的な変動

問題提起
QCD臨界点（CEP）の所在および存在は、QCD相図の研究における未解決かつ重要な課題である。関数的繰り込み群（fRG）、ダイソン・シュウィンガー方程式（DSE）、および格子QCDの補外といった理論的アプローチは、CEPが高バリオン化学ポテンシャル（$\mu_B/T \gtrsim 4-5$）において存在することを示唆しているが、実験による直接的な検証は困難である。揺らぎの観測量、特に正の陽子数（net-proton number）の揺らぎは、CEPを探索するためのバリオン数揺らぎのプロキシ（代理指標）として機能する。しかし、従来の簡略化されたモデルは、衝突エネルギーごとに単一のバリオン化学ポテンシャルと温度を仮定する（化学凍結曲線によってパラメータ化される）ことが多く、重イオン衝突の現実的な時空発展を捉えきれていないため、実験信号に対する信頼性の高い理論的予測を得ることは困難であった。

手法
本研究では、fRG計算と最先端の(3+1)次元粘性流体力学シミュレーションを組み合わせ、9つの衝突エネルギー（$\sqrt{s_{NN}} = 7.7 - 200$ GeV）における正の陽子数揺らぎを計算する。

1. 流体力学的発展： Au+Au衝突の時空発展は、速度依存性を持つ3D初期条件を備えたMUSIC粘性流体力学フレームワークを用いてシミュレートされる。その後、一定のエネルギー密度超曲面（$\epsilon_{fo} = 0.26$ GeV/fm$^3$）上でiS3D粒子化が行われ、UrQMDハドロン後方散乱（afterburner）が適用される。このモデルは、指定されたエネルギーにおけるハドロン観測量（収量、スペクトルなど）に較正されている。
2. 臨界揺らぎのエンコーディング： 凍結超曲面は流体要素の集合として扱われ、各要素は局所的な温度 $T(x_i)$ とバリオン化学ポテンシャル $\mu_B(x_i)$ によって特徴付けられる。fRGを通じて計算されたCEPによる臨界寄与を含む正のバリオン数揺らぎが、これらの流体要素にマッピングされる。
3. 正のバリオン数から正の陽子数の揺らぎへ： fRGは（バリオンと反バリオンを分離したものではなく）正のバリオン数揺らぎ（$\chi_n^B$）を提供するため、著者らはハドロン共鳴ガス（HRG）モデルに基づく関係式（式4）を用いて（反）バリオン揺らぎを推論する。正の陽子数の揺らぎは、確率 $q \approx 1/2$ を仮定したアイソスピン無作為化（二項フィルタリング）を通じて、バリオン揺らぎから導出される。
4. 実験的制約： 計算には以下の実験的実態を実装している：
   • 受容角カット： STAR検出器の運動学（$0.4 < p_T < 2.0$ GeV, $|y| < 0.5$）に一致させるため、二項受容補正が適用される。
   • 全域保存： 全域的なバリオン数保存の効果を補正するために、サブアンサンブル受容法（SAM）が用いられる。
5. 比較ベースライン： 結果は以下と比較される：
   • 同一の流体力学フレームワーク内でのHRG感受性を用いた非臨界ベースライン。
   • 流体力学的発展を考慮せずに、パラメータ化された化学凍結曲線に沿って計算された従来のfRG結果。

主な貢献

• 現実的なフレームワーク： 本論文は、衝突エネルギーごとに単一の凍結点を仮定するという過度に単純化された仮定を超え、fRG由来の臨界揺らぎを現実的な流体力学的凍結超曲面にエンコードする手法を確立した。
• 包括的な実装： 臨界揺らぎを、実験的な受容角カット、アイソスピン無作為化、および全域的なバリオン保存と、動的な背景の中で同時に統合している。
• 系統的なエネルギースキャン： 研究は全ビームエネルギースキャン（BES）範囲（$\sqrt{s_{NN}} = 7.7 - 200$ GeV）をカバーしており、衝突エネルギー依存性の詳細な調査を可能にしている。

結果
本研究では、正の陽子分布の累積モーメント比 $C_2/C_1$, $C_3/C_2$, $C_4/C_2$ を計算している：

• 低次累積モーメント ($C_2/C_1$)： 臨界揺らぎを含む結果（fRG+hydro）と非臨界の結果（HRG+hydro）の差は、全エネルギー範囲にわたって無視できる程度である。
• 高エネルギーにおける高次累積モーメント： $\sqrt{s_{NN} \gtrsim 14.5}$ GeVにおいて、$C_3/C_2$ および $C_4/C_2$ の臨界結果と非臨界結果の差は小さい。これは、これらのエネルギーが臨界揺らぎが抑制される低 $\mu_B$ 領域を探索しているためである。
• 低エネルギーにおける高次累積モーメント： 低衝突エネルギー領域（例：$\sqrt{s_{NN}} = 7.7$ GeV）では、高 $\mu_B$ を探索するため、臨界結果と非臨界結果の差は累積モーメントの次数とともに増大する。
• 非単調挙動： 臨界揺らぎを含めた場合、$C_4/C_2$ において明確な非単調の衝突エネルギー依存性が観察される。この特徴は、非臨界（HRG）ベースラインには存在しない。
• 実験との比較： fRG+hydroの結果は、パラメータ化された凍結曲線上で計算された簡略化されたfRG結果よりも、STAR BES-IIの実験データとよく一致している（簡略化されたモデルは非単調性を過大評価する傾向がある）。

意義と主張
本論文は、fRG計算と流体力学シミュレーションを組み合わせることが、パラメータ化された凍結モデルと比較して、重イオン衝突における揺らぎ観測量を記述するための有望かつより現実的なアプローチであることを主張している。著者らは、低次累積モーメントは非臨界物理学と実験的制約に支配されている一方で、低エネルギー領域の高次累積モーメントはCEPに対して敏感であることを強調している。特に $C_4/C_2$ における非単調挙動の観察は、現実的な流体力学的発展が考慮されている限り、それがCEPのシグネチャーとなり得ることを示唆している。著者らは、CEP領域をさらに探索するために、将来の研究においてこれらの計算を固定ターゲットエネルギー範囲（$\sqrt{s_{NN}} < 7.7$ GeV）へ拡張すべきであると控えめに述べている。