Realization of all-to-all fermion propagator for the first principle high accuracy strong interaction prediction

問題提起

格子量子色力学（LQCD）は、核子の軸ベクトル結合定数（$g_A$）などの第一原理に基づくハドロン構造の計算における主要なツールである。しかし、サブパーセントの精度を達成するには、以下の2つの大きな課題に対処する必要がある：

1. 全対全（all-to-all）フェルミオン・プロパゲーターの計算コスト： QED補正、アイソスピン対称性の破れ（ISB）、および非連結クォークループに不可欠な$N$点相関関数（$N \ge 4$）を計算する場合、全対全フェルミオン・プロパゲーターが必要となる。厳密な計算は格子体積 $V$ に対して $O(V^2)$ でスケールするため、計算コストが極めて高くなる。
2. 励起状態による汚染（ESC）： 基底状態の行列要素を抽出する際、励起状態（例：ローパー様状態、$\pi N$のような多粒子状態）による汚染が問題となる。ディスティレーション法（低運動量モードへの射影）は、ハドロンスペクトルのためのESC抑制には効果的であるが、局所的なカレント挿入に不可欠な高運動量モードを破棄してしまうため、$N$点関数や行列要素への適用には限界がある。

逐次ソース法、低モード置換（LMS）を用いた確率論的手法、あるいはフィールド・スパースニングといった既存の手法は部分的な解決策を提供するものの、膨大な計算コストをかけることなく、高精度な$N$点関数、非連結図、およびESC抑制を同時に扱うことには苦慮している。

手法：「ブレンディング（Blending）」アルゴルズム

著者らは、全対全フェルミオン・プロパゲーターを効率的かつ偏りのない形で近似するための、新しい確率論的手法である「ブレンディング」アルゴリズムを提案している。

• ブレンディング空間の構築： 本手法は、固定された時間スライスに対して、2つの部分空間を組み合わせることで「ブレンディング空間」を構築する。
  1. 低周波数モード（LFM）： ゲージ共変空間ラプラシアンの最低固有値 $N_e$ 個によって張られる部分空間 $L_1$（ディスティレーションと同様）。これは長距離の物理を捉え、空間的に広がったハドロン状態の構築を可能にする。
  2. 高周波数モード（SHFM）： 直交補空間 $L_2$ は、$L_2$ から一様にサンプリングされた $N_{st}$ 個のランダムベクトル（$\eta_j$）を用いて確率論的に近似される。
• 偏りのない推定（Unbiased Estimation）： 全空間上の恒等演算子が、基底ベクトルへの射影の加重和によって近似される。
  • 低モードの射影体は、厳密な重み（1）を保持する。
  • 確率的な高モードの射影体には、補空間上の恒等演算子の偏りのない推定を保証するために、重み $\dim(L_2)/N_{st}$ が割り当てられる。
  • 場の積（例：$N$点関数における）を含む演算子の場合、二重カウントを補正し、特に対角成分の確率的項において推定量が偏らないように、重みテンソル $\Omega^{(2)}_{ij}$ が適用される。
• 圧縮： 全対全クォーク・プロパゲーター $S$ は、このブレンディング空間に射影されることで圧縮される。これにより、ストレージ要件が $O((4N_c N_L^3 N_T)^2)$ から $O((4N_T(N_e + N_{st}))^2)$ へと削減され、任意の補間演算子やカレント挿入を用いた柔軟な相関関数の計算が可能になる。
• 励起状態の制御： 本手法は、追加の反転コストをかけることなく、豊富な基底（$N f_1$ 散乱状態のような「カレントを含む」演算子を含む）を用いた一般化固有値問題（GEVP）解析を容易にする。なぜなら、圧縮されたプロパゲーターは複数の演算子の組み合わせに対して再利用可能だからである。

主な貢献と結果

著者らは、2つのCLQCDアンサンブル、C24P29（$m_\pi \approx 292$ MeV, $a \approx 0.105$ fm）および F64P13（$m_\pi \approx 134$ MeV, $a \approx 0.077$ fm）を用いた応用を通じて、ブレンディング法の効率性と正確性を実証している。

