An expanding spherical fireball model for light hadron production at RHIC ($\sqrt{s_{\rm NN}}=7.7$--$39$ GeV)

技術的要約：RHICにおける軽ハドロン生成のための拡張球状火球モデル

問題提起
相対論的流体力学は、RHIC（相対論的重イオン衝突器）およびLHC（大型ハドロン衝突器）における重イオン衝突（HIC）で生成されるクォーク・グルーオン・プラズマ（QGP）の時空進化を記述する上で非常に成功している。しかし、フル・ハイドロダイナミクス（流体力学）シミュレーションから特定の凍結（freeze-out）特性を抽出することは、計算負荷が高く、初期条件や輸送係数の精密な知識を必要とする。その結果、粒子スペクトルと集団的フローを記述するために「ブラストウェーブ（blast-wave, BW）」パラメトリゼーションが頻繁に採用されている。標準的なBWモデルは多くの場合、円筒形幾何学を仮定しているが、RHICにおける非中央衝突では、初期の空間的異方性が存在し、それが運動量異方性（楕円流）へと進化する。著者らは、方位角方向に積分された横運動量（$p_T$）スペクトルとラピディティ分布は、これらの方位角異方性に対しては概して鈍感であるものの、スペクトルと異方的フローの両方を同時に記述しようとする既存のモデルは、多くの自由パラメータ（例：楕円火筒モデルにおける8つのパラメータ）を必要とすることが多いと指摘している。本研究が取り組む課題は、より少ない自由パラメータを用いる、より単純な球対称火球モデルが、幅広い衝突エネルギー（$\sqrt{s_{NN}} = 7.7–39$ GeV）および中心度において、同定された軽ハドロンの$p_T$スペクトルとラピディティ分布を適切に記述できるかどうかを判断することである。

手法
著者らは、Au+Au衝突によって生成される媒質を記述するために、拡張する球状火球モデルを提案している。その手法は、以下の主要な構成要素からなる。

1. 火球の幾何学と進化： 媒質は、時間依存する半径 $r_B(t)$ を持つ球対称の火球としてモデル化される。半径の進化は、$r_B(t) = r_0 + v_\infty [t - (1 - e^{-At})/A]$ とパラメトリゼーションされる。ここで、$r_0$ は衝突する原子核の初期の楕円状の重なり体積から導かれる初期有効半径であり、$v_\infty$ は漸近的膨張速度、そして $A$ は半径方向のフローの発達率を特徴付ける。
2. フロー・プロファイル： 標準的な円筒形ブラストウェーブとは異なり、流体速度は半径方向に定義される。半径方向の流体ラピディティ $\eta_r$ は、表面における速度が半径の膨張率と一致するように決定され、内部はブラストウェーブに似たプロファイル（$v_r = (r/r_B(t))\dot{r}_B(t)$）に従う。
3. 凍結（Freeze-out）処方： 粒子スペクトルは、Cooper–Fryeの凍結処方を用いて計算される。モデルは、実験室系における一定の時間 $t_f$ で定義される運動学的凍結超曲面上で積分された、局所平衡分布関数（ボース＝アインシュタインまたはフェルミ＝ディラック統計）を仮定している。
4. パラメータ抽出戦略：
   • モデルパラメータである運動学的凍結温度（$T_{kin}$）、漸近速度（$v_\infty$）、フロー発達パラメータ（$A$）、および凍結時間（$t_f$）は、まず様々な中心度におけるパイオン（$\pi^\pm$）の中ラピディティ $p_T$ スペクトルに対してフィッティングされる。
   • これらの集団的フローのパラメータは、他のハドロン種（陽子、反陽子、カオン）に対して固定される。
   • 他の種に対する唯一の追加の自由パラメータは、粒子依存の運動学的化学ポテンシャル（$\mu_{kin}$）であり、これは実験データに適合するように調整される。
   • 2つのフィッティング・スキームがテストされる。Set-1は利用可能な全 $p_T$ 範囲（0.25–2.0 GeV）を使用し、Set-2は共鳴崩壊（resonance feed-down）や硬過程の影響が最小限に抑えられる区間に制限する（標準的なSTAR BW解析と同様）。

主な貢献と結果
本研究は、モデルを6つの中央度ビンと5つの衝突エネルギーにおける $\pi^\pm$, $K^\pm$, $p$, $\bar{p}$ に関するSTARコラボレーションのデータに照らし合わせている。

• スペクトルの記述： モデルは、同定されたすべての軽ハドロンの実験的な $p_T$ スペクトルを良好に記述している。フィッティングにより、抽出されたパラメータに系統的な傾向が得られた：
  • 中心度依存性： 衝突がより中央になるにつれ、凍結時間（$t_f$）が増加し、これはより大きな火球の寿命を反映している。逆に、運動学的凍結温度（$T_{kin}$）は減少しており、これはデカップリング（脱結合）までのより大幅な冷却を示している。
  • フローの発達： 半径方向のフロー発達パラメータ（$A$）は中央衝突においてより大きく、これは強い圧力勾配によるフローの発達がより速いことを示唆している。
  • 化学ポテンシャル： 抽出された有効運動学的化学ポテンシャル（$\mu_{kin}$）は、粒子種および中心度によって変化し、モデルに明示的に含まれていない共鳴崩壊の寄与を効果的に吸収している。
• ラピディティ分布： $p_T$ スペクトルから抽出されたパラメータを用いて、モデルはすべてのハドロン種のラピディティ分布（$dN/dy_p$）を予測する。結果は、検討されたエネルギー範囲および中心度全体にわたって一貫してガウス型に近い分布を示しており、低エネルギー実験（CERN SPS, GSI）の観測結果と一致し、$\sqrt{s_{NN}} = 39$ GeVまで拡張されている。
• パラメータの効率性： モデルは、わずか5つの自由パラメータ（半径進化のための2つ、凍結時間のための1つ、および温度／化学ポテンシャル）でこれらの結果を達成している。これは、スペクトルと異方的フローを同時に記述するために8つのパラメータを必要とする、より精緻な楕円ブラストウェーブ・パラメトリゼーションよりも少ない。

意義と主張
著者らは、彼らの研究が、複雑な重イオン衝突の観測量を記述するための、単純な球対称動的火球モデルの実行可能性を実証していると主張している。具体的には、以下のことを断言している：

1. 物理的に動機付けられた少数のパラメータを持つモデルは、横運動量スペクトルと、それに対応するラピディティ分布を、一貫した枠組みの中で同時に説明できる。
2. 非中央衝突の初期の幾何学的異方性にもかかわらず、方位角方向に積分されたスペクトルおよびラピディティ分布を記述するには、球対称性の仮定は十分である。
3. 抽出された凍結パラメータは、衝突の中心度およびエネルギーに対して物理的に直感的な系統的依存性を示しており、火球のバルクの進化と冷却を捉えるモデルの能力を検証している。
4. このアプローチは、フル・ハイドロダイナミクス・シミュレーションの代わりとなる、凍結特性を抽出するための計算効率の高い代替手段を提供する。これは、特にフル・ハイドロダイナミクスによる記述が困難な可能性があるビームエネルギー・スキャン（BES）のエネルギー領域において重要である。

論文は、展開する球状火球モデルが軽ハドロンに関するSTARのデータをうまく再現したことを結論付けており、このエネルギー範囲における媒質のバルク特性が、簡略化された球対称膨張シナリオによって効果的に捉えられるという仮説を支持している。