In principle, the Luttinger-Ward Green's-function formalism allows one to compute simultaneously the total energy and the quasiparticle band structure of a many-body electronic system from first principles. We present approximate and exact expressions for the correlation energy within the $GW$-random-phase approximation that are more amenable to computation and allow for developing efficient approximations to the self-energy operator and correlation energy. The exact form is a sum over differences between plasmon and interband energies. The approximate forms are based on summing over screened interband transitions. We also demonstrate that blind extremization of such functionals leads to unphysical results: imposing physical constraints on the allowed solutions (Green's functions) is necessary. Finally, we present some relevant numerical results for atomic systems.
In principle, the Luttinger-Ward Green's-function formalism allows one to compute simultaneously the total energy and the quasiparticle band structure of a many-body electronic system from first principles. We present approximate and exact expressions for the correlation energy within the GW-random-phase approximation that are more amenable to computation and allow for developing efficient approximations to the self-energy operator and correlation energy. The exact form is a sum over differences between plasmon and interband energies. The approximate forms are based on summing over screened interband transitions. We also demonstrate that blind extremization of such functionals leads to unphysical results: imposing physical constraints on the allowed solutions (Green's functions) is necessary. Finally, we present some relevant numerical results for atomic systems.
The well-known and well-used Density Functional Theory (DFT) is a variational theory for the ground-state. There are no guarantees that it delivers good results for excited state properties such as quasiparticle energies or wave function. All that is guaranteed is that the total energy is exact, the chemical potential is correct, and that the eective electronic states (Kohn-Sham states) generate the right electron density at the variational minimum. Luttinger-Ward approaches to electronic structure are attractive as beyond-DFT schemes since, in principle, they oer a variational functional of the one-electron Green’s function that yields both the correct total energy and all one-particle properties at the variational extremum. However, compared to DFT, these methods face a number of additional challenges beyond the familiar problem of choosing an exchange-correlation functional: (a) they have a much greater computational cost compared to DFT, (b) one must choose some set of trial Green’s functions, and (c) one must locate an extremum (not necessarily a minimum). This talk aims both to be pedagogical and to provide a summary of our most recent work in this area. We focus mainly on the random-phase approximation (RPA) for the correlation functional which yields the GW self-energy for the Green’s function. We will describe ecient approaches for computing the RPA correlation energy; how various types of approximations to the correlation energy create a ladder of approximate self-energy operators and what type of physics is included in these approximations (exchange, screening, Hubbard U, etc.); and how the extremization is in fact highly problematic. We will end with some ideas for making practical progress with the extremization and their relation to the successful Quasiparticle Self-consistent GW (QSGW) scheme [1].
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 .
In this lecture we present many-body perturbation theory as a method to determine quasiparticle excitations in solids, especially electronic band structures, accurately from first principles.The main ingredient is the electronic self-energy that, in principle, contains all many-body exchange and correlation effects beyond the Hartree potential.As its exact mathematical expression is unknown, approximations must be used in practical calculations.The approximation is obtained using a systematic algebraic approach on the basis of Green function techniques.It constitutes an expansion of the self-energy up to linear order in the screened Coulomb potential, which describes the interaction between the quasiparticles and includes dynamic screening through the creation of exchange-correlation holes around the bare particles.The implementation of the approximation relies on a perturbative treatment starting from density functional theory.Besides a detailed mathematical discussion we focus on the underlying physical concepts and show some illustrative applications.
The total energy and electron addition and removal spectra can, in principle, be obtained exactly from the one-body Green's function (GF). In practice, the GF is obtained from an approximate self-energy. In the framework of many-body perturbation theory, we derive different expressions that are based on an approximate self-energy, but that yield nevertheless, in principle, the exact exchange-correlation contribution to the total energy for any interaction strength. Response functions play a crucial role, which explains why, for example, ingredients from time-dependent density functional theory can be used to build these approximate self-energies. We show that the key requirement for obtaining exact results is the consistent combination of ingredients. Also when further approximations are made, as it is necessary in practice, this consistency remains the key to obtain good results. All findings are illustrated using the exactly solvable symmetric Hubbard dimer.
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.
We derive analytical expressions for chemical potentials within the random phase approximation (RPA), equivalently the GW energy functional evaluated using non interacting Green's functions (Gs). The chemical potential is obtained using two formally equivalent approaches: a direct derivative of the total energy with respect to particle number, and a functional derivative via the chain rule through Gs, both validated with finite difference benchmarks. We show that the functional derivative of the GW correlation energy$\unicode{x2013}$i.e., the GW correlation self energy$\unicode{x2013}$exhibits a discontinuity at integer particle numbers with finite jumps. This resolves the apparent inconsistency between accurate GW quasiparticle energies and the large delocalization errors observed in RPA total energies, as standard GW self energies neglect this nonanalytic behavior. Our results suggest that derivative discontinuities are a fundamental feature of correlation energy functionals, analogous to the known discontinuity in the exact exchange correlation energy.