NRQCD および NRQED における S 波クォーコニウムおよびレプトニウム生成の自動化イベント生成

基本情報

• 論文ID: 2510.26773
• タイトル: Automated event generation for S-wave quarkonium and leptonium production in NRQCD and NRQED
• 著者: Alice Colpani Serri, Chris A. Flett, Jean-Philippe Lansberg, Olivier Mattelaer, Hua-Sheng Shao, Lukas Simon
• 分類: hep-ph（高エネルギー物理-現象論）、hep-ex（高エネルギー物理-実験）、nucl-ex（核物理-実験）、nucl-th（核物理-理論）
• 提出日時: 2025年10月30日
• 論文リンク: https://arxiv.org/abs/2510.26773

要旨

本論文は、MadGraph5_aMC@NLO フレームワークの拡張を提案し、非相対論的 QCD（NRQCD）および非相対論的 QED（NRQED）因子化形式の下で、S 波クォーコニウム（quarkonium）およびレプトニウム（leptonium）生成の領頭次（leading-order）断面積を自動計算できるようにしたものである。本フレームワークは複数のベンチマーク過程を通じて検証され、その堅牢性と現象論的研究における柔軟性が実証されている。本実装の主な利点は、既存の MadGraph5_aMC@NLO 機能とのシームレスな統合にあり、標準模型計算だけでなく、修正された汎用ファインマン出力（UFO）インターフェースを通じた様々な標準模型を超えた（BSM）シナリオや有効場理論（EFT）シナリオにも対応できる。さらに、本フレームワークは標準的なモンテカルロイベント生成器のパートン シャワー（parton shower）および強子化（hadronization）機能との互換性を保持している。多くの例を通じて、著者らはクォーコニウム過程の理論的研究には以下の点の慎重な考慮が必要であることを強調している：次領頭次（NLO）寄与の影響は、結合定数階数および速度スケーリング規則に基づく単純な計数論証のみでは予測困難なことが多い。

研究背景と動機

研究課題

本研究は、特に NRQCD および NRQED 理論フレームワークの下でのクォーコニウムおよびレプトニウム生成過程の自動化イベント生成問題に対処することを目的としている。

問題の重要性

1. 物理的意義：クォーコニウム（J/ψ、Υ など）は QCD の摂動的および非摂動的側面を研究するための重要なプローブである。50年前の J/ψ 発見による「11月革命」以来、クォーコニウム生成メカニズムは依然として完全には理解されていない。
2. 実験的要求：LHC、EIC などの衝突型加速器は膨大なクォーコニウムデータを生成し、比較のための信頼できる理論予測ツールが必要である。
3. 新物理探索：レプトニウム（ポジトロニウム、真ミューオニウム、ダイタウオニウム）は QED の精密検証、新物理の探索、および基本定数の測定に利用できる。

既存手法の限界

1. MadOnia：MadGraph/MadEvent v4 に基づき、単一の S 波または P 波クォーコニウムのみをサポートし、最新の MG5_aMC フレームワークへの移植がなされていない。
2. HELAC-Onia：複数の S 波および P 波クォーコニウムをサポートするが、物理的束縛状態を扱う際に短距離係数の手動結合が必要である。
3. 汎用 MC 生成器：Pythia、Herwig、Sherpa におけるクォーコニウム行列要素およびパートン シャワー機能は完全には互換性がない。
4. レプトニウムツールの欠落：公開利用可能なレプトニウム生成理論ツールは極めて限定的である。

研究動機

以下を可能にする包括的で自動化された、ユーザーフレンドリーなツールの開発：

• NRQCD/NRQED フレームワークの下で任意数の S 波クォーコニウムおよびレプトニウムを処理
• 広く使用されている MG5_aMC フレームワークへのシームレスな統合
• 標準模型および超標準模型シナリオのサポート
• 将来の NLO 精度計算の基礎確立

核心的貢献

