Progress on computing the hadronic vacuum polarization contribution to the muon anomalous magnetic moment with staggered fermions

本論文は、物理的なパイオン質量を持つ新しい HISQ 格子エネブルを用いた計算の更新結果と、ベクトル - ベクトル相関関数の効率的な計算に向けたコードおよびアルゴリズムの改善について報告するものです。

要約

この論文は、**「ミューオン（ミュー粒子）という小さな粒子が、なぜ少しだけ『ふらふら』と揺れているのか？」**という不思議な現象を、スーパーコンピューターを使って解き明かそうとする研究の進捗報告です。

専門用語を避け、わかりやすい比喩を使って説明しましょう。

1. 何をしているのか？（ミューオンの「ふらふら」現象）

ミューオンという粒子は、電子の「お兄さん」のような存在ですが、少し重くて不安定です。この粒子には「磁石としての強さ（g 値）」という性質があり、理論上は「2」になるはずですが、実際には「2 より少しだけ大きい」値を持っています。これを**「異常磁気能率（g-2）」**と呼びます。

この「少しだけ大きい」部分の原因の一つが、真空（何もない空間）の中に一瞬だけ現れて消える**「仮想粒子の嵐」**です。まるで、静かな湖（真空）に突然、小さな波（仮想粒子）が立ち、ミューオンがその波に揺さぶられてふらつくようなイメージです。

この「波の影響」を正確に計算することが、この研究の目的です。もし計算値と実験値（実験室で測った値）がズレれば、それは**「まだ見えない新しい粒子や力」**が隠れている証拠になるかもしれません。

2. 使っている道具（格子 QCD と HISQ）

研究者たちは、この「波の影響」を計算するために、宇宙を小さなタイル（格子）に分割してシミュレーションしています。これを**「格子 QCD」**と呼びます。

今回は、**「HISQ（ハイ・インプロブド・スタッガード・クォーク）」**という、非常に精度の高い「計算のルール（アルゴリズム）」を使っています。

• 昔のルール： 計算が少し雑で、細かい波まで捉えきれなかった。
• 今回のルール： 非常に滑らかで、細かい波まで正確に捉えられるようになった。

さらに、今回は**「物理的な質量を持つパイオン（軽い粒子）」**が使われており、まるで「本物の宇宙」に近い環境で計算が行われています。

3. 何が新しくなったのか？（「ノイズ」を消す工夫）

この計算の最大の問題は**「ノイズ（雑音）」**です。
シミュレーションでは、無数の「波（粒子の動き）」を計算しますが、遠く離れた場所での計算ほど、雑音が大きくなり、真の答えが見えにくくなります。

研究者たちは、この雑音を減らすために、**「低モード平均（LMA）」と「全モード平均（AMA）」という 2 つのテクニックを組み合わせ、さらに「新しい分割方法」**を開発しました。

比喩：大規模な合唱団の練習

この計算を「大規模な合唱団の練習」と想像してみてください。

• 合唱団（計算）： 全員が歌うと、遠くにいる人の声が雑音になって聞こえにくくなります。
• 従来の方法： 指揮者が「全員で歌って、その中から良い音を探しなさい」と言っていたので、雑音に埋もれていました。
• 今回の新手法：
  1. 低音パート（低モード）： 声の低いメンバー（重要な波）だけをまず集めて、全員で歌わせるようにします。これで「基調」を正確に捉えます。
  2. 高音パート（高モード）： 残りのメンバーは、低音パートの歌声を基準にして、効率的に歌います。
  3. 新しい工夫（HL 分離）： 以前は「低音と高音を混ぜて計算」していましたが、今回は**「低音パートが歌っている間、高音パートは静かに待機」させるようにルールを変えました。これにより、低音の雑音が高音の計算に混ざり込むのを防ぎ、結果として「雑音が大幅に減り、歌声（答え）がクリアになりました」**。

4. 具体的な成果（スーパーコンピューターの進化）

