The Journal of Chemical Physics

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

Published Nov 3, 2025• 3 citations

Abstract

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.

DOI: 10.1063/5.0293280

Related papers

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.22990The Journal of Physical Chemistry LettersMar 18, 20259 citations

A Diagrammatic Route to Multi-Reference GW for Strongly Correlated Molecules

Yuqi Wang, Wei‐Hai Fang, Zhendong Li

The GW approximation is a cornerstone of many-body perturbation theory for computing single-particle excitations, yet it can fail qualitatively in strongly correlated systems where the single-reference picture becomes inadequate. To overcome this limitation, we develop a diagrammatic route to multi-reference GW (MR-GW) from an interacting active-space reference, treating strong correlation nonperturbatively and the residual interaction through dynamical screening at the multi-reference random phase approximation level. Because the conventional diagrammatic self-energy expansion no longer applies, we define a practical one-shot MR-GW approximation within a generalized diagrammatic framework, recovering standard GW as a special case. Applications to prototypical correlated atoms and molecules show that the interacting active-space reference captures multiconfigurational satellite structures poorly described by single-reference GW, while dynamical screening generally improves the satellite energies as well as the leading ionization and electron-attachment energies. This work opens a diagrammatic route to Green's-function methods that combine nonperturbative strong correlation with dynamical screening.

arXiv:2604.16013Apr 17, 20260 citations

Communication: Random phase approximation renormalized many-body perturbation theory

Jefferson E. Bates, Filipp Furche

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.

The Journal of Chemical PhysicsNov 6, 2013115 citations

Duality of Particle–Hole and Particle–Particle Theories for Strongly Correlated Electronic Systems

Aleksandra Tucholska, Yang Guo, Katarzyna Pernal

We propose a novel approach to electron correlation for multireference systems. It is based on particle-hole (ph) and particle-particle (pp) theories in the second-order, developed in the random phase approximation (RPA) framework for multireference wave functions. We show a formal correspondence (duality), between contributions to the correlation energy in the ph and pp pictures. It allows us to describe correlation energy by rigorously combining pp and ph terms, avoiding correlation double counting. The multireference ph, pp, and the combined correlation methods are applied to ground and excited states of systems in the intermediate and strong correlation regimes and compared with the multireference second-order perturbation method (MRPT2). It is shown that the pp approximation fails to describe dissociation of multiple bonds. The ph-pp combined method is overall superior to both ph and pp alone. It parallels good accuracy of the second-order perturbation theory for ground states and singlet excitation energies. For the singlet-triplet gaps of biradicals its accuracy is significantly better. This is impressive, taking into account that it relies only on one- and two-body density matrices, while MRPT2 methods typically require density matrices up to the four-body.

arXiv:2409.11284The Journal of Physical Chemistry LettersNov 25, 20246 citations

Particle-particle and quasiparticle random phase approximations: Connections to coupled cluster theory

Gustavo E. Scuseria, Thomas M. Henderson, Ireneusz W. Bulik

We establish a formal connection between the particle-particle (pp) random phase approximation (RPA) and the ladder channel of the coupled cluster doubles (CCD) equations. The relationship between RPA and CCD is best understood within a Bogoliubov quasiparticle (qp) RPA formalism. This work is a follow-up to our previous formal proof on the connection between particle-hole (ph) RPA and ring-CCD. Whereas RPA is a quasibosonic approximation, CC theory is a "correct bosonization" in the sense that the wavefunction and Hilbert space are exactly fermionic, yet the amplitude equations can be interpreted as adding different quasibosonic RPA channels together. Coupled cluster theory achieves this goal by interacting the ph (ring) and pp (ladder) diagrams via a third channel that we here call "crossed-ring" whose presence allows for full fermionic antisymmetry. Additionally, coupled cluster incorporates what we call "mosaic" terms which can be absorbed into defining a new effective one-body Hamiltonian. The inclusion of these mosaic terms seems to be quite important. The pp-RPA and qp-RPA equations are textbook material in nuclear structure physics but are largely unknown in quantum chemistry, where particle number fluctuations and Bogoliubov determinants are rarely used. We believe that the ideas and connections discussed in this paper may help design improved ways of incorporating RPA correlation into density functionals based on a CC perspective.

arXiv:1306.6360The Journal of Chemical PhysicsSep 12, 2013110 citations

RPA as a Hessian Closure: Effective Functionals and Source-Variable Duality Across DFT, LR-TDDFT, 1RDMFT, and MBPT

Nan Sheng

We present a variational formulation of the random phase approximation (RPA) that places density functional theory (DFT), linear-response time-dependent density functional theory (LR-TDDFT), one-body reduced density matrix functional theory (1RDMFT), and Green's function many-body perturbation theory (MBPT) into a common source-variable hierarchy. The central claim is that a broad class of RPA constructions can be organized, independently of any one problem-specific formula, diagrammatic resummation, or small-amplitude equation of motion, as closure approximations to the exact Hessian of an effective functional. In this language, exact linear response is governed by the Hessian of the corresponding effective functional, while RPA is obtained by retaining a reference contribution together with an explicit interaction kernel and discarding the irreducible remainder. The hierarchy has two independent enrichments of the density-level description. One may enlarge the static local density to a time-dependent density, giving the dynamical density channel of LR-TDDFT, or enlarge it to an equal-time bilocal one-body reduced density matrix, giving the static bilocal channel of 1RDMFT. The Green's function level combines both enrichments, since the one-particle Green's function is bilocal in both space and time. This picture clarifies the relation between DFT, LR-TDDFT, 1RDMFT, and MBPT through exact forward reductions and source restrictions, while emphasizing that the corresponding RPA closures need not commute under projection. The hierarchy also distinguishes the local branch-wise Legendre geometry common to all four levels from the stronger global convex duality that may emerge only in source sectors possessing additional positivity and global regularity.

arXiv:2606.09763Jun 8, 20260 citations