Out-of-equilibrium contributions to charm hadrons in a fluid-dynamic approach

技術要約：流体動力学的アプローチによるチャームハドロンの非平衡寄与

問題提起
チャームクォークはその大きな質量により、初期の硬い散乱（hard scattering）で主に生成され、ファイアボールの寿命内に熱的生成を通じて熱平衡に達することはないため、クォーク・グルーオン・プラズマ（QGP）に対するユニークなプローブとして機能する。従来の流体動力学的枠組みは、チャームクォークのダイナミクスを局所熱平衡に近いと仮定することで、チャームクォークの挙動をモデル化することに成功してきた。しかし、これらのモデルは、衝突の初期段階（自由ストリーミング相の後）と凍結（freeze-out）超曲面の両方における重要な段階での非平衡補正を無視していた。凍結補正の欠如は、チャームハドロンの多重度を空間拡散係数（$D_s$）に対して非物理的な依存性を持たせ、チャーム輸送係数の系統的な決定を妨げていた。また、初期分布はしばしば局所熱平衡にあると仮定されており、初期のチャーム生成の非平衡な性質が無視されていた。

手法
著者らは、初期段階と凍結表面の両方において非平衡補正を明示的に計算することにより、以前に開発されたチャームクォークの流体動力学的枠組みを拡張している。

1. 初期段階の初期化:

   • 本モデルは、チャーム生成時刻と流体力学的初期化時刻（$\tau_0$）の間の自由ストリーミング相を導入する。
   • 特性曲線法を用いて、このフェーズにおいて直線軌道で移動する粒子のキネティック方程式を解析的に解く。
   • このアプローチは、初期のチャーム密度を修正し、$\tau_0$ において非ゼロの拡散電流（$\nu^\mu_q$）を生成する。これは、二体衝突密度（$n_{coll}$）の空間依存性と、初期の硬い散乱の運動量分布から導かれる。
   • 全チャーム分布関数の積分を微分生成断面積に関連付けるために、ランダウ・マッチング条件が用いられる。

2. 凍結補正（マルチ種アプローチ）:

   • 著者らは、チャーム・ハドロン・レゾナンス・ガス（HRGC）内でのマルチ種アプローチを採用し、系を、すべてのチャーム・ハドロン種が軽いQGP成分と相互作用する流体として扱う。
   • 各ハドロン種 $i$ の分布関数に対する非平衡補正（$\delta f_{i,p}$）は、モーメント法とナビエ・ストークス近似を用いて導出される。
   • 補正は既約モーメントに展開され、チャーム・ハドロンからの体積圧力および剪断応力の寄与は軽い成分と比較して無視できると仮定して、ランク $l=1$（拡散電流）で打ち切られる。
   • 主要な解析結果（式23）は、補正を次のように表す：
     $$\delta f_{i,p} = f^{(0)}_{i,p} \frac{1}{P_i} \bar{\kappa}_i \nu^\mu_q p_{\langle \mu \rangle}$$
     ここで、$\bar{\kappa}_i$ は種の電荷と密度に依存する正規化された拡散係数であり、$\nu^\mu_q$ は全チャーム拡散電流である。この定式化により、個々の種の電流ではなく、全電流のみの発展が可能となる。

3. 発展とスペクトル:

   • 系は、一定の $D_s T$ を用いた 1+1D 流体動力学フレームワーク（FluiduM）を用いて発展させる。
   • 不変収量および横運動量（$p_T$）スペクトルは、導出された非平衡補正を組み込んだクーパー・フライ（Cooper-Frye）の手順を通じて計算される。
   • モデルは、FONLL pQCD 予測と ALICE 測定のチャーム断面積の両方を入力として使用し、ALICE の実験データ（$\sqrt{s_{NN}} = 5.02$ TeV、0–10% セントラリティの Pb-Pb 衝突）に対してテストされる。

主要な貢献と結果

• 非物理的な依存性の解消: 凍結時における非平衡補正の導入により、チャームハドロンの多重度の空間拡散係数 $D_os$ に対する非物理的な依存性が解消された。これらの補正がない場合、多重度は $D_s$ によって大きく変動したが、補正を加えることで、統合された収量は、テストされた全範囲の $D_s$ において、非拡散ケースおよび初期の総チャームクォーク数と一致した。
• 初期拡散電流: 自由ストリーミングによる初期化は、$\tau_0$ において、後期に生成される電流よりも1桁大きい拡散電流を生成する。しかし、本研究では、この初期条件が、ゼロ初期電流のシナリオと比較して、その後のチャーム場の発展および最終的な不変収量に与える影響は無視できる程度であることを見出した。
• モデルとデータの比較:
  • FONLL 断面積を用いることで、モデルは不確かさの範囲内でチャーム中間子の統合収量を再現する。
  • $\Lambda_c^+$ バリオンの収量については、過小評価されているというテンションが観察される。先行文献に従い、著者らは励起チャームバリオンの統計的重みを3倍にすることで「チャーム・バリオン増強」を適用する。この調整を行い、ALICE 測定の断面積を使用することで、モデルは $p_T \approx 3$ GeV/c までの実験的な収量および $p_T$ スペクトルとの適合性を達成する。
• 適用限界: 本論文は、分布関数の正値性に基づく流体動力学的アプローチの根本的な限界を特定している。
  • 非平衡補正 $\delta f_{i,p}$ は（拡散電流の負のラジアル・プロファイルにより）負になる可能性がある。
  • $|\delta f_{i,p}| \approx f^{(0)}_{i,p}$ となると、全分布関数が負となり、物理的な解釈が無効となる。
  • この $p_T$ スペクトルの「ゼロ交差」は、$D_s T$ が大きくなるにつれてより低い $p_T$ 値で発生する。$D_s T \gtrsim 0.2$ の場合、モデルは $p_T \sim 2$ GeV/c で負のスペクトルを与える。
  • したがって、本フレームワークは、全 $p_T$ 範囲において $D_s T \lesssim 0.1$ である場合、あるいは高い $D_s T$ に対しては制限された低 $p_T$ 領域（例：$D_s T \approx 0.2$ の場合は $p_T \lesssim 2.5$ GeV/c）においてのみ信頼できる。

意義
本論文は、衝突の初期段階と最終段階の両方における非平衡効果を考慮した、チャームハドロンの一貫した流体動力学的フレームワークを確立している。輸送係数に対する収量の非物理的な依存性を補正することで、本研究は、チャーム拡散係数の系統的な抽出（例：ベイズ解析による）のための必要な基礎を築いた。しかし、著者らは、現在の線形化アプローチが、高い拡散係数や横運動量における摂動論的仮定（$|\delta f| \ll f^{(0)}$）の崩壊によって制限されていることを認め、控えめに結論付けている。彼らは、全運動学的範囲において正の分布関数を保証するために、最大エントロピー法のような代替手法を探索する必要があると示唆している。