• より細かいタイル： 以前は「64×64」のタイルでしたが、今回は「144×144」という、非常に細かいタイルで計算を行いました。これは、より微細な「波」を捉えるために必要です。
• メモリ節約の工夫（スパースニング）： 144×144 のタイルは計算量が膨大です。そこで、**「必要な場所だけ計算し、関係ない場所はスキップする」**という知恵を使いました。
  • 例え： 広大な畑の土壌調査をする際、すべての土を掘り起こすのではなく、「一定間隔で掘る」ことで、全体の傾向はわかるのに、作業時間は劇的に短縮できる、という感じです。
• 結果： 新しい方法を使うことで、「計算の誤差（ノイズ）」が約 24% 減りました。特に、遠く離れた場所（長い距離）の計算において、効果が大きかったことが確認できました。

5. 結論と今後の展望

この研究は、「ミューオンのふらつき」を計算するための、より正確で効率的な「計算のレシピ」を完成させたという進捗報告です。

• 現状： 新しい計算方法が有効であることが証明されました。
• 今後： さらに多くの計算を行い、最終的な「理論値」を確定させます。
• 意義： もしこの新しい計算結果と、実験室で測った「ミューオンのふらつき」の値がズレれば、それは**「標準模型（今の物理学の教科書）には載っていない、新しい物理の発見」**につながる可能性があります。

つまり、この論文は**「新しい物理学の扉を開けるための、より鋭い『計算の道具』を作りましたよ」**という報告なのです。

技術概要

この論文は、ミューオンの異常磁気能率（$g-2$）におけるハドロン真空偏極（HVP）への寄与を、特に軽クォーク（アップ、ダウン、ストレンジ）の連結図（connected diagrams）に焦点を当てて計算する、格子 QCD による研究の進捗報告です。MILC コラボレーションが生成した 2+1+1 味（flavor）の HISQ（Highly-Improved Staggered Quark）アンサンブルを用いており、物理的なパイオン質量と非常に細かい格子間隔（0.042 fm）を持つ新しいデータセットでの予備結果を含んでいます。

以下に、問題提起、手法、主要な貢献、結果、および意義について詳細な技術的サマリーを記述します。

1. 問題提起と背景

• 研究の目的: ミューオンの異常磁気能率 $a_\mu = (g-2)/2$ は、実験（フェルミ国立加速器研究所の E989 実験など）と理論の高精度な比較を通じて、標準模型を超える新しい物理（新しい粒子や相互作用）を発見する可能性を秘めています。実験精度が向上するにつれ、理論計算の精度も同等に向上させる必要があります。
• 現在の課題: HVP への寄与は、連結図と非連結図（disconnected diagrams）に分類されますが、本研究では最大の寄与を持つ「軽クォークの連結図」に焦点を当てています。
• 計算上の困難: 格子 QCD 計算において、特に長距離領域（大きな Euclidean 時間 $t$）での統計誤差を低減することは困難です。従来の手法では、低モード（low-modes）と高モード（high-modes）の扱いに課題があり、特に大規模な格子（$1443^3 \times 288$）において計算コストと統計精度のバランスを最適化する必要がありました。

2. 手法とアルゴリズムの改善

本研究では、統計誤差を低減し計算効率を向上させるための 2 つの主要な技術的改善を提案・実装しています。

A. 時間 - 運動量表現とスペクトル分解

$a_\mu$ の計算には、時間 - 運動量表現（time-momentum representation）が用いられます。クォーク伝播関数 $S(x, y)$ をディラック演算子の固有ベクトル（固有モード）を用いて低モード部分（$S_L$）と高モード部分（$S_H$）に分解し、相関関数 $C(t)$ を以下の 4 つの成分に分割します：
$$C(t) = C_{LL} + C_{LH} + C_{HL} + C_{HH}$$

B. 高 - 低（HL）成分の全領域平均化

