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

Communication: Random phase approximation renormalized many-body perturbation theory

Jefferson E. Bates, Filipp Furche

Published Nov 6, 2013• 115 citations
Publisher

Abstract

We derive a renormalized many-body perturbation theory (MBPT) starting from the random phase approximation (RPA). This RPA-renormalized perturbation theory extends the scope of single-reference MBPT methods to small-gap systems without significantly increasing the computational cost. The leading correction to RPA, termed the approximate exchange kernel (AXK), substantially improves upon RPA atomization energies and ionization potentials without affecting other properties such as barrier heights where RPA is already accurate. Thus, AXK is more balanced than second-order screened exchange [A. Grüneis et al., J. Chem. Phys. 131, 154115 (2009)], which tends to overcorrect RPA for systems with stronger static correlation. Similarly, AXK avoids the divergence of second-order Møller-Plesset (MP2) theory for small gap systems and delivers a much more consistent performance than MP2 across the periodic table at comparable cost. RPA+AXK thus is an accurate, non-empirical, and robust tool to assess and improve semi-local density functional theory for a wide range of systems previously inaccessible to first-principles electronic structure calculations.

DOI
10.1063/1.4827254

Related papers

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

Adiabatic Connection without Coupling Constant Integration

Jefferson E. Bates, Niladri Sengupta, J. Sensenig +1 more

Using a second-order approximation to Random Phase Approximation renormalized (RPAr) many-body perturbation theory for the interacting density-density response function, we have developed a so-called higher-order terms (HOT) approximation for the correlation energy. In combination with the first-order RPAr correction, our new method faithfully captures the infinite-order correlation for a given exchange-correlation kernel, yielding errors of the total correlation energy on the order of 1% or less for most systems. For exchange-like kernels, our new method has the further benefit that the coupling-strength integration can be completely eliminated resulting in a modest reduction in computational cost compared to the traditional approach. When the correlation energy is accurately reproduced by the HOT approximation, structural properties and energy differences are also accurately reproduced, as we demonstrate for several periodic solids and some molecular systems. Energy differences involving fragmentation are challenging for the HOT method, however, due to errors that may not cancel between a composite system and its constituent pieces.

Journal of Chemical Theory and Computation•May 7, 2018•14 citations

Random phase approximation-based local natural orbital coupled cluster theory

Ruiheng Song, Xiliang Gong, Aamy Bakry +1 more

Practical applications of fragment embedding and closely related local correlation methods critically depend on a judicious choice of a low-level theory to define the local embedding subspace and to capture long-range electrostatic and correlation effects outside the embedding region. Second-order Møller-Plesset perturbation theory (MP2) is by far the most widely used correlated low-level theory; however, its applicability becomes questionable in systems where MP2 is known to fail either quantitatively or qualitatively. In this work, we present the random phase approximation (RPA) as a promising alternative low-level theory to MP2 within the local natural orbital-based coupled-cluster (LNO-CC) framework. We demonstrate that RPA-based LNO-CC closely matches the performance of its MP2-based counterpart for systems with sizable energy gaps, while delivering significantly faster convergence toward the canonical coupled-cluster limit for metallic systems, particularly as the thermodynamic limit is approached. These results highlight the critical role of the low-level theory in fragment embedding and local correlation methods and identify RPA as a compelling alternative to the commonly used MP2.

arXiv:2601.00131•Dec 31, 2025•0 citations

Generalized Many-Body Perturbation Theory for the Electron Correlation Energy: Multireference Random Phase Approximation via Diagrammatic Resummation

Yuqi Wang, Wei‐Hai Fang, Zhendong Li

