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

Towards numerically accurate many-body perturbation theory: Short-range correlation effects

Andris Guļāns

Published Oct 28, 2014• 28 citations
Publisher

Abstract

The example of the uniform electron gas is used for showing that the short-range electron correlation is difficult to handle numerically, while it noticeably contributes to the self-energy. Nonetheless, in condensed-matter applications studied with advanced methods, such as the GW and random-phase approximations, it is common to neglect contributions due to high-momentum (large q) transfers. Then, the short-range correlation is poorly described, which leads to inaccurate correlation energies and quasiparticle spectra. To circumvent this problem, an accurate extrapolation scheme is proposed. It is based on an analytical derivation for the uniform electron gas presented in this paper, and it provides an explanation why accurate GW quasiparticle spectra are easy to obtain for some compounds and very difficult for others.

DOI
10.1063/1.4900447

Related papers

Many-body perturbation theory: The GW approximation

Christoph Friedrich, Arno Schindlmayr

In this lecture we present many-body perturbation theory as a method to determine quasiparticle excitations in solids, especially electronic band structures, accurately from first principles.The main ingredient is the electronic self-energy that, in principle, contains all many-body exchange and correlation effects beyond the Hartree potential.As its exact mathematical expression is unknown, approximations must be used in practical calculations.The approximation is obtained using a systematic algebraic approach on the basis of Green function techniques.It constitutes an expansion of the self-energy up to linear order in the screened Coulomb potential, which describes the interaction between the quasiparticles and includes dynamic screening through the creation of exchange-correlation holes around the bare particles.The implementation of the approximation relies on a perturbative treatment starting from density functional theory.Besides a detailed mathematical discussion we focus on the underlying physical concepts and show some illustrative applications.

•Jan 1, 2006•19 citations

Correlation beyond the random phase approximation: A consistent many-body perturbation theory approach

F. Bechstedt

The correlation energy of an inhomogeneous electron gas is ruled by the density correlation function $L$. It obeys a Bethe-Salpeter equation (BSE), whose kernel is the functional derivative of the particle self-energy. It is solved on equal footing for two-particle and single-particle Green's functions within Hedin's GW approximation. According to the Hartree and exchange-correlation contributions to the kernel, the correlation expression decays into random-phase-approximation (RPA) -like and screened-exchange-like parts. The actual treatment of the screened potential $W$ in the BSE for $L$ leads to different second-order screened exchange corrections (SOSEX). Including correlation in the kernel even yields approximations beyond SOSEX. The progress beyond RPA is illustrated for spinless and spin-polarized homogeneous electron gases, indicating that the screening treatment is substantial.

Physical review. B./Physical review. B•Jun 15, 2018•5 citations

Perturbation calculation of the uniform electron gas with a transcorrelated Hamiltonian

Hongjun Luo, Ali Alavi

With a transcorrelated Hamiltonian, we perform a many body perturbation calculation on the uniform electron gas in the high density regime. By using a correlation factor optimized for a single determinant Jastrow ansatz, the second order correlation energy is calculated as 1−ln⁡2π2ln(rs)−0.05075. This already reproduces the exact logarithmic term of the random phase approximation (RPA) result, while the constant term is roughly 7% larger than the RPA one. The close agreement with the RPA method demonstrates that the transcorrelated method offers a viable and potentially efficient method for treating metallic systems.

arXiv:2206.03228•The Journal of Chemical Physics•Jul 18, 2022•5 citations

Core electrons and self-consistency in the GW approximation from a PAW perspective

Matteo Giantomassi

Density functional theory (DFT) performs reasonably well for the determination of structural properties of many materials, but fails to predict electronic band gap values accurately. \nSuch a failure of DFT is not unexpected, since there exists no formal justification for interpreting the DFT eigenvalues as addition or removal energies of the many-body system (quasiparticle energies).\nAn alternative approach to the study of exchange-correlation effects in many-particle systems is provided by many-body perturbation theory (MBPT), which defines a rigorous approach to description of excited-state properties based on the Green's function formalism and the concept of quasiparticle electrons. Within MBPT, one can calculate the quasiparticle (QP) energies and the QP amplitudes by approximately solving, within the so-called GW approximation, the set of coupled integro-differential equations proposed by Hedin in 1965.\nThe GW method generally yields significantly better values for QP energies with respect to DFT as it accounts for the dynamic many-body effects in the electron-electron interaction going beyond the mean-field approximation of independent Kohn-Sham particles. The drawback of such an involved approach, however, is its large computational cost, which is mainly due to the evaluation of dielectric matrices, their inversion, and the solution of a non-Hermitian and nonlinear eigenvalue problem. \nThe first part of this PhD research is devoted to the implementation of an efficient and scalable coarse-grained parallel algorithms and the development of state-of-the-art methods to solve the GW equations.\nAll these techniques are then applied to the study of the quasiparticle band structures of SiO2 in the alpha quartz crystalline structure. The effects of the different approximations involved in the theory, the influence of self-consistency, and a systematic analysis of the reliability of the different plasmon-pole models are presented and discussed. \nThe second part of this PhD work is dedicated to the implementation of a methodological approach to the solution of the GW equations based on the so-called projector augmented wave method (PAW) proposed by Blochl in 1994. This new approach permits to remove many of the limitations intrinsic to the use of the pseudopotential technique widely used for ab-initio calculations, allowing one to reach GW results closer to an all-electron implementation while still maintaining computational efficiency, and flexibility. This allows us to achieve a coherent implementation of many different algorithms which enables a detailed comparison of GW calculations not yet performed until now. Results for the quasiparticle band structure and optical spectra of prototype s- and p-compounds are discussed and compared to recent studies reported in the scientific literature (when available).

