The GW compendium: A practical guide to theoretical photoemission spectroscopy Dorothea Golze, Marc Dvorak, and Patrick Rinke Department of Applied Physics, Aalto University School of Science, 00076-Aalto, Finland (Dated: December 11, 2019) The GW approximation in electronic structure theory has become a widespread tool for predicting electronic excitations in chemical compounds and materials. In the realm of theoretical spectroscopy, the GW method provides access to charged excitations as mea- sured in direct or inverse photoemission spectroscopy. The number of GW calculations in the past two decades has exploded with increased computing power and modern codes. The success of GW can be attributed to many factors: favorable scaling with respect to system size, a formal interpretation for charged excitation energies, the importance of dynamical screening in real systems, and its practical combination with other theories. In this review, we provide an overview of these formal and practical considerations. We expand, in detail, on the choices presented to the scientist performing GW calculations for the first time. We also give an introduction to the many-body theory behind GW , a review of modern applications like molecules and surfaces, and a perspective on methods which go beyond conventional GW calculations. This review addresses chemists, physi- cists and material scientists with an interest in theoretical spectroscopy. It is intended for newcomers to GW calculations but can also serve as an alternative perspective for experts and an up-to-date source of computational techniques. CONTENTS B. Eigenvalue self-consistency and level alignment 32 C. Self-consistency via a new ground state 33 I. Introduction 2 VI. Solids 34 II. Theoretical spectroscopy 3 A. Band gaps 35 A. Direct and inverse photoemission spectroscopy 3 B. Band structures and band parameters 36 B. The quasiparticle concept 4 C. Lifetimes 37 C. Comparison to experimental spectra 6 D. More challenging solids 38 E. Defects in solids 39 III. Hedin's GW equations 7 F. Outlook on solids 40 IV. The G0W0 approach: Concept and implementation 11 VII. Surfaces 40 A. The G0W0 equations 11 B. Procedure 12 VIII. Two-dimensional materials 42 C. Frequency treatment 14 1. Plasmon-pole models 14 IX. Molecules 44 2. Contour deformation 15 A. First ionization potentials and electron affinities 44 B. Ionization spectra 45 3. Analytic Continuation 16 C. The GW 100 benchmark set 46 4. Fully Analytic Approach 17 D. Molecular crystals 47 5. Comparison of accuracy and computational cost 17 D. Basis sets 18 X. Total energy and the electronic ground state 48 1. Plane waves 18 A. Electron density 48 2. Localized basis sets 18 B. Total energy 49 3. Augmented basis sets 19 4. Pseudopotentials 20 XI. Current challenges and beyond GW 50 5. Projector augmented-wave method (PAW) 21 A. Challenges 50 E. Basis set convergence 22 B. Quantum chemistry 50 arXiv:1912.04893v1 [cond-mat.mtrl-sci] 10 Dec 2019 F. Elimination of unoccupied state summation 23 C. Non-equilibrium Green's functions 51 G. Starting point dependence and optimization 25 D. Vertex corrections 51 1. How severe is the dependence on the reference E. Optical properties 53 state? 25 F. T -matrix 53 2. Consistent starting point scheme 26 G. Cumulant expansion 54 3. Deviation from the straight line scheme 26 H. Local vertex 55 4. IP-theorem schemes 27 H. Computational complexity and cost 28 XII. Conclusion 55 I. Practical guidelines 29 XIII. Acknowledgements 55 V. Beyond G0W0: self-consistency schemes 30 A. Fully self-consistent GW 30 A. Integral notation 56 2 B. The many-body problem 56 its much higher computational expense. Soon after it was 1. Green's function formalism 56 realized that GW can overcome some of the most noto- 2. Hedin's equations 57 rious deficiencies of common density functionals such as C. Computational details 58 the self-interaction error, the absence of long-range po- larization effects and the Kohn-Sham band-gap problem. References 59 The GW approach is now an integral part of elec- tronic structure theory and readily available in major electronic structure codes. It is taught at summer schools I. INTRODUCTION along side DFT and other electronic structure methods. This review is intended as a tutorial that complements Electronic structure theory derives from the fundamen- showcases of GW 's achievements with a practical guide tal laws of quantum mechanics and describes the behav- through the theory, its implementation and actual use. ior of electrons { the glue that shapes all matter. To For GW novices, the review offers a gentle introduction understand the properties of matter and the behavior of to the GW concept and its application areas. Regular molecules, the quantum mechanical laws must be solved GW users can consult this review as a handbook in their numerically because a pen and paper solution is not pos- day-to-day use of the GW method. For seasoned GW sible. In this context, Hedin's GW method (Hedin, 1965) users and GW experts it might serve as a reference for has become the de facto standard for electronic structure key applications and some of the subtler points of the properties as measured by direct and inverse photoemis- GW framework. sion experiments, such as quasiparticle band structures In this review we take a more practical approach to- and molecular excitations. wards the GW method. We will recap the basic the- Electronic structure theory covers the quantum me- ory starting from theoretical spectroscopic view point chanical spectrum of computational materials science as a probe of the electronic structure. Aiming at GW and quantum chemistry. The fundamental aim of compu- practitioners, we will illustrate how the GW approach tational science is to derive understanding entirely from emerges from the theoretical spectroscopy framework as the basic laws of physics, i.e. quantum mechanical first a practical scheme for electronic structure calculations. principles, and increasingly also to make predictions of A more in-depth discussion of the theoretical founda- new properties or new materials and new molecules for tions of many-body theory can be found in textbooks, specific tasks. The exponential increase in available com- e.g., (Bechstedt, 2014; Fetter and Walecka, 1971; Gross puter power and new methodological developments are et al., 1991; Martin et al., 2016; Szabo and Ostlund, two major factors in the growing impact of this field for 1989), while the GW theory itself is covered in excel- practical applications to real systems. As a result of these lent early reviews (Aryasetiawan and Gunnarsson, 1998; advances, computational science is establishing itself as Hedin, 1999; Hedin and Lundqvist, 1970). GW began a viable complement to the purely experimental and the- to flourish at the beginning of the 21st century and two oretical sciences. reviews succinctly summarized the state of the field until Hedin published the GW method in 1965, in the same then (Aulbur et al., 2000; Onida et al., 2002). Our re- time period as the foundational density-functional theory view bridges the ensuing gap of almost twenty years, af- (DFT) papers (Hohenberg and Kohn, 1964; Kohn and ter which only several more specialized reviews addressed Sham, 1965). While DFT has shaped the realm of first different aspects and applications of GW calculations principles computational science like no other method to- (Bruneval and Gatti, 2014; Faber et al., 2014; Gerosa day, GW 's fame took a little longer to develop1. DFT's et al., 2018a; Giantomassi et al., 2011; Kang et al., 2019; success has been facilitated by the computational effi- Marom, 2017; Ping et al., 2013; Rinke et al., 2008a) and ciency of the local-density (Kohn and Sham, 1965) or complements a recent review (Reining, 2017). generalized gradient approximations (Becke, 1988; Lee A considerable part of our review is devoted to the et al., 1988; Perdew et al., 1996a) (LDA and GGA) of presentation of different GW implementations. We will the exchange-correlation functional that make DFT ap- discuss the practical considerations that GW users have plicable to polyatomic systems containing up to several to make when they decide on a particular GW imple- thousand atoms. GW , however, only saw its first appli- mentation or code for their work. Moreover, we will cations to realistic materials 20 years after its inception guide the reader through computational decisions that (Godby et al., 1986; Hybertsen and Louie, 1985), due to might affect the accuracy of their GW calculations and illustrate them with concrete examples from the GW lit- erature. An important aspect in this regard is the issue of self-consistency in GW , which we cover in detail. 1 The theory directly comparable to DFT is the Luttinger-Ward formalism, not the GW approximation. Here, we focus on the Although the GW method might be best known for chronological development of practical electronic structure cal- its success in predicting the band gaps of solids, we will culations with Kohn-Sham DFT or GW . present its diversity and discuss a range of different quan- 3 tities that can be computed with the GW method. Since a Photoemission b Inverse Photoemission no method is perfect, we will conclude with a critical outlook on the challenges faced by the GW method and - discuss ways to go beyond GW . Ekin E We conclude this introduction with a quote from H. hv kin Evac Evac J. Monkhorst, who wrote in a laudation in 2005: It is therefore my conviction that, rather sooner then later, we hv ϵs - will see a resurgence of the precise many-body approach Φ Φ to solid-state theory as we envisioned. Almost assuredly -EAs EF EF the GW method will be the tool of choice. (Monkhorst, - - - - 2005). In 2019, we can say that Monkhorst was right.