It is shown that the procedure suggested in a previous work [Ⅰ]~[I] provides a method to calculate the off-energy-shell property of all the G-matrices involved in an infinite series of G-matrix diagrams in a closed form. Explicit expressions are derived for (1) the random phase approximation (RPA), (2) the RPA self-scrcening of the particle-hole interaction and (3) the renormalization of the effective partlcle-particle interaction by core polarization processes which include the RPA self-screening of the particle-hole interaction to all orders.
It is shown that the procedure suggested in a previous work [Ⅰ]~[I] provides a method to calculate the off-energy-shell property of all the G-matrices involved in an infinite series of G-matrix diagrams in a closed form. Explicit expressions are derived for (1) the random phase approximation (RPA), (2) the RPA self-scrcening of the particle-hole interaction and (3) the renormalization of the effective partlcle-particle interaction by core polarization processes which include the RPA self-screening of the particle-hole interaction to all orders.
We present a method for the treatment of correlations in finite Fermi systems in a more consistent way than the random phase approximation (RPA). The main differences with respect to previous approaches are underlined and discussed. In particular, no use of renormalized operators is made. By means of the method of linearization of equations of motion, we derive a set of RPA-like equations depending only on the one-body density matrix. The latter is no more assumed to be diagonal and it is expressed in terms of the X and Y amplitudes of the particle-hole phonon operators. This set of nonlinear equations is solved via an iterative procedure, which allows us to calculate the energies and the wave functions of the excited states. After presenting our approach in its general formulation, we test its quality by using metal clusters as a good test laboratory for a generic many-body system. Comparison to RPA shows significant improvements. We discuss also how the present approach can be further extended beyond the particle-hole configuration space, getting an approximation scheme in which energy weighted sum rules are exactly preserved, thus solving a problem common to all extension of RPA proposed until now. The implementation of such approach will be afforded in a future work.
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.
Besides the proof of a modified Luttinger–Ward–Klein variational principle and a related self-consistency idea, there is not much new in principle in this paper. L. Hedin, Phys. Rev . 139, A796–823 (1965) Summary In this chapter a set of equations is formulated that determine the self-energy and the one-body Green's function in terms of the screened Coulomb interaction between classical charges. The equations contain a correction to the classical picture in terms of a vertex function. The physical meaning of the various contributions is discussed. The simplest approximation for the vertex yields the random phase approximation for the polarizability and the GW approximation for the self-energy. Various aspects of the GWA are analyzed, with a focus on the physics that is added beyond Hartree– Fock. A brief summary of model cases illustrates the domain of validity and the limits of the GWA. In this chapter we elaborate in more detail on the question of how to calculate the onebody Green's function from a Dyson equation with a self-energy kernel. In the previous chapter a scheme was introduced to design approximations to the self-energy. However, the question of where to stop, which pieces of physics to include and which to neglect, is not yet settled. Of course, there is no unique answer, besides the exact solution, but different strategies can be more or less advantageous in practice. In a system with a few electrons, for example, different aspects will be important than in a system with many electrons. Here we are mostly interested in solids, or more generally in extended systems. In such systems, screening plays an essential role: the interaction between two charges is strongly modified, in general reduced, by the rearrangement of all the other charges. It is therefore most convenient to reformulate the equations such that screening appears explicitly. Some steps in this direction can be found in earlier chapters, in particular in Sec. 8.3, the formulation of the Ψ [ G , W ]-functional of the screened interaction W instead of the Ф[ G , v c ]-functional of the bare v c . In Sec. 10.5 the screened interaction approximation for the self-energy is derived, with the self-energy as a product of the one-body Green's function and the screened interaction W .
A self-consistent RPA method is presented for N-particle Fermi systems with two- and three-body forces. The RPA equations for the generalised 'particle-hole' and 'particle-particle' problem, respectively, are derived by minimising the mean excitation energies defined analogously to mean single-particle energies. This corresponds to the equation-of-motion method. For the generalised 'one-particle/two-particle-one-hole problem' the Green-function scheme is used in order to establish a connection between two possible ways of treating the problem. It is shown that the RPA Hamiltonians have almost the same structure as in the case of pure two-body forces, and that Rowe's general RPA equations are contained as a special case of the treatment.
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.
A general theory of long range correlations in the ground state of a closed shell nucleus is outlined, using a number of ideas introduced or adumbrated, but not always successfully developed, in previous work of the authors and others. These include the following: (i) Any operator can be formally expressed as a series in particle-hole (p-h) excitation and destruction operators which should converge rapidly for a suitable choice of single-particle basis at the closed shells. By means of these relations and with the help of sum rules all observables can be computed without the explicit use of wave functions. (ii) For a fixed choice of single-particle basis, the ground state energy can be computed from a variational principle in which in lowest nontrivial approximation, in consequence of (i), the parameters are amplitudes of the same type as occur in the random phase approximation (RPA). The precise relation to RPA is fully detailed. (iii) The single particle basis chosen is the natural orbital basis of Lowdin, which presents itself as the obvious choice both for formulation of the dynamics and for the algorithm of solution. (iv) Though the theory deals with pairs of fermion operators, particular attention is paid that all Pauli principle restrictions are satisfied. (v) It is verified that in the perturbation limit, known results are reproduced. The method is generally applicable to any system for which the correlations can be viewed as particle-hole excitations with respect to a reference Slater determinant.