[論文レビュー] The Möbius Domain Wall Fermion Algorithm

Richard C. Brower, H. Neff+1 more|arXiv (Cornell University)|Jun 22, 2012

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

ひとこと要約

本論文は、ウィルソンオーバlaplace核にモビウス変換を施したMöbiusドメインウォールフェルミオン（MDWF）アルゴリズムを導入する。この手法は、計算コストを増加させることなく、カイラル対称性の破れを低減する一般化されたドメインウォールフェルミオン形式である。固定された$L_s$において残留質量（$m_{\text{res}}$）が1桁改善され、$L_s$が大きい場合に$m_{\text{res}} = O(1/L_s^2)$のスケーリングが達成可能となり、最適化された$\alpha$のチューニングにより、$L_s = 32$で$10^3$に達するシャミアの性能を実現する。

ABSTRACT

We present a review of the properties of generalized domain wall Fermions, based on a (real) Möbius transformation on the Wilson overlap kernel, discussing their algorithmic efficiency, the degree of explicit chiral violations measured by the residual mass ($m_{res}$) and the Ward-Takahashi identities. The Möbius class interpolates between Shamir's domain wall operator and Boriçi's domain wall implementation of Neuberger's overlap operator without increasing the number of Dirac applications per conjugate gradient iteration. A new scaling parameter ($α$) reduces chiral violations at finite fifth dimension ($L_s$) but yields exactly the same overlap action in the limit $L_s ightarrow \infty$. Through the use of 4d Red/Black preconditioning and optimal tuning for the scaling $α(L_s)$, we show that chiral symmetry violations are typically reduced by an order of magnitude at fixed $L_s$. At large $L_s$ we argue that the observed scaling for $m_{res} = O(1/L_s)$ for Shamir is replaced by $m_{res} = O(1/L_s^2)$ for the properly tuned Möbius algorithm with $α= O(L_s)$

研究の動機と目的

• 有限の5次元目（$L_s$）における明示的カイラル対称性の破れを計算コストを増加させることなく低減する、より効率的かつ高精度なドメインウォールフェルミオンアルゴリズムの開発。
• オーバラップ核にモビウス変換を施すことにより、シャミアとニューバーグの実装の間を滑らかに補間するドメインウォールフェルミオン作用の一般化。
• カイラル対称性を有限$L_s$で改善するスケーリングパラメータ$\alpha$の導入。$L_s \to \infty$の極限では正確なオーバラップ作用が保持される。
• $m_{\text{res}} = O(1/L_s^2)$のスケーリングを達成することで、$L_s$を小さくしても実用的なシミュレーションが可能となる。
• 5次元ドメインウォール作用を4次元オーバラップ作用素に写像する形式的枠組みの提供。これには、Möbiusフェルミオンのためのウォード＝タカハシ恒等式および保存電流が含まれる。

提案手法

• スケーリングパラメータ$\alpha$でパrameter化されたモビウス変換をウィルソンオーバラップ核に適用し、連続的なドメインウォールフェルミオン作用の族を生成する。
• 計算コストを維持しつつ、共役勾配反復の高速化を図るため、4次元のレッド/ブラックプリコンディショニングを採用する。
• カイラル対称性の破れの指標として残留質量$m_{\text{res}}$を計算し、最適な$\alpha(L_s)$のチューニングにより$m_{\text{res}} = O(1/L_s^2)$を達成する。
• 5次元ドメインウォール作用を有効な4次元オーバラップ作用素に写像し、ギンスパルグ＝ウィルソン関係を保持することで、一貫性のあるカイラル対称性解析を可能にする。
• Möbiusフェルミオンのための保存および部分的に保存される軸ベクトルカレントとそのウォード＝タカハシ恒等式を形式的に導出する。これにより、一貫性のある相関関数解析が可能となる。
• ゲージアンサンブル上での数値的検証により、MDWFとシャミアおよびギャップフェルミオン法を比較。$m_{\text{res}}$は最大で1桁の低減が確認された。

実験結果

リサーチクエスチョン

• RQ1 オーバラップ核にモビウス変換を施すことで、計算コストを増加させることなく、有限$L_s$におけるカイラル対称性の破れを低減できるか？
• RQ2 スケーリングパラメータ$\alpha$の導入により、$L_s$が大きい場合に$m_{\text{res}} = O(1/L_s^2)$のスケーリングが達成されるか？
• RQ3 $L_s \sim 10^3$のシャミアの性能を、$L_s = 32$程度の小さな$L_s$値で達成できるか？
• RQ4 Möbiusフェルミオン形式は、4次元有効作用においてギンスパルグ＝ウィルソン関係とウォード＝タカハシ恒等式をどのように保持するか？
• RQ5 Hasenbuschプリコンディショニングやマルチグリッド技術といった他の手法と効果的に組み合わせることで、シミュレーション効率がさらに向上するか？

主な発見

• 最適なスケーリングパラメータ$\alpha$のチューニングにより、固定された$L_s$においてカイラル対称性の破れが1桁低減される。
• スケーリングパラメータ$\alpha = O(L_s)$の下で、$L_s$が大きい場合に残留質量が$m_{\text{res}} = O(1/L_s^2)$のスケーリングを示し、シャミアの$O(1/L_s)$スケーリングを上回る。
• 数値的結果から、$L_s = 32$におけるMDWFは、$L_s \sim 10^3$におけるシャミアのカイラル性と同等の性能を示す。$L_s$の必要値が著しく削減された。
• モビウス変換により、同じアルゴリズムフレームワーク内でシャミアとニューバーグのドメインウォールフェルミオン作用の間を滑らかに補間可能となり、柔軟性が向上する。
• 形式的枠組みにより、Möbiusフェルミオンのための保存および部分的に保存される軸ベクトルカレントとそのウォード＝タカハシ恒等式の一貫性のある導出が可能となる。
• 本アルゴリズムはQUDAやMADWFといった高度なソルバーと互換性があり、ギャップフェルミオンやマルチグリッドプリコンディショニングと組み合わせることで、さらなる$m_{\text{res}}$低減が期待できる。