We present an approach that uses the doubly occupied configuration interaction (DOCI) wave function as the trial wave function in phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC). DOCI is a seniority-zero method focused on electron pairs. Although DOCI considers much fewer electron configurations than the complete active space (CAS) configuration interaction method, it efficiently captures the static correlation, while the consequent ph-AFQMC recovers the dynamical correlation across all orbitals. We also explore an orbital-optimized version (OO-DOCI) to further improve accuracy. We test this approach on several chemical systems, including single O-H bond breaking in water and polymer additives. In these cases, OO-DOCI-AFQMC closely matches CAS-based ph-AFQMC and even outperforms coupled-cluster singles, doubles, and perturbative triples. However, for strongly correlated systems, such as the carbon dimer and multibond dissociation in hydrogen systems and water, the method's accuracy drops. This suggests that seniority-zero space models may be insufficient as trial wave functions in ph-AFQMC for strongly correlated systems, suggesting the need for trial wave functions in an extended space. Despite such a limitation, our study demonstrates that DOCI- and OO-DOCI-based ph-AFQMC can reduce the steep cost of CAS approaches, offering a path to accurate multireference calculations for larger, more complex systems.
We present an approach that uses the doubly occupied configuration interaction (DOCI) wave function as the trial wave function in phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC). DOCI is a seniority-zero method focused on electron pairs. Although DOCI considers much fewer electron configurations than the complete active space (CAS) configuration interaction method, it efficiently captures the static correlation, while the consequent ph-AFQMC recovers the dynamical correlation across all orbitals. We also explore an orbital-optimized version (OO-DOCI) to further improve accuracy. We test this approach on several chemical systems, including single O-H bond breaking in water and polymer additives. In these cases, OO-DOCI-AFQMC closely matches CAS-based ph-AFQMC and even outperforms coupled-cluster singles, doubles, and perturbative triples. However, for strongly correlated systems, such as the carbon dimer and multibond dissociation in hydrogen systems and water, the method's accuracy drops. This suggests that seniority-zero space models may be insufficient as trial wave functions in ph-AFQMC for strongly correlated systems, suggesting the need for trial wave functions in an extended space. Despite such a limitation, our study demonstrates that DOCI- and OO-DOCI-based ph-AFQMC can reduce the steep cost of CAS approaches, offering a path to accurate multireference calculations for larger, more complex systems.
We present an approach that uses the doubly occupied configuration interaction (DOCI) wave function as the trial wave function in phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC). DOCI is a seniority-zero method focused on electron pairs. Although DOCI considers much fewer electron configurations than the complete active space (CAS) configuration interaction method, it efficiently captures the static correlation, while the consequent ph-AFQMC recovers the dynamical correlation across all orbitals. We also explore an orbital-optimized version (OO-DOCI) to further improve accuracy. We test this approach on several chemical systems, including single O-H bond breaking in water and polymer additives. In these cases, OO-DOCI-AFQMC closely matches CAS-based ph-AFQMC and even outperforms coupled-cluster singles, doubles, and perturbative triples. However, for strongly correlated systems, such as the carbon dimer and multi-bond dissociation in hydrogen systems and water, the method's accuracy drops. This suggests that seniority-zero space models may be insufficient as trial wave functions in ph-AFQMC for strongly correlated systems, suggesting the need for trial wave functions in an extended space. Despite such a limitation, our study demonstrates that DOCI- and OO-DOCI-based ph-AFQMC can reduce the steep cost of CAS approaches, offering a path to accurate multi-reference calculations for larger, more complex systems.
We propose doubly occupied configuration interaction-quantum selected configuration interaction (DOCI-QSCI), which samples from the seniority-zero space. While the use of this space effectively doubles the qubit budget, equaling the number of spatial orbitals, this sector restriction can compromise quantitative accuracy. To compensate for this, we expand sampled bitstrings via their Cartesian product into a larger space that includes seniority-breaking determinants. The resulting wave function is also proposed using the trial state in phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) to recover dynamical correlations across the full orbital space (DOCI-QSCI-AFQMC). We evaluate the proposed methods on the H6 chain, N2 dissociation, and the addition of singlet O2 to a BODIPY dye. For the H6 chain, DOCI-QSCI-AFQMC reproduces the accuracy of the level of the complete-active-space counterpart with the quantum device ibm kobe. For N2 and BODIPY-O2, with (14e, 28o) and up to (20e, 20o) active spaces, it yields reasonable results, whereas single-reference CCSD(T) fails qualitatively. These results demonstrate that the DOCI-QSCI doubles the orbital space accessible to conventional QSCI and subsequent ph-AFQMC post-processing delivers reasonably high accuracy.
We present efficient algorithms for using selected configuration interaction (sCI) trial wave functions in phaseless auxiliary field quantum Monte Carlo (ph-AFQMC). These advances, geared toward optimizing computational performance for longer configuration interaction expansions, allow us to use up to a million configurations in the trial state for ph-AFQMC. In one example, we found the cost of ph-AFQMC per sample to increase only by a factor of about 3 for a calculation with 104 configurations compared to that with a single one, demonstrating the tiny computational overhead due to a longer expansion. This favorable scaling allows us to study the systematic convergence of the phaseless bias in auxiliary field quantum Monte Carlo calculations with an increasing number of configurations and provides a means to gauge the accuracy of ph-AFQMC with other trial states. We also show how the scalability issues of sCI trial states for large system sizes could be mitigated by restricting them to a moderately sized orbital active space and leveraging the near-cancellation of out of active space phaseless errors.
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.
The phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) method is a stochastic imaginary-time projection method for computing ground-state properties of strongly correlated quantum systems, whose accuracy depends critically on the trial wavefunction. We present a benchmarking study of trial wavefunctions prepared using parameterized quantum circuits, spanning unitary coupled cluster (UCC), Hamiltonian-informed, Jastrow-inspired, and adaptively constructed ansatzes, and evaluate their accuracy, expressibility, and resource requirements within the quantum-computing-enhanced AFQMC (QC-AFQMC) framework. Trial states are prepared with GPU-accelerated statevector simulations using CUDA-Q and tested on stretched linear hydrogen chains and nine closed-shell molecules at equilibrium with varying static correlation. We find that at similar parameter counts, ph-AFQMC largely removes differences between families for weakly and moderately correlated systems: every trial yields a ph-AFQMC energy within chemical accuracy despite trial energy errors spanning nearly three orders of magnitude, whereas ansatz choice becomes consequential as static correlation grows. Our results indicate that trial variational energy is not always a reliable indicator of ph-AFQMC quality, and that adding variational parameters can lower the trial energy while degrading the ph-AFQMC energy. Across our benchmark systems, trial fidelity ranks trials more reliably than variational energy, and for spin-contaminated trials, the singlet-projected fidelity is more reliable than unprojected fidelity. The UCCSD ansatz initialized from unrestricted CCSD amplitudes provides a robust and inexpensive-to-optimize baseline across benchmark systems, with circuit depth as its principal limitation, while adaptive and compact layered ansatzes remain competitive in accuracy under strong correlation with substantially shallower circuits.
Generalized Hartree Fock (GHF) is a long-established electronic structure method that can lower the energy (compared to spin-restricted variants) by breaking physical wavefunction symmetries, namely S^2 and Sz. After an exposition of GHF theory, we assess the use of GHF trial wavefunctions in phaseless auxiliary field quantum Monte Carlo (ph-AFQMC-G) calculations of strongly correlated molecular systems including symmetrically stretched hydrogen rings, carbon dioxide, and dioxygen. We find comparable or better accuracy than CCSD(T) with unrestricted HF and GHF references, and consistently smooth dissociation curves -- a remarkable result given the relative scalability of ph-AFQMC-G to larger system sizes. The present exploration of model strongly correlated systems marks a promising starting point for future studies of more chemically-relevant molecules, and demonstrates that ph-AFQMC-G provides a highly accurate (and, in contrast to active-space-based trials, relatively black box and always size-consistent) description of challenging systems exhibiting, e.g., antiferromagnetic coupling and/or geometric spin frustration.