We present a general procedure for obtaining progressively more accurate functional expressions for the electron self-energy by iterative solution of Hedin's coupled equations. The iterative process starting from Hartree theory, which gives rise to the $\mathrm{GW}$ approximation, is continued further, and an explicit formula for the vertex function from the second full cycle is given. Calculated excitation energies for a Hubbard Hamiltonian demonstrate the convergence of the iterative process and provide further strong justification for the $\mathrm{GW}$ approximation.
We present a general procedure for obtaining progressively more accurate functional expressions for the electron self-energy by iterative solution of Hedin's coupled equations. The iterative process starting from Hartree theory, which gives rise to the GW approximation, is continued further, and an explicit formula for the vertex function from the second full cycle is given. Calculated excitation energies for a Hubbard Hamiltonian demonstrate the convergence of the iterative process and provide further strong justification for the GW approximation.
In many ways the solution to the Hedin equations represents an exact solution to the many body problem. However, for most systems of practical interest, the solution to the Hedin equations is rendered nearly numerically intractable because the Hedin equations are of functional derivative form. Integral equations, on the other hand, are much more numerically tractable, than functional derivative equations, as they can often be solved iteratively. In this work we present a systematic set of integral equations (with no functional derivatives) - Hedin approximations I, II, III, IV etc. - whose solutions converge to the solutions of the exact Hedin equations. The Hedin approximations are well suited to iterative numerical solutions (which we also describe). Furthermore Hedin approximation I is just the GW approximation (as such this work may be viewed as a systematic improvement of the GW approximation). We present a systematic study of the Hedin equations for zero dimensional field theory (which, in particular, is a method to enumerate Feynman diagrams for field theories in arbitrary dimensions) and show better and better convergence to the solutions of the Hedin equations for higher and higher Hedin approximations, with Hedin approximations I, II and III being explicitly studied. We, in particular, show that the higher Hedin approximations capture more and and more Feynman diagrams for the self energy. We also show that already Hedin approximation II captures more diagrams than the state of the art diagrammatic vertex corrections approach. Furthermore Hedin approximation III is a near perfect match to the exact solutions of the Hedin equations, at least in the zero dimensional case, and enumerates a large number of Feynman diagrams.
Self-consistent solutions of Hedin's equations (HE) for the two-site Hubbard model (HM) have been studied. They have been found for three-point vertices of increasing complexity (Γ = 1 (GW approximation), Γ1 from the first-order perturbation theory, and the exact vertex Γ(E)). Comparison is made between the cases when an additional quasiparticle (QP) approximation for Green's functions is applied during the self-consistent iterative solving of HE and when QP approximation is not applied. The results obtained with the exact vertex are directly related to the present open question-which approximation is more advantageous for future implementations, GW + DMFT or QPGW + DMFT. It is shown that in a regime of strong correlations only the originally proposed GW + DMFT scheme is able to provide reliable results. Vertex corrections based on perturbation theory (PT) systematically improve the GW results when full self-consistency is applied. The application of QP self-consistency combined with PT vertex corrections shows similar problems to the case when the exact vertex is applied combined with QP sc. An analysis of Ward Identity violation is performed for all studied in this work's approximations and its relation to the general accuracy of the schemes used is provided.
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.
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.
Within many-body perturbation theory, Hedin's formalism offers a systematic way to iteratively compute the self-energy $\mathrm{\ensuremath{\Sigma}}$ of any dynamically correlated interacting system, provided one can evaluate the interaction vertex $\mathrm{\ensuremath{\Gamma}}$ exactly. This is, however, impossible, in general, for it involves the functional derivative of $\mathrm{\ensuremath{\Sigma}}$ with respect to the Green's function. Here, we analyze the structure of this derivative, splitting it into four contributions and outlining the type of quasiparticle interactions that each of them generate. Moreover, we show how, in the implementation of self-consistency, the action of these contributions can be classified into two: A quantitative renormalization of previously included interaction terms and the inclusion of qualitatively distinct interaction terms through successive functional derivatives of $\mathrm{\ensuremath{\Gamma}}$ itself. Implementing this latter type of self-consistency can extend the validity of perturbative approximations based on Hedin's equations toward the high interaction limit, as we show in the example of the Hubbard dimer. Our analysis also provides a unifying perspective on the perturbation theory landscape, showing how the T-matrix approach is completely contained in Hedin's formalism.