[論文レビュー] Algorithm Shootout: R versus RHMC

M. A. Clark, Ph. de Forcrand+1 more|ArXiv.org|Sep 30, 2005

Particle physics theoretical and experimental studies参考文献 3被引用数 4

ひとこと要約

この論文は、カイラルフェルミオンおよびステガードフェルミオンを用いた格子QCDシミュレーションにおいて、カイラル相転移付近で、Rational Hybrid Monte Carlo (RHMC) 法とR法を比較している。RHMCはフェルミオン行列式の有理近似を用いることで、R法に内在するステップサイズ誤差を排除し、擬スピン場を介して正確なHMCダイナミクスを実現する。その結果、特に小さなクォーク質量および相転移付近では、大幅に高い効率性と正確性を達成する。

ABSTRACT

We present initial results comparing the RHMC and R algorithms on large lattices with small quark masses using chiral fermions. We also present results concerning staggered fermions near the deconfinement/chiral phase transition. We find that the RHMC algorithm not only eliminates the step-size error of the R algorithm, but is also considerably more efficient. We discuss several possibilities for further improvement to the RHMC algorithm.

研究の動機と目的

• カイラルフェルミオンを用いた格子QCDシミュレーションにおいて、R法と正確なRHMC法の性能および正確性を評価すること。
• 有限ステップサイズ誤差がカイラル転移点付近の相転移位置および物理量に与える影響を調査すること。
• 非整数フェルミオン行列式を必要とするシミュレーションにおいて、RHMCがR法に置き換え可能かどうかを特定すること。
• 多時スケール積分や緩和されたソルバ収束を含む、RHMCの最適化戦略を検討すること。
• 有限温度および低温領域における両法のアルゴリズム的効率性および系統的誤差制御を比較すること。

提案手法

• RHMCは、$ \mathcal{M}^{-\alpha} $ の有理近似を用いてフェルミオン行列式を再定式化し、擬スピン場を介して正確なHMCダイナミクスを実現する。
• 有理関数は部分分数分解 $ r(x) = \sum_{k=1}^{n} \alpha_k / (x + \beta_k) $ に分解され、マルチシフト共役勾配法の使用が可能になる。
• 近似誤差を補正するための確率的受容/拒否ステップが導入され、低次の有理近似であっても正確性が保証される。
• フェルミオン力はHMCに類似した項の和として計算される：$ S'_{\text{F}} = -\sum_{i=1}^{\bar{n}} \bar{\alpha}_i \phi^\dagger (\mathcal{M} + \bar{\beta}_i)^{-1} \mathcal{M}' (\mathcal{M} + \bar{\beta}_i)^{-1} \phi $。
• ドメインウォールフェルミオンでは、二段階の有理近似が用いられる：$ r_1(x) \approx x^{1/4} $、$ r_2(x) \approx x^{-1/2} $、これによりヒートバースト更新が可能になる。
• 計算コストを削減するために、多時スケール積分および小さなシフトに対する緩和された収束が検討されている。

実験結果

リサーチクエスチョン

• RQ1 R法における有限ステップサイズ誤差は、QCDのカイラル相転移位置にどのように影響を与えるか？
• RQ2 RHMCは、R法と同等の物理的結果を、はるかに大きな積分ステップサイズで達成できるか？
• RQ3 小さなクォーク質量およびカイラル転移付近のシミュレーションにおいて、RHMCとR法の相対的な計算コストおよび効率性はどのようになるか？
• RQ4 多時スケール積分や緩和されたソルバ収束は、バイアスを導入することなくRHMCの性能を向上させられるか？
• RQ5 RHMCアルゴリズムは、現代の格子QCDシミュレーションにおいてR法に置き換えるのに十分な正確性と効率性を有するか？

主な発見

• R法はステップサイズ依存性に起因する $ O(\delta\tau^2) $ の系統的誤差を示すが、RHMCにはそのような誤差が存在しない。
• RHMCとR法が一貫性を持つのは、R法が非常に小さなステップサイズ（$ \delta\tau \leq m_{\text{ud}}/4 $）を使用する場合に限る。これは計算コストが非常に高くなる。
• $ \delta\tau = 0.055 $ のとき、RHMCの受容率は $ \approx 84\% $ であったが、R法は $ \delta\tau = 0.0019 $ でも受容率がなく、依然として顕著な誤差を示した。
• $ m_{\text{ud}} = 0.02 $ のとき、R法は $ \delta\tau = 0.005 $ で $ \langle \bar{\psi}\psi \rangle = 0.60817(3) $ を得たが、RHMCは $ \delta\tau = 0.0185 $ で $ 0.60809(1) $ を得た。これは、RHMCがより大きなステップサイズでも正確であることを示している。
• 多時スケール積分や緩和されたソルバ収束を用いることで、RHMCの性能は少なくとも2倍向上可能であると示唆されている。