[論文レビュー] Improving DWF Simulations: the Force Gradient Integrator and the Möbius Accelerated DWF Solver

Hantao Yin, Robert D. Mawhinney|arXiv (Cornell University)|Nov 21, 2011

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

ひとこと要約

この論文では、格子QCDシミュレーションの効率を向上させるために、力勾配積分子（FGI）とモビウス加速型ドメインウォールフェルミオン（DWF）ソルバーを提案する。FGIは積分誤差を低減し、より大きなステップサイズを可能にする。一方、モビウスソルバーは、より小さい第五次元を用いることでディラック作用素の応用回数を約2倍削減し、32³×64体積、パイオン質量約180 MeVの大きな格子でも2倍の高速化を達成した。

ABSTRACT

We have implemented a variant of the force gradient integrator proposed by Kennedy et.al. and are using it in our production 2+1 flavor DWF simulations with pion masses of 180 MeV in (4.5fm)3 volumes. We find modest speed-ups (\sim 20%) from using the force gradient integrator, compared to our previously used Omelyan integrator. On other ensembles, primarily finite temperature 2+1 flavor DWF QCD, we have extensively tuned the Hasenbusch preconditioning masses and achieved speed-ups of 2-3x. Here we have also switched to the force gradient integrator, but this change has not had any impact on the speed. We also report on an improved solver for DWF, which uses Möbius fermions, with a smaller fifth dimension than the original DWF fermions, as an intermediate step in the generation of solutions of the Dirac equation. This approach cuts the number of effective Dirac applications by approximately a factor of 2 when the conjugate gradient iteration count is large.

研究の動機と目的

• 分子動力学統合スキームの最適化により、2+1フレーバーのドメインウォールフェルミオン（DWF）格子QCDシミュレーションの効率を向上させること。
• HMCアルゴリズムにおけるシミュレーション時間の大部分を占めるDWFディラック方程式の解法にかかる計算コストを低減すること。
• モビウスフェルミオンと標準DWFの数学的類似性を活用して、DWFディラック方程式の解法を高速化すること。
• 対称的かつシンプレクティック性・可逆性を保ちつつ安定性とステップサイズを向上させる新しい力勾配積分子（FGI）の実装と検証。
• FGIとハーゼンバスヒューズ質量分割および商作用素を組み合わせることで、生産スケールのHMCシミュレーションにおける顕著な高速化を達成すること。

提案手法

• まず、力勾配の小さな指数を用いてゲージ場を摂動し、その後更新された場を使って再び力を再計算する2段階の更新を用いて、対称的かつシンプレクティックな力勾配積分子（FGI）を実装する。
• 2次力勾配項をテイラー展開近似により計算し、$ U' = e^{-\frac{\tau^2}{24}F^j T_j}U $ によるゲージ場のシフトに置き換えることで、効率的な実装を可能にする。
• 最初のステップで得た解を2番目のステップの初期推定値として再利用することで「力勾配予測」を実施し、力評価におけるディラックソルブコストを削減する。
• モビウスフェルミオン作用素の第五次元を小さくすることで、完全なDWF解を近似するモビウス加速型DWFソルバーを導入し、4次元ウィルソンディラック作用素の応用回数を削減する。
• 残差誤差を反復的にモビウス作用素を用いて解くことで、モビウス近似の精度を向上させる欠損補正法を用いる。
• FGIとハーゼンバスヒューズ質量分割および商作用素を組み合わせ、力関連のディラック方程式の解法をやや不正確に許容することで、HMCコストを低減する。

実験結果

リサーチクエスチョン

• RQ1 力勾配積分子（FGI）は、大規模なDWFシミュレーションにおいて、オメリアン積分子と比較して測定可能な性能向上を達成できるか？
• RQ2 モビウスフェルミオン形式は、有効な第五次元サイズを小さくすることで、DWFディラック方程式の解法を高速化できるか？
• RQ3 DWFソルバーにおける精度と計算コストのバランスを最適化するためのモビウス近似の最適パラメータセットは何か？
• RQ4 FGIとハーゼンバスヒューズ質量分割および商作用素の組み合わせにより、生産レベルのHMCシミュレーションで顕著な高速化が達成できるか？
• RQ5 力勾配予測技術は、可逆性やシンプレクティック性を損なわずに、HMCにおける力評価コストをどのように低減するか？

主な発見

• 4.5\textrm{fm}^3$ 格子、180 MeVパイオンの条件下で、力勾配積分子はわずかに20％の高速化を達成したが、より大きな体積ではより大きな向上が期待できる。
• $32^3 \times 64$ 格子では、モビウス加速型DWFソルバーにより、4次元ウィルソンディラック作用素の応用回数が3.6×10⁵から2.0×10⁵に低下し、1.8倍の削減が達成された。
• モビウス加速型DWFソルバーは、直接解法と比較して、ウォールクロック時間で2倍の高速化を達成し、ソルバー自体はわずかに42％の時間しか使用しなかった。
• 欠損補正法により、モビウス近似の精度が向上し、追加コストを最小限に抑えつつ高精度な解が得られた。
• ハーゼンバスヒューズ質量分割と商作用素を組み合わせることで、力関連のディラック方程式の解法をやや不正確に許容でき、有限温度系のシミュレーションで2〜3倍の高速化が達成された。
• モビウスパラメータは1回のエンsemblesごとの調整で十分であり、構成やソースベクトルを問わず有効であり、生産シミュレーションへの広範な再利用が可能であった。