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

Journal of Chemical Theory and Computation

Self-Refinement of Auxiliary-Field Quantum Monte Carlo via Non-Orthogonal Configuration Interaction

Zoran Sukurma, Martin Schlipf, Georg Kresse

Published Apr 28, 2025• 7 citations
Publisher

Abstract

For optimal accuracy, auxiliary-field quantum Monte Carlo (AFQMC) requires trial states consisting of multiple Slater determinants. We develop an efficient algorithm to select the determinants from an AFQMC random walk eliminating the need for other methods. When determinants contribute significantly to the nonorthogonal configuration interaction energy, we include them in the trial state. These refined trial wave functions significantly reduce the phaseless bias and sampling variance of the local energy estimator. With 100 to 200 determinants, we lower the error of AFQMC by up to a factor of 10 for second-row elements that are not accurately described with a Hartree-Fock trial wave function. For the HEAT set, we improve the average error to within chemical accuracy. For benzene, the largest studied system, we reduce AFQMC error by 80% with 214 Slater determinants and find a 10-fold increase of the time to solution. We show that phaseless errors prevail in systems with static correlation or strong spin contamination. For such systems, improved trial states enable stable free-projection AFQMC calculations, achieving chemical accuracy even in the strongly correlated regime.

DOI
10.1021/acs.jctc.5c00127

Related papers

Self-Refinement of Auxiliary-Field Quantum Monte Carlo via Non-Orthogonal Configuration Interaction

Zoran Sukurma, Martin Schlipf, Georg Kresse

For optimal accuracy, auxiliary-field quantum Monte Carlo (AFQMC) requires trial states consisting of multiple Slater determinants. We develop an efficient algorithm to select the determinants from an AFQMC random walk eliminating the need for other methods. When determinants contribute significantly to the non-orthogonal configuration interaction energy, we include them in the trial state. These refined trial wave functions significantly reduce the phaseless bias and sampling variance of the local energy estimator. With 100 to 200 determinants, we lower the error of AFQMC by up to a factor of 10 for second row elements that are not accurately described with a Hartree-Fock trial wave function. For the HEAT set, we improve the average error to within the chemical accuracy. For benzene, the largest studied system, we reduce AFQMC error by 80% with 214 Slater determinants and find a 10-fold increase of the time to solution. We show that the remaining error of the method prevails in systems with static correlation or strong spin contamination.

arXiv:2501.12765•Jan 22, 2025•0 citations

Selected configuration interaction wave functions in phaseless auxiliary field quantum Monte Carlo

Ankit Mahajan, Joonho Lee, Sandeep Sharma

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.

arXiv:2201.10458•The Journal of Chemical Physics•Apr 14, 2022•40 citations

Taming the sign problem in auxiliary field quantum Monte Carlo using accurate trial wave functions

Ankit Mahajan, Sandeep Sharma

We explore different ways of incorporating accurate trial wave functions into free projection auxiliary field quantum Monte Carlo (fp-AFQMC). Trial states employed include coupled cluster singles and doubles, multi-Slater, and symmetry projected mean-field wave functions. We adapt a recently proposed fast multi-Slater local energy evaluation algorithm for fp-AFQMC, making the use of long expansions from selected configuration interaction methods feasible. We demonstrate how these trial wave functions serve to mitigate the sign problem and accelerate convergence in quantum chemical problems, allowing the application of fp-AFQMC to systems of substantial sizes. Our calculations on the widely studied $\text{Cu}_2\text{O}_2^{2+}$ model system show that many previously reported isomerization energies differ substantially from the near-exact fp-AFQMC value.

arXiv:2104.06597•Apr 14, 2021•0 citations

Benchmark Phaseless Auxiliary-Field Quantum Monte Carlo Method for Small Molecules

Zoran Sukurma, Martin Schlipf, Moritz Humer +2 more

High Resolution Image Download MS PowerPoint Slide We report a scalable Fortran implementation of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) and demonstrate its excellent performance and beneficial scaling with respect to system size. Furthermore, we investigate modifications of the phaseless approximation that can help to reduce the overcorrelation problems common to the ph-AFQMC. We apply the method to the 26 molecules in the HEAT set, the benzene molecule, and water clusters. We observe a mean absolute deviation of the total energy of 1.15 kcal/mol for the molecules in the HEAT set, close to chemical accuracy. For the benzene molecule, the modified algorithm despite using a single-Slater-determinant trial wavefunction yields the same accuracy as the original phaseless scheme with 400 Slater determinants. Despite these improvements, we find systematic errors for the CN, CO 2, and O 2 molecules that need to be addressed with more accurate trial wavefunctions. For water clusters, we find that the ph-AFQMC yields excellent binding energies that differ from CCSD(T) by typically less than 0.5 kcal/mol.

