Mean-$p_T$ fluctuations in Au+Au collisions at $\sqrt{s_{\rm NN}}=3.0$--$19.6$ GeV within JAM2

技術要約：JAM2における $\sqrt{s_{NN}} = 3.0–19.6$ GeV の Au+Au 衝突における平均 $p_T$ ゆらぎ

問題と動機
平均横運動量（$\langle p_T \rangle$）のイベントごとのゆらぎは、相対論的重イオン衝突における初期状態のエネルギー沈着、横方向のサイズゆらぎ、およびそれらが集団ダイナミクスへと変換される過程を調べるための敏感なプローブとして機能する。RHICのビームエネルギー・スキャン（BES）の文脈において、系が大きな正バリオン密度に達する中では、これらのゆらぎを理解することは極めて重要である。このエネルギー領域における系の進化は、バリオンの停止、ハドロン相互作用、核平均場、および有限密度QCDの状態方程式（EoS）の影響を受ける。流体モデルはこれらのゆらぎを初期幾何学や集団膨張に関連付けるが、HIJINGやAMPTのような他のフレームワークは参照値を提供するものである。一方で、バリオン平均場を組み込んでいるが、明示的なパートン段階を欠いた微視的な輸送計算の必要性が存在する。本研究は、二次的な $p_T$ 相関、特に粒子種依存性が、これらの条件下でどのように振る舞うかを扱うことで、臨界力学を分離する前の、バリオンに富んだ物質のベースラインを確立することを目的としている。

手法とモデルの設定
本研究では、RQMDv 平均場モードおよび MH2 パラメトリゼーションを用いた JAM2 微視的輸送モデルを利用している。主な手法の特徴は以下の通りである：

• ダイナミクス： モデルは、バイナリ散乱、レゾナンスの励起／崩壊、およびストリング断片化を介した非平衡なハドロン進化をシミュレートする。BES エネルギーにおいては、ローレンツ・ベクトル・ポテンシャルとして Skyrme 型の密度および運動量依存性を持つハドロン平均場を実装している。スカラー・ポテンシャルは消失している（$S_i=0$）が、運動量依存性は一粒子分散関係に残存している。
• パラメータ： MH2 パラメトリゼーションは、飽和密度において $K = 380$ MeV の核非圧縮率を採用している。同一のパラメータがすべての衝突エネルギー（$\sqrt{s_{NN}} = 3.0–19.6$ GeV）および中心度に対して適用され、進化は 50 fm/c まで追跡される。
• 観測量： 本研究では、二次相関関数 $\langle c_2 \rangle$ および無次元のスケールされたゆらぎ $k_2 = \sqrt{\langle c_2 \rangle}/\langle \langle p_T \rangle \rangle$ を定義する。これらは $|\eta| < 0.5$ かつ $0.2 < p_T < 2.0$ GeV/c の範囲内の荷電粒子に対して計算される。
• イベント選択： 中心度は、実験的な多重度に基づく選択とは異なり、JAM2 のイベント記録から直接導出された参加陽子数（$\langle N_{part} \rangle$）による真のレベルの分類を用いて定義される。同定粒子スペクトルは、$\pi^\pm$ および $p+\bar{p}$ について個別に分析される。

主要な結果

1. 一粒子ベースライン： JAM2/RQMDv の計算は、$\sqrt{s_{NN}} = 7.7–19.6$ GeV における $\pi^+$, $K^+$, および陽子のミッドラピディティにおける同定粒子 $p_T$ スペクトルを定性的に再現している。モデルは中心度の順序付けとスペクトルの形状を捉えており、ゆらぎ分析のための一粒子ベースラインを提供しているが、二粒子相関や 7.7 GeV 未満のスペクトルを独立して検証するものではない。
2. スケールされたゆらぎ ($k_2$)： スケールされたゆらぎ $k_2$ は、すべてのエネルギーにおいて $\langle N_{part} \rangle$ の増加とともに減少する。これは、ソースの数が増加することで相関のある粒子対が希釈されることと一致している。
   • $\sqrt{s_{NN}} = 7.7–19.6$ GeV において、計算された大きさおよび中心度依存性は、利用可能な STAR の測定値と概ね一致している。
   • $\sqrt{s_{NN}} = 3.0–4.5$ GeV において、計算された中心度依存性は、特に中心的な衝突において、観測されたデータよりも強い傾向を示す。著者らは、この不一致について、中心度選択、受容角、および粒子組成の効果に関するさらなる調査が必要であると述べている。
3. 粒子種依存性 ($\langle c_2 \rangle$)：
   • 陽子 ($p+\bar{p}$): 正規化されていない相関関数 $\langle c_2 \rangle$ は、すべてのエネルギーおよび中心度において正の値を維持している。その大きさは、特に周辺的な衝突において、荷電パイオンよりも大幅に大きい。著者らは、この階層関係は、相対的なゆらぎの強さではなく、陽子のより大きな特徴的な横運動量スケールを反映している可能性があると注意を促している。
   • 荷電パイオン ($\pi^\pm$): 相関関数は明確なエネルギー依存性を示す。3.0 および 3.5 GeV ではほとんどの centrality 間隔で負またはゼロに近い値であるが、4.5 GeV では弱く正となり、より高いエネルギー（7.7–19.6 GeV）では正の値を維持する。
4. 楕円流との比較： 荷電パイオンの $\langle c_2 \rangle$ のエネルギー進化は、反応面に対する楕円流（$v_2$）が、スペクテーターによる squeeze-out（押し出し）から in-plane 膨張へと転移する 4.5 GeV 付近で負から正へと変化する様子を定性的に模している。しかし、本論文は、スカラー $p_T$ 相関と方位角高次成分が異なるダイナミクスをプローブしているため、この比較が共通の微視的な起源を確立するものではないことを明記している。

意義と主張
本論文は、明示的なパートン段階を含まないものの、バリオン平均場を含む輸送フレームワーク内での平均 $p_T$ ゆらぎのモデル参照を確立している。

• ベースラインの確立： 本研究は、バリオン平均場を含む輸送計算内での、種依存の予測を提供し、これらは BES 領域における実験データの解釈のためのベースラインとして機能する。
• 限定的な範囲： 著者らは、本研究が EoS を抽出したり、個々のメカニズム（平均場 vs 再散乱 vs スペクテーター）を分離したり、あるいは支配的なパートン相互作用の出現を確認したりするものではないことを強調している。$v_2$ の傾向との比較は、記述的なものであり因果的なものではない。
• 今後の方向性： 著者らは、特定の寄与を孤立させるには、輸送ダイナミクスの制御された変化が必要であると示唆している。将来の研究には、電荷分解された相関関数、中心度および崩壊の寄与のチェック、および代替の平均場設定との比較が含まれるべきであると提案している。さらに、低エネルギーへの計算の拡張、および HIRFL-CSR の CEE プログラムや HIAF プログラムに関連する定量的予測のための実験的受容角との一致も必要である。

要約すると、本研究は、実験データで見られる定性的な傾向（パイオン相関の符号の変化やエネルギー依存性など）が、パートン相相を持たないハドロン輸送モデルからも生じ得ることを示しており、基礎となる微視的なメカニズムを特定するためには、制御された理論的変分が必要であることを浮き彫りにしている。