[論文レビュー] Dynamical overlap fermions: techniques and results

Stefan Schaefer, Thomas DeGrand|ArXiv.org|Aug 26, 2005

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

ひとこと要約

本稿では、ストラウトスムージング処理を施したゲージリンクと複数時間スケール積分法を用いた、格子QCDにおけるダイナミカルなオーバラップフェルミオンシミュレーションの安定的で効率的なアルゴリズムを提示する。この手法により、ディラック固有モードのランダム行列理論（RMT）解析を用いて、正確なカイラル凝集体の計算が可能となった。本手法は10倍の高速化を達成し、軽いクォーク質量におけるΣ ≈ (280 MeV)³の測定に成功し、RMTの予測と整合することを確認した。

ABSTRACT

We summarize our recent investigations of lattice QCD with dynamical overlap fermions. We sketch algorithmic issues and our approach to solving them. We show our measurement of the topological susceptibility. We describe a computation of the chiral condensate using an analysis of the distribution of eigenmodes of the Dirac operator and Random Matrix Theory.

研究の動機と目的

• 格子QCDにおけるダイナミカルなオーバラップフェルミオンシミュレーションのための安定的で効率的なアルゴリズムの開発。
• 主な2つの課題に対処する：シミュレーション速度の遅さとトポロジカル・チャージの変更の難しさ。
• 低エネルギーのディラック固有モードの分布とランダム行列理論（RMT）を用いて、カイラル凝集体を測定する。
• 異なるトポロジカル・セクターにおけるRMTの予測とシミュレートされた固有モード分布を比較することで、手法の妥当性を検証する。
• 軽いクォーク質量における系の挙動を調査し、トポロジカル・セクターを制限した場合のエルゴディシティを確認する。

提案手法

• 核演算子の条件数を改善するためにストラウトスムージング処理を施したゲージリンクを用いることで、小さな固有値の数を減少させ、オーバラップフェルミオン作用の評価を高速化する。
• ゲージ力とフェルミオン力のスケール差を扱うために、複数時間スケール分子動力学積分法（Sexton-Weingarten）を採用し、ゲージ時間ステップとフェルミオン時間ステップの比率を1/12に設定する。
• 核演算子の固有値がゼロを crosses する際に有効作用に不連続性が生じるのを処理するため、屈折／反射手順を導入し、HMCアルゴリズムにおける詳細平衡を保持する。
• トポロジカル変化面における有効作用とフェルミオン行列式のステップ変化を計算するために、Sherman–Morrisonの公式を適用する。
• オーバラップディラック演算子における符号関数のスペクトル表現を用いて、エルミート的オーバラップ演算子とその二乗形を定義する。
• 次元なし固有モード（ζ = ρλₖΣV）のRMT予測を測定データにフィットさせ、カイラル凝集体Σを抽出する。

実験結果

リサーチクエスチョン

• RQ1 ストラウトスムージング処理を施したゲージリンクの使用は、ダイナミカルなオーバラップフェルミオンシミュレーションの計算コストを著しく低減できるか？
• RQ2 ランダム行列理論（RMT）を用いることで、低エネルギーのディラック固有モードの分布からカイラル凝集体をどれほど正確に抽出できるか？
• RQ3 ダイナミカルなオーバラップフェルミオンを用いた有限体積シミュレーションにおいても、RMTによる固有モード分布の予測は成立するか？
• RQ4 トポロジカル・セクターを制限した場合、系の挙動はいかなるものか？また、トポロジカル変化が禁止された場合、シミュレーションはエルゴディックか？
• RQ5 カイラル凝集体の結果はクォーク質量にどのように依存するか？また、カイラル極限における値は何か？

主な発見

• ストラウトスムージング処理を施したゲージリンクの使用により、1トレジェクトリあたりのディラック演算子の行列・ベクトル乗算回数が約1桁減少した。
• フェルミオン力が約1桁減少し、複数時間スケール統合が効率的に行えるようになった。
• 2番目に低い2つのディラック固有モードの分布について、RMTの予測が測定データとよく一致しており、特にν = 0およびν = ±1セクターで顕著であった。
• 物理的パイオン崩壊定数と有限体積補正係数ρ ≈ 1.4を用いて、カイラル凝集体をΣ ≈ (280 MeV)³として抽出した。
• amₚ = 0.01およびamₚ = 0.005のシミュレーションはν = 0セクターで安定に進行し、エルゴディックでない兆候や系統的バイアスの兆候は観察されなかった。
• 3番目の固有モードの分布はRMTの予測とよく一致したが、|ν| = 1セクターではわずかなずれが観察された。これは、トゥーレスエネルギースケールに近づいている可能性があるためと推測される。