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:cond-mat/0510791 · Physical Review B

Auxiliary-field quantum Monte Carlo study of TiO and MnO molecules

Wissam A. Al-Saidi, Henry Krakauer, Shiwei Zhang

Published Feb 2, 2006• 38 citations
arXivHTMLPDF

Abstract

Calculations of the binding energy of the transition-metal oxide molecules TiO and MnO are presented, using a recently developed phaseless auxiliary-field quantum Monte Carlo approach. This method maps the interacting many-body problem onto a linear combination of noninteracting problems by a complex Hubbard-Stratonovich transformation, and controls the phase and sign problem with a phaseless approximation relying on a trial wave function. It employs random walks in Slater determinant space to project the ground state of the system, and allows use of much of the same machinery as in standard density functional theory calculations using the plane-wave basis and nonlocal pseudopotentials. The calculations used a single Slater determinant trial wave function obtained from a density functional calculation, with no further optimization. The calculated binding energies are in good agreement with experiment and with recent diffusion Monte Carlo results. Together with previous results for $\mathit{sp}$-bonded systems, the present study indicates that the phaseless auxiliary-field method is a robust and promising approach for the study of correlation effects in real materials.

DOI
10.1103/physrevb.73.075103

Related papers

Phaseless auxiliary-field quantum Monte Carlo calculations with plane waves and pseudopotentials: Applications to atoms and molecules

Malliga Suewattana, Wirawan Purwanto, Shiwei Zhang +2 more

The phaseless auxiliary-field quantum Monte Carlo (AF QMC) method [S. Zhang and H. Krakauer, Phys. Rev. Lett. 90, 136401 (2003)] is used to carry out a systematic study of the dissociation and ionization energies of second-row group 3A--7A atoms and dimers: Al, Si, P, S, and Cl. In addition, the ${\mathrm{P}}_{2}$ dimer is compared to the third-row ${\mathrm{As}}_{2}$ dimer, which is also triply bonded. This method projects the many-body ground state by means of importance-sampled random walks in the space of Slater determinants. The Monte Carlo phase problem, due to the electron-electron Coulomb interaction, is controlled via the phaseless approximation, with a trial wave function $\ensuremath{\mid}{\ensuremath{\Psi}}_{T}⟩$. As in previous calculations, a mean-field single Slater determinant is used as $\ensuremath{\mid}{\ensuremath{\Psi}}_{T}⟩$. The method is formulated in the Hilbert space defined by any chosen one-particle basis. The present calculations use a plane wave basis under periodic boundary conditions with norm-conserving pseudopotentials. Computational details of the plane wave AF QMC method are presented. The isolated systems chosen here allow a systematic study of the various algorithmic issues. We show the accuracy of the plane wave method and discuss its convergence with respect to parameters such as the supercell size and plane wave cutoff. The use of standard norm-conserving pseudopotentials in the many-body AF QMC framework is examined.

arXiv:cond-mat/0702085•Physical Review B•Jun 22, 2007•77 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

Bond breaking with auxiliary-field quantum Monte Carlo

Wissam A. Al-Saidi, Shiwei Zhang, Henry Krakauer

Bond stretching mimics different levels of electron correlation and provides a challenging test bed for approximate many-body computational methods. Using the recently developed phaseless auxiliary-field quantum Monte Carlo (AF QMC) method, we examine bond stretching in the well-studied molecules BH and N(2) and in the H(50) chain. To control the sign/phase problem, the phaseless AF QMC method constrains the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. With single Slater determinants from unrestricted Hartree-Fock as trial wave function, the phaseless AF QMC method generally gives better overall accuracy and a more uniform behavior than the coupled cluster CCSD(T) method in mapping the potential-energy curve. In both BH and N(2), we also study the use of multiple-determinant trial wave functions from multiconfiguration self-consistent-field calculations. The increase in computational cost versus the gain in statistical and systematic accuracy are examined. With such trial wave functions, excellent results are obtained across the entire region between equilibrium and the dissociation limit.

arXiv:0705.2827•The Journal of Chemical Physics•Oct 9, 2007•61 citations

Classical and quantum trial wave functions in auxiliary-field quantum Monte Carlo applied to oxygen allotropes and a CuBr2 model system

Maximilian Amsler, Peter Deglmann, Matthias Degroote +11 more

In this work, we test a recently developed method to enhance classical auxiliary-field quantum Monte Carlo (AFQMC) calculations with quantum computers against examples from chemistry and material science, representative of classes of industry-relevant systems. As molecular test cases, we calculate the energy curve of H4 and the relative energies of ozone and singlet molecular oxygen with respect to triplet molecular oxygen, which is industrially relevant in organic oxidation reactions. We find that trial wave functions beyond single Slater determinants improve the performance of AFQMC and allow it to generate energies close to chemical accuracy compared to full configuration interaction or experimental results. In the field of material science, we study the electronic structure properties of cuprates through the quasi-1D Fermi-Hubbard model derived from CuBr2, where we find that trial wave functions with both significantly larger fidelities and lower energies over a mean-field solution do not necessarily lead to AFQMC results closer to the exact ground state energy.

The Journal of Chemical Physics•Jul 28, 2023•16 citations

Some recent developments in auxiliary-field quantum Monte Carlo for real materials

Hao Shi, Shiwei Zhang

