Local Embedding and Effective Downfolding in the Auxiliary-Field Quantum Monte Carlo Method

Toward Linear Scaling Auxiliary-Field Quantum Monte Carlo with Local Natural Orbitals

We develop a local correlation variant of auxiliary-field quantum Monte Carlo (AFQMC) based on local natural orbitals (LNO-AFQMC). In LNO-AFQMC, independent AFQMC calculations are performed for each localized occupied orbital using a truncated set of tailored orbitals. Because the size of this space does not grow with the system size for a target accuracy, the method has linear scaling. Applying LNO-AFQMC to molecular problems containing a few hundred to a thousand orbitals, we demonstrate convergence of total energies with significantly reduced costs. The savings are more significant for larger systems and larger basis sets. However, even for our smallest system studied, we find that LNO-AFQMC is cheaper than canonical AFQMC, in contrast with many other reduced-scaling methods. Perhaps most significantly, we show that energy differences converge much more quickly than total energies, making the method ideal for applications in chemistry and material science. Our work paves the way for linear scaling AFQMC calculations of strongly correlated systems, which would have a transformative effect on ab initio quantum chemistry.

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within "chemical accuracy" (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N⁷) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob's ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N³ - N⁴) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Twenty Years of Auxiliary-Field Quantum Monte Carlo in Quantum Chemistry: An Overview and Assessment on Main Group Chemistry and Bond-Breaking

In this work, we present an overview of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) approach from a computational quantum chemistry perspective and present a numerical assessment of its performance on main group chemistry and bond-breaking problems with a total of 1004 relative energies. While our benchmark study is somewhat limited, we make recommendations for the use of ph-AFQMC for general main-group chemistry applications. For systems where single determinant wave functions are qualitatively accurate, we expect the accuracy of ph-AFQMC in conjunction with a single-determinant trial wave function to be between that of coupled-cluster with singles and doubles (CCSD) and CCSD with perturbative triples (CCSD(T)). For these applications, ph-AFQMC should be a method of choice when canonical CCSD(T) is too expensive to run. For systems where multireference (MR) wave functions are needed for qualitative accuracy, ph-AFQMC is far more accurate than MR perturbation theory methods and competitive with MR configuration interaction (MRCI) methods. Due to the computational efficiency of ph-AFQMC compared to MRCI, we recommended ph-AFQMC as a method of choice for handling dynamic correlation in MR problems. We conclude with a discussion of important directions for future development of the ph-AFQMC approach.

Ab initio calculations in atoms, molecules, and solids, treating spin-orbit coupling and electron interaction on an equal footing

We incorporate explicit, non-perturbative treatment of spin-orbit coupling into ab initio auxiliary-field quantum Monte Carlo (AFQMC) calculations. The approach allows a general computational framework for molecular and bulk systems in which material specificity, electron correlation, and spin-orbit coupling effects can be captured accurately and on an equal footing, with favorable computational scaling vs system size. We adopt relativistic effective-core potentials that have been obtained by fitting to fully relativistic data and that have demonstrated a high degree of reliability and transferability in molecular systems. This results in a two-component spin-coupled Hamiltonian, which is then treated by generalizing the ab initio AFQMC approach. We demonstrate the method by computing the electron affinity in Pb, the bond dissociation energy in Br₂ and I₂, and solid Bi.

Surface science using coupled cluster theory via local Wannier functions and in-RPA-embedding: The case of water on graphitic carbon nitride

A first-principles study of the adsorption of a single water molecule on a layer of graphitic carbon nitride is reported employing an embedding approach for many-electron correlation methods. To this end, a plane-wave based implementation to obtain intrinsic atomic orbitals and Wannier functions for arbitrary localization potentials is presented. In our embedding scheme, the localized occupied orbitals allow for a separate treatment of short-range and long-range correlation contributions to the adsorption energy by a fragmentation of the simulation cell. In combination with unoccupied natural orbitals, the coupled cluster ansatz with single, double, and perturbative triple particle-hole excitation operators is used to capture the correlation in local fragments centered around the adsorption process. For the long-range correlation, a seamless embedding into the random phase approximation yields rapidly convergent adsorption energies with respect to the local fragment size. Convergence of computed binding energies with respect to the virtual orbital basis set is achieved employing a number of recently developed techniques. Moreover, we discuss fragment size convergence for a range of approximate many-electron perturbation theories. The obtained benchmark results are compared to a number of density functional calculations.

Regional Embedding Enables High-Level Quantum Chemistry for Surface Science

Compared to common density functionals, ab initio wave function methods can provide greater reliability and accuracy, which could prove useful when modeling adsorbates or defects of otherwise periodic systems. However, the breaking of translational symmetry necessitates large supercells that are often prohibitive for correlated wave function methods. As an alternative, this paper introduces the regional embedding approach, which enables correlated wave function treatments of only a target fragment of interest through small, fragment-localized orbital spaces constructed using a simple overlap criterion. Applications to the adsorption of water on lithium hydride, hexagonal boron nitride, and graphene substrates show that regional embedding combined with focal-point corrections can provide converged CCSD(T) (coupled-cluster) adsorption energies with very small fragment sizes.

Local embedding of coupled cluster theory into the random phase approximation using plane waves

We present an embedding approach to treat local electron correlation effects in periodic environments. In a single consistent framework, our plane wave based scheme embeds a local high-level correlation calculation [here, Coupled Cluster (CC) theory], employing localized orbitals, into a low-level correlation calculation [here, the direct Random Phase Approximation (RPA)]. This choice allows for an accurate and efficient treatment of long-range dispersion effects. Accelerated convergence with respect to the local fragment size can be observed if the low-level and high-level long-range dispersions are quantitatively similar, as is the case for CC in RPA. To demonstrate the capabilities of the introduced embedding approach, we calculate adsorption energies of molecules on a surface and in a chabazite crystal cage, as well as the formation energy of a lattice impurity in a solid at the level of highly accurate many-electron perturbation theories. The absorption energy of a methane molecule in a zeolite chabazite is converged with an error well below 20 meV at the CC level. As our largest periodic benchmark system, we apply our scheme to the adsorption of a water molecule on titania in a supercell containing more than 1000 electrons.