1. 初の実装：MG5_aMC フレームワークにおける S 波クォーコニウムおよびレプトニウムの自動化イベント生成の初の実装
2. 統一フレームワーク：クォーコニウムおよびレプトニウムを扱うための統一された NRQCD/NRQED 実装の提供
3. 新しい UFO モデル：Fock 状態内容を定義する専用の boundstates.py ファイルを含む sm_onia モデルの開発
4. 広範な検証：HELAC-Onia および解析表現との詳細な比較を通じた実装の正確性の検証
5. 豊富な応用：pp、e⁺e⁻、e⁻p 衝突型加速器における様々な生成過程の実証：
   • 単一クォーコニウム生成およびその電弱ボソン、ジェット との関連生成
   • 二重/三重クォーコニウム生成
   • レプトニウム（ポジトロニウム、真ミューオニウム、ダイタウオニウム）生成
6. 理論的洞察：多くの例を通じて、次領頭次寄与の重要性は速度スケーリング規則のみでは推定不可能であることを実証

方法の詳細

理論フレームワーク

NRQCD/NRQED 因子化

束縛状態 B の微分断面積は以下のように表現できる： $d\sigma(N_1 + N_2 \to B + X) = \sum_{I_1,I_2,n} \int dx_1 dx_2 f_{I_1/N_1}(x_1) f_{I_2/N_2}(x_2) d\hat{\sigma}(I_1I_2 \to (C_1C_2)[n] + X_p) \langle \mathcal{O}^B_n \rangle$

ここで：

• $(C_1C_2)$：重クォーク対 $(Q\bar{Q}')$ または軽子対 $(\ell^-\ell'^+)$
• $n = {}^{2S+1}L_J^{[C]}$：量子状態（スピン S、軌道角運動量 L、全角運動量 J、色状態 C）
• $d\hat{\sigma}$：短距離断面積（摂動 QCD/QED で計算可能）
• $\langle \mathcal{O}^B_n \rangle$：長距離行列要素（LDME）

投影法

共変投影法を通じて振幅を特定の Fock 状態に投影：

1. 色投影（クォーコニウムのみ）：
   • 色一重項：$P_{C=1} = \delta_{c_3c_4}/\sqrt{N_c}$
   • 色八重項：$P_{C=8} = \sqrt{2} t^{c_{34}}_{c_4c_3}$
2. スピン投影： $P_S = \frac{\bar{v}_{\lambda_{C_2}}(k_4) \Gamma_S u_{\lambda_{C_1}}(k_3)}{2\sqrt{2m_{C_1}m_{C_2}}}$ ここで $\Gamma_{S=0} = \gamma_5$（一重項）、$\Gamma_{S=1} = \not{\varepsilon}^*_{\lambda_s}(K)$（三重項）
3. 軌道角運動量投影（S 波：L=0）： $\mathcal{A}^{(n,0)}_{\{[C],S,L=0\}}(r) = \left[\mathcal{A}^{(n,0)}_{\{[C],S\}}(r)\right]_{q=0}$

実装アーキテクチャ

UFO モデル拡張

1. 新しい Boundstate クラス：Particle クラスと同様に、束縛状態の Fock 状態を定義
2. 命名規則：Jpsi(1|3S11) は J/ψ の主量子数 N=1、スピン三重項 ³S₁、色一重項を表す
3. PDG エンコーディング：
   • 色一重項：標準 PDG エンコーディングを使用（例：J/ψ: 443）
   • 色八重項：非標準エンコーディング 99nqnsnrnLnJ
   • レプトニウム：新しいエンコーディング ±990nℓ₁nℓ₂nrnLnJ

コード生成フロー

1. 色投影：Python レイヤーでファインマン図生成時に適用
2. スピン投影：ALOHA モジュールを通じて専用サブルーチンを生成、HELAS ライブラリに依存
3. LDME 処理：振幅二乗計算に直接統合、onia_card.dat を通じて入力
4. 位相空間積分：疑似マルチチャネル戦略を採用、VEGAS インテグレーターで並列化

