Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

Size-extensive auxiliary-field quantum Monte Carlo with perturbative coupled cluster trial wavefunctions

Yichi Zhang, Ankit Mahajan, Yann Damour, Sandeep Sharma

Published Sep 1, 2026• 0 citations
Publisher

Abstract

In this work, we develop a size-extensive Auxiliary-Field Quantum Monte Carlo (AFQMC) approach that scales as O(N5) for local energy evaluation by perturbatively treating the Coupled Cluster Singles and Doubles (CCSD) trial wavefunctions. Comprehensive numerical examinations, spanning from main-group molecules to 3d transition metal complexes, demonstrate that this perturbative treatment introduces negligible bias. For small systems, our method achieves an accuracy and level of noise comparable to AFQMC with configuration interaction singles and doubles trial wavefunctions while outperforming CCSD(T). This size extensivity offers a decisive advantage for large systems, as suggested by the ground state energies of non-interacting monomers and one-dimensional atomic chains. Finally, the numerical simulations of the uniform electron gas provide evidence that, unlike the CCSD(T) method, our new approach does not suffer from infrared divergence in the thermodynamic limit.

DOI
10.1063/5.0341498

Related papers

Size Extensive Auxiliary-Field Quantum Monte Carlo with Perturbative Coupled Cluster Trial Wavefunction

Yichi Zhang, Ankit Mahajan, Yann Damour +1 more

In this work, we develop a size extensive Auxiliary-Field Quantum Monte Carlo (AFQMC) approach that scales as $O(N^5)$ for local energy evaluation by treating the Coupled Cluster Singles and Doubles (CCSD) trial wavefunctions perturbatively. Comprehensive numerical examinations, spanning from main-group molecules to $3d$ transition metal complexes, demonstrate that this perturbative treatment introduces negligible bias. For small systems, our method achieves an accuracy and level of noise comparable to AFQMC with configuration interaction singles and doubles (CISD) trial wavefunctions while outperforming CCSD(T). This size extensivity offers a decisive advantage for large systems, as suggested by the ground state energies of non-interacting monomers and one-dimensional atomic chains. Finally, the numerical simulations of the uniform electron gas (UEG) provide evidence that, unlike the CCSD(T) method, our new approach does not suffer from infrared divergence in the thermodynamic limit (TDL).

arXiv:2605.11257•May 11, 2026•0 citations

Beyond CCSD(T) Accuracy at Lower Scaling with Auxiliary Field Quantum Monte Carlo

Ankit Mahajan, James H. Thorpe, Jo S. Kurian +3 more

We introduce a black-box auxiliary field quantum Monte Carlo (AFQMC) approach to perform highly accurate electronic structure calculations using configuration interaction singles and doubles (CISD) trial states. This method consistently provides more accurate energy estimates than coupled cluster singles and doubles with perturbative triples (CCSD(T)), often regarded as the gold standard in quantum chemistry. This level of precision is achieved at a lower asymptotic computational cost, scaling as O ( N 6 ) compared to the O ( N 7 ) scaling of CCSD(T). We provide numerical evidence supporting these findings through results for challenging main group and transition metal-containing molecules.

arXiv:2410.02885•Journal of Chemical Theory and Computation•Feb 5, 2025•18 citations

Excited states in auxiliary field quantum Monte Carlo

Ankit Mahajan, Sandeep Sharma, Shiwei Zhang +1 more

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.

arXiv:2509.11616•The Journal of Chemical Physics•Jan 23, 2026•0 citations

Systematic improvement of trial states in phaseless auxiliary-field quantum Monte Carlo

Eirik F. Kjønstad, Yann Damour, Sandeep Sharma +1 more

We extend the use of coupled cluster (CC) trial states in the phaseless auxiliary-field quantum Monte Carlo (AFQMC) method beyond single and double excitations to include both triple and quadruple excitations. With this AFQMC/CC hierarchy, we are able to systematically benchmark the method's performance on molecular systems as the quality of the trial is improved. Our results show that the phaseless AFQMC energy improves systematically and is typically significantly more accurate than the energy of the underlying trial state. However, the relative improvement compared to the trial CC energy decreases as we ascend the CC hierarchy. As the CC wavefunction is usually further approximated when used as an AFQMC trial, we also explore the relationship between the components of the CC wavefunction and the resulting AFQMC/CC error. Our results suggest that improving the representation of the CC wave function in the AFQMC trial does not always lower the bias even when it increases the fidelity of the trial with the exact ground state.

