Unbiased Data-Driven Determination of the Nuclear Dipole Amplitude in the Color Glass Condensate

技術要約：カラーグラス・コンデンセートにおける核ディポール振幅の、偏りのないデータ駆動型決定

問題提起
量子色力学（QCD）の高エネルギー極限において、ビョルケン$x$が小さい領域でのパルトン密度の増大はグルオン飽和によって抑制され、カラーグラス・コンデンセート（CGC）領域へと至る。飽和スケール $Q_s(x)$ は原子核質量数 $A$ とともに増大するが（$Q_s^2 \propto A^{1/3}$）、核グルオン・ディポール振幅の定量的な抽出は、歴史的にパラメータ化された初期条件に依存してきた。これらの仮定は、核プロファイルや初期振幅の関数形式に関する制御不能なモデル依存性を導入する。その結果、抽出された核飽和スケールの値は文献間で大きく異なり、真の核効果を不明瞭にし、飽和物理学の精密プローブとしての核データの利用を妨げている。さらに、既存の陽子に関する抽出では、標準的なパラメータ化を用いた際に運動量空間における正値性の破れに苦慮してきた。

手法
著者らは、初期条件の関数形式を仮定することなく、$^{208}\text{Pb}$ のインパクトパラメータ平均化された核ディポール振幅 $N(r, x)$ を決定するための、物理情報ニューラルネットワーク（PINN）フレームワークを導入する。

• 物理情報アーキテクチャ: ディポール振幅は、ニューラルネットワーク $N_\theta(r, Y)$ （ここで $Y = \ln(1/x)$）によって表現される。このネットワークは、フィッティングループ中に方程式を反復的に解くのではなく、損失関数内で直接、高次補正を含むバリツキー・コフチェフ（ciBK）進化方程式を満たすように訓練される。これにより、抽出された振幅がQCD進化ダイナミクスを厳密に遵守することが保証される。
• 訓練制約: 訓練目的関数は、以下の要素からなる総損失関数 $L_{tot}$ を最小化する：
  • ciBK進化損失: 自動微分と微分可能な求積法を用いて、ciBK方程式の残差にペナルティを課す。
  • 観測量損失: 実験データ、具体的には前方ハドロン生成における核修飾因子 $R_{p\text{Pb}}$ および、超フェリカル（UPC）Pb-Pb衝突におけるコヒーレント $J/\psi$ 光生成断面積で制約をかける。本フレームワークは、実験的な不確かさと学習されたデータ依存の分散を考慮するために、ヘテロスセダスティック（不等分散）ガウス損失を利用する。
  • 物理的正則化: 大きなディポールサイズにおけるブラックディスク限界、運動量空間におけるフーリエ正値性、および小さな $r$ における異常次元の境界を強制する。
• 最適化: ネットワークは、異なる損失成分（進化、データ、および物理的制約）の重みをバランスさせるための勾配キャリブレーション手順を採用している。不確かさの定量化は、データ依存の不確かさとモデル依存の不確かさを分離したディープアンサンブルを介して行われる。
• データ入力: 解析には、$\sqrt{s_{NN}} = 5.02$ TeV におけるLHCbの $R_{p\text{Pb}}$ 測定値、およびALICE、CMS、ATLAS、LHCbからのコヒーレント $J/\psi$ 光生成データが使用される。

主要な結果

• Pb振幅の抽出: 本フレームワークは、開始スケール $x_0 = 0.01$ における $^{208}\text{Pb}$ ディポール振幅の抽出に成功した。進化した振幅は、利用可能な運動学的範囲全体にわたって、測定された $R_{p\text{Pb}}$ およびコヒーレント $J/\psi$ 断面積を再現する。
• 飽和スケールの比: 抽出されたインパクトパラメータ平均化された初期飽和スケールは、比 $Q_{s0,\text{Pb}}^2 / Q_{s0,p}^2 = 3.17^{+0.17}_{-0.10}$ を与える。この値は、単純な幾何学的スケーリングの推定値（$A S_{\perp,p} / S_{\perp,A}$）と一致しているものの、それよりも約20%高く、重い原子核におけるカラー電荷密度の高さが予想されることを反映している。
• 初期条件の形式: 陽子の場合、MV型（マクレラン・ヴェヌガポラン）のパラメータ化では不十分であることが判明したが、抽出されたPbの初期条件は、MV型の関数形式によって良好に記述される。これはデータ駆動型のフィットから自然に導かれ、大きな原子核におけるガウス源近似の妥当性を検証している。
• 運動量空間の正値性: 抽出された振幅は、進化の間を通じて運動量空間における正値性を維持しており、従来のパラメータ化されたフィットにおける既知の問題を解決している。
• 予測能: 抽出された同一の振幅を用いて、著者らは $pp$、$p\text{Pb}$、および$\text{Pb}p$衝突における横運動量比$\langle p_T \rangle / \langle p_T \rangle(\eta_0)$ のラピディティ依存性を予測している。これらの予測は、システム依存のパラメータを導入することなく、低マルチプリシティ領域における最新のLHCb測定値と一致している。高マルチプリシティにおける不一致は、初期状態のディポールフレームワークの範囲外である最終状態効果に起因するとされている。
• 汎用性: 著者らは、抽出された陽子およびPbのオパシティに基づき、中間的な原子核（例：$^{12}\text{C}$、$^{16}\text{O}$、$^{197}\text{Au}$）の初期条件を構築するための、物理的に動機付けられた補間スキームを提供しており、将来の小 $x$ フェノメノロジーのためのベンチマーク入力を提供する。

意義および主張
本論文は、飽和領域における核ディポール振幅の、初の偏りのないデータ駆動型の決定を提供すると主張している。非線形なciBK進化方程式を機械学習の訓練目的関数に直接組み込むことにより、本研究は、核の初期条件に関する既定の関数形式に関連する持続的なバイアスの源を除去している。

著者らは、この戦略が、QCDにおける動的に制約された観測量を抽出するための一般的な手法を確立するものであると主張している。結果は、幾何学的スケーリングの議論に対するモデルに依存しない定量的なテストを提供し、LHCおよびRHICにおける前方粒子生成、ならびに将来の電子イオン衝突器（EIC）での測定を解釈するための強固な基礎を提供する。本研究は、QCD進化と実験データの相互作用が、飽和領域における核構造を一意に決定できることを示しており、異なる核種間での飽和の予測を比較する際の歴史的な曖昧さを解消するものである。