[論文レビュー] Delay Update in determinant quantum Monte Carlo

Fanjie Sun, Xiao Yan Xu|arXiv (Cornell University)|Aug 23, 2023

Advanced Chemical Physics Studies参考文献 37被引用数 5

ひとこと要約

本論文は、局所的および拡張された相互作用を有する強い電子相関系のシミュレーションを高速化する、一般化された遅延更新スキームを決定的量子モンテカルロ（DQMC）に導入する。全グリーン関数の更新を延期し、対角成分と中間ベクトルのみを追跡することで、1スイープあたりの計算コストを O(βN³) から O(βN²) に削減し、大幅に効率が向上し、ゼロ温度および有限温度のDQMCにおける大規模系のシミュレーションが可能になる。

ABSTRACT

Determinant quantum Monte Carlo (DQMC) is a widely used unbiased numerical method for simulating strongly correlated electron systems. However, the update process in DQMC is often a bottleneck for its efficiency. To address this issue, we propose a generalized delay update scheme that can handle both onsite and extended interactions. Our delay update scheme can be implemented in both zero-temperature and finite-temperature versions of DQMC. We apply the delay update scheme to various strongly correlated electron models and evaluate its efficiency under different conditions. Our results demonstrate that the proposed delay update scheme significantly improves the efficiency of DQMC simulations, enable it to simulate larger system size.

研究の動機と目的

• 頻繁な全グリーン関数更新に起因する決定的量子モンテカルロ（DQMC）シミュレーションにおける計算ボトルネックを解消すること。
• 従来、局所的相互作用および有限温度に限定されていた遅延更新技術を、拡張（複数サイト）相互作用およびゼロ温度DQMCに対応する一般化すること。
• 更新ステップの計算複雑さを低減することで、符号問題のないモデルにおけるDQMCの効率を向上させること。
• 遅延された全グリーン関数更新を通じて、O(βN³)スケーリングの前因子を低減し、より大きな系サイズのシミュレーションを可能にすること。

提案手法

• 本手法は、局所的更新を扱うための一般化されたkランク更新形式を採用し、局所的および拡張的相互作用の両方に適用可能である。
• 時間発展演算子のSVD分解から得られる対角成分および中間ベクトル（U, D, V）のみを保持することで、全グリーン関数行列の更新を延期する。
• 行列式比は、グリーン関数の対角成分のみを用いて計算され、各局所的更新後に逐次更新される。
• 測定または安定化の必要がある場合にのみ、保存されたSVD成分とSherman-Morrison型更新を用いて全グリーン関数を再構成する。
• Trotter分解および Hubbard-Stratonovich 変換を用いて相互作用を分解することで、有限温度およびゼロ温度DQMCの両方の設定に実装されている。
• 計算精度を維持するため、定期的に全SVD成分を用いた数値的安定化を実施し、複雑さはO(βN³)であるが、頻度は低く抑えられる。

実験結果

リサーチクエスチョン

• RQ1 DQMCにおける遅延更新技術は、局所的相互作用を超えて、拡張（複数サイト）相互作用に対しても一般化可能か？
• RQ2 一般化された遅延更新は、ゼロ温度および有限温度DQMCの両方において、計算コストの低減と並行して数値的安定性および精度を維持できるか？
• RQ3 従来の高速更新と比較して、壁時計時間および系サイズスケーラビリティの観点で、遅延更新はどの程度の性能向上を達成するか？
• RQ4 拡張相互作用を有するモデルにおいて、本手法は系サイズNおよび逆温度βに対してどのようにスケーリングするか？

主な発見

• 一般化された遅延更新により、1モンテカルロスイープあたりの計算複雑さがO(βN³)からO(βN²)に低減され、著しく効率が向上した。
• 複雑さスケーリングの前因子が低減されたため、従来の高速更新では不可能だったより大きな系サイズのシミュレーションが可能になった。
• 定期的な全グリーン関数再構成とSVD安定化を組み合わせることで、数値的安定性および精度を維持した。
• 本手法はゼロ温度および有限温度DQMCの両方に対応可能であり、遅延更新の適用範囲を従来の範囲を超えて拡張した。
• 特に、全更新のO(βN⁴)コストが非現実的になるようなモデルでは、実用的な高速化が達成された。