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:2103.04419 · Physical review. B./Physical review. B

Random phase approximation for gapped systems: Role of vertex corrections and applicability of the constrained random phase approximation

Erik G. C. P. van Loon, Malte Rösner, Mikhail I. Katsnelson, Tim Oliver Wehling

Published Jul 22, 2021• 37 citations
arXivHTMLPDF

Abstract

The many-body theory of interacting electrons poses an intrinsically difficult problem that requires simplifying assumptions. For the determination of electronic screening properties of the Coulomb interaction, the random phase approximation (RPA) provides such a simplification. Here we explicitly show that this approximation is justified for band structures with sizable band gaps. This is when the electronic states responsible for the screening are energetically far away from the Fermi level, which is equivalent to a short electronic propagation length of these states. The RPA contains exactly those diagrams in which the classical Coulomb interaction covers all distances, whereas neglected vertex corrections involve quantum tunneling through the barrier formed by the band gap. Our analysis of electron-electron interactions provides a real-space analogy to Migdal's theorem on the smallness of vertex corrections in electron-phonon problems. An important application is the increasing use of constrained RPA calculations of effective interactions. We find that their usage of Kohn-Sham energies accounts for the leading local (excitonic) vertex correction in insulators.

DOI
10.1103/physrevb.104.045134

Related papers

Band-gap energy in the random-phase approximation to density-functional theory

Yann‐Michel Niquet, Xavier Gonze

We calculate the interacting bandgap energy of a solid within the random-phase approximation (RPA) to density functional theory (DFT). The interacting bandgap energy is defined as ${E}_{g}={E}^{\mathrm{RPA}}(N+1)+{E}^{\mathrm{RPA}}(N\ensuremath{-}1)\ensuremath{-}2{E}^{\mathrm{RPA}}(N)$, where ${E}^{\mathrm{RPA}}(N)$ is the total DFT-RPA energy of the $N$-electron system. We compare the interacting bandgap energy with the Kohn-Sham bandgap energy, which is the difference between the conduction and valence band edges in the Kohn-Sham band structure. We show that they differ by an unrenormalized ``${G}_{0}{W}_{0}$'' self-energy correction (i.e., a $GW$ self-energy correction computed using Kohn-Sham orbitals and energies as input). This provides a well-defined and meaningful interpretation to ${G}_{0}{W}_{0}$ quasiparticle bandgap calculations, but questions the physics behind the renormalization factors in the expression of the bandgap energy. We also separate the kinetic from the Coulomb contributions to the DFT-RPA bandgap energy, and discuss the related problem of the derivative discontinuity in the DFT-RPA functional. Last we discuss the applicability of our results to other functionals based on many-body perturbation theory.

Physical Review B•Dec 16, 2004•55 citations

Effective interactions in multiband systems from constrained summations

Carsten Honerkamp A

In correlated electron materials, the application of many-body techniques for the study of interaction effects or unconventional superconductivity often requires the formulation of an effective low-energy model that contains only the relevant bands near the Fermi level. However, the bands away from the Fermi level are known to renormalize the low-energy interactions substantially. Here we compare different schemes to derive low-energy effective theories for interacting electrons in solids. The frequently used constrained random phase approximation (cRPA) is identified as a particular resummation of higher-order interaction terms that includes important virtual corrections. We then propose an adapted functional renormalization group scheme that includes the cPRA, but also allows one to go beyond the cRPA approximation. We study a simple two-band model in order to explore the differences between the approximations.

arXiv:1112.5143•Physical Review B•May 15, 2012•25 citations

Beyond the random phase approximation with a local exchange vertex

Maria Hellgren, Nicola Colonna, Stefano de Gironcoli

With the aim of constructing an electronic structure approach that systematically goes beyond the GW and random phase approximation (RPA) we introduce a vertex correction based on the exact-exchange (EXX) potential of time-dependent density functional theory. The EXX vertex function is constrained to be local but is expected to capture similar physics as the Hartree-Fock vertex. With the EXX vertex, we then unify different beyond-RPA approaches such as the various resummations of RPA with exchange and the second-order screened exchange approximation. The theoretical analysis is supported by numerical studies on the hydrogen dimer and the electron gas, and we discuss the role of including the vertex correction in both the screened interaction and the self-energy. Finally, we give details on our implementation within the plane-wave pseudo potential framework and demonstrate the excellent performance of the different RPA with exchange methods in describing the energetics of hydrogen and van der Waals bonds.

arXiv:1804.10566•Physical review. B./Physical review. B•Jul 12, 2018•42 citations

Convergence behavior of the random phase approximation renormalized correlation energy

Jefferson E. Bates, J. Sensenig, Adrienn Ruzsinszky

