In this thesis, the GW approximation (GWA, Green's function G times screened interaction W) and the Hubbard I approximation (HIA) are combined in a non-self-consistent one-shot calculation to determine the electronic structure of a one-dimensional strongly-correlated model. The scheme was chosen to incorporate both screening effects through the GWA and strong on-site correlations through the atomic HIA. The resulting self-energies of both methods are summed together and another term is subtracted to correct double-counting. The examined double-counting terms were the local self-energy of the GWA result, the GWA self-energy of the atomic model used in the HIA and lastly the impurity self-energy within the GWA. The self-energy of the atomic model in the HIA was solved analytically. It was found that all three approaches yield non-causal features, which increase with the on-site Coulomb repulsion, in the resulting spectral function. The local GWA self-energy was observed to perform best in terms of causality and computational effort. Changes in the resulting quasi-particle structure in the spectral function showed that screening effects and local correlations were included successfully.
In this thesis, the GW approximation (GWA, Green's function G times screened interaction W) and the Hubbard I approximation (HIA) are combined in a non-self-consistent one-shot calculation to determine the electronic structure of a one-dimensional strongly-correlated model. The scheme was chosen to incorporate both screening effects through the GWA and strong on-site correlations through the atomic HIA. The resulting self-energies of both methods are summed together and another term is subtracted to correct double-counting. The examined double-counting terms were the local self-energy of the GWA result, the GWA self-energy of the atomic model used in the HIA and lastly the impurity self-energy within the GWA. The self-energy of the atomic model in the HIA was solved analytically. It was found that all three approaches yield non-causal features, which increase with the on-site Coulomb repulsion, in the resulting spectral function. The local GWA self-energy was observed to perform best in terms of causality and computational effort. Changes in the resulting quasi-particle structure in the spectral function showed that screening effects and local correlations were included successfully.
The main focus of this thesis is the detailed investigation of computational methods to tackle strongly correlated materials in which a rich variety of exotic phenomena are found. A many-body problem with sizable electronic correlations can no longer be explained by independent-particle approximations such as density functional theory (DFT) or tight-binding approaches. The influence of an electron to the others is too strong for each electron to be treated as an independent quasiparticle and consequently those standard band-structure methods fail even at a qualitative level. One of the most powerful approaches for strong correlations is the dynamical mean-field theory (DMFT), which has enlightened the understanding of the Mott transition based on the Hubbard model. For realistic applications, the dynamical mean-field theory is combined with various independent-particles approaches. The most widely used one is the DMFT combined with the DFT in the local density approximation (LDA), so-called LDA+DMFT. In this approach, the electrons in the weakly correlated orbitals are calculated by LDA while others in the strongly correlated orbitals are treated by DMFT. Recently, the method combining DMFT with Hedin's GW approximation was also developed, in which the momentum-dependent self-energy is also added. In this thesis, we discuss the application of those methodologies based on DMFT. First, we apply the dynamical mean-field theory to solve the 3-dimensional Hubbard model in Chap. 3. In this application, we model the interface between the thermodynamically coexisting metal and Mott insulator. We show how to model the required slab geometry and extract the electronic spectra. We construct an effective Landau free energy and compute the variation of its parameters across the phase diagram. Finally, using a linear mixture of the density and double-occupancy, we identify a natural Ising order parameter which unifies the treatment of the bandwidth and filling controlled Mott transitions. Secondly, we study the double-counting problem, a subtle issue that arises in LDA+DMFT. We propose a highly precise double-counting functional, in which the intersection of LDA and DMFT is calculated exactly, and implement a parameter-free version of the LDA+DMFT that is tested on one of the simplest strongly correlated systems, the H2 molecule. We show that the exact double-counting treatment along with a good DMFT projector leads to very accurate and total energy and excitation spectrum of H2 molecule. Finally, we implement various versions of GW+DMFT, in its fully self-consistent way, one shot GW approximation, and quasiparticle self-consistent scheme, and studied how well these combined methods perform on H2 molecule as compared to more established methods such as LDA+DMFT. We found that most flavors of GW+DMFT break down in strongly correlated regime due to causality violation. Among GW+DMFT methods, only the self-consistent quasiparticle GW+DMFT with static double-counting, and a new method with causal double-counting, correctly recover the atomic limit at large H-atom separation. While some flavors of GW+DMFT improve the single-electron spectra of LDA+DMFT, the total energy is best predicted by LDA+DMFT, for which the exact double-counting is known, and is static.
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.
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 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 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.
Many observables such as the density, total energy, or electric current, can be expressed explicitly in terms of the one‐body Green's function, which describes electron addition or removal to or from a system. An efficient way to determine such a Green's function is to introduce a self‐energy, which is a nonlocal and dynamic effective potential that influences the propagation of particles in an interacting system. The state‐of‐the art approximation for the self‐energy is the GW approximation, where the system to (or from) which the electron is added (or removed) is described as a polarizable, screening, medium. This is expressed by the name of the approximation: ‘GW’ stands for the one‐body Green's function G and for W , the dynamically screened Coulomb interaction. The GW approximation is very popular for the calculation of band structures in solids, and increasingly used also to describe nanostructures, clusters, and molecules. As compared to static mean‐field approximations for the effective potential, the dynamical screening of the Coulomb interaction in GW leads to a renormalization of energies, to broadening and/or to the observation of additional excitations. An analysis of the approximations that lead to the GW self‐energy, and of the underlying picture, explains the successes and the limitations of the approach. This article is categorized under: Electronic Structure Theory > Density Functional Theory Electronic Structure Theory > Ab Initio Electronic Structure Methods Theoretical and Physical Chemistry > Spectroscopy Structure and Mechanism > Computational Materials Science