• 従来の課題: 以前の手法（AMA: All-Mode Averaging）では、HL 成分と HH 成分をまとめて計算していました。これにより、HL 成分の統計誤差が十分に低減されず、長距離領域でのノイズが残っていました。
• 新しいアプローチ: 低モード固有ベクトルをソースとして用い、高モード部分の伝播関数を個別に計算することで、HL 成分に対して**全領域平均（full-volume average）**を実現しました。
• コスト削減の工夫: 各時間スライスで全ての固有モードに対して伝播関数を反転させるのは計算コストが高すぎます。そこで、各時間スライスで低モードソースの線形結合をランダムな重み（ヒット）を用いて作成し、シンク（sink）で同じランダム重みで縮約することで、不要な交差項を平均的に打ち消す手法を採用しました。これにより、計算コストを大幅に抑えつつ、統計誤差を低減しています。

C. 低 - 低（LL）成分のスパース化（Sparsening）

• 課題: LL 成分の計算には「メソン場（meson field）」の構築が必要ですが、格子サイズと固有ベクトル数が増えると計算コストとメモリ使用量が爆発的に増加します（特に $1443^3 \times 288$ のような大規模格子）。
• 解決策: 固有ベクトルの空間的なサンプリングを「スパース化」しました。規則的なパターン（例えば、$s$ 個の点をスキップして次の点を取る）で格子点を省略し、メソン場を構築します。
  • 隣接点は相関が高いため、全ての点を含めなくても統計誤差への影響は限定的です。
  • この手法により、必要なメモリフットプリントと計算量が劇的に減少し、特に細かい格子間隔（$a \to 0$）での計算を可能にします。

D. 保存電流の採用

局所的な電磁流は保存されず、格子間隔依存の補正が必要になるため、厳密に保存される点分割電流（point-split current）を使用しています。

3. 主要な結果

• 統計誤差の低減:
  • 既存の $64^3 \times 96$ アンサンブル（格子間隔 0.08787 fm）を用いた比較において、新しい手法（HL 成分の分離計算＋10 ヒット）を採用することで、長距離領域（2.3〜3.3 fm）での統計誤差が**23.7%**改善されました。
  • 全体としての誤差改善率は6.78%（1 ヒットの場合）でした。
  • 誤差の主要な寄与源は LL 成分から HL 成分へシフトしており、HL 成分の全領域平均化が有効であることが確認されました。
• 新しいアンサンブル（$1443^3 \times 288$）での予備結果:
  • 格子間隔 0.042 fm、物理パイオン質量を持つ新しい MILC 生成データセットでの計算を開始しました。
  • 中間ウィンドウ（intermediate window）と長距離ウィンドウ（long-distance window）に対する $a_\mu^{\text{HVP, win}}$ の予備値を算出しました（統計誤差のみ）。
  • 現時点では LL 成分のスパース化は適用していませんが、HL 成分の新しい手法を適用しています。

4. 意義と結論

• 計算効率の劇的な向上: 低モードの全領域平均化とスパース化の組み合わせにより、大規模かつ微細な格子（$1443^3 \times 288$）上での高精度計算が現実的な計算コストで可能になりました。
• 理論精度の向上: 実験精度の向上（2025 年頃の E989 最終結果）に追随するため、理論誤差をさらに低減する必要があります。本研究で提案された手法は、特に長距離領域のノイズを抑制し、HVP 寄与の計算精度を高める鍵となります。
• 将来の展望: 現在、より多くの構成（configurations）での計算が進められており、新しい手法が長距離ウィンドウの誤差をさらに低減できるか検証中です。また、LL 成分へのスパース化の適用も進行しており、これにより計算速度がさらに向上することが期待されています。

総じて、この論文は、格子 QCD による $g-2$ 計算において、統計誤差低減と計算コスト削減を両立させるための革新的なアルゴリズム改良を提示し、次世代の高精度計算の基盤を築く重要な進歩を示しています。