Based on the random phase approximation (RPA), RPA renormalization [J. E. Bates and F. Furche, J. Chem. Phys. 139, 171103 (2013)] is a robust many-body perturbation theory that works for molecules and materials because it does not diverge as the Kohn-Sham gap approaches zero. Additionally, RPA renormalization enables the simultaneous calculation of RPA and beyond-RPA correlation energies since the total correlation energy is the sum of a series of independent contributions. The first-order approximation (RPAr1) yields the dominant beyond-RPA contribution to the correlation energy for a given exchange-correlation kernel, but systematically underestimates the total beyond-RPA correction. For both the homogeneous electron gas model and real systems, we demonstrate numerically that RPA renormalization beyond first order converges monotonically to the infinite-order beyond-RPA correlation energy for several model exchange-correlation kernels and that the rate of convergence is principally determined by the choice of the kernel and spin polarization of the ground state. The monotonic convergence is rationalized from an analysis of the RPA renormalized correlation energy corrections, assuming the exchange-correlation kernel and response functions satisfy some reasonable conditions. For spin-unpolarized atoms, molecules, and bulk solids, we find that RPA renormalization is typically converged to 1 meV error or less by fourth order regardless of the band gap or dimensionality. Most spin-polarized systems converge at a slightly slower rate, with errors on the order of 10 meV at fourth order and typically requiring up to sixth order to reach 1 meV error or less. Slowest to converge, however, open-shell atoms present the most challenging case and require many higher orders to converge.

Physical review. B./Physical review. B•May 26, 2017•17 citations

Limitations of constrained random phase approximation downfolding

Carsten Honerkamp A, Hiroshi Shinaoka, Fakher F. Assaad +1 more

We check the accuracy of the constrained random phase approximation (cRPA) downfolding scheme by considering one-dimensional two- and three-orbital Hubbard models with a target band at the Fermi level and one or two screening bands away from the Fermi level. Using numerically exact quantum Monte Carlo simulations of the full and downfolded model, we demonstrate that depending on filling, the effective interaction in the low-energy theory is either barely screened or antiscreened, in contrast to the cRPA prediction. This observation is explained by a functional renormalization group analysis, which shows that the cRPA contribution to the screening is to a large extent canceled by other diagrams in the direct particle-hole channel. We comment on the implications of this finding for the ab initio estimation of interaction parameters in low-energy descriptions of solids.

arXiv:1809.03742•Physical review. B./Physical review. B•Dec 26, 2018•57 citations

Efficient vertex parametrization for the constrained functional renormalization group for effective low-energy interactions in multiband systems

Carsten Honerkamp A

We describe an efficient approximation for the electron-electron interaction in the determination of the low-energy effective interaction in multiband lattice systems. By using ideas for channel decomposition, form-factor expansion, and the truncated-unity technique we describe the interaction as arising from the nonlocal and orbital-dependent coupling of particle-hole and particle-particle bilinears formed by fields residing in the same one or two orbitals. This allows us to employ the constrained functional renormalization group (cfRG) with a suitable momentum and frequency discretization. The approach gives insights into the nonlocal screening of spin and charge interactions when bands away from the Fermi level are integrated out. Specifically, we compute the effective low-energy interactions in the low-energy target band of a three-band model with on-site and nonlocal bare interactions. We show that the cfRG adds important features to the effective target-band interaction that cannot be found using the constrained random phase approximation (cRPA).

