Extraction of baryon number susceptibilities at finite density from heavy-ion collisions

技術要約：重イオン衝突におけるバリオン数感受率のベイズ抽出

問題提起
保存電荷、特にバリオン数の揺らぎは、量子色力学（QCD）の相構造、とりわけ予想される臨界点（CP）付近の性質を解明するための敏感なプローブとして機能する。第一原理に基づく格子QCD計算は、零バリオン化学ポテンシャル（$\mu_B = 0$）においてバリオン数感受率（$\chi_B^n$）を正確に決定できるが、有限密度におけるそれらの振る舞いは、外挿または代替的な理論的アプローチに依存している。実験的には、相対論的重イオン衝突器（RHIC）のビームエネルギー・スキャン（BES）プログラムにおいて、Au-Au衝突におけるネット・プロトン数のイベントごとの揺らぎが、$\sqrt{s_{NN}} = 7.7–27$ GeVの範囲で測定されている。しかし、これらのデータの定量的な解釈は、いくつかの非平衡および実験的効果によって複雑化している。具体的には、測定された陽子は全バリオン数の一部しか担っていないこと、実験的な受容領域が運動量空間において限定的であること、そして厳密な全域的バリオン数保存が、グランドカノニカルな期待値に対して非自明に揺らぎを修飾することなどが挙げられる。従来の戦略は、通常、基礎となる感受率を直接推論するのではなく、特定のモデル計算（例：ハドロン共鳴ガス、排除体積を持つ流体力学、あるいは輸送モデル）をデータに照らし合わせる手法をとってきた。

手法
本研究は、著者らの知る限り、重イオン衝突データからバリオン数感受率（$\chi_B^2, \chi_B^3, \chi_B^4$）を直接抽出する初のベイズ的抽出法を提示するものである。このフレームワークは、理論的な感受率から測定されたプロトンの高次モーメントへの順方向写像を、以下の4つの異なるステップを通じて構築する：

1. 流体背景（Hydrodynamic Background）： 3次元初期条件と一定のエネルギー密度（$\epsilon_{sw} = 0.26$ GeV/fm$^3$）での粒子化を用いたMUSICコードを使用して、イベント平均された流体的な超曲面を生成する。
2. 最大エントロピー凍結（MaxEnt）： 超曲面上の局所的なバリオン感受率を、結合バリオン・反バリオンの累積量へと変換する。本フレームワークでは、局所的な感受率の入力として、理想的なハドロン共鳴ガス（iHRG）のベースラインに対する比率 $\gamma = (\chi_B^2/\bar{\chi}_B^2, \chi_B^3/\chi_B^1, \chi_B^4/\chi_B^2)$ をパラメータ化する。MaxEnt処方箋は、流体場とその揺らぎを、ハドロン多重度の既約相対累積量へと写像する。
3. 受容性とフィルタリング： クーパー・フライ（Cooper-Frye）受容およびバリオン・プロトン・フィルタリングにより、独立した二項サンプリングを通じて、結合累積量を受理された陽子および反陽子へと写像し、STAR検出器の運動学的受容領域を考慮する。
4. 厳密な保存： サブエンサンブル受容法（SAM-3.0）を用いて厳密な全域的バリオン数保存を課し、グランドカノニカルな累積量を、実験データと比較可能なカノニカルな累積量へと変換する。

ベイズ推論は、BES-IIの範囲内（具体的には $\sqrt{s_{NN}} = 7.7, 9.2, 11.5, 14.6, 17.3, 19.6, 27$ GeV）の特定の衝突エネルギーにおける、0–5%中心度Au-Au衝突のSTARデータに対して行われる。補完的な2つの抽出が、(i) プロトン階乗累積量比（$\hat{C}_n/\hat{C}_1$）および (ii) ネット・プロトン累積量比（$\kappa_n/\kappa_1$）を用いて実施される。感受率の比に対して平坦な事前分布を適用し、統計的および系統的な不確かさを考慮したガウス型尤度を構築する。

主要な結果

• 二次の感受率 ($\chi_B^2$)： 解析により、iHRGの値で正規化された二次の感受率（$\chi_B^2/\bar{\chi}_B^2$）に対するタイトな制約が得られ、相対的な不確かさは5–10%であった。抽出された値は、測定されたエネルギー全域にわたって非単調な衝突エネルギー依存性を示した：最低エネルギー（$\sqrt{s_{NN}} = 7.7$ GeV）での増強（$\sim 1.06–1.07$）、中間エネルギー（$11.5–19.6$ GeV）での極小（$\sim 0.73–0.82$）、および高エネルギーでの平坦な抑制（$\sim 0.8–0.9$）。
• 格子QCDとの比較： 化学的凍結線にマッピングすると、抽出された $\chi_B^2$ の値は、$\mu_B \lesssim 300$ MeV において格子QCDに基づく推定値（具体的には4D-TExS状態方程式）と定量的に一致している。しかし、より大きな $\mu_B$（$\sqrt{s_{NN}} = 7.7$ GeV に対応）においては、抽出された値は格子QCDのバンドから上方へ逸脱し、最大 $\sim 50\%$ の増強を示している。
• 高次の感受率： 三次および四次の感受率（$\chi_B^3, \chi_B^4$）は、バリオンからプロトンへの写像における効率の損失と、バリオン保存による歪みのために、制約が弱い。四次の比については、$\sqrt{s_{NN}} = 19.6$ GeV において負の値への局所的な選好（$\chi_B^4/\chi_B^2 = -33 \pm 14$）が見られるものの、大きな不確かさのため、これらの高次モーメントのみからQCD相構造について決定的な結論を下すことはできない。
• 相関の支配性： 二バリオン間の既約相関のみを仮定する最小限のシナリオ（$\Delta\chi_B^3 = \Delta\chi_B^4 = 0$ とする）は、$\sqrt{s_{NN}} \approx 11$ GeV 付近におけるプロトン階乗累積量比 $\hat{C}_3/\hat{C}_1$ の非単調なピークをうまく再現している。これは、観測された構造が、三バリオンまたは四バリオンの既約相関を必要とするのではなく、$\chi_B^2$ のエネルギー依存性と厳密なバリオン数保存の相互作用から生じていることを示唆している。

意義と主張
本論文は、特定の状態方程式（EoS）の仮定を回避し、化学的凍結線に沿ったバリオン数感受率を決定する、初のデータ駆動型のベイズ的決定法を提供すると主張している。著者らは、本フレームワークが二次の感受率を効果的に分離することに成功し、低バリオン密度では格子QCDと一致するものの、高密度（$\mu_B \gtrsim 300$ MeV）においては、非臨界的な外挿と比較してバリオン数揺らぎが増強されている可能性を示していると断言している。

本研究は、観測された現象は二次の相関を超えた既約多体相関なしでも概ね記述可能であるが、これはそれらの存在を否定するものではなく、むしろ現在の測定は主に二バリオンの寄与に対して感度を持っているということを控えめに結論付けている。著者らは、低衝突エネルギーで見られる増強は将来の固定ターゲット実験による精査に値すると述べており、現在の抽出は固定された流体および凍結の設定に依存しており、これはグローバル解析における不確かさの下限値を示すものであるとしている。本フレームワークは、有限密度のQCD状態方程式を制約するための一般的なツールとして提示されており、流体および凍結のモデル化に関する不確かさがより制御可能になった後、固定ターゲットエネルギー、体積揺らぎの補正、および他の保存電荷へと拡張が可能である。