[論文レビュー] Disconnected contributions to hadronic structure: a new method for stochastic noise reduction

Sara Collins, Gunnar Bali+1 more|arXiv (Cornell University)|Sep 20, 2007

Particle physics theoretical and experimental studies参考文献 9被引用数 18

ひとこと要約

本稿では、収束する前に共役勾配反復を打ち切ることで、ラティスQCDにおける分散クォークループ計算における確率的ノイズを偏りのない方法で低減する新規手法を提案する。大多数のノイズ低減に安価な部分的逆行列を用い、誤差補正に僅か数回の完全な逆行列計算を行うことで、計算時間で4–12倍の高速化を達成し、はしごパラメータ展開と組み合わせると最大30倍まで向上する。これにより、ハドロン構造観測量の計算の実現可能性が著しく向上する。

ABSTRACT

We present a new method for reducing the stochastic noise of all-to-all propagators based on stopping the inversion of the propagator before convergence. The method is easy to implement, unbiased and independent of the quark action. Applying this method to the calculation of disconnected loops needed for hadronic structure observables we find savings in computer time of factors of 4-12 depending on the operator inserted in the loop. When combined with a hopping parameter expansion technique we obtain combined gains of up to factors of 30 for some operators.

研究の動機と目的

• ハドロン構造観測量の非連結クォークループ計算における高い計算コストと確率的ノイズを解消すること。
• 偏りのない、クォーク作用法に普遍的に適用可能で、容易に実装可能なノイズ低減技術を開発すること。
• ヌクレオンフォームファクターやパートン分布のモーメントへの非連結寄与の計算可能性を向上させること。
• 既存の手法（例：はしごパラメータ展開）と組み合わせて、さらなる性能向上を達成すること。

提案手法

• 共役勾配ソルバーを収束する以前の少数の反復数 $ n_t \ll n_c $ で打ち切ることで、近似逆行列 $ \mathrm{M}^{-1}_{n_t} $ を得る。これにより計算コストを著しく低減する。
• 全逆行列の不偏推定値は、$ L_1 $ 個の部分的逆行列と $ L_2 $ 個の補正項 $ \mathrm{M}^{-1}_{n_c} - \mathrm{M}^{-1}_{n_t} $ の完全な逆行列を組み合わせることで構築される。
• 2つのノイズベクトルは互いに独立であるため、全体の推定値は不偏であり、確率的誤差は $ \mathcal{O}(1/\sqrt{L}) $ のスケーリングを示す。
• 固定された計算コストのもとで、合計分散を最小化することにより、最適な $ n_t $ と比 $ L_1/L_2 $ を決定する。
• はしごパラメータ展開（HPE）と組み合わせることで、分散をさらに低減し、特に重クォーク領域で顕著な効果を示す。
• CGおよびBiCGStabソルバーの両方でテストした結果、両者で一貫した性能向上が得られた。

実験結果

リサーチクエスチョン

• RQ1 収束を早期に打ち切ったソルバー法は、すべての全結合クォーク伝播関数において確率的ノイズを低減しつつ、偏りを生じさせないか？
• RQ2 与えられた分散に対して計算コストを最小化するための、早期打ち切りと補正逆行列の最適なトレードオフは何か？
• RQ3 $ \mathbb{1} $, $ \gamma_5 $, $ \gamma_3\gamma_5 $ などの異なるディラック構造 $ \Gamma $ において、この手法はどのように性能を示すか？
• RQ4 はしごパラメータ展開と組み合わせることで、分散をどの程度さらに低減できるか？
• RQ5 クォーク質量が減少するに従い、特に strange クォークおよび軽いクォーク領域に近づくと、この手法の性能はどのように変化するか？

主な発見

• 収束打ち切りソルバー法により、ストレンジクォーク質量付近のクォーク質量において、非連結ループ計算の計算時間が4–12倍短縮された。これは、演算子 $ \Gamma $ に依存する。
• $ \Gamma = \gamma_3\gamma_5 $ の場合、はしごパラメータ展開と組み合わせることで、合計で最大30倍の高速化が達成された。
• 最適な打ち切り点 $ n_t^{\text{opt}} $ はHPEの導入に伴い増加するが、依然として完全収束の $ n_c $ よりも著しく小さいため、効果的なノイズ低減が可能である。
• この手法は $ m_{PS} \approx 300 \, \text{MeV} $ まで有効であり、いくつかの $ \Gamma $ 演算子について、標準的な確率的手法よりも少なくとも2倍の高速化を維持した。
• CGの代わりにBiCGStabを用いた場合でも、同程度の高速化（例：$ \Gamma = \gamma_3\gamma_5 $ で9倍および24倍）が得られ、異なるソルバー間でも安定した性能を示した。
• 計算コストは非常に効率的であり、$ \not\!D^4 $ の適用コストは完全な逆行列計算の5%にとどまり、HPEに基づく補正は容易に実装可能である。