Many-body perturbation theory (MBPT) based on Green’s functions and Feynman diagrams provides a fundamental theoretical framework for various ab initio computational approaches in molecular and materials science, including the random phase approximation (RPA) and GW approximation. Unfortunately, this perturbation expansion often fails in systems with strong multireference characters. Extending diagrammatic MBPT to the multireference case is highly nontrivial and remains largely unexplored, primarily due to the breakdown of Wick’s theorem. In this work, we develop a diagrammatic multireference generalization of MBPT for computing correlation energies of strongly correlated systems, by using the cumulant expansion of many-body Green’s functions in place of Wick’s theorem. This theoretical framework bridges the gap between MBPT in condensed matter physics and multireference perturbation theories (MRPT) in quantum chemistry, which had been almost exclusively formulated within time-independent wave function frameworks prior to this work. Our formulation enables the explicit incorporation of strong correlation effects from the outset as in MRPT, while treating residual weak interactions through a generalized diagrammatic perturbation expansion as in MBPT. As a concrete demonstration, we formulate a multireference (MR) extension of the standard single-reference (SR) RPA by systematically resumming generalized ring diagrams, which naturally leads to a unified set of equations applicable to both SR and MR cases. Benchmark calculations on prototypical molecular systems reveal that MR-RPA successfully resolves the well-known failure of SR-RPA in strongly correlated systems. This theoretical advancement paves the way for advancing ab initio computational methods through diagrammatic resummation techniques in future.

arXiv:2410.22990•The Journal of Physical Chemistry Letters•Mar 18, 2025•9 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

A unified diagrammatic formulation of single-reference and multi-reference random phase approximations: The particle–hole and particle–particle channels

Yuqi Wang, Wei‐Hai Fang, Zhendong Li

A diagrammatic multi-reference generalization of many-body perturbation theory was recently introduced [Wang et al., J. Phys. Chem. Lett. 16, 3047 (2025)]. This framework allows us to extend single-reference (SR) Green's function methods defined at the diagrammatic level naturally to the multi-reference case, as previously exemplified by the formulation of multi-reference direct random phase approximation (MR-dRPA) and the multi-reference second-order screened exchange approximation. In this work, we further elaborate on this framework and use it to develop MR generalizations of two other RPA variants, namely, particle-hole (ph) RPA with exchange (MR-RPAx) and particle-particle RPA (MR-ppRPA). We define these two MR generalizations by infinite-order resummations of the generalized "ring" and "ladder" diagrams with antisymmetrized interaction vertices, respectively, which incorporate the contributions from the active-space connected two-body Green's functions. As for MR-dRPA, we derive unified sets of equations that hold at both SR and MR levels for RPAx and ppRPA, respectively. We perform numerical studies of prototypical systems using the three MR-RPA methods and carry out a perturbative analysis to gain a deeper understanding of their behaviors. We find that error cancellation between the second and third orders is a key factor for both SR-RPA and MR-RPA. In addition, we observe that MR-phRPA (MR-dRPA and MR-RPAx) and MR-ppRPA tend to overestimate and underestimate correlation energies, respectively, suggesting that better accuracy can be achieved by further combining these two channels in the future.