arXiv:2303.04256•Journal of Chemical Theory and Computation•Jul 20, 2023•17 citations

Implementing advanced trial wave functions in fermion quantum Monte Carlo via stochastic sampling

Zhi-Yu Xiao, Tao Xiang, Zixiang Lu +2 more

We introduce an efficient approach to implement correlated many-body trial wave functions in auxiliary-field quantum Monte Carlo (AFQMC). To control the sign/phase problem in AFQMC, a constraint is derived from an exact gauge condition but is typically imposed approximately through a trial wave function or trial density matrix, whose quality can affect the accuracy of the method. Furthermore, the trial wave function can also affect the efficiency through importance sampling. The most natural form of the trial wave function has been single Slater determinants or their linear combinations. More sophisticated forms, for example, with the inclusion of a Jastrow factor or other explicit correlations, have been challenging to use, and their implementation is often assumed to require a quantum computer. In this work, we demonstrate that a large class of correlated wave functions, written in the general form of multi-dimensional integrals over hidden or auxiliary variables times Slater determinants, can be implemented as a trial wave function by coupling the random walkers to a generalized Metropolis sampling. We discuss the fidelity of AFQMC with stochastically sampled trial wave functions, which are relevant to both quantum and classical algorithms. We illustrate the method and show that an efficient implementation can be achieved, which preserves the computational scaling of AFQMC. We test our method in molecules under bond stretching and in transition metal diatomics. Significant improvements are seen in both accuracy and efficiency over typical trial wave functions, and the method yields total ground-state energies systematically within chemical accuracy. The method can be useful for incorporating other advanced wave functions, for example, neural quantum state wave functions optimized from machine learning techniques, or for other forms of fermion quantum Monte Carlo.

arXiv:2505.18519•The Journal of Chemical Physics•Oct 22, 2025•5 citations

Selecting Optimal Unrestricted Hartree–Fock Trial Wave Functions for Phaseless Auxiliary-Field Quantum Monte Carlo: Accuracy and Limitations in Modeling Three Iron–Sulfur Clusters

Don Danilov, Brad Ganoe, Leon Otis +3 more

Phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) has emerged as a promising electronic structure method for correlated electronic systems. However, the quality of its predictions depends critically on the choice of trial wave function, and it is not obvious how to make an optimal choice, especially for strongly correlated states of large systems. Mean-field wave functions are compelling trial wave function candidates, as they map directly to chemical concepts and can be obtained withO(N4) cost. Yet in the strongly correlated regime, one faces a symmetry dilemma and the existence of multiple nearly degenerate solutions. In this work, we investigate active space models of [2Fe-2S]2+, mixed-valent [4Fe-4S]2+, and [4Fe-4S]4+ and explore the sensitivity of ph-AFQMC to the choice of unrestricted Hartree-Fock trial wave function. We find that surprisingly accurate ground-state energies for these systems can be obtained when trial selection is guided by chemical properties and physical symmetries, rather than the variational energy. However, in all cases, we find a rapidly decaying overlap between the stochastic wave function and the UHF trial, indicating that the trials are acting as suboptimal importance functions. By analogy to a similar situation in the stretched helium dimer cation, we show how this sampling bias pushes ph-AFQMC with UHF trials toward artificially negative energies, which evidently can be compensated for by the phaseless bias in certain cases.

