Light-by-light scattering at three loops in massless QCD and QED: amplitudes and cross sections

本論文では、't Hooft-Veltman 次元正則化法を用いたローレンツテンソル分解により計算を簡素化し、質量ゼロの QCD および QED における 3 ループ光 - 光散乱のヘリチティ振幅を導出するとともに、それらを用いて超重心イオン衝突実験の ATLAS データと一致する NNLO 微分断面積を予測した。

要約

1. 物語の舞台：「光の衝突」という不思議な現象

普段、私たちが光（電球の光や太陽光）を見ていても、2 つの光線がぶつかり合うと、お互いにすり抜けていくだけです。まるで幽霊同士がぶつかり合っても何の干渉も起こらないように見えます。これは古典物理学の常識でした。

しかし、量子力学の世界では話は違います。
真空中には「見えない粒子の海（仮想粒子）」が常に湧き上がっています。光が通ると、この海がざわつき、光同士が「見えない壁」を介して相互作用します。これを**「光 - 光散乱（Light-by-Light scattering）」**と呼びます。

• 日常の例え：
  静かな湖（真空）に、2 人の人がボート（光）で進んでいるとします。通常、ボート同士はぶつかりません。しかし、湖の底に巨大なクジラ（仮想粒子）が潜んでいて、そのクジラが波を起こすことで、2 つのボートが互いに揺らされ、結果として「ぶつかった」ように見える現象です。

この論文は、その「クジラが波を起こす仕組み」を、これまで誰も成し遂げなかった**「超・高解像度」**で計算し、実測データと完璧に一致させました。

2. 研究者たちの挑戦：「3 回ループ」という難関

物理学では、現象の計算を「ループ（輪っか）」の数で精度を表します。

• 1 ループ： 基本的な計算（低解像度）。
• 2 ループ： さらに詳しく計算（中解像度）。
• 3 ループ： 究極の計算（超・高解像度）。

これまでの研究は 2 ループまででしたが、この論文は**「3 ループ」の計算に挑みました。
これは、「1 枚の絵を描く作業」**で例えると、以下のような違いがあります。

• 1 ループ： 輪郭線だけを描く。
• 2 ループ： 輪郭線に影をつけて立体感を出す。
• 3 ループ： 影、質感、光の反射、背景の空気感まで全て描き込み、写真と見分けがつかないレベルにする。

この「3 ループ」の計算は、式があまりにも複雑で、従来のスーパーコンピュータでも処理しきれないほどの膨大なデータ（中間段階で 500GB 以上の式）が生まれるため、**「計算の地獄」**と呼ばれていました。

3. 解決策：「整理整頓」と「魔法の道具」

研究者たちは、この「計算の地獄」をどう乗り越えたのでしょうか？

① 複雑な箱を分解する（テンソル分解）

計算対象は、光の振動方向（ヘリシティ）ごとに異なる 8 種類の「箱（テンソル）」が入っていました。これをバラバラに分解し、最もシンプルな形に整理する「『t Hooft-Veltman 法」というテクニックを使いました。

• 例え： 巨大で複雑なパズルを、まずは「赤いピース」「青いピース」のように色ごとに分類し、さらに「角のピース」「真ん中のピース」に分けてから組み立てるようなものです。

② 数学の「魔法の辞書」を使う（積分の削減）

計算の核心部分である「積分（面積や体積を求める計算）」は、何百万通りもありました。しかし、研究者たちは「IBP 法」という数学的なルールを使って、何百万通りもの計算を「486 個の基本的な計算（マスター積分）」にまで減らしました。

• 例え： 世界中のすべての料理のレシピを計算しようとしたとき、「卵料理」「肉料理」「野菜料理」の 3 つの基本パターンさえ分かれば、残りはそれらを組み合わせるだけでいいと気づき、作業量を劇的に減らしたようなものです。

③ 巨大な式を「圧縮」する（部分分数分解）

計算途中に式が巨大化しないよう、独自のプログラムで式を常に整理し、無駄な部分を削除しました。

• 例え： 膨大な荷物を運ぶ際、箱を詰め替えて隙間をなくし、トラックの荷台を最大限に有効活用する「パッキング術」です。