実行モード

1. Leading-order モード：完全な LO イベントシミュレーション、LHEF 形式でイベント生成
2. スタンドアロン モード：行列要素計算コード（Fortran）のみ生成、外部フレームワーク統合用

技術的革新点

1. 自動 Fock 状態結合：ユーザーが物理的束縛状態（例：Jpsi）を指定すると、コードが自動的にすべての関連 Fock 状態寄与を含める
2. 柔軟な LDME 設定：実行時に LDME 値を調整可能、再コンパイル不要
3. BSM 互換性：UFO インターフェースを通じて任意の新物理モデル（例：HEFT 例）をサポート
4. パートン シャワー インターフェース：Pythia などの MC 生成器とのシームレスな対接
5. レプトニウム サポート：レプトニウムの自動化計算能力を初めて提供

実験設定

計算設定

基本パラメータ

• クォーク質量：$m_c = 1.55$ GeV（3FS）、$m_b = 4.7$ GeV、$m_t = 173$ GeV
• 軽子質量：$m_e = 0.511$ MeV、$m_\mu = 105.66$ MeV、$m_\tau = 1.777$ GeV
• 電弱パラメータ：$m_Z = 91.188$ GeV、$m_H = 125$ GeV、$\alpha^{-1}_{G_\mu} = 132.507$
• 繰り込み/因子化スケール：$\mu_R = \mu_F = H_T/2 = \frac{1}{2}\sum_i \sqrt{k^2_{i,T} + m^2_i}$

衝突型加速器設定

1. 陽子-陽子衝突（LHC）：$\sqrt{s} = 13$ TeV、PDF4LHC21_40 NNLO PDF
2. 電子-陽電子衝突（Belle/BaBar）：$\sqrt{s} = 10.58$ GeV（Υ(4S) 共鳴）
3. 電子-陽子衝突（EIC）：$\sqrt{s} = 140.7$ GeV（$E_{e^-} = 18$ GeV、$E_p = 275$ GeV）

LDME 値

表 2 のデフォルト値を採用：

• 色一重項：波動関数 $R_Q(0)$ から計算、例えば $\langle \mathcal{O}^{J/\psi}_{{}^3S_1^{[1]}} \rangle = 1.16$ GeV³
• 色八重項：実験データからのフィッティング、例えば $\langle \mathcal{O}^{J/\psi}_{{}^1S_0^{[8]}} \rangle = 0.0146$ GeV³

検証戦略

行列要素レベルの検証

• 42 の異なる過程の振幅二乗を HELAC-Onia と比較
• 相対偏差 $\Delta_{\text{rel}} < 10^{-6}$（浮動小数点精度レベル）

断面積レベルの検証

• 全断面積および微分断面積を HELAC-Onia と比較
• プル値 $< 3$、統計的期待値と一致

解析表現との検証

• レプトニウム過程を文献の解析公式と比較
• Z および H ボソンのデカップリング極限で千分の一精度で一致

実験結果

主要結果の提示

1. 単一クォーコニウム生成（pp 衝突、√s = 13 TeV）

2→1 過程の全断面積：

• $\sigma(pp \to \eta_c + X) = 2.9366(5)$ μb
• $\sigma(pp \to J/\psi + X) = 536.14(6)$ nb
• $\sigma(pp \to \eta_b + X) = 5.4935(7)$ μb
• $\sigma(pp \to \Upsilon + X) = 6.0655(4)$ nb

主要な発見：

• 疑似スカラー状態の断面積はベクトル状態よりはるかに大きい。これはベクトル状態の 2→1 QCD 過程が C パリティ保存により禁止されるためである
• $\eta_c$ 断面積は低 x 領域の PDF に極めて敏感であり、異なる PDF セットは $\eta_b$ に対する $\eta_c$ 断面積の大きさの反転をもたらす可能性がある

2. 関連生成（pp 衝突）

電弱ボソンとの関連：

