[論文レビュー] Fluctuations and reweighting of the quark determinant on large lattices

Martin Lüscher, Filippo Palombi|arXiv (Cornell University)|Oct 6, 2008

Quantum Chromodynamics and Particle Interactions参考文献 14被引用数 10

ひとこと要約

本論文は、非常に軽量なウィルソンクォーカーを伴う格子QCDのハイブリッドモンテカルロ（HMC）シミュレーションを安定化するための行列式因数分解手法を提案する。この手法では、クォーカー行列式に寄与する低モード成分を再重み付け要因として分離する。行列式のフラクチュエーションの体積に比例する指数的増加にもかかわらず、低固有値のスペクトル剛性のおかげで効率が保たれ、大規模な格子上でも非常に小さなクォーク質量での安定シミュレーションが可能になる。

ABSTRACT

We propose to stabilise HMC simulations of lattice QCD with very light Wilson quarks by splitting the quark determinant into two factors and by treating the factor that includes the contribution of the low modes of the Dirac operator as a reweighting factor. In general, determinant reweighting becomes inefficient on large lattices, because the statistical fluctuations of quark determinants increase exponentially with the lattice volume. Random matrix theory and some numerical studies now suggest that the low-mode contribution to the determinant behaves differently, which allows factorisations to be devised that preserve the efficiency of the simulation on large lattices.

研究の動機と目的

• 非常に軽量なウィルソンクォークを伴う格子QCDのHMCシミュレーションにおけるアルゴリズム的不安定性とサンプリング非効率性を解消すること。
• 通常は再重み付け効率を損なう、格子体積に比例して指数的に増加するクォーク行列式のフラクチュエーションを克服すること。
• クォーク行列式の低モード寄与を分離し、再重み付け要因として扱う行列式因数分解戦略を構築し、シミュレーションの安定性を維持すること。
• 低ディラック固有値のスペクトル剛性が再重み付け要因のフラクチュエーションを抑制することを示し、大規模体積での効率的シミュレーションを可能にすること。
• ステオキスティック推定子と擬スフェルミオン場を用いた実用的な再重み付けフレームワークを提供し、大規模格子上での実行可能性を数値的証明すること。

提案手法

• クォーク行列式を $ \det(D^\dagger D) = W \det(\tilde{D}^\dagger \tilde{D}) $ と因数分解し、$ W $ を再重み付け要因、$ \tilde{D} $ をスペクトルをシフトさせた修正ディラック作用素とする。
• 2つの特定の因数分解を用いる：$ \tilde{D}_1 = D + i\mu\gamma_5 $ および $ \tilde{D}_2 = (D + i\mu\gamma_5) \frac{\gamma_5 D - i\mu}{\gamma_5 D - i\sqrt{2}\mu} $、両者とも $ \mu $ のオーダーのスペクトルギャップを保証する。
• 近似ゼロモードによる不安定性を避けるために、HMCアルゴリズム内での修正行列式 $ \det(\tilde{D}^\dagger \tilde{D}) $ の取り扱いを実施する。
• 再重み付けを $ \langle \mathcal{O} \rangle = \frac{\langle \mathcal{O} W_l \rangle_{\rm m}}{\langle W_l \rangle_{\rm m}} $ で行い、$ \langle \cdot \rangle_{\rm m} $ は修正理論における期待値を表す。
• $ N $ 個の擬スフェルミオン場を用いて再重み付け要因 $ W_l $ を推定し、$ W_{l,N} = \frac{1}{N} \sum_{k=1}^N \exp\left( (\eta_k, [1 - w_l(D^\dagger D)^{-1}] \eta_k ) \right) $ とすることで、$ N $ の増加に伴い分散が低下する。
• ランダム行列理論の知見を活用し、低固有値のフラクチュエーションが $ (\Sigma V)^{-1} $ のオーダーでスケーリングされることを主張し、$ W_l $ の体積依存性が抑制されることを示す。

実験結果

リサーチクエスチョン

• RQ1 行列式再重み付けは、体積に比例して指数的に増加するクォーク行列式フラクチュエーションの影響を受けるが、大規模格子上でも効率的に実行可能だろうか？
• RQ2 ウィルソンフェルミオン形式における低ディラック固有値のスペクトル剛性が、再重み付け要因 $ W_l $ のフラクチュエーションを抑制するのか？ その結果、安定なシミュレーションが可能になるのか？
• RQ3 因数分解の選択（例：$ \tilde{D}_1 $ と $ \tilde{D}_2 $）が、再重み付け要因の統計的フラクチュエーションおよびシミュレーション効率に与える影響は何か？
• RQ4 再重み付け要因のフラクチュエーションは比 $ \mu / m $ にどのように依存するのか？ また、フラクチュエーションを最小化するためにはこの比をどのように調整すべきか？
• RQ5 複数の擬スフェルミオン場を用いたステオキスティック推定子は、実際のところ再重み付け要因の分散を効果的に低減できるのか？

主な発見

• 再重み付け要因 $ W_l $ の統計的フラクチュエーションは、体積に比例して指数的に増加しない。低ディラック固有値の剛性のおかげで、そのフラクチュエーションは $ \mathcal{O}(1/V) $ のオーダーに抑えられる。
• $ a \approx 0.08 $ fm の $ 48 \times 24^3 $ および $ 64 \times 32^3 $ 格子上での数値的結果から、再重み付け要因のフラクチュエーションは体積に伴い増加するが、$ \mu / m $ の調整によって制御可能であることが示された。
• $ \mu_{\rm R}/m_{\rm R} = 1.0 $ の場合、$ W_{2,24} $ のフラクチュエーションは $ \mu_{\rm R}/m_{\rm R} = 1.5 $ 時の $ W_{1,24} $ と同等の大きさであり、より大きな $ \mu $ が体積の増大を補償できることを示している。
• ステオキスティック再重み付け推定子 $ W_{l,N} $ の分散は $ N $ の増加に伴い減少し、真の $ W_l $ の分散に収束するため、$ N $ を大きくすることは精度向上に有効である。
• スペクトル剛性によるフラクチュエーション抑制はウィルソンフェルミオン形式でも継続的に成立するため、非常に軽量なクォーク質量に対しても本手法は有効である。
• 著者らは、非常に軽量なウィルソンクォークを伴うHMCシミュレーションが、この再重み付けアプローチによって安定化可能であると結論づけるが、完全な妥当性の確認には実際の大規模シミュレーションが必要であると述べている。