Quantum-enhanced quantum Monte Carlo: an industrial view
Maximilian Amsler, Peter Deglmann, Matthias Degroote, Michael Kaicher, Kiser, Matthew, Kühn, Michael, Chandan Kumar, Andreas Maier, Ge. G. Samsonidze, Anna K. Schroeder, Michael Streif, Davide Vodola, Christopher Wever
In this work, we test a recently developed method to enhance classical auxiliary-field quantum Monte Carlo (AFQMC) calculations with quantum computers against examples from chemistry and material science, representatives of classes of industry-relevant systems. As molecular test cases, we calculate the energy curve of H4 and relative energies of ozone and singlet molecular oxygen with respect to triplet molecular oxygen, which are industrially relevant in organic oxidation reactions. We find that trial wave functions beyond single Slater determinants improve the performance of AFQMC and allow to generate energies close to chemical accuracy compared to full configuration interaction (FCI) or experimental results. As a representative for material science we study a quasi-1D Fermi-Hubbard model derived from CuBr2, a compound displaying electronic structure properties analogous to cuprates. We find that trial wave functions with both, significantly larger fidelities and lower energies over a Hartree-Fock solution, do not necessarily lead to better AFQMC results.
In this work, we test a recently developed method to enhance classical auxiliary-field quantum Monte Carlo (AFQMC) calculations with quantum computers against examples from chemistry and material science, representatives of classes of industry-relevant systems. As molecular test cases, we calculate the energy curve of H4 and relative energies of ozone and singlet molecular oxygen with respect to triplet molecular oxygen, which are industrially relevant in organic oxidation reactions. We find that trial wave functions beyond single Slater determinants improve the performance of AFQMC and allow to generate energies close to chemical accuracy compared to full configuration interaction (FCI) or experimental results. As a representative for material science we study a quasi-1D Fermi-Hubbard model derived from CuBr2, a compound displaying electronic structure properties analogous to cuprates. We find that trial wave functions with both, significantly larger fidelities and lower energies over a Hartree-Fock solution, do not necessarily lead to better AFQMC results.
In this work, we test a recently developed method to enhance classical auxiliary-field quantum Monte Carlo (AFQMC) calculations with quantum computers against examples from chemistry and material science, representative of classes of industry-relevant systems. As molecular test cases, we calculate the energy curve of H4 and the relative energies of ozone and singlet molecular oxygen with respect to triplet molecular oxygen, which is industrially relevant in organic oxidation reactions. We find that trial wave functions beyond single Slater determinants improve the performance of AFQMC and allow it to generate energies close to chemical accuracy compared to full configuration interaction or experimental results. In the field of material science, we study the electronic structure properties of cuprates through the quasi-1D Fermi-Hubbard model derived from CuBr2, where we find that trial wave functions with both significantly larger fidelities and lower energies over a mean-field solution do not necessarily lead to AFQMC results closer to the exact ground state energy.
to develop quantum Monte Carlo (QMC) approaches for excited states and for more systematic and accurate treatment of ground states of strongly correlated sytems and heavier elements. We use the auxiliary-field quantum Monte Carlo (AFQMC) framework, which is a many-body total energy method that looks like many coupled mean-field calculations. I will discuss this method’s relations and dierences with traditional density-functional and quantum chemistry methods, and with other forms of QMC. The AFQMC method uses importance-sampled random walks in the space of nonorthogonal Slater determinants to project out the many-body ground state, or excited state. The random walks are constrained by a trial wave function (TWF) to remove the so-called sign problem, giving approximate (but non-perturbative) results. With singledeterminant TWFs from Hartree-Fock or density-functional theory (DFT), this approach can be thought of as a post-DFT method with no parameter. In this form, the accuracy is typically comparable to the quantum chemistry coupled-cluster CCSD(T) method near equilibrium geometries, and better than CCSD(T) when bonds are stretched. The computational cost of AFQMC scales as N 3 with characteristic size of the system. Its accuracy can be further improved at the cost of increased computing (releasing the constraint) or better TWF. With AFQMC, we can systematically target selected orbitals to correlate. This allows many-body calculations in solids with a frozen core, with down-folded Hamiltonians, or possibly embeded in a DFT calculations to extend scales. Results from a few on-going applications will be presented, including band structures in simple solids, magnetic phases in MnO, and Cobalt adsoption on graphene.
High Resolution Image Download MS PowerPoint Slide Quantum Monte Carlo (QMC) methods are powerful approaches for solving electronic structure problems. Although they often provide high-accuracy solutions, the precision of most QMC methods is ultimately limited by the trial wave function that must be used. Recently, an approach has been demonstrated to allow the use of trial wave functions prepared on a quantum computer [Huggins et al., Unbiasing fermionic quantum Monte Carlo with a quantum computer. Nature 2022, 603, 416] in the auxiliary-field QMC (AFQMC) method using classical shadows to estimate the required overlaps. However, this approach has an exponential post-processing step to construct these overlaps when performing classical shadows obtained using random Clifford circuits. Here, we study an approach to avoid this exponential scaling step by using a fixed-node Monte Carlo method based on full configuration interaction quantum Monte Carlo. This method is applied to the local unitary cluster Jastrow ansatz. We consider H 4, ferrocene, and benzene molecules using up to 12 qubits as examples. Circuits are compiled to native gates for typical near-term architectures, and we assess the impact of circuit-level depolarizing noise on the method. We also provide a comparison of AFQMC and fixed-node approaches, demonstrating that AFQMC is more robust to errors, although extrapolations of the fixed-node energy reduce this discrepancy. Although the method can be used to reach chemical accuracy, the sampling cost to achieve this is high even for small active spaces, suggesting caution about the prospect of outperforming conventional QMC approaches.
We systematically investigate the calculation of excited states in quantum chemistry using auxiliary field quantum Monte Carlo (AFQMC). Symmetry allows targeting of the lowest triplet excited states in AFQMC based on restricted open-shell determinants, effectively as a ground-state calculation. For open-shell singlet states, excited-state calculations can be stabilized with the appropriate trial states, but their quality can have a larger effect on the accuracy in AFQMC. We find that active space-based configuration interaction trial states are often not sufficient to obtain accurate results. We instead use truncated equation of motion coupled cluster with single and double excitations (EOM-CCSD) as trial states. We benchmark the performance of these calculations on a set of small and medium molecules and polyacenes, focusing on predominantly single excitations. We find that the AFQMC results, obtained at a per-sample cost scaling of O(N6), are systematically more accurate than those obtained using EOM-CCSD, reducing excitation energy errors by approximately half for open-shell singlets. In regimes where EOM-CC triples are impractical, these results position AFQMC as a scalable, higher-accuracy complement for low-lying excited states.
The auxiliary‐field quantum Monte Carlo (AFQMC) method provides a computational framework for solving the time‐independent Schrödinger equation in atoms, molecules, solids, and a variety of model systems. AFQMC has recently witnessed remarkable growth, especially as a tool for electronic structure computations in real materials. The method has demonstrated excellent accuracy across a variety of correlated electron systems. Taking the form of stochastic evolution in a manifold of nonorthogonal Slater determinants, the method resembles an ensemble of density‐functional theory (DFT) calculations in the presence of fluctuating external potentials. Its computational cost scales as a low‐power of system size, similar to the corresponding independent‐electron calculations. Highly efficient and intrinsically parallel, AFQMC is able to take full advantage of contemporary high‐performance computing platforms and numerical libraries. In this review, we provide a self‐contained introduction to the exact and constrained variants of AFQMC, with emphasis on its applications to the electronic structure of molecular systems. Representative results are presented, and theoretical foundations and implementation details of the method are discussed. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Structure and Mechanism > Computational Materials Science Computer and Information Science > Computer Algorithms and Programming
We introduce an efficient approach to implement neural network quantum states (NNQS) as trial wave functions in auxiliary-field quantum Monte Carlo (AFQMC). NNQS are a recently developed class of variational ansätze capable of flexibly representing many-body wave functions, though they often incur a high computational cost during optimization. AFQMC, on the other hand, is a powerful stochastic projector approach for ground-state calculations, but it normally requires an approximate constraint via a trial wave function or trial density matrix, whose quality affects the accuracy. Recently, it has been shown (Xiao et al., arXiv2505.18519) that a broad class of highly correlated wave functions can be integrated into AFQMC through stochastic sampling techniques. In this work, we apply this approach and present a direct integration of NNQS with AFQMC, allowing NNQS to serve as high-quality trial wave functions for AFQMC with manageable computational cost. We test the NNQS-AFQMC method on the challenging nitrogen molecule (N 2 ) at stretched geometries. Our results demonstrate that AFQMC with an NNQS trial wave function can attain near-exact total energies, highlighting the potential of AFQMC with NNQS to overcome longstanding challenges in strongly correlated electronic structure calculations. We also outline future research directions for improving this promising methodology.