4. 結果：実験データとの「完璧な一致」

計算が終わった後、彼らはこの結果を**「LHC（大型ハドロン衝突型加速器）」**という巨大な実験施設で観測されたデータと比べました。
特に、重い原子核（鉛）を衝突させた際、その周囲の強い電磁気力で光が光とぶつかる現象（ATLAS 実験）のデータを対象にしました。

• 結果：
  彼らが計算した「3 ループ（超・高解像度）」の予測値は、実験で観測されたデータと驚くほど一致しました。
  以前の計算（1 ループや 2 ループ）では、実験データと少しズレがあったり、誤差の範囲が広かったりしましたが、今回の計算ではそのズレが解消され、理論と実験が見事に握手を交わしました。

5. この研究が意味するもの

この研究の成功は、単に「光の計算ができた」だけではありません。

1. 標準模型の勝利： 私たちが知る物理の法則（標準模型）が、極めて高い精度で正しいことを証明しました。
2. 新物理への扉： 「もし実験データと理論がズレたら、それは『未知の新しい粒子』の発見につながる」と言われます。逆に、「ズレがないこと」は、私たちがまだ知らない「新しい物理」を探すための、より厳密な基準線（ものさし）を確立したことを意味します。
3. 計算技術の進化： 複雑すぎる数学の問題を、最新のアルゴリズムと計算機で解き明かす技術が、物理学の最前線でどう使われているかを示す好例です。

まとめ

この論文は、「光が光とぶつかる」という不思議な現象を、人類最高レベルの計算技術で「3 回ループ」の超・高解像度でシミュレーションし、実験結果と完璧に一致させたという、物理学における「大成功」の報告書です。

まるで、「見えないクジラが湖に作る波の形」を、微細な波紋一つまで計算し尽くして、実際に観測された波の形と見事に一致させたような、壮大で美しい研究なのです。

技術概要

この論文は、質量ゼロの量子色力学（QCD）および量子電磁力学（QED）における、光 - 光散乱（Light-by-Light scattering: LbL）の 3 ループ振幅と、それに対応する次々次-leading order（NNLO）の微分断面積の計算を報告したものです。以下に、問題設定、手法、主要な貢献、結果、および意義について詳細な技術的サマリーを記述します。

1. 問題設定と背景

• 物理的意義: 古典的なマクスウェル方程式では光の自己相互作用は禁止されていますが、QED では仮想の真空レプトン - 反レプトン対との相互作用により、光 - 光散乱（LbL）が可能になります。4 光子散乱は、ファリー（Furry）の定理により 3 光子過程が消滅するため、光の非線形性の最低次数の例となります。
• 実験的状況: 大型ハドロン衝突型加速器（LHC）における重イオンの超周辺衝突（UPC）において、ATLAS および CMS 実験によって LbL 散乱が直接観測されています。特に、鉛（Pb）イオン衝突では電荷 $Z=82$ に比例して断面積が $Z^4$ 倍増大するため、高感度で測定可能です。
• 理論的課題: これまでの研究では、1 ループ（LO）および 2 ループ（NLO）の計算は行われていましたが、3 ループ（NNLO）の計算は未完了でした。標準模型の精密検証や、超対称性や余剰次元などの新物理探索のためには、NNLO 精度での理論予測が不可欠です。

2. 計算手法と技術的アプローチ

本研究では、3 ループ質量ゼロの QCD/QED ヘリシティ振幅を計算するために、以下の高度な数値・解析的手法を組み合わせました。

• テンソル分解と正則化:
  • 't Hooft-Veltman (tHV) 次元正則化（$d = 4 - 2\epsilon$）を採用し、ループ運動量を $d$ 次元、外部運動量を 4 次元に保つことで計算を簡略化しました。
  • 外部粒子がカラーを持たないため、振幅は SU($N_c$) のカシミール不変量（$C_F, C_A$）と電荷の関数として多項式で表されます。
  • 8 つのヘリシティ配置に対応する 8 つのローレンツテンソル基底（$T_i$）と、それに対応する形式因子（$F_i$）に振幅を分解しました。
