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:1202.4990 · Physical Review B

Convergence of many-body wave-function expansions using a plane-wave basis: From homogeneous electron gas to solid state systems

James J. Shepherd, Andreas Grüneis, George H. Booth, Georg Kresse, Ali Alavi

Published Jul 9, 2012• 129 citations
arXivHTMLPDF

Abstract

Using the finite simulation-cell homogeneous electron gas (HEG) as a model, we investigate the convergence of the correlation energy to the complete-basis-set (CBS) limit in methods utilizing plane-wave wave-function expansions. Simple analytic and numerical results from second-order M\o{}ller-Plesset theory (MP2) suggest a $1/M$ decay of the basis-set incompleteness error where $M$ is the number of plane waves used in the calculation, allowing for straightforward extrapolation to the CBS limit. As we shall show, the choice of basis-set truncation when constructing many-electron wave functions is far from obvious, and here we propose several alternatives based on the momentum transfer vector, which greatly improve the rate of convergence. This is demonstrated for a variety of wave-function methods, from MP2 to coupled-cluster doubles theory and the random-phase approximation plus second-order screened exchange. Finite basis-set energies are presented for these methods and compared with exact benchmarks. A transformation can map the orbitals of a general solid state system onto the HEG plane-wave basis and thereby allow application of these methods to more realistic physical problems. We demonstrate this explicitly for solid and molecular lithium hydride.

DOI
10.1103/physrevb.86.035111

Related papers

Efficient Explicitly Correlated Many-Electron Perturbation Theory for Solids: Application to the Schottky Defect in MgO

Andreas Grüneis

We introduce a novel and efficient explicitly correlated implementation of second-order perturbation theory for solids. The required three-electron integrals are computed directly using a plane wave basis set. We parametrize the employed correlation factors using results previously obtained for a finite uniform electron gas simulation cell. We demonstrate for a range of solids that basis set converged correlation energies, equilibrium volumes, and bulk moduli can be obtained efficiently in this theory using a few ten orbitals per atom. To stretch the capabilities of this novel method we compute the Schottky defect formation energy in MgO, studying systems with 54 atoms in the supercell. We verify the accuracy of the calculated formation energies using the more accurate coupled cluster singles and doubles theory. Furthermore, we discuss other potential applications for the derived and implemented expressions such as an occupied orbital only correlation energy functional.

Physical Review Letters•Aug 6, 2015•46 citations

Coupled Cluster Theory in Materials Science

Igor Ying Zhang, Andreas Grüneis

The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.

arXiv:2004.06424•Frontiers in Materials•Jun 5, 2019•138 citations

Full configuration interaction perspective on the homogeneous electron gas

James J. Shepherd, George H. Booth, Andreas Grüneis +1 more

Highly accurate results for the homogeneous electron gas (HEG) have only been achieved to date within a diffusion Monte Carlo (DMC) framework. Here, we introduce a recently developed stochastic technique, full configuration interaction quantum Monte Carlo (FCIQMC), which samples the exact wave function expanded in plane-wave Slater determinants. Despite the introduction of a basis-set incompleteness error, we obtain a finite-basis energy, which is significantly and variationally lower than any previously published work for the 54-electron HEG at ${r}_{s}$ $=$ 0.5 a.u., in a Hilbert space of ${10}^{108}$ Slater determinants. At this value of ${r}_{s}$, as well as of 1.0 a.u., we remove the remaining basis-set incompleteness error by extrapolation, yielding results comparable to state-of-the-art DMC backflow energies. In doing so, we demonstrate that it is possible to yield highly accurate results with the FCIQMC method in periodic systems.

arXiv:1109.2635•Physical Review B•Feb 7, 2012•147 citations

Combining the Transcorrelated Method with Full Configuration Interaction Quantum Monte Carlo: Application to the Homogeneous Electron Gas

Hongjun Luo, Ali Alavi

We suggest an efficient method to resolve electronic cusps in electronic structure calculations through the use of an effective transcorrelated Hamiltonian. This effective Hamiltonian takes a simple form for plane wave bases, containing up to two-body operators only, and its use incurs almost no additional computational overhead compared to that of the original Hamiltonian. We apply this method in combination with the full configuration interaction quantum Monte Carlo (FCIQMC) method to the homogeneous electron gas. As a projection technique, the non-Hermitian nature of the transcorrelated Hamiltonian does not cause complications or numerical difficulties for FCIQMC. The rate of convergence of the total energy to the complete basis set limit is improved from [Formula: see text] to [Formula: see text], where M is the total number of orbital basis functions.