•Jan 1, 2009•2 citations

Configuration interaction approach for the computation of the electronic self-energy

Y. Pavlyukh, Wolfgang Hübner

In order to properly understand the utility of many-body perturbation theory as applied to finite systems, we use the configuration interaction approach to compute the electronic self-energy. The validity of the commonly used $GW$ approximation from many-body perturbation theory is tested by comparing the self-energy explicitly computed according to the diagrammatic expansions to that from the inversion of the Dyson equation. It is constructed as a functional of the interacting Green function $G$ and screened Coulomb interaction $W$. $W$ is explicitly computed through the diagonalization of the many-body Hamiltonian and, thus, takes polarization effects into account beyond the random phase approximation. The $\mathrm{Na}_{9}{}^{+}$ cluster and the dissociation of the ${\mathrm{C}}_{2}$ molecule are discussed as examples. We find that the $GW$ approximation yields accurate results for weakly correlated systems with a predominantly single-determinant ground state such as $\mathrm{Na}_{9}{}^{+}$ clusters. The ${\mathrm{C}}_{2}$ molecule is a pathological system with multideterminantal ground state where the $GW$ approach ceases to be valid.

Physical Review B•May 30, 2007•17 citations

Short-range correlations in the electron gas

Kevin S. Bedell, G. E. Brown

