[論文レビュー] Improvement of algorithms for dynamical overlap fermions

JLQCD Collaboration, Hideo Matsufuru+8 more|ArXiv.org|Oct 3, 2006

Quantum Chromodynamics and Particle Interactions参考文献 7被引用数 6

ひとこと要約

この論文は、ブルー・ジェン・システム上で標準的手法と比較して3倍の性能向上を達成するため、動的オーバラップフェルミオンのハイブリッドモンテカルロ（HMC）アルゴリズムを、5次元共役勾配法ソルバと多重時間ステップ統合、ノイズありメトロポリステストを組み合わせることで改善した。この手法は、追加のウィルソンフェルミオンによる近いゼロ固有値の抑制とソルバ精度の最適化により、計算コストを低減しながらも精度を維持する。

ABSTRACT

We investigate the algorithms for dynamical overlap fermions aiming at improving the performance for large-scale simulations. We look for the best combination of Hybrid Monte Carlo options and iterative quark solvers with respect to the numerical costs. Our main target is a $N_f=2$ simulation with overlap fermion on a $16^3 imes 32$ lattice at lattice spacing around 0.12 fm.

研究の動機と目的

• 大規模な格子QCDにおける動的オーバラップフェルミオンシミュレーションの性能向上。
• HMCシミュレーションで最も計算コストが高くなる、オーバラップディラック作用素の逆行列計算の計算コスト低減。
• 時間ステップやソルバ精度など、HMCアルゴリズムのパラメータ最適化により、効率性の向上。
• 現代のスーパーコンピュータ上で、2つの軽いクォークと小さなクォーク質量を有する大スケールのシミュレーションを可能にすること。
• 追加のウィルソンフェルミオンによるオーバラップ作用素の近いゼロ固有値の抑制により、安定性と精度を確保すること。

提案手法

• オーバラップディラック作用素における符号関数の計算に、部分分数近似を用いた5次元共役勾配法ソルバを採用。
• ゲージ項、$S_{PF1}$、$S_{PF2}$、$S_E$ 項に対して、力の大きさの階層に応じて異なる時間ステップを用いる多重時間ステップ統合を実装。
• フェルミオンソルバの収束を加速するため、より重いクォーク質量 $m'$ を用いた質量前処理を適用。
• $H_W$ の近いゼロ固有値を抑制するために、$\mu=0.2$ のねじれ質量を持つ追加のウィルソンフェルミオンを導入し、反射/屈折のチェック処理の必要性を排除。
• 5次元ソルバの不正確さを補正するため、確率的メトロポリステストを導入し、新しい配置を確率 $P = \min\{1, e^{-dS}\}$ で受け入れる。
• 精度と性能のバランスを取るために、極の数 $N$、$m'$、ステップサイズ比などのソルバパラメータを最適化。

実験結果

リサーチクエスチョン

• RQ1 動的オーバラップフェルミオンシミュレーションにおける計算コストを最小化するための、HMCパラメータと反復的ソルバの最適な組み合わせは何か？
• RQ2 追加のウィルソンフェルミオンの導入は、HMCアルゴリズムの安定性と効率性にどのように影響するか？
• RQ3 5次元共役勾配法ソルバは、標準的および緩和型CGソルバーよりも速度と精度の面で優れているか？
• RQ4 階層的な力項を有するHMCにおいて、多重時間ステップ統合は安定性と性能にどのような影響を与えるか？
• RQ5 ノイズありメトロポリステストは、5次元ソルバの不正確さをどの程度補償できるか、物理的精度を維持したままか？

主な発見

• $N=20$ 極を用いた5次元共役勾配法ソルバは、同じ精度の下で緩和型CG法と比較して2〜3倍速い。
• $\Delta\tau_{(PF2)}/\Delta\tau_{(PF1)}=5$ および $\Delta\tau_{(PF1)}/\Delta\tau_{(G,E)}=6$ の多重時間ステップアプローチは、著しく効率性を向上させる。
• ノイズありメトロポリステストにより、より精度の低い5次元ソルバの使用が可能になり、物理的精度を維持しながら高精度の逆行列計算の必要性が減少する。
• 5次元ソルバ、多重時間ステップ、ノイズありメトロポリステストの組み合わせにより、ブルー・ジェン（512ノード）における標準的なネストドCGと比較して3倍の高速化が達成された。
• 質量 $\mu=0.2$ の追加ウィルソンフェルミオン項は、$H_W$ の近いゼロ固有値を効果的に抑制し、反射/屈折処理の必要性を排除した。
• ブルー・ジェン（512ノード）における性能評価では、改善されたアルゴリズムにより1トラジェクトリの時間が22分に短縮され、標準的なネストドCGの112分から大幅に改善された。