1. 偏りのなさの検証：

   • 本手法は、パイオン（$\pi^+$）の保存ベクトルカレント行列要素を用いてテストされた。低モードのみ（$N_{st}=0$）を用いた場合は偏った結果（$\approx 0.64$）となったが、確率的モード（$N_{st} \ge 50$）を含めることで、理論的期待値である 1.000(29) に復元された。
   • パイオンの電気的フォームファクター $f_{\pi\pi}(Q^2)$ および電荷半径についても、3点および4点相関関数の両方を用いた一貫性のチェックが行われた。結果は不確かさの範囲内で一致しており、4点関数の計算の妥当性が示された。

2. 核子の軸ベクトル結合定数（$g_A$）：

   • アイソベクター（$g_A^{u-d}$）： 様々な低モード数（$\bar{N}_e$）およびソース・シンク間隔にわたるジョイント3状態フィットを用いることで、物理的パイオン質量において $g_A^{u-d} = 1.2337(84)$ を得た。この結果は、特定の $\bar{N}_e$ の選択に伴う曖昧さを排除しており、他の高精度研究と極めて良く一致している。
   • シングレットおよびフレーバー分解： 本手法により、シングレット軸チャージ（$g_A^{u+d+s}$）および個々のフレーバー寄与のための非連結クォークループの計算が可能となった。
     • $N, N f_1, N f'_1$ 補間子を含む基底を用いたGEVP解析により、励起状態の汚染を効果的に抑制した。
     • 物理点（$a=0.077$ fm）における最終結果は以下の通りである：
       • $g_A^u = 0.8408(86)$
       • $g_A^d = -0.3929(86)$
       • $g_A^s = -0.0381(57)$
       • $g_A^{u+d+s} = 0.410(20)$

3. 計算効率：

   • ブレンディング法は、従来の逐次ソース法や確率論的手法と比較して、ディラック演算子の反転回数を大幅に削減する。
   • アイソベクター $g_A$ について、本手法によるコストは、同様の精度を目標とする最近の高精度研究（例：ETM, RQCD, PNDME）と比較して、約**1桁（one order of magnitude）**削減された。
   • 非連結ループを含むシングレット $g_A$ について、本手法は 0.70 百万回の反転で 4% の精度を達成した。これは、比較対象となる先行研究の 1.37～34.1 百万回の反転と比較して極めて高い効率である。この効率性は、圧縮されたプロパゲーターを、新たな反転を行うことなく複数の補間子やカレント挿入に対して再利用できることに由来する。

意義と主張

本論文は、ブレンディング法が高精度なハドロン構造のLQCD計算のための堅牢な枠組みを提供し、計算コストと励起状態による汚染という二重の課題に対処していると主張している。

• 系統的な制御： 同じ量子数を共有する複数の補間演算子を用いて相関関数を構築できることにより、変分解析（GEVP）を通じた基底状態の行列要素の堅牢な抽出を可能にし、効果的にESCを抑制する。
• スケーラビリティ： 物理的パイオン質量において、本手法は、大幅に削減された計算リソースを用いて、核子行列要素のサブパーセント精度の予測を可能にする。
• 汎用性： このアプローチは、非連結図、QED補正、およびISB効果を含む様々な $N$ 点関数に適用可能であり、これらは標準模型の検証において極めて重要である。
• 限界： 著者らは、オーバーラップ・フェルミオン（Overlap fermions）は反転コストが高いため、計算コストが依然として高い課題として残っていることを指摘しており、より広い適用性のためのさらなる最適化の必要性を示唆している。

要約すれば、ブレンディング法は、高精度な第一原理に基づく複雑なハドロン観測量の計算を可能にするための重要な進展であり、減少した系統誤差および統計誤差を伴う、多体散乱やパートン物理学への将来の研究を可能にするものである。