Path-length dependence of parton energy loss across collision systems: a Bayesian analysis of charged-particle RAA, consistent with a universal exponent from O+O to Pb+Pb

技術要約：衝突系におけるパートンのエネルギー損失の経路長依存性

問題提起
クォーク・グルーオン・プラズマ（QGP）におけるパートンのエネルギー損失を支配するメカニズムは、平均エネルギー損失 $\langle \Delta E \rangle$ がイン・メディウムの経路長 $L$ に対してどのようにスケーリングするかという点に刻み込まれている。理論的枠組みは、異なるスケーリング挙動を予測している：衝突型エネルギー損失は線形（$\langle \Delta E \rangle \propto L$）、媒質誘起放射エネルギー損失は二次関数的（$\langle \Delta E \rangle \propto L^2$）、そして強結合（ホログラフィック）モデルは三次関数的（$\langle \Delta E \rangle \propto L^3$）なスケーリングを予測する。歴史的に、このスケーリング指数を抽出することは、単一の大型系（例：Pb+Pb）内での中心度変化を分析することによって困難であった。これは、中心度選択によるバイアスや、周辺事象におけるグローバーに基づく規格化（$T_{AA}$）を導入してしまうためである。LHCにおける軽イオン衝突データ（O+O, Ne+Ne）の登場は、補完的なレバーを提供している。これらは、周辺的な重イオン衝突に匹敵する経路長を最小バイアスで探索することで、システムサイズ依存性を孤立させることが可能である。

手法
著者らは、4つの衝突系（O+O, Ne+Ne, Xe+Xe, Pb+Pb、質量数 $A=16$ から $208$まで）にわたる荷電粒子核修飾因子（$R_{AA}$）のベイズ解析を行う。

1. フォワードモデル： パラメトリックモデルが、分画エネルギー損失 $\Delta p_T$ と $R_{AA}$ を関連付ける。モデルは $\Delta p_T \propto \rho G^n$ と仮定しており、ここで $\rho$ は相対的な媒質密度、$G$ は無次元の幾何学的比、$n$ は決定すべきシステムサイズ指数である。スペクトル指数 $a(p_T)$ は、$pp$衝突からのデータ駆動型の測定を用いて固定され、$x_T$ スケーリングが異なる衝突エネルギーに対して適用される。
2. 幾何学： 衝突幾何学は、2つの独立したグローバー実装（光学グローバーモデル、および硬い芯の排除を伴う原子核の位置をサンプリングするカスタム・モンテカルロ・グローバーモデル）を用いて計算される。主要な幾何学的プロキシには、$A^{1/3}$、$\langle N_{part} \rangle^{1/3}$、および平均脱出経路長 $\langle L \rangle$ が含まれる。
3. 統計的枠組み：
   • 尤度： ビン間の共分散行列を組み込んだ相関ガウス尤度が使用される。系統誤差は、完全な相関による過剰制約を避けるため、部分的な相関を考慮して相関長 $\xi = 4$ ビンを用いてモデリングされる。
   • 推論： パラメータ（$\kappa, n, \beta$）の事後分布は、アフィン不変アンサンブルサンプラー（EMCEE）を用いてサンプリングされる。
   • モデル選択： ネステッドサンプリング（DYNESTY）により、固定された指数（$n=1, 2, 3$）を持つモデルを比較するためにベイズ証拠（Bayesian evidence）を計算する。
   • 検証： フレームワークは、正規化フローを用いたニューラル事後推定器による「シミュレーションベース推論（SBI）」の柱を含み、シミュレーションベースのキャリブレーション（SBC）および一系統抜き出し（LOSO）交差検証によって検証されている。また、保持された系を予測するためにガウス過程エミュレータも訓練される。
4. 堅牢性： 解析は、異なる幾何学的プロキシ、共分散長、スペクトル基底、および事前分布を含む160のバリアントに対してテストされている。

主な結果

• 有効指数： 解析により、$\langle N_{part} \rangle^{1/3}$ 幾何学に対して、有効なシステムサイズ指数 $n = 1.78 \pm 0.15 \text{ (stat)} \pm 0.05 \text{ (syst)}$ が抽出された。
• モデル選択： ベイズ証拠は、衝突型（$n=1$）および強結合（$n=3$）の仮説よりも、放射型スケーリング（$n=2$）を決定的に支持している。対数証拠の差は $2\Delta \ln Z \approx -29$ （$n=1$ の場合）および $-48$（$n=3$ の場合）であり、いずれも「決定的（decisive）」な閾値を超えている。
• 微視的境界： エネルギー損失のゆらぎが存在する中での微視的指数と有効指数の間のマッピングを分析することにより、著者らは、ゆらぎは有効指数を微視的指数よりも低くすることしかできないことを確立した。したがって、この測定は、固定された幾何学において、微視的指数（$n_{micro} \ge 1.78$）の下限を課しており、純粋な衝突型エネルギー損失（$n_{micro}=1$）を事実上除外している。
• 普遍性： ベイズ因子テストにより、小規模系（O+O, Ne+Ne）と大規模系（Xe+Xe, Pb+Pb）の間でエネルギー損失のレジームに変化がある証拠は見つからず、全範囲にわたる普遍的な指数が支持された。
• O+Oの抑制： 本フレームワークは、O+O単独においても非ゼロのパートンエネルギー損失があることに対し、決定的なベイズ証拠（$2\Delta \ln Z \approx 110$）を提供し、最小の系における抑制の兆候を定量化している。
• 輸送係数： 抽出されたエネルギー損失の大きさは、有効なジェット輸送係数 $\hat{q}/T^3 \approx 2\text{--}5$ に対応しており、これは独立したJETSCAPEの決定値と一致している。

意義と主張
本論文は、パス長指数を抽出するための、統計的に較正された、システムサイズに基づく最初の手法を提供すると主張している。その主な貢献は以下の通りである：

1. 堅牢な有効指数： 特定の輸送実装に依存しない、幾何学的レベルの観測量としての $n \approx 1.78$ を確立し、将来のマルチステージ・エネルギー損失計算のベンチマークとして機能させる。
2. 衝突メカニズムの排除： ゆらぎの境界を通じて、純粋な衝突型エネルギー損失は、ゆらぎを考慮しても観測されたシステムサイズ・トレンドと矛盾することを実証している。
3. 手法の厳密性： 相関のある不確実性を伝播させ、閉鎖テストを通じて区間被覆を検証し、結果のクロスチェックのためにシミュレーションベース推論を採用した包括的なベイズフレームワークを導入し、抽出の統計的防御力を確保している。
4. 予測能力： 本フレームワークは、将来の軽イオン測定（Ar+Ar および Kr+Kr）に対する検証可能な予測を与え、$R_{AA}(10 \text{ GeV})$ がそれぞれ $0.55 \pm 0.03$ および $0.41 \pm 0.03$ であると予測している。

限界と範囲
著者らは、この結果が直接的な微視的パス長依存性の測定ではなく、あくまで「有効な」指数であることを明示している。解析は、最小バイアス・システムサイズ・スキャンにおける媒質密度とパス長間の縮退（両者が $\langle N_{part} \rangle^{1/3}$ としてスケールするため）によって制限されており、抽出された $n$ は両方の効果を組み合わせたものである。本論文は、中心度微分データや方位角異方性（$v_2$）の測定のような追加の入力なしには、微視的なメカニズム（例：BDMPS-Z と有限不透明度補正の区別）を解決するものではないと主張している。結論は、解析された4つの系および採用された特定のフォワードモデルの仮定に厳密に限定される。