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:2409.11284 · The Journal of Physical Chemistry Letters

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

Aleksandra Tucholska, Yang Guo, Katarzyna Pernal

Published Nov 25, 2024• 6 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1021/acs.jpclett.4c02788

Related papers

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

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

Communication: Random-phase approximation excitation energies from approximate equation-of-motion coupled-cluster doubles

Timothy C. Berkelbach

The ground-state correlation energy calculated in the random-phase approximation (RPA) is known to be identical to that calculated using a subset of terms appearing in coupled-cluster theory with double excitations (CCD). In particular, for particle-hole (ph) RPA this equivalence requires keeping only those terms that generate time-independent ring diagrams, and for particle-particle (pp) RPA it requires keeping only those terms that generate ladder diagrams. Here I show that these identities extend to excitation energies, for which those calculated in each RPA are identical to those calculated using approximations to equation-of-motion coupled-cluster theory with double excitations (EOM-CCD). The equivalence requires three approximations to EOM-CCD: first, the ground-state CCD amplitudes are obtained from the ring-CCD or ladder-CCD equations (the same as for the correlation energy); second, the EOM eigenvalue problem is truncated to the minimal subspace, which is one particle + one hole for ph-RPA and two particles or two holes for pp-RPA; third, the similarity transformation of the Fock operator must be neglected, as it corresponds to a Brueckner-like dressing of the single-particle propagator, which is not present in the conventional RPA.

arXiv:1804.01037•The Journal of Chemical Physics•Jul 28, 2018•32 citations

Exchange-correlation energy from pairing matrix fluctuation and the particle-particle random-phase approximation

Helen van Aggelen, Yang Yang, Weitao Yang

We formulate an adiabatic connection for the exchange-correlation energy in terms of pairing matrix fluctuation. This connection opens new channels for density functional approximations based on pairing interactions. Even the simplest approximation to the pairing matrix fluctuation, the particle-particle random phase approximation (pp-RPA), has some highly desirable properties. It has no delocalization error with a nearly linear energy behavior for systems with fractional charges, describes van der Waals interactions similarly and thermodynamic properties significantly better than particle-hole RPA, and eliminates static correlation error for single-bond systems. Most significantly, the pp-RPA is the first known functional that has an explicit and closed-form dependence on the occupied and unoccupied orbitals and captures the energy derivative discontinuity in strongly correlated systems. These findings illustrate the potential of including pairing interactions within a density functional framework.

arXiv:1306.4957•Physical Review A•Sep 24, 2013•105 citations

Electron correlations: I. Ground-state results in the high-density regime

R. F. Bishop, K.H. Lührmann

The coupled-cluster or $\mathrm{exp}(S)$ formalism is applied to the problem of correlations in a many-electron system. In this first work in a series of papers we concentrate on exact results and restrict ourselves to various ring approximations applicable essentially to the high-density regime. It is shown that the well-known random-phase approximation (RPA) can be formulated as a nonlinear three-dimensional integral equation for the four-point function ${S}_{2}$ which provides a measure of the two-particle---hole pair component in the true ground-state wave function. Exact analytic solutions are presented for both ${S}_{2}$ and the three-point particle-hole vertex function in the RPA, and their properties are discussed. The Tamm-Dancoff approximation to the ring summation is also formulated, and the analagous exact solutions are presented for the first time. It is demonstrated that both solutions are unique, and a comparison of the two approximations is presented in terms of the correlation energy.

Physical review. B, Condensed matter•May 15, 1978•168 citations

Second-Order Multireference Algebraic Diagrammatic Construction Theory for Photoelectron Spectra of Strongly Correlated Systems

Koushik Chatterjee, Alexander Yu. Sokolov

We present a second-order formulation of multireference algebraic diagrammatic construction theory [ Sokolov, A. Yu. J. Chem. Phys. 2018, 149, 204113] for simulating photoelectron spectra of strongly correlated systems (MR-ADC(2)). The MR-ADC(2) method uses second-order multireference perturbation theory (MRPT2) to efficiently obtain ionization energies and intensities for many photoelectron transitions in a single computation. In contrast to conventional MRPT2 methods, MR-ADC(2) provides information about ionization of electrons in all orbitals (i.e., core and active) and allows computation of transition intensities in a straightforward and efficient way. Although equations of MR-ADC(2) depend on four-particle reduced density matrices, we demonstrate that computation of these large matrices can be completely avoided without introducing any approximations. The resulting MR-ADC(2) implementation has a lower computational scaling compared to conventional MRPT2 methods. We present results of MR-ADC(2) for photoelectron spectra of small molecules, a carbon dimer, and equally spaced hydrogen chains (H 10 and H 30 ) and outline directions for future developments.

