[論文レビュー] A (P)HMC algorithm for N_f = 2 + 1 + 1 flavours of twisted mass fermions

T. Chiarappa, R. Frezzotti+1 more|ArXiv.org|Sep 27, 2005

Quantum Chromodynamics and Particle Interactions被引用数 4

ひとこと要約

本稿では、$N_f = 2+1+1$ ねじれ質量フェルミオンを用いた格子QCDシミュレーションのための (P)HMC アルゴリズムを提示する。偶数奇数前処理、質量シフト前処理、複数時間スケールを採用し、重いクォーク系におけるディラック行列の逆平方根に対して多項式近似を導入することで、擬fermion場と再重み付けを用いた非 degenerate な奇妙クォークと charm クォークの効率的かつ安定なシミュレーションを可能にする。

ABSTRACT

We present a detailed design of a (P)HMC simulation algorithm for N_f = 2 + 1 + 1 maximally twisted Wilson quark flavours. The algorithm retains even/odd and mass-shift preconditionings combined with multiple Molecular Dynamics time scales for both the light mass degenerate, u and d, quarks and the heavy mass non-degenerate, s and c, quarks. Various non-standard aspects of the algorithm are discussed, among which those connected to the use of a polynomial approximation for the inverse (square root of the squared) Dirac matrix in the s and c quark sector.

研究の動機と目的

• 非デジェネレートな重いクォーク系（奇妙と charm）を含む $N_f = 2+1+1$ ねじれ質量フェルミオンのための安定的かつ効率的な (P)HMC アルゴリズムの開発。
• 2 フレーバーの擬フェルミオン場表現を用いた、非デジェネレートな重いクォークのシミュレーションに向けた挑戦への対処。
• ディラック行列の逆平方根に対する多項式近似を用いて、フェルミオン行列式の評価における数値的安定性と正確性を確保。
• 複数時間スケールと前処理技術を統合し、分子動力学更新の性能を向上。

提案手法

• アルゴリズムは、$h$-クォークのディラック作用素 $\tilde{Q}_{2\times 2}$ の行列式を表現するため、2 フレーバーの擬フェルミオン場 $\Phi$ を用い、PHMC のアプローチを可能にする。
• 行列 $\hat{S} = \tilde{Q}_h \tilde{Q}_h^\dagger$ の逆平方根を近似する多項式 $P_n(s) \approx 1/\sqrt{s}$ を用い、正定値かつ安定した評価を保証する。
• 多項式近似は $P_{n,s_L}(\hat{S}) = \mathcal{N} \prod_{k=1}^n (\hat{S} - z_k)$ の形で構築され、複素共役に関して対称な根を持つことで正の性質を維持する。
• 高い精度を求めるための第二の多項式 $\tilde{P}_{\tilde{n},s_L}(\hat{S})$ を用い、擬フェルミオン場の生成と再重み付けに使用。$\tilde{n} = 620$ および $s_L = 10^{-4}$ の場合、相対誤差 $\tilde{\delta}_{IR} \leq 10^{-6}$ を達成。
• 偶数奇数前処理と質量シフト前処理を組み合わせ、軽いクォーク系と重いクォーク系の両方に異なる分子動力学時間ステップを適用。
• 64ビット演算を用い、Clenshaw の再帰的評価（例：Clenshaw の再帰）による安定な多項式評価により、丸め誤差を制御。

実験結果

リサーチクエスチョン

• RQ1 非デジェネレートな重いクォーク（奇妙と charm）を含む $N_f = 2+1+1$ ねじれ質量フェルミオンのための (P)HMC アルゴリズムを、どのように効率的に設計できるか？
• RQ2 重いクォーク系におけるディラック行列の逆平方根に対する最適な多項式近似戦略は何か？
• RQ3 行列 $\hat{S}$ のスペクトルに低固有値が存在する状況下でも、数値的安定性と正確性をどのように維持できるか？
• RQ4 複数時間スケールと前処理の影響は、HMC アルゴリズムの性能と受容率にどのように現れるか？
• RQ5 多項式近似を用いて、必要な擬フェルミオン場をどのように高精度に生成できるか？

主な発見

• 次数 $n = 162$ の多項式近似は、$s \in [10^{-4}, 1]$ の範囲で相対フィット誤差 $|R_{n,s_L}(s)| \leq 0.03$ を達成し、$\delta_{IR} = |R_{n,s_L}(10^{-4})|$ を保守的な誤差上限として用いる。
• 第二の高精度多項式（次数 $\tilde{n} = 620$）は、同じ範囲で相対誤差 $\tilde{\delta}_{IR} \leq 10^{-6}$ を達成し、正確な擬フェルミオン生成と再重み付けを可能にする。
• 2 フレーバーの擬フェルミオン場 $\Phi$ の使用により、非デジェネレートな $h$-クォークのディラック作用素の行列式が $\det[\tilde{Q}_{2\times 2}] \propto \int \mathcal{D}\Phi \, e^{-\Phi^\dagger P(\hat{S}) \Phi}$ の形で表現され、PHMC が可能になる。
• Clenshaw の再帰による安定な多項式評価を備えた 64 ビット算術の下で、数値的安定性が確保されている。
• 複数時間スケールと前処理（偶数奇数および質量シフト）の統合により、分子動力学更新の効率が向上している。
• 低固有値領域における相対フィット誤差が小さいことから、誤差を制御したまま $N_f = 2+1+1$ ねじれ質量フェルミオンの正確なシミュレーションが可能であることが、検証済み。