The effects of short-range correlations on the properties of the electron gas are studied using an approach phrased in the polarization potential language of Pines and Nozi\`eres. It was argued by Lowy and Brown that two-body terms make up the dominant contribution to the effective interaction for large momentum transfers. In the region of small-momentum transfers it is known that the random-phase approximation provides an adequate picture of the electron gas. With these two regions of large- and small-momentum transfers in mind, a polarization propagator is constructed to interpolate between them. A number of sum rules, among them the $f$-sum rule, involving the imaginary part of the polarization propagator, were found to be satisfied. The failure of our theory to satisfy the compressibility sum rule is discussed in detail. We present here the results for the local $t$-matrix, the dynamic and static form factors, and the pair correlation function. That our pair-correlation function remains positive and is essentially the same as in the paper of Lowy and Brown indicates that the physics and not the formalism is responsible for the results we obtain.

Physical review. B, Condensed matter•Jun 15, 1978•38 citations
The example of the uniform electron gas is used for showing that the short-range electron correlation is difficult to handle numerically, while it noticeably contributes to the self-energy. Nonetheless, in condensed-matter applications studied with advanced methods, such as the GW and random-phase approximations, it is common to neglect contributions due to high-momentum (large q) transfers. Then, the short-range correlation is poorly described, which leads to inaccurate correlation energies and quasiparticle spectra. To circumvent this problem, an accurate extrapolation scheme is proposed. It is based on an analytical derivation for the uniform electron gas presented in this paper, and it provides an explanation why accurate GW quasiparticle spectra are easy to obtain for some compounds and very difficult for others.
In this lecture we present many-body perturbation theory as a method to determine quasiparticle excitations in solids, especially electronic band structures, accurately from first principles.The main ingredient is the electronic self-energy that, in principle, contains all many-body exchange and correlation effects beyond the Hartree potential.As its exact mathematical expression is unknown, approximations must be used in practical calculations.The approximation is obtained using a systematic algebraic approach on the basis of Green function techniques.It constitutes an expansion of the self-energy up to linear order in the screened Coulomb potential, which describes the interaction between the quasiparticles and includes dynamic screening through the creation of exchange-correlation holes around the bare particles.The implementation of the approximation relies on a perturbative treatment starting from density functional theory.Besides a detailed mathematical discussion we focus on the underlying physical concepts and show some illustrative applications.
The correlation energy of an inhomogeneous electron gas is ruled by the density correlation function LLL. It obeys a Bethe-Salpeter equation (BSE), whose kernel is the functional derivative of the particle self-energy. It is solved on equal footing for two-particle and single-particle Green's functions within Hedin's GW approximation. According to the Hartree and exchange-correlation contributions to the kernel, the correlation expression decays into random-phase-approximation (RPA) -like and screened-exchange-like parts. The actual treatment of the screened potential WWW in the BSE for LLL leads to different second-order screened exchange corrections (SOSEX). Including correlation in the kernel even yields approximations beyond SOSEX. The progress beyond RPA is illustrated for spinless and spin-polarized homogeneous electron gases, indicating that the screening treatment is substantial.
With a transcorrelated Hamiltonian, we perform a many body perturbation calculation on the uniform electron gas in the high density regime. By using a correlation factor optimized for a single determinant Jastrow ansatz, the second order correlation energy is calculated as 1−ln⁡2π2ln(rs)−0.05075. This already reproduces the exact logarithmic term of the random phase approximation (RPA) result, while the constant term is roughly 7% larger than the RPA one. The close agreement with the RPA method demonstrates that the transcorrelated method offers a viable and potentially efficient method for treating metallic systems.
Density functional theory (DFT) performs reasonably well for the determination of structural properties of many materials, but fails to predict electronic band gap values accurately. \nSuch a failure of DFT is not unexpected, since there exists no formal justification for interpreting the DFT eigenvalues as addition or removal energies of the many-body system (quasiparticle energies).\nAn alternative approach to the study of exchange-correlation effects in many-particle systems is provided by many-body perturbation theory (MBPT), which defines a rigorous approach to description of excited-state properties based on the Green's function formalism and the concept of quasiparticle electrons. Within MBPT, one can calculate the quasiparticle (QP) energies and the QP amplitudes by approximately solving, within the so-called GW approximation, the set of coupled integro-differential equations proposed by Hedin in 1965.\nThe GW method generally yields significantly better values for QP energies with respect to DFT as it accounts for the dynamic many-body effects in the electron-electron interaction going beyond the mean-field approximation of independent Kohn-Sham particles. The drawback of such an involved approach, however, is its large computational cost, which is mainly due to the evaluation of dielectric matrices, their inversion, and the solution of a non-Hermitian and nonlinear eigenvalue problem. \nThe first part of this PhD research is devoted to the implementation of an efficient and scalable coarse-grained parallel algorithms and the development of state-of-the-art methods to solve the GW equations.\nAll these techniques are then applied to the study of the quasiparticle band structures of SiO2 in the alpha quartz crystalline structure. The effects of the different approximations involved in the theory, the influence of self-consistency, and a systematic analysis of the reliability of the different plasmon-pole models are presented and discussed. \nThe second part of this PhD work is dedicated to the implementation of a methodological approach to the solution of the GW equations based on the so-called projector augmented wave method (PAW) proposed by Blochl in 1994. This new approach permits to remove many of the limitations intrinsic to the use of the pseudopotential technique widely used for ab-initio calculations, allowing one to reach GW results closer to an all-electron implementation while still maintaining computational efficiency, and flexibility. This allows us to achieve a coherent implementation of many different algorithms which enables a detailed comparison of GW calculations not yet performed until now. Results for the quasiparticle band structure and optical spectra of prototype s- and p-compounds are discussed and compared to recent studies reported in the scientific literature (when available).
In order to properly understand the utility of many-body perturbation theory as applied to finite systems, we use the configuration interaction approach to compute the electronic self-energy. The validity of the commonly used GWGWGW approximation from many-body perturbation theory is tested by comparing the self-energy explicitly computed according to the diagrammatic expansions to that from the inversion of the Dyson equation. It is constructed as a functional of the interacting Green function GGG and screened Coulomb interaction WWW. WWW is explicitly computed through the diagonalization of the many-body Hamiltonian and, thus, takes polarization effects into account beyond the random phase approximation. The Na9+\mathrm{Na}_{9}{}^{+}Na9​+ cluster and the dissociation of the C2{\mathrm{C}}_{2}C2​ molecule are discussed as examples. We find that the GWGWGW approximation yields accurate results for weakly correlated systems with a predominantly single-determinant ground state such as Na9+\mathrm{Na}_{9}{}^{+}Na9​+ clusters. The C2{\mathrm{C}}_{2}C2​ molecule is a pathological system with multideterminantal ground state where the GWGWGW approach ceases to be valid.
The effects of short-range correlations on the properties of the electron gas are studied using an approach phrased in the polarization potential language of Pines and Nozi\`eres. It was argued by Lowy and Brown that two-body terms make up the dominant contribution to the effective interaction for large momentum transfers. In the region of small-momentum transfers it is known that the random-phase approximation provides an adequate picture of the electron gas. With these two regions of large- and small-momentum transfers in mind, a polarization propagator is constructed to interpolate between them. A number of sum rules, among them the fff-sum rule, involving the imaginary part of the polarization propagator, were found to be satisfied. The failure of our theory to satisfy the compressibility sum rule is discussed in detail. We present here the results for the local ttt-matrix, the dynamic and static form factors, and the pair correlation function. That our pair-correlation function remains positive and is essentially the same as in the paper of Lowy and Brown indicates that the physics and not the formalism is responsible for the results we obtain.