Lattice determination of the higher-order hadronic vacuum polarization contribution to the muon $g-2$

この論文は、CL S エンサンブルを用いた格子 QCD 計算により、ミューオンの異常磁気能率に対するハドロン真空分極の次々項寄与を 0.6% の精度で初めて決定し、その結果がデータ駆動型の評価と統計的に有意な乖離を示していることを報告しています。

要約

この論文は、物理学の大きな謎を解くための「新しい計算方法」を提案した画期的な研究です。専門用語を避け、身近な例えを使って説明します。

1. 何について話しているの？「ミューオン」という不思議なコマ

まず、ミューオンという粒子について考えてみましょう。これは電子に似た粒子ですが、少し重たいです。
このミューオンは、まるで**「コマ」のように回転しています。しかし、このコマの回転の仕方は、私たちが知っている物理の法則（標準模型）で予測されるよりも、少しだけ「ゆがんでいる」ことが分かっています。この「ゆがみ」を「異常磁気能率（g-2）」**と呼びます。

この「ゆがみ」の大きさを正確に測ることは、**「見えない新しい物理（標準模型の外の何か）」**を見つけるための最も鋭い探知機になります。もし実験値と理論値がズレていれば、そこにはまだ見つかっていない新しい粒子や力が隠れているはずです。

2. 問題点：理論計算の「壁」と「揺らぎ」

この「ゆがみ」を理論で計算する際、最大の難関は**「ハドロン真空分極（HVP）」という現象です。
これを「真空の泡」**に例えてみましょう。

• 真空は、何もない空間ではなく、常に**「粒子と反粒子のペア」**が生まれては消える「泡」で満たされています。
• ミューオンがその中を回転すると、この「泡」の影響を受けて、回転のスピードが微妙に変化します。

この「泡」の計算は非常に複雑で、実験データ（過去の測定値）を頼りに推測するしかありませんでした。しかし、最近の実験データ（CMD-3 という実験など）同士で**「値がバラバラ」になってしまい、理論計算がどこを基準にすればいいか迷走していました。まるで、「地図が 2 種類あり、どちらが正しいか分からない」**ような状態です。

3. この論文の解決策：「計算機」でゼロから作り直す

そこで、この論文のチームは**「格子 QCD（格子量子色力学）」という手法を使いました。
これは、「宇宙の最小単位（格子）をコンピュータ上で作り上げ、ミューンの動きを最初からシミュレーションする」**という方法です。実験データに頼らず、物理の法則そのものから計算し直そうという試みです。

最大の工夫：「足してゼロにする」魔法

通常、この「泡」の計算は、遠くまで計算するほどノイズ（誤差）が溜まり、正確になりにくいという弱点がありました。

しかし、この研究チームは**「ある 2 つの効果を足し合わせると、お互いが打ち消し合って、ノイズが劇的に減る」**という不思議な現象を見つけました。

• A という効果（大きなマイナスの値）
• B という効果（大きなプラスの値）

これらを個別に計算すると誤差が大きいのですが、**「A と B をセットにして計算する」と、遠くの方で「ほぼゼロ」になってしまいます。
これは、「大きな波と大きな逆波が重なり合って、海面が静かになる」**ようなものです。この「静かになった部分」を重点的に見ることで、誤差を極限まで小さくすることに成功しました。

4. 結果：驚異的な精度と新しい発見

この新しい計算方法で得られた結果は以下の通りです。

• 驚異的な精度： 誤差が0.6%という、これまでにない高精度になりました。これまでの計算方法（実験データに頼る方法）の精度を2 倍上回っています。
• 大きなズレ： この計算結果は、実験データ（CMD-3 を含まない古いデータ）と比べると、**4.6 倍もの標準偏差（σ）もズレていました。これは「偶然の誤差」の範囲を遥かに超える、「何か新しい物理が起きている可能性が高い」**という強力な証拠です。
• 最新の予測との比較： 2025 年の最新の理論予測（White Paper）とは少しズレていますが、その差は統計的に許容される範囲（1.4σ）です。

5. まとめ：なぜこれが重要なのか？

この研究は、**「実験データに頼らず、純粋な計算だけで、ミューオンの不思議な動きを極めて高い精度で再現できた」**という歴史的な成果です。

• 従来： 「実験データのバラつき」に悩まされていた。
• 今回： 「計算機シミュレーション」でそのバラつきを乗り越え、**「新しい物理の兆候」**をより鮮明に浮かび上がらせた。

まるで、**「霧がかかった山道で、誰かが新しい道（計算機シミュレーション）を開拓し、目的地（新しい物理）への道筋がはっきり見えてきた」**ようなものです。

今後の物理学は、この「計算結果」と「実験結果」のズレが、**「標準模型の壁を越える、新しい物理の発見」**につながるかどうかの決定的な証拠になるでしょう。

技術概要

以下は、提示された論文「Lattice determination of the higher-order hadronic vacuum polarization contribution to the muon 𝑔−2（ミューオン異常磁気能率に対する高次ハドロン真空偏極項の格子 QCD による決定）」の技術的な要約です。

1. 研究の背景と課題

ミューオンの異常磁気能率（$a_\mu$）は、標準模型（SM）を超える物理に対する最も敏感なプローブの一つです。フェルミラボの E989 実験により実験精度は 124 ppb まで向上しましたが、理論予測の主要な不確かさは、$\alpha^2$ オーダーで寄与する一次（LO）ハドロン真空偏極（HVP）項に起因しています。