Journal of Chemical Theory and Computation•Feb 12, 2018•83 citations

Combining the Transcorrelated Method with Full Configuration Interaction Quantum Monte Carlo: Application to the Homogeneous Electron Gas.

Hongjun Luo, Ali Alavi

We suggest an efficient method to resolve electronic cusps in electronic structure calculations by using an effective transcorrelated Hamiltonian. This effective Hamiltonian takes a simple form for plane wave bases, containing up to two-body operators only, and its use incurs almost no additional computational overhead compared to that of the original Hamiltonian. We apply this method in combination with the full configuration interaction quantum Monte Carlo (FCIQMC) method to the homogeneous electron gas. As a projection technique, the non-Hermitian nature of the transcorrelated Hamiltonian does not cause complications or numerical difficulties for FCIQMC. The rate of convergence of the total energy to the complete basis set limit is improved from ${\cal O}(M^{-1})$ to ${\cal O}\left({M^{-5/3}}\right)$, where $M$ is the total number of orbital basis functions.

arXiv:1712.07524•Dec 20, 2017•0 citations

Quantum Solvers for Plane-Wave Hamiltonians: Abridging Virtual Spaces Through the Optimization of Pairwise Correlations

Eric J. Bylaska, Duo Song, Nicholas P. Bauman +3 more

For many-body methods such as MCSCF and CASSCF, in which the number of one-electron orbitals is optimized and independent of the basis set used, there are no problems with using plane-wave basis sets. However, for methods currently used in quantum computing such as select configuration interaction (CI) and coupled cluster (CC) methods, it is necessary to have a virtual space that is able to capture a significant amount of electron-electron correlation in the system. The virtual orbitals in a pseudopotential plane-wave Hartree–Fock calculation, because of Coulomb repulsion, are often scattering states that interact very weakly with the filled orbitals. As a result, very little correlation energy is captured from them. The use of virtual spaces derived from the one-electron operators has also been tried, and while some correlations are captured, the amount is quite low. To overcome these limitations, we have been developing new classes of algorithms to define virtual spaces by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. With these procedures, we have been able to derive virtual spaces, containing only a few orbitals, which are able to capture a significant amount of correlation. The focus in this manuscript is on using these derived basis sets to target full CI (FCI) quality results for H2 on near-term quantum computers. However, the initial results for this approach were promising. We were able to obtain good agreement with FCI/cc-pVTZ results for this system with just 4 virtual orbitals, using both FCI and quantum simulations. The quality of the results using COVOs suggests that it may be possible to use them in other many-body approaches, including coupled cluster and Møller–Plesset perturbation theories, and open up the door to many-body calculations for pseudopotential plane-wave basis set methods.

