[論文レビュー] Systematic lowering of the scaling of Monte Carlo calculations by partitioning andsubsampling

Antoine Bienvenu, Jonas Feldt+2 more|arXiv (Cornell University)|May 2, 2022

Theoretical and Computational Physics参考文献 24被引用数 2

ひとこと要約

本稿では、系を断片に分割し、環境を固定した上でサイドウォークを用いて各断片を部分的にサンプリングすることで、広範な物理量の計算スケーリングを低減する新しいモンテカルロ手法を提案する。条件付き期待値を活用することで、バイアスを導入せずにO(N)の計算効率向上を達成しており、金属系における変分モンテカルロ計算により実証されている。

ABSTRACT

We propose to compute physical properties by Monte Carlo calculations using conditional expectation values. The latter are obtained on top of the usual Monte Carlo sampling by partitioning the physical space in several subspaces or fragments, and subsampling each fragment (i.e., performing side walks) while freezing the environment. No bias is introduced and a zero-variance principle holds in the limit of separability, i.e. when the fragments are independent. In practice, the usual bottleneck of Monte Carlo calculations -- the scaling of the statistical fluctuations as a function of the number of particles N -- is relieved for extensive observables. We illustrate the method in variational Monte Carlo on the 2D Hubbard model and on metallic hydrogen chains using Jastrow-Slater wave functions. A factor O(N) is gained in numerical efficiency.

研究の動機と目的

• フェルミ粒子系における変分モンテカルロで生じるO(N^4)のスケーリングボトルネックを解消すること。
• 系の系に伴う分散スケーリングをバイアスを導入せずに低減すること。
• 空間的局所性と断片間の弱い相関を活用することで、数値的効率を向上させること。
• 明示的な積分を含む任意のマルコフ連鎖モンテカルロ応用に適用可能な一般化されたバイアスのない手法を開発すること。

提案手法

• 物理的性質は、配置空間の断片に分割した上での条件付き期待値を用いて計算する。
• 各断片は環境を固定した状態でサイドウォークを用いて部分的にサンプリングされ、エルゴード性とバイアスのなさが保たれる。
• 改善された推定子 ˜O は、その分散が系のサイズNに対して超線形にスケーリングするように構築される。
• 断片が独立する分離極限におけるゼロ分散原理に依拠している。
• 観測量を断片固有の寄与に分解することで、局所的な分散低減が可能になる。
• 積分が期待値として表現可能な任意のモンテカルロ手法に適用可能である。

実験結果

リサーチクエスチョン

• RQ1 バイアスを導入せずに、広範な物理量のモンテカルロ計算における分散スケーリングをO(N)未満に低下させることは可能か？
• RQ2 追加のサンプリングオーバーヘッドがあるにもかかわらず、断片ベースの部分サンプリング（サイドウォークを用いて）は、計算コストのネット低下をもたらすか？
• RQ3 相関が長距離に及ぶが、大きな断片間では弱い金属系において、この手法はどのように動作するか？
• RQ4 大規模系において効率を最大化するための最適な断片サイズとサイドウォーク長は何か？
• RQ5 理論的に予測されたO(N)の効率向上が、実際の計算でも達成可能か？

主な発見

• 2次元 Hubbardモデルでは実効率向上が約0.075N、金属的水素鎖では0.025Nに達し、系のサイズNに比例して線形に増加する。
• 局所的エネルギーの分散はO(N)からほぼO(1)に低下しており、これが全体の効率向上の要因である。
• Hubbardモデルでは最適なサブシステムサイズがl ≈ √Lに比例し、水素鎖ではn ≈ √(N/2)に比例しており、非分離性が示唆される。
• 最適なサイドウォーク長では相関係数cが最大40%まで低下し、これにより非相関推定値よりも40%高い実効率向上が達成された。
• あらゆるテスト対象系においてE(O)とE(˜O)が誤差範囲内で一致しており、バイアスのなさが維持されている。
• 相関係数c > 1の場合でも、効率向上は頑健に観察され、非相関推定値を40%上回る最適な効率向上が達成された。