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.
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.
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 discuss diagrammatic modifications to the coupled cluster doubles (CCD) equations, wherein different groups of terms out of rings, ladders, crossed-rings, and mosaics can be removed to form approximations to the coupled cluster method, of interest due to their similarity with various types of random phase approximations. The finite uniform electron gas (UEG) is benchmarked for 14- and 54-electron systems at the complete basis set limit over a wide density range and performance of different flavours of CCD is determined. These results confirm that rings generally overcorrelate and ladders generally undercorrelate; mosaics-only CCD yields a result surprisingly close to CCD. We use a recently developed numerical analysis [J. J. Shepherd and A. Grüneis, Phys. Rev. Lett. 110, 226401 (2013)] to study the behaviours of these methods in the thermodynamic limit. We determine that the mosaics, on forming the Brueckner one-body Hamiltonian, open a gap in the effective one-particle eigenvalues at the Fermi energy. Numerical evidence is presented which shows that methods based on this renormalisation have convergent energies in the thermodynamic limit including mosaic-only CCD, which is just a renormalised MP2. All other methods including only a single channel, namely, ladder-only CCD, ring-only CCD, and crossed-ring-only CCD, appear to yield divergent energies; incorporation of mosaic terms prevents this from happening.
This paper is the first of two in which the coupled cluster method (CCM) or exp(S) formalism is applied to two-component Fermi systems, the aim being to describe real metals and superconductivity. In this paper the authors concentrate on exact results and restrict ourselves to a ring approximation, applicable essentially to a high-density regime. They show that in the ground-state formalism the random phase approximation (RPA) can be formulated as a system of coupled, bilinear integral equations satisfied by functions associated with the so-called four-point functions of the system which provide a measure of the two-particle two-hole components in the true ground-state wavefunction. These equations are analysed in the dimensionless parameter formed by the ratio of the species masses and exact analytic solutions obtained. For Coulombic potentials (V 11 V 22 =V 2 12 ) they show that the exact analytic solution is unique and obtain an expression for the correlation energy. For non-Coulombic potentials (V 11 V 22 not=V 2 12 ) they indicate how to obtain a possible analytic solution. An RPA-like treatment of the one- and two-body equations in the excited-state formalism is provided for completeness.
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 correlation energy in the direct random phase approximation (dRPA) can be written, among other possibilities, either in terms of the interaction strength averaged correlation density matrix, or in terms of the coupled cluster doubles amplitudes obtained in the direct ring approximation (drCCD). Although the corresponding dRPA correlation density matrix on the one hand, and the drCCD amplitude matrix on the other hand, differ significantly, they yield identical energies. Similarly, the analogous RPA and rCCD correlation energies calculated from antisymmetrized two-electron integrals are identical to each other despite very different underlying working equations. In the present communication, a direct correspondence between amplitudes and densities is established and investigated with perturbation theory arguments. Our analysis also sheds some light on the properties of recently proposed RPA/rCCD variants which use antisymmetrized integrals in part of the equations and nonantisymmetrized integrals in others.
Coupled cluster theory provides hierarchical many-particle models and is presently considered as the ultimate benchmark in quantum chemistry. Despite is practical significance, a rigorous mathematical analysis of its properties is still in its infancy. The present work focuses on nonlinear models within the random phase approximation (RPA). Solutions of these models are commonly represented by series of a particular class of Goldstone diagrams so-called RPA diagrams. We present a detailed asymptotic analysis of these RPA diagrams using techniques from singular analysis and discuss their computational complexity within adaptive approximation schemes. In particular, we provide a connection between RPA diagrams and classical pseudo-differential operators which enables an efficient treatment of the linear and nonlinear interactions in these models. Finally, we discuss a best N-term approximation scheme for RPA-diagrams and provide the corresponding convergence rates.