arXiv:2009.00080•Frontiers in Chemistry•Mar 18, 2021•19 citations
Using the finite simulation-cell homogeneous electron gas (HEG) as a model, we investigate the convergence of the correlation energy to the complete-basis-set (CBS) limit in methods utilizing plane-wave wave-function expansions. Simple analytic and numerical results from second-order M\o{}ller-Plesset theory (MP2) suggest a 1/M1/M1/M decay of the basis-set incompleteness error where MMM is the number of plane waves used in the calculation, allowing for straightforward extrapolation to the CBS limit. As we shall show, the choice of basis-set truncation when constructing many-electron wave functions is far from obvious, and here we propose several alternatives based on the momentum transfer vector, which greatly improve the rate of convergence. This is demonstrated for a variety of wave-function methods, from MP2 to coupled-cluster doubles theory and the random-phase approximation plus second-order screened exchange. Finite basis-set energies are presented for these methods and compared with exact benchmarks. A transformation can map the orbitals of a general solid state system onto the HEG plane-wave basis and thereby allow application of these methods to more realistic physical problems. We demonstrate this explicitly for solid and molecular lithium hydride.
We introduce a novel and efficient explicitly correlated implementation of second-order perturbation theory for solids. The required three-electron integrals are computed directly using a plane wave basis set. We parametrize the employed correlation factors using results previously obtained for a finite uniform electron gas simulation cell. We demonstrate for a range of solids that basis set converged correlation energies, equilibrium volumes, and bulk moduli can be obtained efficiently in this theory using a few ten orbitals per atom. To stretch the capabilities of this novel method we compute the Schottky defect formation energy in MgO, studying systems with 54 atoms in the supercell. We verify the accuracy of the calculated formation energies using the more accurate coupled cluster singles and doubles theory. Furthermore, we discuss other potential applications for the derived and implemented expressions such as an occupied orbital only correlation energy functional.
The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.
Highly accurate results for the homogeneous electron gas (HEG) have only been achieved to date within a diffusion Monte Carlo (DMC) framework. Here, we introduce a recently developed stochastic technique, full configuration interaction quantum Monte Carlo (FCIQMC), which samples the exact wave function expanded in plane-wave Slater determinants. Despite the introduction of a basis-set incompleteness error, we obtain a finite-basis energy, which is significantly and variationally lower than any previously published work for the 54-electron HEG at rs{r}_{s}rs​ === 0.5 a.u., in a Hilbert space of 10108{10}^{108}10108 Slater determinants. At this value of rs{r}_{s}rs​, as well as of 1.0 a.u., we remove the remaining basis-set incompleteness error by extrapolation, yielding results comparable to state-of-the-art DMC backflow energies. In doing so, we demonstrate that it is possible to yield highly accurate results with the FCIQMC method in periodic systems.
We suggest an efficient method to resolve electronic cusps in electronic structure calculations through the use of an effective transcorrelated Hamiltonian. This effective Hamiltonian takes a simple form for plane wave bases, containing up to two-body operators only, and its use incurs almost no additional computational overhead compared to that of the original Hamiltonian. We apply this method in combination with the full configuration interaction quantum Monte Carlo (FCIQMC) method to the homogeneous electron gas. As a projection technique, the non-Hermitian nature of the transcorrelated Hamiltonian does not cause complications or numerical difficulties for FCIQMC. The rate of convergence of the total energy to the complete basis set limit is improved from [Formula: see text] to [Formula: see text], where M is the total number of orbital basis functions.
We suggest an efficient method to resolve electronic cusps in electronic structure calculations by using an effective transcorrelated Hamiltonian. This effective Hamiltonian takes a simple form for plane wave bases, containing up to two-body operators only, and its use incurs almost no additional computational overhead compared to that of the original Hamiltonian. We apply this method in combination with the full configuration interaction quantum Monte Carlo (FCIQMC) method to the homogeneous electron gas. As a projection technique, the non-Hermitian nature of the transcorrelated Hamiltonian does not cause complications or numerical difficulties for FCIQMC. The rate of convergence of the total energy to the complete basis set limit is improved from O(M−1){\cal O}(M^{-1})O(M−1) to O(M−5/3){\cal O}\left({M^{-5/3}}\right)O(M−5/3), where MMM is the total number of orbital basis functions.
For many-body methods such as MCSCF and CASSCF, in which the number of one-electron orbitals is optimized and independent of the basis set used, there are no problems with using plane-wave basis sets. However, for methods currently used in quantum computing such as select configuration interaction (CI) and coupled cluster (CC) methods, it is necessary to have a virtual space that is able to capture a significant amount of electron-electron correlation in the system. The virtual orbitals in a pseudopotential plane-wave Hartree–Fock calculation, because of Coulomb repulsion, are often scattering states that interact very weakly with the filled orbitals. As a result, very little correlation energy is captured from them. The use of virtual spaces derived from the one-electron operators has also been tried, and while some correlations are captured, the amount is quite low. To overcome these limitations, we have been developing new classes of algorithms to define virtual spaces by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. With these procedures, we have been able to derive virtual spaces, containing only a few orbitals, which are able to capture a significant amount of correlation. The focus in this manuscript is on using these derived basis sets to target full CI (FCI) quality results for H2 on near-term quantum computers. However, the initial results for this approach were promising. We were able to obtain good agreement with FCI/cc-pVTZ results for this system with just 4 virtual orbitals, using both FCI and quantum simulations. The quality of the results using COVOs suggests that it may be possible to use them in other many-body approaches, including coupled cluster and Møller–Plesset perturbation theories, and open up the door to many-body calculations for pseudopotential plane-wave basis set methods.