We study the effects of local vertex corrections to the self energy of the electron gas. We find that a vertex derived from time-dependent density-functional theory can give accurate self energies without including the explicit time dependence of the exchange-correlation potential provided, however, that a proper decay at large momentum transfer (large q) is built into the vertex function. (The local-density approximation for the vertex fails badly.) Total energies are calculated from the Galitskii-Migdal formula and it is shown that a proper large-q behavior, results in a close consistency between the chemical potentials derived from these energies and those obtained directly from the self energy. We show that this internal consistency depends critically on including the same vertex correction in both the self-energy and the screening function. In addition the total energies become almost as accurate as those from elaborate quantum Monte-Carlo (QMC) calculations. We also study the accuracy and utility of the functional for the total energy proposed by Luttinger and Ward and a generalization by Almbladh, von Barth, and van Leeuwen. For the electron gas, even the simplest and readily evaluated approximations to these functionals yield total energies of similar quality as those of QMC calculations. The functionals depend on the one-electron Green's function and the screened Coulomb interaction and already rather crude approximations to these quantities produce accurate energies thus demonstrating the insensitivity of the functionals to their arguments. Different ways of incorporating vertex corrections beyond the $GW$ level are studied in simple, exactly soluble polaron-like models. We study models of a structureless core electron coupled to valence electrons and a local polaron model by Cini, Hewson and Newns. Our model results indicate that the first vertex correction alone will in general not suffice to improve the spectrum away from the quasi-particle peak. By including a subsequence of Mahan's fractal vertex series, however, we obtain results with correct physical properties which agree better with exact model results.
We study the effects of local vertex corrections to the self energy of the electron gas. We find that a vertex derived from time-dependent density-functional theory can give accurate self energies without including the explicit time dependence of the exchange-correlation potential provided, however, that a proper decay at large momentum transfer (large q) is built into the vertex function. (The local-density approximation for the vertex fails badly.) Total energies are calculated from the Galitskii-Migdal formula and it is shown that a proper large-q behavior, results in a close consistency between the chemical potentials derived from these energies and those obtained directly from the self energy. We show that this internal consistency depends critically on including the same vertex correction in both the self-energy and the screening function. In addition the total energies become almost as accurate as those from elaborate quantum Monte-Carlo (QMC) calculations. We also study the accuracy and utility of the functional for the total energy proposed by Luttinger and Ward and a generalization by Almbladh, von Barth, and van Leeuwen. For the electron gas, even the simplest and readily evaluated approximations to these functionals yield total energies of similar quality as those of QMC calculations. The functionals depend on the one-electron Green's function and the screened Coulomb interaction and already rather crude approximations to these quantities produce accurate energies thus demonstrating the insensitivity of the functionals to their arguments. Different ways of incorporating vertex corrections beyond the GW level are studied in simple, exactly soluble polaron-like models. We study models of a structureless core electron coupled to valence electrons and a local polaron model by Cini, Hewson and Newns. Our model results indicate that the first vertex correction alone will in general not suffice to improve the spectrum away from the quasi-particle peak. By including a subsequence of Mahan's fractal vertex series, however, we obtain results with correct physical properties which agree better with exact model results.
It is commonly accepted that the GW approximation for the electron self-energy is successful for the description of the band structure of weakly to moderately correlated systems, whereas it will fail for strongly correlated materials. In the present work, we discuss two important aspects of this approximation: first, the "self-screening error," which is due to an incorrect treatment of induced exchange, and second, the atomic limit, in which, instead, correlation is directly responsible for the observed problem. Using the example of the removal of a particle from a box, we show that the self-screening error stems from the use of test charge-test charge screening and that it can be corrected by a two-point vertex contribution to the self-energy derived from time-dependent density functional theory (TDDFT). We explain why the addition of a particle, instead, requires the use of a different approximate vertex. This illustrates why the general vertex function, valid both for valence and conduction states, must be a three-point function. Moreover, we show that also the bad performance of GW in the atomic limit is due to the neglect of the vertex in the self-energy; in that case, the TDDFT-derived vertex correction is not sufficient in order to remove the error even qualitatively. We discuss the effects of the self-screening error as well as the atomic limit using GW for the exactly solvable two-site Hubbard model.
Within many-body perturbation theory, we apply vertex corrections to various closed-shell atoms and to jellium, using a local approximation for the vertex consistent with starting the many-body perturbation theory from a Kohn-Sham Green's function constructed from density-functional theory in the local-density approximation. The vertex appears in two places---in the screened Coulomb interaction W and in the self-energy $\ensuremath{\Sigma}$---and we obtain a systematic discrimination of these two effects by turning the vertex in $\ensuremath{\Sigma}$ on and off. We also make comparisons to standard GW results within the usual random-phase approximation, which omits the vertex from both. When a vertex is included for closed-shell atoms, both ground-state and excited-state properties demonstrate little improvement over standard GW. For jellium, we observe marked improvement in the quasiparticle bandwidth when the vertex is included only in W, whereas turning on the vertex in $\ensuremath{\Sigma}$ leads to an unphysical quasiparticle dispersion and work function. A simple analysis suggests why implementation of the vertex only in W is a valid way to improve quasiparticle energy calculations, while the vertex in $\ensuremath{\Sigma}$ is unphysical, and points the way to the development of improved vertices for ab initio electronic structure calculations.
We introduce a scheme through which the ground-state electronic density n and total energy Ev [n] of inhomogeneous interacting systems can be determined in a self-consistent (SC) way. It is based on a SC perturbation expansion for the single-particle Green function G. Our employed expansion is argued to be valid under all conditions provided that the ground-state charge density of the interacting system is pure-state non-interacting v representable. The formalism discussed in this paper has been applied to a quasi-one-dimensional periodic system. It has even been possible to carry out SC calculations with the dynamically screened exchange diagram (the so-called GW diagram, where W is the dynamically screened electron-electron interaction function) for the self-energy operator Σ. A coherent mode of interpretation of a number of our results entails that the overwhelming success of the GW framework in relatively accurately reproducing the experimental energy gaps of a number of semiconductors and insulators may be a direct consequence of non-self-consistency effects; whereas a non-self-consistent calculation in our tests increases the value of the local-density approximation (LDA) gap in the range of 40–65% of the ‘experimental’ gaps, self-consistency effects reduce this range to 14–17%. The same line of argument leads us to conclude that, in our model, 26–35% of the LDA error in the gap energies is due to the absence of a well known discontinuity in the LDA exchange-correlation potential, and the remaining part is due to the Kohn-Sham eigenvalues not being identifiable with one-particle excitation energies. A natural conclusion that may be drawn from these observations is that a correct description of excitation energies (e.g. of the gap energy) cannot be achieved without including, at least approximately, the vertex part of Σ, which is entirely absent in the GW scheme.
With the aim of identifying universal trends, we compare fully self-consistent electronic spectra and total energies obtained from the GW approximation with those from an extended $\mathrm{GW}\ensuremath{\Gamma}$ scheme that includes a nontrivial vertex function and the fundamentally distinct Bethe-Goldstone approach based on the T matrix. The self-consistent Green's function G, as derived from Dyson's equation, is used not only in the self-energy but also to construct the screened interaction W for a model system. For all approximations we observe a similar deterioration of the spectrum, which is not removed by vertex corrections. In particular, satellite peaks are systematically broadened and move closer to the chemical potential. The corresponding total energies are universally raised, independent of the system parameters. Our results, therefore, suggest that any improvement in total energy due to self-consistency, such as for the electron gas in the GW approximation, may be fortuitous.
With the aim of constructing an electronic structure approach that systematically goes beyond the GW and random phase approximation (RPA) we introduce a vertex correction based on the exact-exchange (EXX) potential of time-dependent density functional theory. The EXX vertex function is constrained to be local but is expected to capture similar physics as the Hartree-Fock vertex. With the EXX vertex, we then unify different beyond-RPA approaches such as the various resummations of RPA with exchange and the second-order screened exchange approximation. The theoretical analysis is supported by numerical studies on the hydrogen dimer and the electron gas, and we discuss the role of including the vertex correction in both the screened interaction and the self-energy. Finally, we give details on our implementation within the plane-wave pseudo potential framework and demonstrate the excellent performance of the different RPA with exchange methods in describing the energetics of hydrogen and van der Waals bonds.
In this chapter, shortcomings of the GW approximation are analyzed, and different ways are proposed to go beyond. The approaches fall into three broad classes: first, one can iterate Hedin's equations that link the self-energy, the polarizability, the screened Coulomb interaction, and the vertex function with the Dyson equation for the one-body Green's function. Second, one can use other Dyson equations like in the T-matrix approach, and third, one can make an ansatz for the one-body Green's function inspired by physical considerations or insight from models, in particular the idea of electron–boson coupling that leads to the cumulant expansion. Applications are given for illustration. Some more ideas and links to subsequent chapters provide an outlook. The GW approximation is a widely used approach in computational electronic structure. However, the previous chapters show some severe limitations. One example is systems with low-energy excitations, as illustrated by the Hubbard dimer in the dissociation limit, where bonding and antibonding orbitals become degenerate. This can be seen in Sec. 11.7. Another example is satellite series, for example in the photoemission spectrum of silicon in Sec. 13.1. Moreover, even in a situation where the approximation is well justified, it induces some error, though maybe small, and one eventually wants to do better. The full self-energy functional is in principle known: it is the Luttinger–Ward sum of skeleton diagrams derived from Eq. (11.35). However, this is an infinite sum, and no closed form of the result is known. Hence, the task is not to invent new diagrams: it is to make a choice of diagrams or combinations of diagrams, which are evaluated fully or in an approximate way. This is the framework of the present chapter. In some cases we will add diagrams to the GWA ones, which gives rise to vertex corrections ; in other cases some of the GWA diagrams are neglected while different ones are taken into account. An example is the T-matrix approximation, which sums ladder diagrams instead of bubbles. In general one tries to combine this with the GWA, especially in extended systems where screening is important. Therefore the GWA plays a prominent role, and the whole chapter is called beyond GW . The chapter contains advanced material, that is in part exploratory and refers to current research. It is intended to make links to communities outside the GW world, and to other chapters.