Systematic lowering of the scaling of Monte Carlo calculations by partitioning andsubsampling

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.

References in corpus (8)

• Density matrix embedding: A simple alternative to dynamical mean-field theory
• Taming the dynamical sign problem in real-time evolution of quantum many-body problems
• Zero-variance principle for Monte Carlo algorithms
• Quantum Monte Carlo Methods in Nuclear Physics: Recent Advances
• Zero Variance Markov Chain Monte Carlo for Bayesian Estimators
• Simple formalism for efficient derivatives and multi-determinant expansions in quantum Monte Carlo
• Linear scaling DFT calculations for large Tungsten systems using an optimized local basis
• Efficient Monte Carlo sampling by parallel marginalization