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

arXiv:1307.2147 · Physical Review B

Symmetry in auxiliary-field quantum Monte Carlo calculations

Hao Shi, Shiwei Zhang

Published Sep 23, 2013• 93 citations
arXivHTMLPDF

Abstract

We show how symmetry properties can be used to greatly increase the accuracy and efficiency in auxiliary-field quantum Monte Carlo (AFQMC) calculations of electronic systems. With the Hubbard model as an example, we study symmetry preservation in two aspects of ground-state AFQMC calculations, the Hubbard-Stratonovich transformation and the form of the trial wave function. It is shown that significant improvement over state-of-the-art calculations can be achieved. In unconstrained calculations, the implementation of symmetry often leads to shorter convergence time and much smaller statistical errors and thereby a substantial reduction of the sign problem. Moreover, certain excited states become possible to calculate which are otherwise beyond reach. In calculations with constraints, the use of symmetry can reduce the systematic error from the constraint. It also allows release-constraint calculations, leading to essentially exact results in many cases. Detailed comparisons are made with exact diagonalization results. Accurate ground-state energies are then presented for larger system sizes in the two-dimensional repulsive Hubbard model.

DOI
10.1103/physrevb.88.125132

Related papers

Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method

Mingpu Qin, Hao Shi, Shiwei Zhang

Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of $U/t$ from 2 to 8, for a range of densities including half-filling and $n=0.3,0.5,0.6,\phantom{\rule{0.16em}{0ex}}0.75$, and $0.875$. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.

arXiv:1605.09421•Physical review. B./Physical review. B•Aug 4, 2016•100 citations

Auxiliary-Field Quantum Monte Carlo Calculations of Excited States and Strongly Correlated Systems

S. Zhang

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.

•Jan 1, 2012•0 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

Taming the Sign Problem in Auxiliary-Field Quantum Monte Carlo Using Accurate 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). 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 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 model Cu 2 O 2 2+ system show that many previously reported isomerization energies differ substantially from the near-exact fp-AFQMC value.

Journal of Chemical Theory and Computation•Jul 7, 2021•43 citations

Symmetry-projected wave functions in quantum Monte Carlo calculations

Hao Shi, Carlos A. Jiménez-Hoyos, R. R. Rodríguez-Guzmán +2 more

We consider symmetry-projected Hartree-Fock trial wave functions in constrained-path Monte Carlo (CPMC) calculations. Previous CPMC calculations have mostly employed Hartree-Fock (HF) trial wave functions, restricted or unrestricted. The symmetry-projected HF approach results in a hierarchy of wave functions with increasing quality: the more symmetries that are broken and restored in a self-consistent manner, the higher the quality of the trial wave function. This hierarchy is approximately maintained in CPMC calculations: the accuracy in the energy increases and the statistical variance decreases when further symmetries are broken and restored. Significant improvement is achieved in CPMC with the best symmetry-projected trial wave functions over those from simple HF. We analyze and quantify the behavior using the two-dimensional repulsive Hubbard model as an example. In the sign-problem-free region, where CPMC can be made exact but a constraint is deliberately imposed here, spin-projected wave functions remove the constraint bias. Away from half filling, spatial symmetry restoration in addition to that of the spin leads to highly accurate results from CPMC. Since the computational cost of symmetry-projected HF trial wave functions in CPMC can be made to scale algebraically with system size, this provides a potentially general approach for accurate calculations in many-fermion systems.

arXiv:1402.0018•Physical Review B•Mar 31, 2014•60 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
We show how symmetry properties can be used to greatly increase the accuracy and efficiency in auxiliary-field quantum Monte Carlo (AFQMC) calculations of electronic systems. With the Hubbard model as an example, we study symmetry preservation in two aspects of ground-state AFQMC calculations, the Hubbard-Stratonovich transformation and the form of the trial wave function. It is shown that significant improvement over state-of-the-art calculations can be achieved. In unconstrained calculations, the implementation of symmetry often leads to shorter convergence time and much smaller statistical errors and thereby a substantial reduction of the sign problem. Moreover, certain excited states become possible to calculate which are otherwise beyond reach. In calculations with constraints, the use of symmetry can reduce the systematic error from the constraint. It also allows release-constraint calculations, leading to essentially exact results in many cases. Detailed comparisons are made with exact diagonalization results. Accurate ground-state energies are then presented for larger system sizes in the two-dimensional repulsive Hubbard model.
Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of U/tU/tU/t from 2 to 8, for a range of densities including half-filling and n=0.3,0.5,0.6,0.75n=0.3,0.5,0.6,\phantom{\rule{0.16em}{0ex}}0.75n=0.3,0.5,0.6,0.75, and 0.8750.8750.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.
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.
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.
We explore different ways of incorporating accurate trial wave functions into free projection auxiliary-field quantum Monte Carlo (fp-AFQMC). 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 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 model Cu 2 O 2 2+ system show that many previously reported isomerization energies differ substantially from the near-exact fp-AFQMC value.
We consider symmetry-projected Hartree-Fock trial wave functions in constrained-path Monte Carlo (CPMC) calculations. Previous CPMC calculations have mostly employed Hartree-Fock (HF) trial wave functions, restricted or unrestricted. The symmetry-projected HF approach results in a hierarchy of wave functions with increasing quality: the more symmetries that are broken and restored in a self-consistent manner, the higher the quality of the trial wave function. This hierarchy is approximately maintained in CPMC calculations: the accuracy in the energy increases and the statistical variance decreases when further symmetries are broken and restored. Significant improvement is achieved in CPMC with the best symmetry-projected trial wave functions over those from simple HF. We analyze and quantify the behavior using the two-dimensional repulsive Hubbard model as an example. In the sign-problem-free region, where CPMC can be made exact but a constraint is deliberately imposed here, spin-projected wave functions remove the constraint bias. Away from half filling, spatial symmetry restoration in addition to that of the spin leads to highly accurate results from CPMC. Since the computational cost of symmetry-projected HF trial wave functions in CPMC can be made to scale algebraically with system size, this provides a potentially general approach for accurate calculations in many-fermion systems.
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.