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

Shifted-contour auxiliary-field Monte Carlo for molecular electronic structure

Naomi Rom, Eyal Fattal, Ashish Kumar Gupta, Emily A. Carter, Daniel Neuhauser

Published Nov 15, 1998• 26 citations
Publisher

Abstract

The shifted-contour auxiliary-field Monte Carlo (SCAFMC) approach has been recently developed by Rom, Charutz and Neuhauser as an extension of the auxiliary-field Monte Carlo (AFMC) method. AFMC replaces the difficult fully interacting electrons problem by an ensemble of simpler problems where the electrons interact with a fluctuating electric field but not with each other. SCAFMC is based on shifting the auxiliary-field contour of integration to pass through the (imaginary) stationary point, leading to numerical stability at long propagation times. The new approach converges to the full CI energy in electronic structure calculations (both ground and low-lying excited states). Here we expand the application of SCAFMC from atomic to molecular problems. First, we calculate ground-state energies of a highly correlated transition-metal system (Cr2) with a moderate (12 orbitals) active space size, and demonstrate that SCAFMC is able to extract the energies accurately. In addition, we use SCAFMC to calculate a C–C bond-stretch energy in ethane with complete active spaces of up to 28 orbitals.

DOI
10.1063/1.477486

Related papers

SHIFTED CONTOUR AUXILIARY FIELD MONTE CARLO

Roi Baer, Daniel Neuhauser

Shifted contour (SC) auxiliary field Monte Carlo (AFMC) for electronic structure is described in considerable algorithmic detail. Several applications demonstrate the high accuracy achievable in molecular electronic structure, including excited states, using the techniques of correlated sampling and variational AFMC. SC-AFMC is also applicable to large Hubbard lattices of strongly correlated electrons. Applications are shown, up to lattices of size 16 ×16 with 202 electrons. Finally, we present a method for exact coarse graining of the AFMC integrand.

Recent advances in computational•Feb 1, 2002•0 citations

Shifted-contour auxiliary field Monte Carlo for ab initio electronic structure: Straddling the sign problem

Roi Baer, Martin Head‐Gordon, Daniel Neuhauser

The auxiliary field Monte Carlo (AFMC) technique has advantages over other ab initio quantum Monte Carlo methods for fermions, as it does not seem to require approximations for alleviating the sign problem and is directly applicable to excited states. Yet, the method is severely limited by a numerical instability, a numerical sign problem, prohibiting application to realistic electronic structure systems. Recently, the shifted contour auxiliary field method (SC-AFMC) was proposed for overcoming this instability. Here we develop a theory for the AFMC stabilization, explaining the success of SC-AFMC. We show that the auxiliary fields can be shifted into the complex plane in a manner that considerably stabilizes the Monte Carlo integration using the exact one-electron density. Practical stabilization can be achieved when an approximate Hartree–Fock density is used, proposing that an overwhelming part of the sign problem is removed by taking proper account of the Fermion mean-field contribution. The theory is demonstrated by application to H2.

The Journal of Chemical Physics•Oct 15, 1998•56 citations

Excited state calculations using phaseless auxiliary-field quantum Monte Carlo: Potential energy curves of low-lying C2 singlet states

Wirawan Purwanto, Shiwei Zhang, Henry Krakauer

We show that the recently developed phaseless auxiliary-field quantum Monte Carlo (AFQMC) method can be used to study excited states, providing an alternative to standard quantum chemistry methods. The phaseless AFQMC approach, whose computational cost scales as M(3)-M(4) with system size M, has been shown to be among the most accurate many-body methods in ground state calculations. For excited states, prevention of collapse into the ground state and control of the Fermion sign/phase problem are accomplished by the approximate phaseless constraint with a trial wave function. Using the challenging C(2) molecule as a test case, we calculate the potential energy curves of the ground and two low-lying singlet excited states. The trial wave function is obtained by truncating complete active space wave functions, with no further optimization. The phaseless AFQMC results using a small basis set are in good agreement with exact full configuration-interaction calculations, while those using large basis sets are in good agreement with experimental spectroscopic constants.

arXiv:0811.4179•The Journal of Chemical Physics•Mar 6, 2009•71 citations

Shell model Monte Carlo approach: the heavy nuclei

Y. Alhassid, L. Fang, Hitoshi Nakada

The auxiliary-field Monte Carlo (AFMC) method, also known as the shell model Monte Carlo (SMMC), enables us to calculate microscopically statistical and collective properties of nuclei in the presence of strong correlations. These calculations can be carried out in shell model spaces that are many orders of magnitude larger than spaces that can be treated with conventional diagonalization methods. A recent major development has been the extension of the AFMC approach to heavy nuclei. Such applications to heavy nuclei have been a major challenge. On the conceptual level, a crucial question is whether a truncated spherical shell model Hamiltonian can describe the proper collectivity observed in heavy nuclei and, in particular, the rotational character of strongly deformed nuclei. On the technical level, the low excitation energies make it necessary to perform calculations down to much lower temperatures. At such low temperatures the propagator becomes ill-conditioned and requires the introduction of stabilization methods.

Journal of Physics Conference Series•Jan 1, 2011•1 citation

Ab initio computations of molecular systems by the auxiliary‐field quantum Monte Carlo method

Mario Motta, Shiwei Zhang

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

arXiv:1711.02242•Wiley Interdisciplinary Reviews Computational Molecular Science•May 25, 2018•173 citations

The performance of phaseless auxiliary-field quantum Monte Carlo on the ground state electronic energy of benzene

Joonho Lee, Fionn D. Malone, David R. Reichman