The auxiliary-field quantum Monte Carlo (AFQMC) method is a general numerical method for correlated many-electron systems, which is being increasingly applied in lattice models, atoms, molecules, and solids. Here, we introduce the theory and algorithm of the method specialized for real materials and present several recent developments. We give a systematic exposition of the key steps of AFQMC, closely tracking the framework of a modern software library we are developing. The building of a Monte Carlo Hamiltonian, projecting to the ground state, sampling two-body operators, phaseless approximation, and measuring ground state properties are discussed in detail. An advanced implementation for multi-determinant trial wave functions is described, which dramatically speeds up the algorithm and reduces the memory cost. We propose a self-consistent constraint for real materials, and discuss two flavors for its realization, either by coupling the AFQMC calculation to an effective independent-electron calculation or via the natural orbitals of the computed one-body density matrix.

arXiv:2009.13585•The Journal of Chemical Physics•Jan 11, 2021•61 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
Calculations of the binding energy of the transition-metal oxide molecules TiO and MnO are presented, using a recently developed phaseless auxiliary-field quantum Monte Carlo approach. This method maps the interacting many-body problem onto a linear combination of noninteracting problems by a complex Hubbard-Stratonovich transformation, and controls the phase and sign problem with a phaseless approximation relying on a trial wave function. It employs random walks in Slater determinant space to project the ground state of the system, and allows use of much of the same machinery as in standard density functional theory calculations using the plane-wave basis and nonlocal pseudopotentials. The calculations used a single Slater determinant trial wave function obtained from a density functional calculation, with no further optimization. The calculated binding energies are in good agreement with experiment and with recent diffusion Monte Carlo results. Together with previous results for sp\mathit{sp}sp-bonded systems, the present study indicates that the phaseless auxiliary-field method is a robust and promising approach for the study of correlation effects in real materials.
The phaseless auxiliary-field quantum Monte Carlo (AF QMC) method [S. Zhang and H. Krakauer, Phys. Rev. Lett. 90, 136401 (2003)] is used to carry out a systematic study of the dissociation and ionization energies of second-row group 3A--7A atoms and dimers: Al, Si, P, S, and Cl. In addition, the P2{\mathrm{P}}_{2}P2​ dimer is compared to the third-row As2{\mathrm{As}}_{2}As2​ dimer, which is also triply bonded. This method projects the many-body ground state by means of importance-sampled random walks in the space of Slater determinants. The Monte Carlo phase problem, due to the electron-electron Coulomb interaction, is controlled via the phaseless approximation, with a trial wave function $\ensuremath{\mid}{\ensuremath{\Psi}}_{T}⟩$. As in previous calculations, a mean-field single Slater determinant is used as $\ensuremath{\mid}{\ensuremath{\Psi}}_{T}⟩$. The method is formulated in the Hilbert space defined by any chosen one-particle basis. The present calculations use a plane wave basis under periodic boundary conditions with norm-conserving pseudopotentials. Computational details of the plane wave AF QMC method are presented. The isolated systems chosen here allow a systematic study of the various algorithmic issues. We show the accuracy of the plane wave method and discuss its convergence with respect to parameters such as the supercell size and plane wave cutoff. The use of standard norm-conserving pseudopotentials in the many-body AF QMC framework is examined.
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.
Bond stretching mimics different levels of electron correlation and provides a challenging test bed for approximate many-body computational methods. Using the recently developed phaseless auxiliary-field quantum Monte Carlo (AF QMC) method, we examine bond stretching in the well-studied molecules BH and N(2) and in the H(50) chain. To control the sign/phase problem, the phaseless AF QMC method constrains the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. With single Slater determinants from unrestricted Hartree-Fock as trial wave function, the phaseless AF QMC method generally gives better overall accuracy and a more uniform behavior than the coupled cluster CCSD(T) method in mapping the potential-energy curve. In both BH and N(2), we also study the use of multiple-determinant trial wave functions from multiconfiguration self-consistent-field calculations. The increase in computational cost versus the gain in statistical and systematic accuracy are examined. With such trial wave functions, excellent results are obtained across the entire region between equilibrium and the dissociation limit.
In this work, we test a recently developed method to enhance classical auxiliary-field quantum Monte Carlo (AFQMC) calculations with quantum computers against examples from chemistry and material science, representative of classes of industry-relevant systems. As molecular test cases, we calculate the energy curve of H4 and the relative energies of ozone and singlet molecular oxygen with respect to triplet molecular oxygen, which is industrially relevant in organic oxidation reactions. We find that trial wave functions beyond single Slater determinants improve the performance of AFQMC and allow it to generate energies close to chemical accuracy compared to full configuration interaction or experimental results. In the field of material science, we study the electronic structure properties of cuprates through the quasi-1D Fermi-Hubbard model derived from CuBr2, where we find that trial wave functions with both significantly larger fidelities and lower energies over a mean-field solution do not necessarily lead to AFQMC results closer to the exact ground state energy.
The auxiliary-field quantum Monte Carlo (AFQMC) method is a general numerical method for correlated many-electron systems, which is being increasingly applied in lattice models, atoms, molecules, and solids. Here, we introduce the theory and algorithm of the method specialized for real materials and present several recent developments. We give a systematic exposition of the key steps of AFQMC, closely tracking the framework of a modern software library we are developing. The building of a Monte Carlo Hamiltonian, projecting to the ground state, sampling two-body operators, phaseless approximation, and measuring ground state properties are discussed in detail. An advanced implementation for multi-determinant trial wave functions is described, which dramatically speeds up the algorithm and reduces the memory cost. We propose a self-consistent constraint for real materials, and discuss two flavors for its realization, either by coupling the AFQMC calculation to an effective independent-electron calculation or via the natural orbitals of the computed one-body density matrix.
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.