Geometric bias and centrality dependence of jet quenching in high-energy nuclear collisions

以下は、提供された論文「Geometric bias and centrality dependence of jet quenching in high-energy nuclear collisions（高エネルギー核衝突におけるジェットクエンチングの幾何学的バイアスと中心性依存性）」の技術的な要約です。

1. 研究の背景と課題

• 背景: 相対論的重イオン衝突（RHIC, LHC）では、クォーク・グルーオン・プラズマ（QGP）が生成され、その存在の証拠として「ジェットクエンチング（高エネルギー粒子の減衰）」が知られている。
• 課題: 従来の理論モデル（標準的な MC-Glauber モデルなど）では、中心衝突（central collisions）のジェット減衰を説明するためにパラメータを調整すると、周縁衝突（peripheral collisions）において実験データ（核変換係数 $R_{AA}$）を過大評価してしまう傾向があった。
• 核心的な問題: 周縁衝突では、核子間の重なりが希薄であるため、QGP との相互作用が弱いはずだが、実験的には依然として $R_{AA} < 1$（減衰）が観測されている。この矛盾の原因として、QGP によるエネルギー損失だけでなく、**「幾何学的バイアス（Geometric Bias）」**と呼ばれる初期状態の効果が無視されている可能性が指摘されていた。

2. 提案された手法とモデル

本研究では、幾何学的バイアスとジェット-QGP 相互作用の両方を組み合わせて、$R_{AA}$ の中心性依存性を再構築する。

A. HIJING ベースの初期条件モデルの開発

• 標準モデルの限界: 従来の MC-Glauber モデルでは、核子 - 核子（NN）衝突が「硬いカットオフ（$b_{NN} < b_0$）」内で起こり、1 回の NN 非弾性衝突につき必ず 1 つのハード散乱（ジェット生成）が起こると仮定している。
• HIJING モデルの改良:
  • NN 衝突のインパクトパラメータ（$b_{NN}$）依存性を明示的に取り入れる。
  • NN 非弾性衝突の確率と、1 回の非弾性衝突あたりのハード部分子散乱の数が、$b_{NN}$ に依存して変化するモデルを構築（HIJING 生成器に基づき、核子形状因子や断面積を考慮）。
  • 周縁衝突（大きな $b_{NN}$）では、非弾性衝突が起こってもハード散乱を伴わない（ソフト散乱のみ）確率が高くなることを実装。
• 幾何学的バイアス因子 ($R^{bias}_{AA}$):
  • 定義: $R^{bias}_{AA} = \langle N^{hard}_{AA} \rangle / (\langle N_{coll} \rangle \langle N^{hard}_{NN} \rangle)$
  • 意味: 特定の中心性クラスにおける AA 衝突のハード散乱数と、偏りのない pp 衝突の期待値との比。周縁衝突では $R^{bias}_{AA} < 1$ となり、これが $R_{AA}$ の減少に寄与する。

B. Linear Boltzmann Transport (LBT) モデルの改良

• 問題点: 従来の LBT モデルでは、周縁衝突のように QGP が小さい場合、ジェットと QGP の相互作用が弱く、エネルギー損失がほとんどないにもかかわらず、カラーフローの取り扱い（「フェイク」部分子の扱い）により、物理的に不自然な $R_{AA} > 1$ の増大が生じる場合があった。
• 改良策:
  • 「フェイク」部分子への大きな運動量付与: 負の部分子（エネルギー枯渇を表す）をストリング接続する際、その運動量（特に $p_z$）を非常に大きく（$10^4$ GeV）設定する。これにより、pp 衝突におけるストリング配置を再現し、周縁衝突での $R_{AA} \approx 1$ という基底線を正しく回復させる。
  • 負の部分子のハドロン化: 負の部分子を QGP 内の熱的分布からサンプリングし、エネルギー保存則を満たすようにハドロン化を行うことで、$p_T \approx 10$ GeV 付近での過剰な減衰（負の値）を解消する。

3. 主要な結果

• 幾何学的バイアスの定量化:
  • HIJING モデルを用いた計算により、周縁衝突（70-90% 中心性）において幾何学的バイアス因子 $R^{bias}_{AA}$ が 1 よりも有意に小さい（約 0.7-0.8 程度）ことが確認された。
  • これは、周縁衝突では平均的な NN インパクトパラメータが大きいため、1 回の NN 衝突あたりのハード散乱確率が低下することに起因する。
• $R_{AA}$ の中心性依存性の説明:
  • 中央〜半周縁（0-50%）: 標準的な Glauber モデルと HIJING モデルの差は小さく、LBT モデルによるエネルギー損失の説明が支配的。
  • 周縁（50-90%）: 幾何学的バイアス因子を考慮することで、実験データ（CMS データ）との一致が劇的に改善された。
  • 特に 70-90% 中心性では、QGP によるエネルギー損失がほとんどなくても、幾何学的バイアスだけで観測されるようなジェット減衰が説明可能となった。
• モデルの統合:
  • 改良された初期条件モデル（HIJING ベース）と改良された LBT モデルを組み合わせることで、5.02 TeV の Pb+Pb 衝突において、0-10% から 90-100% までの全中心性範囲で、荷電ハドロンの $R_{AA}$ を統一的に記述することに成功した。

4. 結論と意義

• 結論: 高エネルギー核衝突、特に周縁衝突におけるジェット減衰の観測は、QGP によるエネルギー損失だけでなく、初期状態の幾何学的バイアス（核子重なりとハード散乱確率の関係）によって大きく支配されている。
• 科学的意義:
  • 従来の「周縁衝突での減衰＝QGP 効果」という解釈は誤りであり、実際には初期状態の幾何学的効果によるものであることを定量的に示した。
  • これにより、小さなシステム（小規模な核衝突や pA 衝突）における QGP 形成の閾値を正しく評価するための基盤が整備された。
  • 将来的には、冷たい核物質効果や中心性選別バイアスなどの追加的な依存性を考慮することで、さらに精度の高い QGP 特性の抽出が可能になる。

この研究は、ジェットクエンチングの理論的記述において、初期状態の幾何学的効果を体系的に組み込む必要性を強く示唆し、実験データと理論の整合性を高める重要なステップとなった。