The ground state electronic energy of benzene is the focus of a recent blind test by Eriksen and co-workers [arXiv:2008.02678]. In this note, we report the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) results on the identical problem. We examine trial wavefunctions of the spin-restricted Hartree-Fock (RHF), and complete active space self-consistent field (CASSCF) with an active space of 6-electron and 6-orbital (i.e., CAS(6,6)) form for the cc-pVDZ basis set. ph-AFQMC+RHF deviates from the value where many methods agreed on (i.e., -863 m$E_h$) by -3.1(3) m$E_h$ whereas ph-AFQMC+CAS(6,6) deviates from the same value by -1.3(4) m$E_h$. In addition to this, we report the frozen core correlation energy of ph-AFQMC+RHF in the cc-pVTZ and cc-pVQZ basis sets as well as their complete basis set limit. Our findings highlight the accuracy, flexibility, and scalability of ph-AFQMC with simple trial wavefunctions.

arXiv:2008.04736•The Journal of Chemical Physics•Sep 22, 2020•32 citations
The shifted-contour auxiliary-field Monte Carlo (SCAFMC) approach has been recently developed by Rom, Charutz and Neuhauser as an extension of the auxiliary-field Monte Carlo (AFMC) method. AFMC replaces the difficult fully interacting electrons problem by an ensemble of simpler problems where the electrons interact with a fluctuating electric field but not with each other. SCAFMC is based on shifting the auxiliary-field contour of integration to pass through the (imaginary) stationary point, leading to numerical stability at long propagation times. The new approach converges to the full CI energy in electronic structure calculations (both ground and low-lying excited states). Here we expand the application of SCAFMC from atomic to molecular problems. First, we calculate ground-state energies of a highly correlated transition-metal system (Cr2) with a moderate (12 orbitals) active space size, and demonstrate that SCAFMC is able to extract the energies accurately. In addition, we use SCAFMC to calculate a C–C bond-stretch energy in ethane with complete active spaces of up to 28 orbitals.
Shifted contour (SC) auxiliary field Monte Carlo (AFMC) for electronic structure is described in considerable algorithmic detail. Several applications demonstrate the high accuracy achievable in molecular electronic structure, including excited states, using the techniques of correlated sampling and variational AFMC. SC-AFMC is also applicable to large Hubbard lattices of strongly correlated electrons. Applications are shown, up to lattices of size 16 ×16 with 202 electrons. Finally, we present a method for exact coarse graining of the AFMC integrand.
The auxiliary field Monte Carlo (AFMC) technique has advantages over other ab initio quantum Monte Carlo methods for fermions, as it does not seem to require approximations for alleviating the sign problem and is directly applicable to excited states. Yet, the method is severely limited by a numerical instability, a numerical sign problem, prohibiting application to realistic electronic structure systems. Recently, the shifted contour auxiliary field method (SC-AFMC) was proposed for overcoming this instability. Here we develop a theory for the AFMC stabilization, explaining the success of SC-AFMC. We show that the auxiliary fields can be shifted into the complex plane in a manner that considerably stabilizes the Monte Carlo integration using the exact one-electron density. Practical stabilization can be achieved when an approximate Hartree–Fock density is used, proposing that an overwhelming part of the sign problem is removed by taking proper account of the Fermion mean-field contribution. The theory is demonstrated by application to H2.
We show that the recently developed phaseless auxiliary-field quantum Monte Carlo (AFQMC) method can be used to study excited states, providing an alternative to standard quantum chemistry methods. The phaseless AFQMC approach, whose computational cost scales as M(3)-M(4) with system size M, has been shown to be among the most accurate many-body methods in ground state calculations. For excited states, prevention of collapse into the ground state and control of the Fermion sign/phase problem are accomplished by the approximate phaseless constraint with a trial wave function. Using the challenging C(2) molecule as a test case, we calculate the potential energy curves of the ground and two low-lying singlet excited states. The trial wave function is obtained by truncating complete active space wave functions, with no further optimization. The phaseless AFQMC results using a small basis set are in good agreement with exact full configuration-interaction calculations, while those using large basis sets are in good agreement with experimental spectroscopic constants.
The auxiliary-field Monte Carlo (AFMC) method, also known as the shell model Monte Carlo (SMMC), enables us to calculate microscopically statistical and collective properties of nuclei in the presence of strong correlations. These calculations can be carried out in shell model spaces that are many orders of magnitude larger than spaces that can be treated with conventional diagonalization methods. A recent major development has been the extension of the AFMC approach to heavy nuclei. Such applications to heavy nuclei have been a major challenge. On the conceptual level, a crucial question is whether a truncated spherical shell model Hamiltonian can describe the proper collectivity observed in heavy nuclei and, in particular, the rotational character of strongly deformed nuclei. On the technical level, the low excitation energies make it necessary to perform calculations down to much lower temperatures. At such low temperatures the propagator becomes ill-conditioned and requires the introduction of stabilization methods.
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
The ground state electronic energy of benzene is the focus of a recent blind test by Eriksen and co-workers [arXiv:2008.02678]. In this note, we report the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) results on the identical problem. We examine trial wavefunctions of the spin-restricted Hartree-Fock (RHF), and complete active space self-consistent field (CASSCF) with an active space of 6-electron and 6-orbital (i.e., CAS(6,6)) form for the cc-pVDZ basis set. ph-AFQMC+RHF deviates from the value where many methods agreed on (i.e., -863 mEhE_hEh​) by -3.1(3) mEhE_hEh​ whereas ph-AFQMC+CAS(6,6) deviates from the same value by -1.3(4) mEhE_hEh​. In addition to this, we report the frozen core correlation energy of ph-AFQMC+RHF in the cc-pVTZ and cc-pVQZ basis sets as well as their complete basis set limit. Our findings highlight the accuracy, flexibility, and scalability of ph-AFQMC with simple trial wavefunctions.