arXiv:1805.01669•Physical review. B./Physical review. B•Oct 17, 2018•25 citations
The many-body theory of interacting electrons poses an intrinsically difficult problem that requires simplifying assumptions. For the determination of electronic screening properties of the Coulomb interaction, the random phase approximation (RPA) provides such a simplification. Here we explicitly show that this approximation is justified for band structures with sizable band gaps. This is when the electronic states responsible for the screening are energetically far away from the Fermi level, which is equivalent to a short electronic propagation length of these states. The RPA contains exactly those diagrams in which the classical Coulomb interaction covers all distances, whereas neglected vertex corrections involve quantum tunneling through the barrier formed by the band gap. Our analysis of electron-electron interactions provides a real-space analogy to Migdal's theorem on the smallness of vertex corrections in electron-phonon problems. An important application is the increasing use of constrained RPA calculations of effective interactions. We find that their usage of Kohn-Sham energies accounts for the leading local (excitonic) vertex correction in insulators.
We calculate the interacting bandgap energy of a solid within the random-phase approximation (RPA) to density functional theory (DFT). The interacting bandgap energy is defined as ${E}_{g}={E}^{\mathrm{RPA}}(N+1)+{E}^{\mathrm{RPA}}(N\ensuremath{-}1)\ensuremath{-}2{E}^{\mathrm{RPA}}(N)$, where ERPA(N){E}^{\mathrm{RPA}}(N)ERPA(N) is the total DFT-RPA energy of the NNN-electron system. We compare the interacting bandgap energy with the Kohn-Sham bandgap energy, which is the difference between the conduction and valence band edges in the Kohn-Sham band structure. We show that they differ by an unrenormalized ``G0W0{G}_{0}{W}_{0}G0​W0​'' self-energy correction (i.e., a GWGWGW self-energy correction computed using Kohn-Sham orbitals and energies as input). This provides a well-defined and meaningful interpretation to G0W0{G}_{0}{W}_{0}G0​W0​ quasiparticle bandgap calculations, but questions the physics behind the renormalization factors in the expression of the bandgap energy. We also separate the kinetic from the Coulomb contributions to the DFT-RPA bandgap energy, and discuss the related problem of the derivative discontinuity in the DFT-RPA functional. Last we discuss the applicability of our results to other functionals based on many-body perturbation theory.
In correlated electron materials, the application of many-body techniques for the study of interaction effects or unconventional superconductivity often requires the formulation of an effective low-energy model that contains only the relevant bands near the Fermi level. However, the bands away from the Fermi level are known to renormalize the low-energy interactions substantially. Here we compare different schemes to derive low-energy effective theories for interacting electrons in solids. The frequently used constrained random phase approximation (cRPA) is identified as a particular resummation of higher-order interaction terms that includes important virtual corrections. We then propose an adapted functional renormalization group scheme that includes the cPRA, but also allows one to go beyond the cRPA approximation. We study a simple two-band model in order to explore the differences between the approximations.
With the aim of constructing an electronic structure approach that systematically goes beyond the GW and random phase approximation (RPA) we introduce a vertex correction based on the exact-exchange (EXX) potential of time-dependent density functional theory. The EXX vertex function is constrained to be local but is expected to capture similar physics as the Hartree-Fock vertex. With the EXX vertex, we then unify different beyond-RPA approaches such as the various resummations of RPA with exchange and the second-order screened exchange approximation. The theoretical analysis is supported by numerical studies on the hydrogen dimer and the electron gas, and we discuss the role of including the vertex correction in both the screened interaction and the self-energy. Finally, we give details on our implementation within the plane-wave pseudo potential framework and demonstrate the excellent performance of the different RPA with exchange methods in describing the energetics of hydrogen and van der Waals bonds.
Based on the random phase approximation (RPA), RPA renormalization [J. E. Bates and F. Furche, J. Chem. Phys. 139, 171103 (2013)] is a robust many-body perturbation theory that works for molecules and materials because it does not diverge as the Kohn-Sham gap approaches zero. Additionally, RPA renormalization enables the simultaneous calculation of RPA and beyond-RPA correlation energies since the total correlation energy is the sum of a series of independent contributions. The first-order approximation (RPAr1) yields the dominant beyond-RPA contribution to the correlation energy for a given exchange-correlation kernel, but systematically underestimates the total beyond-RPA correction. For both the homogeneous electron gas model and real systems, we demonstrate numerically that RPA renormalization beyond first order converges monotonically to the infinite-order beyond-RPA correlation energy for several model exchange-correlation kernels and that the rate of convergence is principally determined by the choice of the kernel and spin polarization of the ground state. The monotonic convergence is rationalized from an analysis of the RPA renormalized correlation energy corrections, assuming the exchange-correlation kernel and response functions satisfy some reasonable conditions. For spin-unpolarized atoms, molecules, and bulk solids, we find that RPA renormalization is typically converged to 1 meV error or less by fourth order regardless of the band gap or dimensionality. Most spin-polarized systems converge at a slightly slower rate, with errors on the order of 10 meV at fourth order and typically requiring up to sixth order to reach 1 meV error or less. Slowest to converge, however, open-shell atoms present the most challenging case and require many higher orders to converge.
We check the accuracy of the constrained random phase approximation (cRPA) downfolding scheme by considering one-dimensional two- and three-orbital Hubbard models with a target band at the Fermi level and one or two screening bands away from the Fermi level. Using numerically exact quantum Monte Carlo simulations of the full and downfolded model, we demonstrate that depending on filling, the effective interaction in the low-energy theory is either barely screened or antiscreened, in contrast to the cRPA prediction. This observation is explained by a functional renormalization group analysis, which shows that the cRPA contribution to the screening is to a large extent canceled by other diagrams in the direct particle-hole channel. We comment on the implications of this finding for the ab initio estimation of interaction parameters in low-energy descriptions of solids.
We describe an efficient approximation for the electron-electron interaction in the determination of the low-energy effective interaction in multiband lattice systems. By using ideas for channel decomposition, form-factor expansion, and the truncated-unity technique we describe the interaction as arising from the nonlocal and orbital-dependent coupling of particle-hole and particle-particle bilinears formed by fields residing in the same one or two orbitals. This allows us to employ the constrained functional renormalization group (cfRG) with a suitable momentum and frequency discretization. The approach gives insights into the nonlocal screening of spin and charge interactions when bands away from the Fermi level are integrated out. Specifically, we compute the effective low-energy interactions in the low-energy target band of a three-band model with on-site and nonlocal bare interactions. We show that the cfRG adds important features to the effective target-band interaction that cannot be found using the constrained random phase approximation (cRPA).