CGC-py: A Monte Carlo Event Generator for Gluon Saturation Physics

技術要約：「CGC-py：グルオン飽和物理のためのモンテカルロ・イベント生成器」

問題提起
高エネルギーかつ小さなビョルケン-$x$において、線形QCD進化はグルオン密度の無制限な増大を予測するが、これはユニタリティを破綻させる。カラーグラス凝縮（CGC）有効場理論は、動的な飽和スケール $Q_s(x)$ によって特徴付けられる、結果として生じる非線形飽和領域を記述する。深非弾性散乱（DIS）における単一ハドロン、二ハドロン、およびダイジェット生成は、この領域の確立されたプローブであるが、既存の理論的ツールには限界がある。主要次数（LO）のCGC計算だけでは、ソフトおよびコリニアな放射がサチュレーションによる広がりと競合するスダコフ・ログ（Sudakov logarithms）を生成するため、方位角相関を定量的に記述することができない。逆に、既存のイベント生成器は、ハード過程をバック・トゥ・バックの相関限界付近に限定している（例：MCDijet）、あるいは完全なBK/JIMWLK階層ではなくGLRのような簡略化された分岐方程式に依存している（例：Cascade3）、あるいはCGCの完全なワイルソン線相関構造を欠いている。その結果、フルな小-$x$ CGC断面積を現実的なパートンシャワーおよびハドロン化に結合し、単一包含および二ハドロン観測量の両方に対して一貫した予測を提供できる生成器が必要とされている。

手法
著者らは、DIS過程（$\gamma^* p(A) \to q\bar{q} + X$）のために設計されたモンテカルロ・イベント生成器である CGC-py を提示している。そのワークフローは、3つの異なるコンポーネントを統合している：

1. フルCGC断面積： この生成器は、バック・トゥ・バックの相関限界（$q_\perp \ll P_\perp$）を仮定することなく、CGCの枠組みにおけるクォーク・反クォーク生成断面積を計算する。断面積を、ターゲット弾性的（コヒーレント）寄与とターゲット非弾性的寄与に分離する。計算にはダイポールおよび四重極相関を用い、四重極については大きな$N_c$ガウス近似を採用している。非弾性的積分の数値評価を加速させるため、著者らはフィロン型アルゴリズムを用いて機械学習エミュレータの訓練データを生成している。プロトンおよび金原子核のダイポール振幅は、HERAおよびLHCのデータによって制約された、以前の物理情報ニューラルネットワーク（PINN）によるグローバルフィットから導出されている。
2. 初期状態放射（ISR）： $q\bar{q}$ ペアは、パートン分岐（PB）横運動量依存（TMD）後方初期状態シャワーを受ける。このシャワーは、ハードスケールからカットオフ $Q_0$（2.5 GeVに設定）まで進化し、小-$x$ ターゲットからのグルオン放出をシミュレートする。これらの放出による反跳は、不変質量保存処方を用いて $q\bar{q}$ ペアに吸収され、システムが最終状態生成器に渡される前に運動量保存を保証する。
3. 最終状態放射（FSR）およびハドロン化： 反跳したクォークと反クォークは、最終状態放射およびラント・ストリング・フラグメンテーションのために Pythia 8 に供給される。独立したフラグメンテーションのためのカラー中性を維持するために、ゼロ運動量を持つソフトなコンパニオン・パートンが各ハード・パートンに付着される。

主な貢献

• 統一されたフレームワーク： バック・トゥ・バックの限界に制限されている従来の生成器とは異なり、CGC-pyはフルなCGC断面積の全運動学を評価する。これにより、単一包含ハドロン（またはジェット）スペクトルと二ハドロン相関を同一のイベントサンプルから同時に生成することが可能となり、一貫した正規化を確保し、個別の計算を回避できる。
• 検証： 生成器は、観測されない反クォークに関するフルな $q\bar{q}$ サンプルのモンテカルロ積分を、解析的な単一包含クォークスペクトルと比較するという厳格な「クロージャ・テスト」を受ける。結果は、数桁のオーダーにわたって $\sim 5\%$ 以内の一致を示している。
• 現象論的ベンチマーク： 著者らは、HERAエネルギーにおける $ep$ 衝突における荷電ハドロンスペクトルに対するH1実験データを用いて、生成器を検証している。

結果

• 単一包含スペクトル： CGC-pyは、Pythia 8のパラメータを再調整することなく、横運動量（$p_T^*$）において5桁近い範囲にわたってH1の $ep$衝突における荷電ハドロンデータを再現する。$eAu$衝突については、モデルはより平坦なスペクトルと、中間および高$p_T^*$ ビンにおける系統的な増強を予測する。これは核的な広がりによるものである。
• 核修飾因子（$R_{eAu}^h$）： 予測される核修飾因子は、期待される飽和パターンを示す：低 $p_T^*$ での抑制（$\sim 0.4$ に達する）と、高 $p_T^*$ での 1 への上昇である。これは、より大きな核飽和スケールによる、より高い横運動量への強度の再分配を反映している。
• 二ハドロン方位角相関： （DGLAP進化を含むが小-$x$ 飽和を含まない）Pythia 6ベースラインとの比較により、二ハドロン相関のアウェイサイド（away-side）の幅（$\sigma_{RMS}$）は、両モデルにおいて衝突エネルギーとともに増大することが明らかになった。
  • 低エネルギー（$E_p = 250$ GeV および $1000$ GeV）では、CGC-pyとPythia 6の結果は密接に一致している。
  • 考慮された最高エネルギー（$E_p = 25000$ GeV、$x_g \sim 4 \times 10^{-4}$ に相当）において、顕著な差が生じ、CGC-pyはより広いアウェイサイド・ピークを示す。
  • $x_g$ 再スケーリング研究（高エネルギーサンプルの平均 $x_g$ を低エネルギーのものと一致させる強制的な操作）は、幅の増大が、コリニアなDGLAP進化と小-$x$ CGC進化の両方によって駆動されていることを示している。真の飽和駆動による広がりは、最高エネルギーにおいてのみ実質的となる。
• $eAu$ における核効果： CGC-py内では、$eAu$ 衝突は $ep$ と比較して増大したアウェイサイド幅と抑制されたバック・トゥ・バック収率を示す。しかし、これらの核効果は最高エネルギーのEICキネマティクスにおいても控えめであり、二ハドロン相関における飽和信号は微妙であることを示唆している。

意義および主張
本論文は、CGC-pyが、フルな小-$x$ CGC断面積を現実的なパートンシャワーに結合した、単一包含および二ハドロン観測量の両方を一貫して記述できる最初のイベント生成器であることを主張している。その主な意義は、以下の点にある：

1. レジューメーション・スキームと飽和現象論を、排他的なハドロンレベルの観測量に対してテストするための実験室を提供すること。
2. 次世代の電子・イオン衝突器（EIC）の $ep$および$eA$ プログラムに対する、現実的な擬似データおよび検出器レベルの予測を提供すること。
3. 飽和効果は存在するものの、EICのキネマティクスにおける二ハドロン相関におけるそれは控えめであることを示し、グルオン飽和への感度を最大化するために、最も前方（forward）でアクセス可能な運動学や補完的な観測量（ダイジェット、光子・ハドロン相関、重いフレーバー生成など）を使用することを動機付けること。

著者らは、このフレームワークが、ターゲット残骸の組み込み、次数の拡張（NLO）、およびエネルギー・エネルギー相関体（EEC）のような他の観測量への適用を含む、将来の研究への方向性を切り開くと結論づけている。