• 積分の削減（IBP 還元）:
  • 3 ループで約 100 万個（1,046,489 個）のファインマン積分を生成し、Reduze 2、Kira 3、および独自開発の Finred フレームワークを用いて、486 個のマスター積分（MIs）に還元しました。
  • 中間式の爆発的な成長（最大 500GB 規模）を抑制するため、段階的な IBP 関係式の適用と、MultivariateApart に基づく多変数部分分数分解を適用しました。
  • 非平面トポロジーが初めて寄与する 3 ループ計算において、既存の還元テーブルを再利用しつつ、新たに 1 万 2 千以上の積分の還元を生成しました。
• 関数基底と解析接続:
  • 有限部分の振幅は、調和多対数関数（Harmonic Polylogarithms: HPLs）の線形結合として表現されました。
  • 物理領域外への解析接続（外部運動量の入れ替え）には、多重対数関数のシャッフル代数関係を利用し、対数関数の変換則を適用しました。
  • 最終的な結果は、最大重み 6 の 23 個の HPL の積で表される極めてコンパクトな形式に整理されました。
• QED への Abel 化（Abelianisation）:
  • 純粋な QCD 振幅を計算した後、グルーオンを光子に置き換える「Abel 化」手順により、QED および混合 QCD-QED 補正を導出しました。
  • 色構造とフレーバー構造に基づき、ゲージ不変な部分振幅（$f^{(i,j)}$）に分解して処理しました。

3. 主要な貢献と結果

• 3 ループ振幅の導出:
  • 質量ゼロの QCD と QED における、すべての独立したヘリシティ振幅（全プラス、単一マイナス、二重マイナスなど）の 3 ループ寄与を初めて解析的に計算しました。
  • 結果は非常にコンパクトであり、数値評価が効率的に行えることが確認されました。
• NNLO 微分断面積の予測:
  • 得られた振幅を用いて、二光子系の不変質量（$m_{\gamma\gamma}$）と rapidity（$y_{\gamma\gamma}$）分布における NNLO 微分断面積を計算しました。
  • QCD と QED の両方の補正、およびそれらの混合項をすべて含めました。
• 実験データとの比較:
  • ATLAS 実験の重イオン UPC データ（$\sqrt{s_{NN}} = 5.02$ TeV）と比較しました。
  • 不変質量分布: LO 予測は低質量領域では実験値と一致しますが、高質量領域でわずかに乖離します。一方、NLO および NNLO 予測はすべてのデータビンで実験値と良好に一致します。
  • K ファクター: NNLO 補正は NLO 補正よりも大きく、特に不変質量分布の端点付近では最大で約 3.5% の補正効果を示しました。これは、3 ループで初めて現れる新しいトポロジー（$N^{(3,4)}$）による寄与が顕著であることが原因です。
  • 全断面積: 計算された NNLO 全断面積は $\sigma_{NNLO} = 102.10^{+0.44}_{-0.25}$ nb であり、ATLAS の測定値 $120 \pm 22$ nb と矛盾しません。

4. 意義と将来展望

• 標準模型の精密検証: 本研究は、LbL 散乱に対するこれまでに最も高精度な理論予測を提供し、標準模型の枠組み内での理解を深めました。
• 新物理探索への基盤: 高精度な理論予測は、クォークやレプトンの束縛状態、あるいは axion のような新粒子探索における背景事象の理解に不可欠です。
• 計算手法の進展: 3 ループレベルでの複雑な積分削減と、巨大な中間式を管理するための数値的・代数的手法（syzygy 関係や有限体復元など）の適用は、他の高次ループ計算への道筋を示しています。
• 今後の課題: 本研究では質量ゼロの近似を用いましたが、低エネルギー領域（チャーム、ボトムクォーク、タウレプトンの閾値付近）では質量効果を考慮する必要があります。また、より高エネルギー領域でのトップクォークループの影響や、重質量粒子を含む計算への拡張が今後の課題として挙げられています。

総じて、この論文は高エネルギー物理学における高次ループ計算の最先端を示すものであり、LHC における重イオン衝突実験のデータ解釈と、標準模型を超える物理の探索に重要な役割を果たす成果です。