• $\sigma(pp \to J/\psi + Z) = 1.3425(4)$ pb（色八重項が主導）
• $\sigma(pp \to \Upsilon + Z) = 340.7(1)$ fb（色一重項が主導）

ジェット/光子との関連（運動学的カット適用）：

• $\sigma(pp \to J/\psi + j) = 329.8(2)$ nb（$p_{T,j} > 10$ GeV）
• $\sigma(pp \to J/\psi + \gamma) = 19.13(1)$ nb（$p_{T,\gamma} > 2$ GeV）

3. 二重クォーコニウム生成

同種クォーコニウム対：

• $\sigma(pp \to J/\psi + J/\psi) = 10.756(3)$ nb（色一重項が主導）
• $\sigma(pp \to \Upsilon + \Upsilon) = 44.63(1)$ pb

混合生成：

• $\sigma(pp \to J/\psi + \Upsilon) = 10.503(3)$ pb（色一重項 QCD 過程は $\mathcal{O}(\alpha^4_s)$ で消失）

$B_c$ 対生成：

• $\sigma(pp \to B^+_c + B^-_c) = 472.5(1)$ pb
• $\sigma(pp \to B^{\ast+}_c + B^{\ast-}_c) = 940.9(2)$ pb

4. 三重 J/ψ 生成

$\sigma(pp \to J/\psi + J/\psi + J/\psi) = 1.038(3)$ pb

三色一重項寄与はわずか 54.57 fb のみである。これは $gg \to 3J/\psi[^3S_1^{[1]}]$ が C パリティ保存により禁止されるためであり、主要寄与は $q\bar{q}$ 初期状態から来ている。

レプトニウム生成結果

ポジトロニウム（Ps）

pp 衝突（$\sqrt{s} = 13$ TeV）：

• $\sigma(pp \to \text{Ps}_1 + j) = 9.483(2)$ fb

e⁺e⁻ 衝突（$\sqrt{s} = 10.58$ GeV）：

• $\sigma(e^+e^- \to \text{Ps}_0 + \gamma) = 3.300(1)$ ab
• $\sigma(e^+e^- \to \text{Ps}_1 + \gamma) = 8.420(2)$ ab

真ミューオニウム（TM）

pp 衝突：

• $\sigma(pp \to \text{TM}_1 + j) = 9.460(2)$ fb（ポジトロニウムとほぼ同一）

e⁺e⁻ 衝突：

• $\sigma(e^+e^- \to \text{TM}_0 + \gamma) = 438.16(1)$ yb（極めて小さい、$4m^2_\mu/s$ による抑制のため）
• $\sigma(e^+e^- \to \text{TM}_1 + \gamma) = 8.422(2)$ ab

ダイタウオニウム（T）

pp 衝突：

• $\sigma(pp \to \text{T}_1) = 9.1906(4)$ fb（直接生成、初の信頼できる計算）
• $\sigma(pp \to \text{T}_1 + j) = 6.055(1)$ fb
• $\sigma(pp \xrightarrow{\gamma\gamma} p\text{T}_0p) = 108.055(7)$ ab（UPC 過程）

e⁺e⁻ 衝突：

• $\sigma(e^+e^- \to \text{T}_0 + \gamma) = 110.001(3)$ zb
• $\sigma(e^+e^- \to \text{T}_1 + \gamma) = 9.576(2)$ ab

特殊な場合：J/ψ + H 生成（HEFT フレームワーク）

ヒッグス有効場理論において、有効 $Hgg$ 頂点を通じて：

• $\sigma(pp \to J/\psi + H) = 1.53^{+0.40}_{-0.29}$ fb（色八重項寄与のみ）
• $\sigma(pp \to \Upsilon + H) = 91^{+23}_{-17}$ ab

ラピディティ分布の特性：

