[論文レビュー] Frequency-splitting estimators for single-propagator traces

Leonardo Giusti, T. H. Harris+2 more|arXiv (Cornell University)|Jan 1, 2020

Particle physics theoretical and experimental studies被引用数 1

ひとこと要約

この論文は、大クォーク質量におけるホッピング展開を用いた分割偶数法と組み合わせることで、格子QCDにおける単一伝播関数トレース計算における統計的ノイズを低減する周波数分割推定器を導入する。この手法は、紫外フラクチュエーションを分離し、効率的な確率的推定を可能にすることで、特にベクトルカレントにおいて1〜2桁の分散低減を達成し、ハドロン的真空偏移およびミュオンのg−2計算における計算コストを顕著に低減する。

ABSTRACT

In these proceedings we address the computation of quark-line disconnected diagrams in lattice QCD. The evaluation of these diagrams is required for many phenomenologically interesting observables, but suffers from large statistical errors due to the vacuum and random-noise contributions to their variances. Motivated by a theoretical analysis of the variances, we introduce a new family of stochastic estimators of single-propagator traces built upon a frequency splitting combined with a hopping expansion of the quark propagator, and test their efficiency in two-flavour QCD with pions as light as 190 MeV. The use of these estimators reduces the cost of the computation by one to two orders of magnitude over standard estimators depending on the fermion bilinear. As a concrete application, we show the impact of these findings on the computation of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment.

研究の動機と目的

• クォーク線の非連結図における大きな統計的分散、特にランダムノイズおよび真空寄与の低減を目的とする。
• 標準的なモンテカルロシミュレーションと互換性を保ちつつノイズを抑える新しいクラスの確率的推定器を、単一伝播関数トレースのために開発すること。
• ミュオンの異常磁気モーメント（aμ）へのハドロン的真空偏移寄与の計算効率を向上させること。
• 電磁的カレント相関関数やアイソスピン破れ補正などの観測量における非連結寄与を、コスト効率よく評価可能にする。

提案手法

• 異なるクォーク質量における伝播関数の差を扱う分割偶数法を用いて、クォーク伝播関数トレースを分解する周波数分割推定器を提案する。
• 系列における最大クォーク質量に対してホッピング展開（最大2n次まで）を適用し、トレースを計算可能な部分（M2n,mr）と剰余項（R）に分離し、確率的推定を可能にする。
• 連続する質量間のK−1個の差に分割偶数推定器を適用し、最大質量においてホッピング展開を適用してUVフラクチュエーションを制御する。
• 全周波数分割推定器を、(1) 質量差の分割偶数推定器、(2) M2n,mr 項の正確な計算、および (3) 剰余項Rの確率的推定器の和として構築する。
• 剰余項を推定するために、単位分散の独立したガウスノイズ場を用い、次式で表す： ̄τRΓ,r(x0) = −1/(aL3Ns) Σi,x Re[aΓ (η†i Hn mr)(x)Γ (D−1mr Hn mr ηi)(x)]。
• 最終的な推定器は、 ̄τfsΓ,r1(x0) = Σk=1 to K−1 ̄τΓ,rkrk+1(x0) + ̄tMΓ,rK(x0) + ̄τRΓ,rK(x0) で定義され、すべての構成要素を統合する。

実験結果

リサーチクエスチョン

• RQ1 周波数（質量）スケールに基づいて寄与を分離することで、単一伝播関数トレースの分散を低減できるか？
• RQ2 大クォーク質量におけるホッピング展開と分割偶数推定器を組み合わせることで、ランダムノイズの抑制にどの程度効果的か？
• RQ3 物理的点付近のベクトルカレントチャネルにおいて、この手法が計算コストをどの程度低減できるか？
• RQ4 この手法を、より高い効率でミュオンの異常磁気モーメントへの非連結寄与を計算するために応用可能か？

主な発見

• 周波数分割推定器は、F7およびG8アンサンブルにおいて、ランダムノイズ寄与の分散を約2桁低減した。
• ベクトルカレントにおいては、1回の反復あたりのコストが3.3〜6倍に増加するものの、クォーク質量に依存して10〜30倍の計算コスト低減が達成された。
• ゲージノイズ寄与は大N s において飽和しており、小N s ではランダムノイズが支配的であることが確認され、新しい手法がこれを効果的に抑制している。
• 特にホッピング展開を高品質で適用した場合、標準推定器よりも厳密に小さい分散を達成した。
• 数値的テストでは、1.0 fmまで非連結カレント相関関数に良好な信号が得られ、分割偶数推定器を用いてaμ寄与が効率的に計算された。
• この技術は、低モード平均化や希釈と同様に他の分散低減手法とも互換性があり、今後の応用においてさらなるコスト低減が可能である。