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

Physical Review B

Brueckner-Goldstone quantum Monte Carlo for correlation energies and quasiparticle energy bands of one-dimensional solids

Soohaeng Yoo Willow, Kwang S. Kim, So Hirata

Published Nov 18, 2014• 20 citations
Publisher

Abstract

A quantum Monte Carlo method that combines the second-order many-body perturbation theory and Monte Carlo (MC) integration has been developed for correlation and correlation-corrected (quasiparticle) energy bands of one-dimensional solids. The sum-of-product expressions of correlation energy and self-energy are transformed, with the aid of a Laplace transform, into high-dimensional integrals, which are subject to a highly scalable MC integration with the Metropolis algorithm for importance sampling. The method can compute correlation energies of polyacetylene and polyethylene within a few $\text{m}{E}_{\text{h}}$ and quasiparticle energy bands within a few tenths of an eV. It does not suffer from the fermion sign problem and its description can be systematically improved by raising the perturbation order.

DOI
10.1103/physrevb.90.201110

Related papers

Fast second-order many-body perturbation method for extended systems

So Hirata, Tomomi Shimazaki

An application of second-order many-body perturbation theory to energies and energy bands of polymers is often hindered by the steep polynomial dependence of its computational cost on the number of wave vector sampling points $(K)$ in the Brillouin zone (BZ). It is shown that a Hartree-Fock (HF) calculation with a large value of $K$ (120 in the first BZ) followed by a second-order many-body perturbation calculation with a much smaller value $(K=6)$ can lead to reliable, interpolated correlated energy bands and density of states of a polymer at less than 1% of the computational cost of the conventional approach. Quantitative simulations on photoelectron spectra of trans- and cis-polyacetylenes and polyethylene show that the correlated energy bands and densities of states thus obtained agree quantitatively with the observed and are significant (sometimes qualitative) improvements over the HF results. The energy bands and photoelectron spectra of polydiacetylene are predicted by this method to assist in the interpretation of future high-resolution measurements.

Physical Review B•Aug 26, 2009•45 citations

Stochastic many-body perturbation theory for electron correlation energies

Zhendong Li

Treating electron correlation more accurately and efficiently is at the heart of the development of electronic structure methods. In the present work, we explore the use of stochastic approaches to evaluate high-order electron correlation energies, whose conventional computational scaling is unpleasantly steep, being O(Nn+3) with respect to the system size N and the perturbation order n for the Møller-Plesset (MP) series. To this end, starting from Goldstone’s time-dependent formulation of ab initio many-body perturbation theory (MBPT), we present a reformulation of MBPT, which naturally leads to a Monte Carlo scheme with O(nN2 + n2N + f(n)) scaling at each step, where f(n) is a function of n depending on the specific numerical scheme. Proof-of-concept calculations demonstrate that the proposed quantum Monte Carlo algorithm successfully extends the previous Monte Carlo approaches for MP2 and MP3 to higher orders by overcoming the factorial scaling problem. For the first time, Goldstone’s time-dependent formulation is made useful numerically for electron correlation energies, not only being purely as a theoretical tool.

arXiv:1909.00307•The Journal of Chemical Physics•Dec 27, 2019•12 citations

Electronic structure quantum Monte Carlo

Michal Bajdich, Luboš Mitáš

Electronic structure quantum Monte Carlo Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. The QMC approaches combine analytical insights with stochastic computational techniques for efficient solution of several classes of important many-body problems such as the stationary Schrödinger equation. QMC methods of various flavors have been applied to a great variety of systems spanning continuous and lattice quantum models, molecular and condensed systems, BEC-BCS ultracold condensates, nuclei, etc. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.

arXiv:1008.2369•Acta Physica Slovaca Reviews and Tutorials•Apr 1, 2009•33 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

On the Brillouin‐zone integrations in second‐order many‐body perturbation calculations for extended systems of one‐dimensional periodicity

Tomomi Shimazaki, So Hirata

Abstract The validity and accuracy of various ways of drastically reducing the number of k‐points in the Brillouin zone integrations occurring in second‐order many‐body perturbation calculations of one‐dimensional solids has been investigated. The most promising approximation can recover correlation energies of polyethylene and polyacetylene within 1% of converged values at less than a tenth of usual computational cost. The quasi‐particle energy bands have also been reproduced quantitatively with the same approximation. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009

International Journal of Quantum Chemistry•Apr 30, 2009•26 citations

Stochastic, real-space, imaginary-time evaluation of third-order Feynman–Goldstone diagrams

Soohaeng Yoo Willow, So Hirata