• ${}^1S_0^{[8]}$ 状態は中心ラピディティ領域に凹陥を示す（$|y| \approx 4$ 付近にピーク）
• ${}^3S_1^{[8]}$ 状態は中心ラピディティ領域にピークを示す
• 凹陥は有効 $Hgg$ 頂点図に由来し、ユカワ結合図により部分的に相殺される

パートン シャワー効果

$pp \to J/\psi + c\bar{c}$ 過程について、固定次数（fLO）と パートン シャワー マッチング（LO+PS）を比較：

• ${}^3S_1^{[1]}$ および ${}^1S_0^{[8]}$：PS が $p_T$ 分布を拡大
• ${}^3S_1^{[8]}$：PS が $p_T$ 分布を狭小化
• 総合効果：部分的な相殺、適度な歪み

アブレーション分析：Fock 状態寄与

$pp \to \eta_c$ の例として、6 つの部分子チャネルに分解：

• $gg \to \eta_c[^1S_0^{[1]}]$：2.2995 μb（主導）
• $gg \to \eta_c[^1S_0^{[8]}]$：33.43 nb
• $q\bar{q} \to \eta_c[^3S_1^{[8]}]$：603.7 nb
• その他のチャネル：$\sim 0$（量子数保存により禁止）

重要な洞察：単純な速度スケーリング推定は頻繁に失敗する。理由は以下の通り：

1. 量子数保存により特定のチャネルが消失
2. プロパゲーター増強効果（例：$g^* \to (Q\bar{Q})[^3S_1^{[8]}]$）
3. 運動学的/動力学的抑制または増強

関連研究

汎用 MC 生成器

• Pythia v8.3/Herwig v7.4：クォーコニウム特定のパートン シャワーを導入したが、行列要素機能との完全な互換性はまだない
• 限界：クォーコニウム生成は主に強子化またはデケイを通じて行われ、行列要素実装は限定的

過程非依存のパートン レベル イベント生成器

• MadOnia：MadGraph v4 に基づき、単一の S/P 波クォーコニウムをサポートするが、MG5_aMC への移植がなされていない
• HELAC-Onia：複数の S/P 波クォーコニウムをサポートするが、Fock 状態の手動結合が必要であり、LDME は実行時調整不可
• 本研究の利点：
  1. 最新の MG5_aMC フレームワーク、積極的なメンテナンス
  2. 自動 Fock 状態結合
  3. 実行時に調整可能な LDME
  4. レプトニウムのサポートが初めて

過程特定生成器

• SuperChic/STARlight/eSTARlight：排他過程に焦点
• EPOS4：重イオン衝突におけるクォーコニウム、Wigner 密度行列形式を採用
• BCVEGPY：$B_c$ 中間子強子生成に特化
• 本研究の利点：汎用性、任意の S 波束縛状態の組み合わせをサポート

断面積計算器

• FDCHQHP：S/P 波クォーコニウム強子生成の NLO 計算
• FONLL：非瞬発 J/ψ の NLO+NLL 計算
• 本研究の位置付け：LO 自動化、NLO 拡張の基礎確立

記号計算ツールキット

• FeynOnium/AmpRed：半自動記号計算
• 本研究との違い：完全自動化された数値計算とイベント生成

結論と議論

主要な結論

1. 技術的成果：MG5_aMC フレームワークにおける S 波クォーコニウムおよびレプトニウムの自動化イベント生成の実装に成功
2. 十分な検証：HELAC-Onia および解析表現との詳細な比較を通じて、実装の正確性を確認
3. 広範な応用：pp、e⁺e⁻、e⁻p 衝突型加速器における様々な生成過程をカバー
4. 理論的洞察：次領頭次寄与の重要性は速度スケーリング規則のみでは推定不可能であり、過程ごとの具体的分析が必要
5. レプトニウム：レプトニウム生成の初の体系的予測を提供し、将来の実験観測に指針を与える

限界

