Radiative return at NLOPS accuracy

技術要約：NLOPS精度における放射性リターン

問題提起
フレーバーファクトリー（$e^+e^-$衝突型加速器）における放射性リターン法は、パイオン形式因子 $F_\pi(Q^2)$ を測定するための主要な手法であり、これは、ミューオン異常磁気モーメント $a_\mu$ への主要なオーダーのハドロン真空偏極（HVP）寄与に対するデータ駆動型の分散的計算において、極めて重要な入力値となる。格子QCDは進展しているものの、格子計算と分散的計算の結果との間の不一致（「新しいミューオンg-2の謎」）により、実験的な測定および理論的予測の精度向上が求められている。

Phokhara のような既存のモンテカルロ生成器は、$e^+e^- \to \mu^+\mu^-\gamma$ および $e^+e^- \to \pi^+\pi^-\gamma$ の過程に対して完全な次世代（NLO）放射補正を組み込んでいるが、多重光子放出による高次寄与が含まれていないため、その理論的精度は約0.5%に制限されている。この不確かさは、精密測定における系統誤差の重要な要因となっている。非放射過程（例：BabaYaga@NLO, KKMC）のための既存の生成器は、NLO補正とパートンシャワー（PS）の指数関数化を組み合わせることに成功しているが、放射性リターン過程（硬い光子を伴う $2 \to 3$ 過程）に対応する同様のツールは欠落していた。

手法
著者らは、$e^+e^- \to X^+X^-\gamma$ （ここで $X = \{\mu, \pi\}$）の過程に対する、パートンシャワー（PS）とマッチングされた正確なNLO補正を計算した。この計算は、更新されたバージョンのモンテカルロイベント生成器 BabaYaga@NLO に実装されている。

1. 固定次（Fixed-Order）計算:

   • NLO断面積は、ソフト/仮想（SV）成分と硬い実光子成分を分離するために、位相空間スライシングを用いて計算される。
   • 初期状態放射（ISR）、終状態放射（FSR）、および初期・終状態干渉（IFI）を含むすべての放射補正の源が組み込まれている。
   • パイオンチャネル（$e^+e^- \to \pi^+\pi^-\gamma$）については、スカラーQEDの振幅に適切な仮想性におけるパイオン形式因子 $F_\pi(Q^2)$ を乗じる、因子分解されたスカラーQED（F$\times$sQED）アプローチを採用している。
   • 仮想補正には、一ループおよび二ループ（ペンタゴン）図が含まれ、これらは Colier, LoopTools, Recola を用いて評価され、オンシェル・プロパゲーター付近での数値的不安定性が慎重に処理されている。
   • 非摂動効果（実行電磁結合 $\alpha(Q^2)$ や、$\rho, \omega, \phi$ 共鳴のためのGounaris-SakuraiおよびBreit-Wigner関数によってパラメータ化されたパイオン形式因子を含む）は、光子プロパゲーターに対する乗法的修正として扱われる。

2. パートンシャワーのマッチング:

   • 著者らは、多重光子放出に対する対数近似（LL）補正を指数関数化する、新しいPSアルゴリズムを開発した。
   • $2 \to 2$ 過程のための従来の定式化とは異なり、この手法は、二重硬いブレムストラールング過程（$e^+e^- \to X^+X^-\gamma\gamma$）の正確な行列要素の上にLL近似を構築する。この選択により、排他的イベントのダイナミクスが保持され、三光子終状態の正確なシミュレーションが可能となる。
   • 硬い光子を特定し、$n$ 光子の運動学を正確な行列要素の評価のための縮退した $2 \to 4$ 運動学へと写像するために、CKKW風のクラスタリングアルゴリズムが使用される。
   • LL指数関数化を保持しつつ、NLOPS断面積が正確なNLO SV補正を再現することを保証するために、マッチング係数 $F_{SV}$ が導入される。この係数は、最も硬い光子の運動量を保存する運動学を用いて評価される。

3. 位相空間サンプリング:

   • ソフトおよび衝突特異点、ならびに非摂動的な形式因子の存在下での共鳴構造（例：$\rho$ メソン）を効率的に扱うために、マルチチャネル・インポータンス・サンプリング技術が実装されている。

主な貢献

• 初のNLOPS実装: 本研究は、放射性リターン過程（$e^+e^- \to X^+X^-\gamma$）に対する、PSとマッチングされた正確なNLO補正の初の実装を提示しており、固定次NLOまたはLO+LLツールしか存在していなかった文献上の空白を埋めるものである。
• 新規のマッチングスキーム: 著者らは、$2 \to 2$ の定式化を放射性リターンに適用した際に発生する、共鳴付近での偽の増大という問題を回避する、$2 \to 3$ 過程のための特定のマッチング定式化を導入した。
• 検証: 実装は、独立した固定次コード（McMule, Phokhara, Recola, Alpha）に対して検証されており、ソフト・ハード光子分離パラメータ $\epsilon$ に対して安定していることが示されている。
• ゲージ不変解析: 様々な実験カット（KLOE Large Angle, KLOE Small Angle, BES3, B-factory）の下での、LOおよびNLOにおけるISR、FSR、およびIFIの相対的な大きさの寄与を定量化し、不変質量分布と前方後方非対称性が提供する相補的な情報を実証している。

結果

• 断面積の予測: 四つの実験シナリオについて、積分および微分断面積の数値結果が示されている。NLO補正は有意であり（特定の運動学的領域では数十パーセントに達することもあるが、全域的には通常 $\sim 1\%$）、イベント選択に依存して変化する。
• 高次効果: 多重光子放出（NLOPS）の包含は、NLOを超える補正をもたらす。これらの高次補正は非自明であり、高不変質量領域や共鳴付近でパーセントレベルに達する。例えば、KLOE Large Angleシナリオでは、高次補正が高不変質量領域をパーセントレベルで修正する。
• 比率の予測: パイオンとミューオンの微分断面積の比 $R(s')$ が計算されている。この比に対する高次補正は、Small Angleシナリオでは小さい（$\sim 0.1\%$）が、Large Angleの設定では $\sim 2\%$ に達する場合がある。
• $2 \to 2$ との比較: 本論文は、$2 \to 2$ のNLOPS定式化を放射性リターン（$2 \to 3$）シナリオに適用すると、パーセントレベルの偏差が生じることを示しており、専用の $2 \to 3$ 定式化が必要であることを裏付けている。

意義
著者らは、本計算が放射性リターン実験に対して改善された精度を持つ最先端の理論的予測を提供すると主張している。正確なNLO補正と排他的な多重光子放出の両方を包含することで、新しい BabaYaga@NLO は放射補正に関連する理論的不確かさを低減する。著者らは、観察された高次効果が現在のフレーバーファクトリーの測定における系統誤差と同等か、あるいはそれよりも大きいことを指摘している。したがって、データ解析におけるこの生成器の使用は、パイオン形式因子の抽出および $a_\mu$ へのHVP寄与の決定における不一致の解決に寄与することが期待される。コードは、完全な排他的シミュレーションおよびデータ解析を容易にするために公開されている。