The Journal of Chemical Physics•Nov 3, 2025•3 citations
We derive a renormalized many-body perturbation theory (MBPT) starting from the random phase approximation (RPA). This RPA-renormalized perturbation theory extends the scope of single-reference MBPT methods to small-gap systems without significantly increasing the computational cost. The leading correction to RPA, termed the approximate exchange kernel (AXK), substantially improves upon RPA atomization energies and ionization potentials without affecting other properties such as barrier heights where RPA is already accurate. Thus, AXK is more balanced than second-order screened exchange [A. Grüneis et al., J. Chem. Phys. 131, 154115 (2009)], which tends to overcorrect RPA for systems with stronger static correlation. Similarly, AXK avoids the divergence of second-order Møller-Plesset (MP2) theory for small gap systems and delivers a much more consistent performance than MP2 across the periodic table at comparable cost. RPA+AXK thus is an accurate, non-empirical, and robust tool to assess and improve semi-local density functional theory for a wide range of systems previously inaccessible to first-principles electronic structure calculations.
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.
Using a second-order approximation to Random Phase Approximation renormalized (RPAr) many-body perturbation theory for the interacting density-density response function, we have developed a so-called higher-order terms (HOT) approximation for the correlation energy. In combination with the first-order RPAr correction, our new method faithfully captures the infinite-order correlation for a given exchange-correlation kernel, yielding errors of the total correlation energy on the order of 1% or less for most systems. For exchange-like kernels, our new method has the further benefit that the coupling-strength integration can be completely eliminated resulting in a modest reduction in computational cost compared to the traditional approach. When the correlation energy is accurately reproduced by the HOT approximation, structural properties and energy differences are also accurately reproduced, as we demonstrate for several periodic solids and some molecular systems. Energy differences involving fragmentation are challenging for the HOT method, however, due to errors that may not cancel between a composite system and its constituent pieces.
Practical applications of fragment embedding and closely related local correlation methods critically depend on a judicious choice of a low-level theory to define the local embedding subspace and to capture long-range electrostatic and correlation effects outside the embedding region. Second-order Møller-Plesset perturbation theory (MP2) is by far the most widely used correlated low-level theory; however, its applicability becomes questionable in systems where MP2 is known to fail either quantitatively or qualitatively. In this work, we present the random phase approximation (RPA) as a promising alternative low-level theory to MP2 within the local natural orbital-based coupled-cluster (LNO-CC) framework. We demonstrate that RPA-based LNO-CC closely matches the performance of its MP2-based counterpart for systems with sizable energy gaps, while delivering significantly faster convergence toward the canonical coupled-cluster limit for metallic systems, particularly as the thermodynamic limit is approached. These results highlight the critical role of the low-level theory in fragment embedding and local correlation methods and identify RPA as a compelling alternative to the commonly used MP2.
Many-body perturbation theory (MBPT) based on Green’s functions and Feynman diagrams provides a fundamental theoretical framework for various ab initio computational approaches in molecular and materials science, including the random phase approximation (RPA) and GW approximation. Unfortunately, this perturbation expansion often fails in systems with strong multireference characters. Extending diagrammatic MBPT to the multireference case is highly nontrivial and remains largely unexplored, primarily due to the breakdown of Wick’s theorem. In this work, we develop a diagrammatic multireference generalization of MBPT for computing correlation energies of strongly correlated systems, by using the cumulant expansion of many-body Green’s functions in place of Wick’s theorem. This theoretical framework bridges the gap between MBPT in condensed matter physics and multireference perturbation theories (MRPT) in quantum chemistry, which had been almost exclusively formulated within time-independent wave function frameworks prior to this work. Our formulation enables the explicit incorporation of strong correlation effects from the outset as in MRPT, while treating residual weak interactions through a generalized diagrammatic perturbation expansion as in MBPT. As a concrete demonstration, we formulate a multireference (MR) extension of the standard single-reference (SR) RPA by systematically resumming generalized ring diagrams, which naturally leads to a unified set of equations applicable to both SR and MR cases. Benchmark calculations on prototypical molecular systems reveal that MR-RPA successfully resolves the well-known failure of SR-RPA in strongly correlated systems. This theoretical advancement paves the way for advancing ab initio computational methods through diagrammatic resummation techniques in future.
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.
A diagrammatic multi-reference generalization of many-body perturbation theory was recently introduced [Wang et al., J. Phys. Chem. Lett. 16, 3047 (2025)]. This framework allows us to extend single-reference (SR) Green's function methods defined at the diagrammatic level naturally to the multi-reference case, as previously exemplified by the formulation of multi-reference direct random phase approximation (MR-dRPA) and the multi-reference second-order screened exchange approximation. In this work, we further elaborate on this framework and use it to develop MR generalizations of two other RPA variants, namely, particle-hole (ph) RPA with exchange (MR-RPAx) and particle-particle RPA (MR-ppRPA). We define these two MR generalizations by infinite-order resummations of the generalized "ring" and "ladder" diagrams with antisymmetrized interaction vertices, respectively, which incorporate the contributions from the active-space connected two-body Green's functions. As for MR-dRPA, we derive unified sets of equations that hold at both SR and MR levels for RPAx and ppRPA, respectively. We perform numerical studies of prototypical systems using the three MR-RPA methods and carry out a perturbative analysis to gain a deeper understanding of their behaviors. We find that error cancellation between the second and third orders is a key factor for both SR-RPA and MR-RPA. In addition, we observe that MR-phRPA (MR-dRPA and MR-RPAx) and MR-ppRPA tend to overestimate and underestimate correlation energies, respectively, suggesting that better accuracy can be achieved by further combining these two channels in the future.