arXiv:1905.13301•Journal of Chemical Theory and Computation•Sep 11, 2019•35 citations
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.
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.
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 ground-state correlation energy calculated in the random-phase approximation (RPA) is known to be identical to that calculated using a subset of terms appearing in coupled-cluster theory with double excitations (CCD). In particular, for particle-hole (ph) RPA this equivalence requires keeping only those terms that generate time-independent ring diagrams, and for particle-particle (pp) RPA it requires keeping only those terms that generate ladder diagrams. Here I show that these identities extend to excitation energies, for which those calculated in each RPA are identical to those calculated using approximations to equation-of-motion coupled-cluster theory with double excitations (EOM-CCD). The equivalence requires three approximations to EOM-CCD: first, the ground-state CCD amplitudes are obtained from the ring-CCD or ladder-CCD equations (the same as for the correlation energy); second, the EOM eigenvalue problem is truncated to the minimal subspace, which is one particle + one hole for ph-RPA and two particles or two holes for pp-RPA; third, the similarity transformation of the Fock operator must be neglected, as it corresponds to a Brueckner-like dressing of the single-particle propagator, which is not present in the conventional RPA.
We formulate an adiabatic connection for the exchange-correlation energy in terms of pairing matrix fluctuation. This connection opens new channels for density functional approximations based on pairing interactions. Even the simplest approximation to the pairing matrix fluctuation, the particle-particle random phase approximation (pp-RPA), has some highly desirable properties. It has no delocalization error with a nearly linear energy behavior for systems with fractional charges, describes van der Waals interactions similarly and thermodynamic properties significantly better than particle-hole RPA, and eliminates static correlation error for single-bond systems. Most significantly, the pp-RPA is the first known functional that has an explicit and closed-form dependence on the occupied and unoccupied orbitals and captures the energy derivative discontinuity in strongly correlated systems. These findings illustrate the potential of including pairing interactions within a density functional framework.
The coupled-cluster or exp(S)\mathrm{exp}(S)exp(S) formalism is applied to the problem of correlations in a many-electron system. In this first work in a series of papers we concentrate on exact results and restrict ourselves to various ring approximations applicable essentially to the high-density regime. It is shown that the well-known random-phase approximation (RPA) can be formulated as a nonlinear three-dimensional integral equation for the four-point function S2{S}_{2}S2​ which provides a measure of the two-particle---hole pair component in the true ground-state wave function. Exact analytic solutions are presented for both S2{S}_{2}S2​ and the three-point particle-hole vertex function in the RPA, and their properties are discussed. The Tamm-Dancoff approximation to the ring summation is also formulated, and the analagous exact solutions are presented for the first time. It is demonstrated that both solutions are unique, and a comparison of the two approximations is presented in terms of the correlation energy.
We present a second-order formulation of multireference algebraic diagrammatic construction theory [ Sokolov, A. Yu. J. Chem. Phys. 2018, 149, 204113] for simulating photoelectron spectra of strongly correlated systems (MR-ADC(2)). The MR-ADC(2) method uses second-order multireference perturbation theory (MRPT2) to efficiently obtain ionization energies and intensities for many photoelectron transitions in a single computation. In contrast to conventional MRPT2 methods, MR-ADC(2) provides information about ionization of electrons in all orbitals (i.e., core and active) and allows computation of transition intensities in a straightforward and efficient way. Although equations of MR-ADC(2) depend on four-particle reduced density matrices, we demonstrate that computation of these large matrices can be completely avoided without introducing any approximations. The resulting MR-ADC(2) implementation has a lower computational scaling compared to conventional MRPT2 methods. We present results of MR-ADC(2) for photoelectron spectra of small molecules, a carbon dimer, and equally spaced hydrogen chains (H 10 and H 30 ) and outline directions for future developments.