Stochastic Resolution-of-the-Identity Auxiliary-Field Quantum Monte Carlo: Scaling Reduction without Overhead

Abstract

We explore the use of the stochastic resolution-of-the-identity (sRI) with the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) method. sRI is combined with four existing local energy evaluation strategies in ph-AFQMC, namely (1) the half-rotated electron repulsion integral tensor (HR), (2) Cholesky decomposition (CD), (3) tensor hypercontraction (THC), or (4) low-rank factorization (LR). We demonstrate that HR-sRI achieves no scaling reduction, CD-sRI scales as $\mathcal{O}(N^3)$, and THC-sRI and LR-sRI scale as $\mathcal{O}(N^2)$, albeit with a potentially large prefactor. Furthermore, the walker-specific extra memory requirement in CD is reduced from $\mathcal{O}(N^3)$ to $\mathcal{O}(N^2)$ with sRI, while sRI-based THC and LR algorithms lead to a reduction from $\mathcal{O}(N^2)$ extra memory to $\mathcal{O}(N)$. Based on numerical results for one-dimensional hydrogen chains and water clusters, we demonstrated that, along with the use of a variance reduction technique, CD-sRI achieves cubic-scaling {\it without overhead}. In particular, we find for the systems studied the observed scaling of standard CD is $\mathcal{O}(N^{3-4})$ while for CD-sRI it is reduced to $\mathcal{O}(N^{2-3})$. Once a memory bottleneck is reached, we expect THC-sRI and LR-sRI to be preferred methods due to their quadratic-scaling memory requirements and their quadratic-scaling of the local energy evaluation (with a potentially large prefactor). The theoretical framework developed here should facilitate large-scale ph-AFQMC applications that were previously difficult or impossible to carry out with standard computational resources.

References in corpus (10)

• Quantum embedding theories
• Auxiliary-field quantum Monte Carlo calculations of molecular systems with a Gaussian basis
• On Achieving High Accuracy in Quantum Chemical Calculations of 3d Transition Metal Systems: A Comparison of Auxiliary-Field Quantum Monte Carlo with Coupled Cluster, Density Functional Theory, and Experiment for Diatomic Molecules
• Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- and Third-Order Møller-Plesset Perturbation Theory
• Phaseless auxiliary-field quantum Monte Carlo calculations with planewaves and pseudopotentials--applications to atoms and molecules
• Embedded fragment stochastic density functional theory
• Eliminating spin contamination in auxiliary-field quantum Monte Carlo: realistic potential energy curve of F2
• An Auxiliary-Field Quantum Monte Carlo Study of the Chromium Dimer
• Computation of ground-state properties in molecular systems: back-propagation with auxiliary-field quantum Monte Carlo
• Hamiltonian symmetries in auxiliary-field quantum Monte Carlo calculations for electronic structure

Cited by in corpus (6)

• A Phaseless Auxiliary-Field Quantum Monte Carlo Perspective on the Uniform Electron Gas at Finite Temperatures: Issues, Observations, and Benchmark Study
• The performance of phaseless auxiliary-field quantum Monte Carlo on the ground state electronic energy of benzene
• Unbiasing Fermionic Auxiliary-Field Quantum Monte Carlo with Matrix Product State Trial Wavefunctions
• Benchmark Phaseless Auxiliary-Field Quantum Monte Carlo Method for Small Molecules
• Improved modularity and new features in ipie: Toward even larger AFQMC calculations on CPUs and GPUs at zero and finite temperatures
• Selected configuration interaction wave functions in phaseless auxiliary field quantum Monte Carlo