さらに、標準模型予測の精度を高めるためには、LO だけでなく二次（NLO）HVP 項も同様の精度で評価する必要があります。従来のデータ駆動型（分散関係を用いた）アプローチでは、LO 項と同様に $e^+e^- \to \pi^+\pi^-$ 断面積の測定値に依存しており、CMD-3 などの最近の測定結果と以前のデータ間に矛盾が生じています。このため、第一原理に基づく格子 QCD 計算による NLO 項の独立した検証が急務でした。

2. 手法と計算枠組み

本研究では、初めてサブパーセント精度（相対誤差 0.6%）で NLO HVP 項を格子 QCD で計算しました。主な手法は以下の通りです。

• 計算対象: ミューオン $g-2$ に対する NLO HVP 寄与を構成する 3 つの図（NLOa: 光子線とミューオンループ、NLOb: 電子・タウループ、NLOc: 2 つの QCD 挿入）を網羅的に計算。
• 時空表現: 時間 - 運動量表現（TMR）を採用し、空間的に総和されたベクトル相関関数 $G(t)$ を利用。
• 格子設定: CLS（Coordinated Lattice Simulations）プロジェクトによって生成された 35 個のゲージ・アンサンブルを使用。
  • 味数：$N_f = 2 + 1$（O($a$) 改善されたウィルソン・フェルミオン）。
  • 格子間隔：0.039 fm から 0.097 fm の 6 種類。
  • パイオン質量：物理値を含む広範囲（約 430 MeV から物理点まで）。
• 解析戦略:
  • 時間領域の分解: 短距離（SD）、中間距離（ID）、長距離（LD）に領域を分割して解析。
  • 核関数の相殺効果の活用: NLOa と NLOb の時間核関数は、Euclidean 時間 $t \approx 1.5$ fm 以降で強く相殺し合う（和がゼロに近づく）という構造的な特徴を利用。これにより、統計的ノイズや有限体積効果、スケール設定の不確かさなどの長距離系が支配的な誤差源の影響を大幅に低減。
  • 短距離処理: 対数的に増幅される格子アーティファクト（O($a^2 \ln^2 a$)）を、摂動 QCD（O($\alpha_s^4$)）を用いた核関数の引き算と再付加によって除去。
  • 有限体積補正: 2 段階の手順（Hansen-Patella 法と Meyer-Lellouch-Lüscher 法の組み合わせ）を用いて補正。
  • アイソスピン破れ補正: 電磁相互作用と強い相互作用による破れを、時空表現のスカラー関数を用いた現象論的見積もりで補正（NLOa と NLOb で相殺するため誤差は小さい）。

3. 主要な結果

計算結果は以下の通りです（単位は $10^{-11}$）。

• 個々の寄与:
  • NLOa: $-216.83(81)(87)... \times 10^{-11}$
  • NLOb: $+111.51(55)(64)... \times 10^{-11}$
  • NLOc: $+3.75(6)(6)... \times 10^{-11}$
• 最終的な NLO HVP 寄与（連続極限）:
  $$a_\mu^{\text{hvp, nlo}} = (-101.57 \pm 0.26_{\text{stat}} \pm 0.54_{\text{syst}}) \times 10^{-11}$$
  • 総誤差: $0.59 \times 10^{-11}$（相対誤差 0.6%）。
  • これは、Muon $g-2$ 理論イニシアチブの 2025 年ホワイトペーパー（WP25）のデータ駆動型見積もりに比べて2 倍の精度を達成したことを意味します。

4. 他研究との比較と意義

• WP25 評価との比較: 本研究の結果は、WP25 の平均値 $(-99.6 \pm 1.3) \times 10^{-11}$ よりも 1.4$\sigma$ 低い値ですが、統計的に矛盾するものではありません。
• KNT19 評価との比較: CMD-3 測定値を含まない従来のデータ駆動型評価（KNT19）とは、4.6$\sigma$ の大きな緊張関係（tension）を示しています。これは、LO HVP 項において格子 QCD と CMD-3 以前のデータ間に観測された傾向と完全に一致しています。
• 科学的意義:
  1. 初の高精度格子計算: NLO HVP 項に対する初のサブパーセント精度の格子 QCD 計算であり、LO と NLO の両方を同じ第一原理枠組みで評価する一貫性を初めて実現しました。
  2. 精度向上の鍵: NLOa と NLOb の時間核関数間の強力な相殺効果により、長距離領域の系統誤差が抑制され、LO 項よりも高精度な計算が可能になったことを実証しました。
  3. 標準模型予測の確立: 実験値と理論値の不一致（$g-2$ 異常）の解釈において、LO 項だけでなく高次項も格子 QCD で確実に見積もる必要性を裏付け、データ駆動型アプローチの矛盾を解決する重要なマイルストーンとなりました。

結論

本研究は、ミューオン $g-2$ 問題における理論予測の精度を飛躍的に向上させ、標準模型を超える物理の探索において、格子 QCD がデータ駆動型アプローチを凌駕する可能性を示唆する重要な成果です。特に、NLOa と NLOb の相殺効果を利用した手法は、将来の高次項計算や他の精密観測量への応用においても極めて重要です。