A new, alternative set of interpretation rules of Feynman-Goldstone diagrams for many-body perturbation theory is proposed, which translates diagrams into algebraic expressions suitable for direct Monte Carlo integrations. A vertex of a diagram is associated with a Coulomb interaction (rather than a two-electron integral) and an edge with the trace of a Green's function in real space and imaginary time. With these, 12 diagrams of third-order many-body perturbation (MP3) theory are converted into 20-dimensional integrals, which are then evaluated by a Monte Carlo method. It uses redundant walkers for convergence acceleration and a weight function for importance sampling in conjunction with the Metropolis algorithm. The resulting Monte Carlo MP3 method has low-rank polynomial size dependence of the operation cost, a negligible memory cost, and a naturally parallel computational kernel, while reproducing the correct correlation energies of small molecules within a few mEh after 10(6) Monte Carlo steps.

The Journal of Chemical Physics•Jan 13, 2014•22 citations
A quantum Monte Carlo method that combines the second-order many-body perturbation theory and Monte Carlo (MC) integration has been developed for correlation and correlation-corrected (quasiparticle) energy bands of one-dimensional solids. The sum-of-product expressions of correlation energy and self-energy are transformed, with the aid of a Laplace transform, into high-dimensional integrals, which are subject to a highly scalable MC integration with the Metropolis algorithm for importance sampling. The method can compute correlation energies of polyacetylene and polyethylene within a few mEh\text{m}{E}_{\text{h}}mEh​ and quasiparticle energy bands within a few tenths of an eV. It does not suffer from the fermion sign problem and its description can be systematically improved by raising the perturbation order.
An application of second-order many-body perturbation theory to energies and energy bands of polymers is often hindered by the steep polynomial dependence of its computational cost on the number of wave vector sampling points (K)(K)(K) in the Brillouin zone (BZ). It is shown that a Hartree-Fock (HF) calculation with a large value of KKK (120 in the first BZ) followed by a second-order many-body perturbation calculation with a much smaller value (K=6)(K=6)(K=6) can lead to reliable, interpolated correlated energy bands and density of states of a polymer at less than 1% of the computational cost of the conventional approach. Quantitative simulations on photoelectron spectra of trans- and cis-polyacetylenes and polyethylene show that the correlated energy bands and densities of states thus obtained agree quantitatively with the observed and are significant (sometimes qualitative) improvements over the HF results. The energy bands and photoelectron spectra of polydiacetylene are predicted by this method to assist in the interpretation of future high-resolution measurements.
Treating electron correlation more accurately and efficiently is at the heart of the development of electronic structure methods. In the present work, we explore the use of stochastic approaches to evaluate high-order electron correlation energies, whose conventional computational scaling is unpleasantly steep, being O(Nn+3) with respect to the system size N and the perturbation order n for the Møller-Plesset (MP) series. To this end, starting from Goldstone’s time-dependent formulation of ab initio many-body perturbation theory (MBPT), we present a reformulation of MBPT, which naturally leads to a Monte Carlo scheme with O(nN2 + n2N + f(n)) scaling at each step, where f(n) is a function of n depending on the specific numerical scheme. Proof-of-concept calculations demonstrate that the proposed quantum Monte Carlo algorithm successfully extends the previous Monte Carlo approaches for MP2 and MP3 to higher orders by overcoming the factorial scaling problem. For the first time, Goldstone’s time-dependent formulation is made useful numerically for electron correlation energies, not only being purely as a theoretical tool.
Electronic structure quantum Monte Carlo Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. The QMC approaches combine analytical insights with stochastic computational techniques for efficient solution of several classes of important many-body problems such as the stationary Schrödinger equation. QMC methods of various flavors have been applied to a great variety of systems spanning continuous and lattice quantum models, molecular and condensed systems, BEC-BCS ultracold condensates, nuclei, etc. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.
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.
Abstract The validity and accuracy of various ways of drastically reducing the number of k‐points in the Brillouin zone integrations occurring in second‐order many‐body perturbation calculations of one‐dimensional solids has been investigated. The most promising approximation can recover correlation energies of polyethylene and polyacetylene within 1% of converged values at less than a tenth of usual computational cost. The quasi‐particle energy bands have also been reproduced quantitatively with the same approximation. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009
A new, alternative set of interpretation rules of Feynman-Goldstone diagrams for many-body perturbation theory is proposed, which translates diagrams into algebraic expressions suitable for direct Monte Carlo integrations. A vertex of a diagram is associated with a Coulomb interaction (rather than a two-electron integral) and an edge with the trace of a Green's function in real space and imaginary time. With these, 12 diagrams of third-order many-body perturbation (MP3) theory are converted into 20-dimensional integrals, which are then evaluated by a Monte Carlo method. It uses redundant walkers for convergence acceleration and a weight function for importance sampling in conjunction with the Metropolis algorithm. The resulting Monte Carlo MP3 method has low-rank polynomial size dependence of the operation cost, a negligible memory cost, and a naturally parallel computational kernel, while reproducing the correct correlation energies of small molecules within a few mEh after 10(6) Monte Carlo steps.