arXiv:2510.06486•Oct 7, 2025•1 citation

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

Joonho Lee, Hung Q. Pham, David R. Reichman

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.

arXiv:2208.01280•Journal of Chemical Theory and Computation•Oct 18, 2022•76 citations

Benchmarking quantum trial wavefunctions for phaseless auxiliary-field quantum Monte Carlo

Rod Rofougaran, Neil Mehta, Katherine Klymko +5 more

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.

quant-pharXiv:2605.02056•May 5, 2026•0 citations
In this work, we develop a size-extensive Auxiliary-Field Quantum Monte Carlo (AFQMC) approach that scales as O(N5) for local energy evaluation by perturbatively treating the Coupled Cluster Singles and Doubles (CCSD) trial wavefunctions. Comprehensive numerical examinations, spanning from main-group molecules to 3d transition metal complexes, demonstrate that this perturbative treatment introduces negligible bias. For small systems, our method achieves an accuracy and level of noise comparable to AFQMC with configuration interaction singles and doubles trial wavefunctions while outperforming CCSD(T). This size extensivity offers a decisive advantage for large systems, as suggested by the ground state energies of non-interacting monomers and one-dimensional atomic chains. Finally, the numerical simulations of the uniform electron gas provide evidence that, unlike the CCSD(T) method, our new approach does not suffer from infrared divergence in the thermodynamic limit.
In this work, we develop a size extensive Auxiliary-Field Quantum Monte Carlo (AFQMC) approach that scales as O(N5)O(N^5)O(N5) for local energy evaluation by treating the Coupled Cluster Singles and Doubles (CCSD) trial wavefunctions perturbatively. Comprehensive numerical examinations, spanning from main-group molecules to 3d3d3d transition metal complexes, demonstrate that this perturbative treatment introduces negligible bias. For small systems, our method achieves an accuracy and level of noise comparable to AFQMC with configuration interaction singles and doubles (CISD) trial wavefunctions while outperforming CCSD(T). This size extensivity offers a decisive advantage for large systems, as suggested by the ground state energies of non-interacting monomers and one-dimensional atomic chains. Finally, the numerical simulations of the uniform electron gas (UEG) provide evidence that, unlike the CCSD(T) method, our new approach does not suffer from infrared divergence in the thermodynamic limit (TDL).
We introduce a black-box auxiliary field quantum Monte Carlo (AFQMC) approach to perform highly accurate electronic structure calculations using configuration interaction singles and doubles (CISD) trial states. This method consistently provides more accurate energy estimates than coupled cluster singles and doubles with perturbative triples (CCSD(T)), often regarded as the gold standard in quantum chemistry. This level of precision is achieved at a lower asymptotic computational cost, scaling as O ( N 6 ) compared to the O ( N 7 ) scaling of CCSD(T). We provide numerical evidence supporting these findings through results for challenging main group and transition metal-containing molecules.
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.
We extend the use of coupled cluster (CC) trial states in the phaseless auxiliary-field quantum Monte Carlo (AFQMC) method beyond single and double excitations to include both triple and quadruple excitations. With this AFQMC/CC hierarchy, we are able to systematically benchmark the method's performance on molecular systems as the quality of the trial is improved. Our results show that the phaseless AFQMC energy improves systematically and is typically significantly more accurate than the energy of the underlying trial state. However, the relative improvement compared to the trial CC energy decreases as we ascend the CC hierarchy. As the CC wavefunction is usually further approximated when used as an AFQMC trial, we also explore the relationship between the components of the CC wavefunction and the resulting AFQMC/CC error. Our results suggest that improving the representation of the CC wave function in the AFQMC trial does not always lower the bias even when it increases the fidelity of the trial with the exact ground state.
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.