1. LO 精度：現在は領頭次のみを実装、NLO 修正は大きい可能性がある（例：三 J/ψ 生成）
2. S 波限定：S 波状態のみをサポート、P 波状態の実装は将来の課題
3. DPS 欠落：二重/三重部分子散射寄与を含まない、これが特定の過程で主導する可能性がある
4. 前馈効果：励起状態のデケイからの前馈寄与を考慮していない
5. 位相空間積分：簡略化された疑似マルチチャネル戦略を採用、極めて複雑な過程では効率が不十分な可能性
6. 低エネルギー領域：$pp \to \text{Ps}_1$ などの極低エネルギー過程では、因子化仮説が失効する可能性がある

将来の方向

1. NLO 実装：拡張 FKS 減除スキームに基づく NLO 精度の実装（理論的基礎は既に存在67）
2. P 波サポート：軌道および全角運動量投影演算子の実装
3. マルチチャネル最適化：束縛状態に特化した高効率位相空間パラメータ化の開発
4. グローバルフィッティング：実験データを組み合わせた LDME グローバルフィッティング
5. DPS 統合：DPS フレームワークとの結合、完全な予測の提供
6. 重イオン衝突：核 PDF および冷核物質効果への拡張

物理的影響

1. 実験指導：LHC、EIC、Belle II などの実験によるクォーコニウム/レプトニウム研究に信頼できる予測を提供
2. メカニズム理解：色一重項/色八重項メカニズム、SPS/DPS 寄与の区別を支援
3. 新物理探索：BSM シナリオにおける束縛状態生成研究をサポート
4. 精密検証：QCD/QED 精密検証のための理論ツールを提供

深い評価

利点

技術的革新性

1. 優れた設計：Boundstate クラスの導入と UFO 形式の修正を通じて、束縛状態を MG5_aMC フレームワークに優雅に統合
2. ユーザーフレンドリー：命名規則が明確（例：Jpsi(1|3S11)）、LDME は実行時調整可能、学習曲線が緩い
3. モジュール化実装：色/スピン/軌道投影の分離により、拡張と保守が容易
4. 多層検証：行列要素、断面積、解析公式の三層検証により、信頼性を確保

科学的貢献

1. 先駆性：主流 MC フレームワークにおけるレプトニウム自動化計算の初の実装
2. 体系性：3 種類の衝突型加速器、40 以上の過程をカバーし、広範な適用性を実証
3. 深い洞察：多くの例を通じて速度スケーリング規則の限界を説明し、誤解を正す
4. 実用的価値：実験物理学者に使いやすいツールを提供し、理論計算の敷居を低下

論文の質

1. 構造が明確：理論-実装-検証-応用のロジックチェーンが完全
2. 詳細が充分：完全なパラメータ設定と再現可能なコマンド例を提供
3. 図表が豊富：68 ページに 5 つの図と 10 の表を含み、情報密度が高い
4. 文献が包括的：297 の参考文献により、関連研究を十分に整理

不足と改善提案

方法的限界

1. 位相空間効率：疑似マルチチャネル戦略は高多重度過程の効率が限定的。提案：
   • 束縛状態特定の重要度サンプリング開発
   • 機械学習による位相空間マッピング最適化
2. 低エネルギー挙動：レプトニウムの極低エネルギー領域での因子化仮説の慎重な論証が必要
3. LDME 不確定性：固定 LDME 値を採用、不確定性の影響を体系的に評価していない

実験設定

1. カット選択：特定の運動学的カット（例：$p_{T,j} > 2$ GeV）は恣意的であり、体系的研究が不足
2. PDF 依存性：単一 PDF セットのみを使用、PDF 不確定性を評価していない
3. スケール依存性：特定の過程ではスケール変化不確定性が示されていない

物理的議論

1. メカニズム区別：異なる Fock 状態寄与を示しているが、実験的区別方法の深い議論が不足
2. 可観測量：主に全断面積と単粒子分布に焦点、偏極や相関などの可観測量が不足
3. 背景推定：実験背景と検出効率の議論がない

影響力評価

学術的影響

