[論文レビュー] Large-scale simulations with chiral symmetry

The JLQCD Collaboration, T. Kaneko+5 more|arXiv (Cornell University)|Nov 27, 2013

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

ひとこと要約

本稿では、正確なチャーミカル対称性を備えた格子QCDシミュレーションにおけるオーバlapの代替として、計算的に効率的な選択肢を特定するための5次元フェルミオン形式の比較的な研究を提示している。カーネル作用素の最適化、符号関数の近似、およびスミアリングにより、2.4 GeVで約0.5 MeV、3.6 GeVで約0.1 MeVの残差クォーク質量が達成され、高い受容率と低減された自己相関性を有する大規模シミュレーションが可能となった。

ABSTRACT

We carry out a comparative study among five-dimensional formulations of chirally symmetric fermions about the algorithmic performance, chiral symmetry violation and topological tunneling to find a computationally inexpensive formulation with good chiral symmetry. With our choice of the lattice action, we have launched large-scale simulations on fine lattices aiming at a precision study of light and heavy quark physics. We report on the comparative study, current status of the large-scale simulations, and preliminary results on the residual quark mass and auto-correlation.

研究の動機と目的

• 正確なチャーミカル対称性を保存するが、計算的に効率的な格子フェルミオン形式を特定し、軽いおよび重いクォーク物理学の大規模シミュレーションを可能にする。
• 物理的クォーク質量未満の残差クォーク質量を低下させつつ、高いアルゴリズム的性能とトポロジカル安定性を維持する。
• オーバラップ作用素と比較して20倍の速度向上を達成し、重いフレーバー物理学の高精度な研究を可能にする。
• 最適化されたカーネル作用素、符号関数の近似、およびスミアリングにより、ハイブリッドモンテカルロ（HMC）シミュレーションにおける高い受容率を実現する。
• 細かい格子（a⁻¹ ≈ 2.4–3.6 GeV）上で高統計的シミュレーションを可能にし、軽いハドロンおよび重いフレーバー物理学の高精度な研究を実現する。

提案手法

• 可変パラメータを有する5次元ドメインウォールフェルミオン形式を用い、スティュートスミアリングステップ数（N_smr）、カーネル作用素（H_W、H_T、2H_T）、符号関数近似（ε_Z、ε_p）を調整した。
• 木レベルのシマンジックゲージ作用素を用い、分子動力学統合を用いたハイブリッドモンテカルロ（HMC）アルゴリズムを適用し、80％の受容率を維持するようにステップサイズを調整した。
• 軽いクォークに対してはハーゼンブルシュプリコンディショニング（質量パラメータam′ = 0.150）を、ストレンジクォークに対しては有理的HMCを用いて効率性を向上させた。
• HMCシミュレーションにおいて、リープ・フロッグ統合法に比べ2倍の高速化を達成するため、Omelian統合法を採用した。
• アルゴリズム的性能を測定するため、MDステップ数（N_MD,P=0.80）およびCG反復回数（N_inv）を測定し、CPUコストの代理指標とした。
• 物理的ストレンジクォーク質量との一致と、混合因子の推定を用いて、2 GeVにおけるMS̄スケールでの残差クォーク質量（m_res）を計算した。

実験結果

リサーチクエスチョン

• RQ1 カーネル作用素、符号関数近似、およびスミアリングのどの組み合わせが、高いHMC効率を維持しながら最低の残差クォーク質量を達成するか？

主な発見

• 2H_Tカーネル、ε_p近似、6段階スミアリング（N_smr = 6）を用いた形式が、最も低いCPUコストを達成し、オーバラップ作用素と比較してN_MD,P=0.80 × N_invが約20倍低減された。
• a⁻¹ ≈ 2.4 GeVにおけるMS̄スケールでの残差クォーク質量は約0.5 MeVであり、N₅ = 8を用いることでa⁻¹ ≈ 3.6 GeVでさらに0.1 MeVにまで低下した。
• トポロジカルトンネルリングは、軌道長τの増加に伴い改善され、a⁻¹ ≈ 3.6 GeVでτを1から2に増加させることで、トポロジカルセクターのサンプリングを加速した。
• すべてのシミュレーションで高い受容率（P ≈ 0.8–0.9）を維持し、⟨e^{-ΔH}⟩ ≈ 1.0であったため、面積保存性およびHMC実装の正しさが確認された。
• β = 4.17およびβ = 4.35の32³×64×12および48³×96×8の格子上で、それぞれa⁻¹ ≈ 2.4 GeVおよび3.6 GeVのシミュレーションが実施され、BlueGene/Q上で約6か月で10,000 MDタイムが予想された。
• Omelian統合法の使用により、リープ・フロッグ法と比較してMDステップ数が2倍減少し、性能向上が顕著に見られた。