Journal of Chemical Theory and Computation•Aug 6, 2026•1 citation
For optimal accuracy, auxiliary-field quantum Monte Carlo (AFQMC) requires trial states consisting of multiple Slater determinants. We develop an efficient algorithm to select the determinants from an AFQMC random walk eliminating the need for other methods. When determinants contribute significantly to the nonorthogonal configuration interaction energy, we include them in the trial state. These refined trial wave functions significantly reduce the phaseless bias and sampling variance of the local energy estimator. With 100 to 200 determinants, we lower the error of AFQMC by up to a factor of 10 for second-row elements that are not accurately described with a Hartree-Fock trial wave function. For the HEAT set, we improve the average error to within chemical accuracy. For benzene, the largest studied system, we reduce AFQMC error by 80% with 214 Slater determinants and find a 10-fold increase of the time to solution. We show that phaseless errors prevail in systems with static correlation or strong spin contamination. For such systems, improved trial states enable stable free-projection AFQMC calculations, achieving chemical accuracy even in the strongly correlated regime.
For optimal accuracy, auxiliary-field quantum Monte Carlo (AFQMC) requires trial states consisting of multiple Slater determinants. We develop an efficient algorithm to select the determinants from an AFQMC random walk eliminating the need for other methods. When determinants contribute significantly to the non-orthogonal configuration interaction energy, we include them in the trial state. These refined trial wave functions significantly reduce the phaseless bias and sampling variance of the local energy estimator. With 100 to 200 determinants, we lower the error of AFQMC by up to a factor of 10 for second row elements that are not accurately described with a Hartree-Fock trial wave function. For the HEAT set, we improve the average error to within the chemical accuracy. For benzene, the largest studied system, we reduce AFQMC error by 80% with 214 Slater determinants and find a 10-fold increase of the time to solution. We show that the remaining error of the method prevails in systems with static correlation or strong spin contamination.
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.
We explore different ways of incorporating accurate trial wave functions into free projection auxiliary field quantum Monte Carlo (fp-AFQMC). Trial states employed include coupled cluster singles and doubles, multi-Slater, and symmetry projected mean-field wave functions. We adapt a recently proposed fast multi-Slater local energy evaluation algorithm for fp-AFQMC, making the use of long expansions from selected configuration interaction methods feasible. We demonstrate how these trial wave functions serve to mitigate the sign problem and accelerate convergence in quantum chemical problems, allowing the application of fp-AFQMC to systems of substantial sizes. Our calculations on the widely studied Cu2O22+\text{Cu}_2\text{O}_2^{2+}Cu2​O22+​ model system show that many previously reported isomerization energies differ substantially from the near-exact fp-AFQMC value.
High Resolution Image Download MS PowerPoint Slide We report a scalable Fortran implementation of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) and demonstrate its excellent performance and beneficial scaling with respect to system size. Furthermore, we investigate modifications of the phaseless approximation that can help to reduce the overcorrelation problems common to the ph-AFQMC. We apply the method to the 26 molecules in the HEAT set, the benzene molecule, and water clusters. We observe a mean absolute deviation of the total energy of 1.15 kcal/mol for the molecules in the HEAT set, close to chemical accuracy. For the benzene molecule, the modified algorithm despite using a single-Slater-determinant trial wavefunction yields the same accuracy as the original phaseless scheme with 400 Slater determinants. Despite these improvements, we find systematic errors for the CN, CO 2, and O 2 molecules that need to be addressed with more accurate trial wavefunctions. For water clusters, we find that the ph-AFQMC yields excellent binding energies that differ from CCSD(T) by typically less than 0.5 kcal/mol.
We introduce an efficient approach to implement correlated many-body trial wave functions in auxiliary-field quantum Monte Carlo (AFQMC). To control the sign/phase problem in AFQMC, a constraint is derived from an exact gauge condition but is typically imposed approximately through a trial wave function or trial density matrix, whose quality can affect the accuracy of the method. Furthermore, the trial wave function can also affect the efficiency through importance sampling. The most natural form of the trial wave function has been single Slater determinants or their linear combinations. More sophisticated forms, for example, with the inclusion of a Jastrow factor or other explicit correlations, have been challenging to use, and their implementation is often assumed to require a quantum computer. In this work, we demonstrate that a large class of correlated wave functions, written in the general form of multi-dimensional integrals over hidden or auxiliary variables times Slater determinants, can be implemented as a trial wave function by coupling the random walkers to a generalized Metropolis sampling. We discuss the fidelity of AFQMC with stochastically sampled trial wave functions, which are relevant to both quantum and classical algorithms. We illustrate the method and show that an efficient implementation can be achieved, which preserves the computational scaling of AFQMC. We test our method in molecules under bond stretching and in transition metal diatomics. Significant improvements are seen in both accuracy and efficiency over typical trial wave functions, and the method yields total ground-state energies systematically within chemical accuracy. The method can be useful for incorporating other advanced wave functions, for example, neural quantum state wave functions optimized from machine learning techniques, or for other forms of fermion quantum Monte Carlo.
Phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) has emerged as a promising electronic structure method for correlated electronic systems. However, the quality of its predictions depends critically on the choice of trial wave function, and it is not obvious how to make an optimal choice, especially for strongly correlated states of large systems. Mean-field wave functions are compelling trial wave function candidates, as they map directly to chemical concepts and can be obtained withO(N4) cost. Yet in the strongly correlated regime, one faces a symmetry dilemma and the existence of multiple nearly degenerate solutions. In this work, we investigate active space models of [2Fe-2S]2+, mixed-valent [4Fe-4S]2+, and [4Fe-4S]4+ and explore the sensitivity of ph-AFQMC to the choice of unrestricted Hartree-Fock trial wave function. We find that surprisingly accurate ground-state energies for these systems can be obtained when trial selection is guided by chemical properties and physical symmetries, rather than the variational energy. However, in all cases, we find a rapidly decaying overlap between the stochastic wave function and the UHF trial, indicating that the trials are acting as suboptimal importance functions. By analogy to a similar situation in the stretched helium dimer cation, we show how this sampling bias pushes ph-AFQMC with UHF trials toward artificially negative energies, which evidently can be compensated for by the phaseless bias in certain cases.