• 引用可能性：高い。MG5_aMC 束縛状態実装の初例として、この分野の標準ツールになる見込み
• 後続研究：多くの応用研究（LHC データ分析、EIC 物理予測など）を促進する見込み
• 理論発展：NLO 自動化、DPS 統合などの後続研究の基礎を確立

実用的価値

• 実験応用：ATLAS、CMS、LHCb、Belle II などの実験が直接使用可能
• 教育的価値：束縛状態物理研究の敷居を低下させ、新人育成に有利
• オープンソース貢献：コードは Launchpad と NLOAccess を通じて公開予定、コミュニティ発展を促進

再現性

• コード利用可能：公開発表予定、補助資料を提供
• ドキュメント詳細：論文に完全な使用例とパラメータ設定を含む
• 検証充分：複数の交差検証により結果の信頼性を確保

適用シナリオ

理想的な応用

1. LHC 物理：クォーコニウム生成断面積予測、メカニズム研究、新物理探索
2. EIC 物理：深非弾性散乱におけるクォーコニウム生成、グルオン TMD 抽出
3. B ファクトリー：精密 QED 検証、レプトニウム探索
4. 理論研究：LDME フィッティング、色八重項メカニズム検証、DPS 研究

不適用なシナリオ

1. 重イオン衝突：熱/冷核物質効果が不足、EPOS4 などのツールとの結合が必要
2. 超高精度予測：LO 精度では不十分、NLO 実装を待つ必要
3. P 波主導過程：χc 生成など、P 波実装を待つ必要
4. 排他過程：中心排他生成など、SuperChic などの専用ツールとの結合が必要

提案と展望

著者への提案

1. NLO 精度の実装を急ぐ。これが最も緊急な要求
2. P 波サポートを開発し、Fock 状態カバレッジを完善
3. LDME グローバルフィッティングフレームワークを提供
4. 実験グループとの協力を強化し、ツール応用を推進

ユーザーへの提案

1. LO 精度の限界を理解し、絶対予測の解釈に慎重
2. 体系的な不確定性評価（スケール、PDF、LDME）を実施
3. 他のツール（DPS 生成器など）と結合して完全な画像を取得
4. 使用体験をフィードバックし、ツール改善を促進

分野への提案

1. 標準化された LDME データベースの構築
2. 束縛状態物理の NLO 自動化技術の発展
3. より多くの過程と可観測量の実験測定カバレッジ拡大
4. 理論-実験協力の強化、精密 QCD/QED 検証の実現

参考文献（精選）

核心的理論基礎

• 9 Bodwin, Braaten, Lepage (1995)：NRQCD 因子化形式の基礎的研究
• 64 Caswell, Lepage (1986)：NRQED 有効ラグランジアン
• 67 A A H, Shao, Simon (2024)：クォーコニウム生成の NLO FKS 減除

関連ツール

• 25 Artoisenet et al. (2008)：MadOnia - 前世代 NRQCD 自動化ツール
• 27,28 Shao (2013, 2016)：HELAC-Onia - 主要な比較基準
• 6,7 Alwall et al. (2014, 2018)：MadGraph5_aMC@NLO フレームワーク

重要な応用

• 112 Lansberg, Shao (2013)：J/ψ 対生成の α⁵s 修正の重要性
• 194 Shao, Zhang (2019)：三 J/ψ 生成を多部分子散射プローブとして
• 4 d'Enterria, Shao (2023)：ダイタウオニウム発見の見通し

総合評価：これは高品質な方法論論文であり、技術実装が堅牢で、検証が充分で、応用が広範である。束縛状態物理の自動化ツール、特にレプトニウム分野における重要な空白を埋めている。LO 精度と S 波限定などの限界があるが、シリーズ研究の第一歩として、将来の発展のための堅実な基礎を確立している。この分野の標準ツールになることが予想され、広範な影響を生じる見込みである。実験および理論研究者は、特に NLO 実装と P 波拡張の後続発展に注視することを推奨する。