We investigate the use of interpolative separable density fitting (ISDF) as a means to reduce the memory bottleneck in auxiliary field quantum Monte Carlo (AFQMC) simulations of real materials in Gaussian basis sets. We find that ISDF can reduce the memory scaling of AFQMC simulations from [Formula: see text] to [Formula: see text]. We test these developments by computing the structural properties of carbon in the diamond phase, comparing to results from existing computational methods and experiment.
We investigate the use of interpolative separable density fitting (ISDF) as a means to reduce the memory bottleneck in auxiliary field quantum Monte Carlo (AFQMC) simulations of real materials in Gaussian basis sets. We find that ISDF can reduce the memory scaling of AFQMC simulations from [Formula: see text] to [Formula: see text]. We test these developments by computing the structural properties of carbon in the diamond phase, comparing to results from existing computational methods and experiment.
auxiliary-field quantum Monte Carlo (AFQMC) is a systematically improvable many-body method, but its application to extended solids has been severely limited by unfavorable computational scaling and memory requirements that obstruct rigorous access to the thermodynamic and complete-basis-set limits. By combining tensor hypercontraction with k -point symmetry, we reduce the computational and memory scaling of AFQMC for solids to O ( N 3 ) and O ( N 2 ) , respectively, with an arbitrary basis, comparable to diffusion Monte Carlo. This enables simultaneous, controlled extrapolation of full-system AFQMC calculations to the thermodynamic and complete-basis-set limits across insulating, metallic, and strongly correlated solids, without embedding, local approximations, empirical finite-size corrections, or composite schemes. Our results establish AFQMC as a general-purpose, systematically improvable alternative to diffusion Monte Carlo and coupled-cluster methods for predictive simulations of solids, enabling accurate energies and magnetic observables within a unified framework.
The auxiliary-field quantum Monte Carlo (AFQMC) method is a general numerical method for correlated many-electron systems, which is being increasingly applied in lattice models, atoms, molecules, and solids. Here, we introduce the theory and algorithm of the method specialized for real materials and present several recent developments. We give a systematic exposition of the key steps of AFQMC, closely tracking the framework of a modern software library we are developing. The building of a Monte Carlo Hamiltonian, projecting to the ground state, sampling two-body operators, phaseless approximation, and measuring ground state properties are discussed in detail. An advanced implementation for multi-determinant trial wave functions is described, which dramatically speeds up the algorithm and reduces the memory cost. We propose a self-consistent constraint for real materials, and discuss two flavors for its realization, either by coupling the AFQMC calculation to an effective independent-electron calculation or via the natural orbitals of the computed one-body density matrix.
Auxiliary-field quantum Monte Carlo (AFQMC) has repeatedly demonstrated itself as one of the most accurate quantum many-body methods, capable of simulating both real and model systems. In this article, we investigate the application of AFQMC to realistic strongly correlated materials in periodic Gaussian basis sets. Using nickel oxide (NiO) as an example, we investigate the importance of finite size effects and basis set errors on the structural properties of the correlated solid. We provide benchmark calculations for NiO and compare our results to both experimental measurements and existing theoretical methods.
We describe an algorithm to reduce the cost of auxiliary-field quantum Monte Carlo (AFQMC) calculations for the electronic structure problem. The technique uses a nested low-rank factorization of the electron repulsion integral (ERI). While the cost of conventional AFQMC calculations in Gaussian bases scales as, where N is the size of the basis, we show that ground-state energies can be computed through tensor decomposition with reduced memory requirements and subquartic scaling. The algorithm is applied to hydrogen chains and square grids, water clusters, and hexagonal BN. In all cases, we observe significant memory savings and, for larger systems, reduced, subquartic simulation time.
We investigate the use of optimized correlation-consistent Gaussian basis sets for the study of insulating solids with auxiliary-field quantum Monte Carlo (AFQMC). The exponents of the basis set are optimized through the minimization of the second-order Møller-Plesset perturbation theory (MP2) energy in a small unit cell of the solid. We compare against other alternative basis sets proposed in the literature, namely, calculations in the Kohn-Sham basis and in the natural orbitals of an MP2 calculation. We find that our optimized basis sets accelerate the convergence of the AFQMC correlation energy compared to a Kohn-Sham basis and offer similar convergence to MP2 natural orbitals at a fraction of the cost needed to generate them. We also suggest the use of an improved, method independent, MP2-based basis set correction that significantly reduces the required basis set sizes needed to converge the correlation energy. With these developments, we study the relative performance of these basis sets in LiH, Si, and MgO and determine that our optimized basis sets yield the most consistent results as a function of volume. Using these optimized basis sets, we systematically converge the AFQMC calculations to the complete basis set and thermodynamic limit and find excellent agreement with experiment for the systems studied. Although we focus on AFQMC, our basis set generation procedure is independent of the subsequent correlated wavefunction method used.
Phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) has recently emerged as a promising method for the production of benchmark-level simulations of medium- to large-sized molecules because of its accuracy and favorable polynomial scaling with system size. Unfortunately, the memory footprints of standard energy evaluation algorithms are nontrivial, which can significantly impact timings on graphical processing units (GPUs) where memory is limited. Previous attempts to reduce scaling by taking advantage of the low-rank structure of the Coulombic integrals have been successful but exhibit high prefactors, making their utility limited to very large systems. Here we present a complementary cubic-scaling route to reduce memory and computational scaling based on the low rank of the Coulombic interactions between localized orbitals, focusing on the application to ph-AFQMC. We show that the error due to this approximation, which we term localized-orbital AFQMC (LO-AFQMC), is systematic and controllable via a single variable and that the method is computationally favorable even for small systems. We present results demonstrating robust retention of accuracy versus both experiment and full ph-AFQMC for a variety of test cases chosen for their potential difficulty for localized-orbital-based methods, including the singlet–triplet gaps of the polyacenes benzene through pentacene, the heats of formation for a set of Platonic hydrocarbon cages, and the total energy of ferrocene, Fe(Cp) 2 . Finally, we reproduce our previous result for the gas-phase ionization energy of Ni(Cp) 2, agreeing with full ph-AFQMC to within statistical error while using less than 1/15th of the computer time.