Hadronic rescattering effects on net-proton cumulants from functional renormalization group calculations

技術要約：ネット・プロトン・キュムラントに対するハドロン再散乱効果

問題提起
重イオン衝突のビーム・エネルギー・スキャン（BES）において測定されるネット・プロトン・キュムラントは、量子色力学（QCD）相図における予想される臨界点（CEP）を探索するための主要な観測量である。関数的繰り込み群（fRG）計算などの理論的枠組みは、化学的凍結（chemical freeze-out）における平衡状態のキュムラントを提供するが、その後のハドロン再散乱がこれらの観測量に与える影響については、依然として十分に定量化されていない。既存の研究は、初期状態や相転移の寄与に焦点を当てることが多く、運動学的凍結（kinetic freeze-out）に至るまでのハドロン・カスケードがいかにゆらぎを修飾するかについては軽視されている。さらに、厳密な保存則（例：正味のバリオン数）や火の玉（fireball）の有限な寿命は、非自明な背景効果をもたらし、臨界信号の抽出を困難にする。

手法
著者らは、ハドロン再散乱の効果を定量化するために、平衡状態のfRGキュムラントをハドロン輸送モデルへと伝播させる枠組みを構築している。その手法は主に以下の3つのステップで構成される。

1. 初期分布の再構成:

   • 入力: 本研究では、fRGを介してグランドカノニカルアンサンブル（GCE）で計算された正味のバリオン数の平衡キュムラントを利用する [10]。
   • マッピング: 正味のバリオン・キュムラントをネット・プロトン・キュムラント（$C_n^B \approx C_n^P$）として近似する。これは標準的な仮定であるが、衝突エネルギーが低くなるにつれて不確実性が増す。
   • 保存則: 有限系における厳密な正味バリオン数保存を考慮するため、サブアンサンブル受容法（SAM）を適用し、GCE感受率をカノニカルアンサンブル（CE）のキュムラントへと変換する。
   • 最大エントロピー原理: 利用可能なのは最初の4つのキュムラントのみであるため、バイアスを最小限に抑えるために最大エントロピー原理を用いてネット・プロトンの多重度に関する確率密度関数を再構成する。
   • イベントごとの生成: 各イベントに対して、プロトン（$N_p$）と反プロトン（$N_{\bar{p}}$）の多重度を個別に生成する。$N_p$をサンプリングし、$N_{\bar{p}} = N_p - N_{net}^p$ を通じて $N_{\bar{p}}$ を導出することで、$N_{\bar{p}} \geq 0$ を保証する。運動学は、半径 $R=5.89$ fm の球状体内に初期化された、径方向流（radial flow）を伴う熱分布からサンプリングされる。

2. ハドロン輸送シミュレーション:

   • 生成された粒子リストは、SMASH（Simulating Many Accelerated Strongly-interacting Hadrons）輸送モデル（バージョン3.2）を通じて伝播される。
   • システムは、散乱率が無視できる程度になるまで（運動学的凍結）、$t = 101$ fm/c まで進化する。
   • シミュレーションは、衝突エネルギー $\sqrt{s_{NN}} = 3.0, 3.9, 4.9, 7.2,$ および $7.7$ GeV をカバーする。
   • 球状の設定では実験的な受容角カット（例：$p_T$ やラピディティ）は適用せず、代わりに全運動量カット（$0.4 < p < 2$ GeV/c）を使用する。

3. 解析:

   • ハドロン進化後、ネット・プロトン・キュムラント（$C_1$ から $C_4$）およびその比（特に $C_4/C_2$）を再計算する。
   • 弾性散乱、共鳴、および非弾性チャネルの衝突率と重心エネルギー分布を解析し、キュムラント修飾の微視的な駆動要因を特定する。

主な結果

• 高次キュムラントの抑制: ハドロン再散乱は、初期の平衡信号を著しく修飾する。厳密な保存を課したカノニカルアンサンブルにおいて、$\sqrt{s_{NN}} = 4.9$ GeV での比 $C_4/C_2$ は、カスケードの初期段階で約 20% 抑制される。
• 非単調性の生存: $C_4/C_2$ の非単調なエネルギー依存性（fRG入力および想定されるCEPに関連するピーク構造）は、GCEおよびCEの両方の枠組みにおいて、ハドロン進化後も生存する。しかし、このピークの大きさは大幅に減少する。
• 時間進化ダイナミクス:
  • $C_4$ および $C_2$ の抑制は、ハドロンカスケードの初期の密度の高い段階（最初の数 fm/c）で最も顕著である。
  • 初期の減少に続き、システムが膨張・希薄化するにつれて部分的な回復が見られる。これは、ハドロン進化が単純な単調減衰メカニズムとして機能しているわけではないことを示している。
  • $C_4$ の抑制は $C_2$ よりも深刻であり、これが $C_4/C_2$ 比の純減少につながっている。
• 微視的メカニズム: 修飾は主に共鳴形成によって駆動されており、これが衝突率の大部分を占め、次いで弾性散乱が続く。この効果は 4.9 GeV で最も強く、初期のゆらぎが最大であることから、fRG入力から継承された高次のゆらぎは、ハドロン相互作用の影響を特に受けやすいことが示唆される。
• カノニカル vs グランドカノニカル: CEの結果はGCEの結果よりも強い抑制を示しており、中間エネルギー領域における厳密な正味バリオン数保存の重要な役割を浮き彫りにしている。

意義と主張
本論文は、ハドロン再散乱が、実験データからQCD臨界点信号を抽出する際に考慮すべき無視できない背景効果を構成すると主張している。具体的には：

• 臨界信号の大きさ（$C_4/C_2$ のピークの高さ）は、ハドロン段階によって著しく修飾（減少）されるが、定性的な非単調構造は維持される。
• 本研究は、直接的なfRG入力を用いることで（モデルによるパラメータ化ではなく）、BESエネルギー範囲におけるこれらの効果を初めて定量的に評価している。
• 結果は、CEPの位置は（生存する）ピーク位置にエンコードされている可能性があるものの、ピークの高さ（抑制される）については、正しく解釈するためにハドロン・アフターバーナー効果に対する注意深い補正が必要であることを示唆している。

限界と今後の課題
著者らは以下の系統的な不確実性を指摘している：

• バリオン・プロトン写像: $C_n^B \approx C_n^P$ という近似は、ハイパーオンの寄与が無視できなくなる低エネルギー領域において不確実性を導入する。
• 受容分率: SAMにおける受容分率 $\alpha$ の現象論的な固定には、さらなる専用の研究が必要である。
• 幾何学: 一様な球状膨張の使用は、凍結が表面から内側に向かって段階的に起こる現実的な流体動力学的進化を簡略化している。今後の課題は、この枠組みを、現実的なゆらぎを持つ初期状態やハイブリッド流体動力学シミュレーションと結合させることである。