The concepts of quasiparticle and collective excitation are considered in the context of continued-fraction expansions. This leads to a practical scheme allowing a highly accurate description of the screened interaction function W and the electron self-energy \ensuremath{\Sigma}, which can be achieved at a reasonable computational cost. The usefulness of this scheme is demonstrated on a simple two-band model of a semiconductor, where it leads to an unambiguous calculation of the electron self-energy in the GW approximation. Conceptually, the theory clarifies some of the limitations of the standard quasiparticle concept in many-body theory. Furthermore, it allows a proper assessment of the quality of the standard approximations employed in GW calculations.
The concepts of quasiparticle and collective excitation are considered in the context of continued-fraction expansions. This leads to a practical scheme allowing a highly accurate description of the screened interaction function W and the electron self-energy \ensuremath{\Sigma}, which can be achieved at a reasonable computational cost. The usefulness of this scheme is demonstrated on a simple two-band model of a semiconductor, where it leads to an unambiguous calculation of the electron self-energy in the GW approximation. Conceptually, the theory clarifies some of the limitations of the standard quasiparticle concept in many-body theory. Furthermore, it allows a proper assessment of the quality of the standard approximations employed in GW calculations.
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.
We present self-consistent GW calculations for a quasi-one-dimensional model semiconductor, using multipole representations for W, the screened interaction, and G, the electron Green function. In the case of G, we distinguish between the quasiparticle contribution and the so-called incoherent contribution. We consider various strengths of the external potential. For two different starting points of the self-consistency cycle, local-density approximation and Hartree-Fock, the band gaps converge to the same values. The self-consistent GW band gaps differ very little from the Hartree-Fock band gaps. There is considerable disagreement with quantum Monte Carlo calculations for the same model, indicating the importance of vertex corrections.
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
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.
Using the simple (symmetric) Hubbard dimer, we analyze some important features of the GW approximation. We show that the problem of the existence of multiple quasiparticle solutions in the (perturbative) one-shot GW method and its partially self-consistent version is solved by full self-consistency. We also analyze the neutral excitation spectrum using the Bethe-Salpeter equation (BSE) formalism within the standard GW approximation and find, in particular, that 1) some neutral excitation energies become complex when the electron-electron interaction U increases, which can be traced back to the approximate nature of the GW quasiparticle energies; 2) the BSE formalism yields accurate correlation energies over a wide range of U when the trace (or plasmon) formula is employed; 3) the trace formula is sensitive to the occurrence of complex excitation energies (especially singlet), while the expression obtained from the adiabatic-connection fluctuation-dissipation theorem (ACFDT) is more stable (yet less accurate); 4) the trace formula has the correct behavior for weak (i.e., small U) interaction, unlike the ACFDT expression.
Density functional theory (DFT) performs reasonably well for the determination of structural properties of many materials, but fails to predict electronic band gap values accurately. \nSuch a failure of DFT is not unexpected, since there exists no formal justification for interpreting the DFT eigenvalues as addition or removal energies of the many-body system (quasiparticle energies).\nAn alternative approach to the study of exchange-correlation effects in many-particle systems is provided by many-body perturbation theory (MBPT), which defines a rigorous approach to description of excited-state properties based on the Green's function formalism and the concept of quasiparticle electrons. Within MBPT, one can calculate the quasiparticle (QP) energies and the QP amplitudes by approximately solving, within the so-called GW approximation, the set of coupled integro-differential equations proposed by Hedin in 1965.\nThe GW method generally yields significantly better values for QP energies with respect to DFT as it accounts for the dynamic many-body effects in the electron-electron interaction going beyond the mean-field approximation of independent Kohn-Sham particles. The drawback of such an involved approach, however, is its large computational cost, which is mainly due to the evaluation of dielectric matrices, their inversion, and the solution of a non-Hermitian and nonlinear eigenvalue problem. \nThe first part of this PhD research is devoted to the implementation of an efficient and scalable coarse-grained parallel algorithms and the development of state-of-the-art methods to solve the GW equations.\nAll these techniques are then applied to the study of the quasiparticle band structures of SiO2 in the alpha quartz crystalline structure. The effects of the different approximations involved in the theory, the influence of self-consistency, and a systematic analysis of the reliability of the different plasmon-pole models are presented and discussed. \nThe second part of this PhD work is dedicated to the implementation of a methodological approach to the solution of the GW equations based on the so-called projector augmented wave method (PAW) proposed by Blochl in 1994. This new approach permits to remove many of the limitations intrinsic to the use of the pseudopotential technique widely used for ab-initio calculations, allowing one to reach GW results closer to an all-electron implementation while still maintaining computational efficiency, and flexibility. This allows us to achieve a coherent implementation of many different algorithms which enables a detailed comparison of GW calculations not yet performed until now. Results for the quasiparticle band structure and optical spectra of prototype s- and p-compounds are discussed and compared to recent studies reported in the scientific literature (when available).