Self-consistent GW and higher-order calculations of electron states in metals

Reference Energies for Valence Ionizations and Satellite Transitions

Upon ionization of an atom or a molecule, another electron (or more) can be simultaneously excited. These concurrently generated states are called "satellites" (or shakeup transitions) as they appear in ionization spectra as higher-energy peaks with weaker intensity and larger width than the main peaks associated with single-particle ionizations. Satellites, which correspond to electronically excited states of the cationic species, are notoriously challenging to model using conventional single-reference methods due to their high excitation degree compared to the neutral reference state. This work reports 42 satellite transition energies and 58 valence ionization potentials (IPs) of full configuration interaction quality computed in small molecular systems. Following the protocol developed for the quest database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; and Loos, P.-F. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], these reference energies are computed using the configuration interaction using a perturbative selection made iteratively (CIPSI) method. In addition, the accuracy of the well-known coupled-cluster (CC) hierarchy (CC2, CCSD, CC3, CCSDT, CC4, and CCSDTQ) is gauged against these new accurate references. The performances of various approximations based on many-body Green's functions (GW, GF2, and T-matrix) for IPs are also analyzed. Their limitations in correctly modeling satellite transitions are discussed.

Going Beyond the GW Approximation Using the Time-Dependent Hartree-Fock Vertex

The time-dependent Hartree-Fock (TDHF) vertex of many-body perturbation theory (MBPT) makes it possible to extend TDHF theory to charged excitations. Here we assess its performance by applying it to spherical atoms in their neutral electronic configuration. On a theoretical level, we recast the TDHF vertex as a reducible vertex, highlighting the emergence of a self-energy expansion purely in orders of the bare Coulomb interaction; then, on a numerical level, we present results for polarizabilities, ionization energies (IEs), and photoemission satellites. We confirm the superiority of THDF over simpler methods such as the random phase approximation for the prediction of atomic polarizabilities. We then find that the TDHF vertex reliably provides better IEs than GW and low-order self-energies do in the light-atom, few-electron regime; its performance degrades in heavier, many-electron atoms instead, where an expansion in orders of an unscreened Coulomb interaction becomes less justified. New relevant features are introduced in the satellite spectrum by the TDHF vertex, but the experimental spectra are not fully reproduced due to a missing account of nonlinear effects connected to hole relaxation. We also explore various truncations of the self-energy given by the TDHF vertex, but do not find them to be more convenient than low-order approximations such as GW and second Born (2B), suggesting that vertex corrections should be carried out consistently both in the self-energy and in the polarizability.

Total Energy beyond GW: Exact Results and Guidelines for Approximations

The total energy and electron addition and removal spectra can, in principle, be obtained exactly from the one-body Green's function (GF). In practice, the GF is obtained from an approximate self-energy. In the framework of many-body perturbation theory, we derive different expressions that are based on an approximate self-energy, but that yield nevertheless, in principle, the exact exchange-correlation contribution to the total energy for any interaction strength. Response functions play a crucial role, which explains why, for example, ingredients from time-dependent density functional theory can be used to build these approximate self-energies. We show that the key requirement for obtaining exact results is the consistent combination of ingredients. Also when further approximations are made, as it is necessary in practice, this consistency remains the key to obtain good results. All findings are illustrated using the exactly solvable symmetric Hubbard dimer.

The three channels of many-body perturbation theory: GW, particle-particle, and electron-hole T-matrix self-energies

We derive the explicit expression of the three self-energies that one encounters in many-body perturbation theory: the well-known GW self-energy, as well as the particle-particle and electron-hole T-matrix self-energies. Each of these can be easily computed via the eigenvalues and eigenvectors of a different random-phase approximation linear eigenvalue problem that completely defines their corresponding response function. For illustrative and comparative purposes, we report the principal ionization potentials of a set of small molecules computed at each level of theory. The performance of these schemes on strongly correlated systems (B2 and C2) is also discussed.

A Similarity Renormalization Group Approach to Green's Function Methods

The family of Green's function methods based on the GW approximation has gained popularity in the electronic structure theory thanks to its accuracy in weakly correlated systems combined with its cost-effectiveness. Despite this, self-consistent versions still pose challenges in terms of convergence. A recent study [Monino and Loos J. Chem. Phys. 2022, 156, 231101.] has linked these convergence issues to the intruder-state problem. In this work, a perturbative analysis of the similarity renormalization group (SRG) approach is performed on Green's function methods. The SRG formalism enables us to derive, from first-principles, the expression of a naturally static and Hermitian form of the self-energy that can be employed in quasiparticle self-consistent GW (qsGW) calculations. The resulting SRG-based regularized self-energy significantly accelerates the convergence of qsGW calculations, slightly improves the overall accuracy, and is straightforward to implement in existing code.

Numerically Precise Benchmark of Many-Body Self-Energies on Spherical Atoms

We investigate the performance of beyond-GW approaches in many-body perturbation theory by addressing atoms described within the spherical approximation via a dedicated numerical treatment based on B-splines and spherical harmonics. We consider the GW, second Born (2B), and GW + second order screened exchange (GW+SOSEX) self-energies and use them to obtain ionization potentials from the quasi-particle equation (QPE) solved perturbatively on top of independent-particle calculations. We also solve the linearized Sham–Schlüter equation (LSSE) and compare the resulting xc potentials against exact data. We find that the LSSE provides consistent starting points for the QPE but does not present any practical advantage in the present context. Still, the features of the xc potentials obtained with it shed light on possible strategies for the inclusion of beyond-GW diagrams in the many-body self-energy. Our findings show that solving the QPE with the GW+SOSEX self-energy on top of a PBE or PBE0 solution is a viable scheme to go beyond GW in finite systems, even in the atomic limit. However, GW shows a comparable performance if one agrees to use a hybrid starting point. We also obtain promising results with the 2B self-energy on top of Hartree–Fock, suggesting that the full time-dependent Hartree–Fock vertex may be another viable beyond-GW scheme for finite systems.

Assessing the G0W0Γ0(1) Approach: Beyond G0W0 with Hedin's Full Second-Order Self-Energy Contribution

We present and benchmark a self-energy approach for quasiparticle energy calculations that goes beyond Hedin's GW approximation by adding the full second-order self-energy (FSOS-W) contribution. The FSOS-W diagram involves two screened Coulomb interaction (W) lines, and adding the FSOS-W to the GW self-energy can be interpreted as first-order vertex correction to GW (GWΓ⁽¹⁾). Our FSOS-W implementation is based on the resolution-of-identity technique and exhibits better than O(N⁵) scaling with system size for small- to medium-sized molecules. We then present one-shot GWΓ⁽¹⁾ (G₀W₀Γ₀⁽¹⁾) benchmarks for the GW100 test set and a set of 24 acceptor molecules. For semilocal or hybrid density functional theory starting points, G₀W₀Γ₀⁽¹⁾ systematically outperforms G₀W₀ for the first vertical ionization potentials and electron affinities of both test sets. Finally, we demonstrate that a static FSOS-W self-energy significantly underestimates the quasiparticle energies.

Tuplewise Material Representation Based Machine Learning for Accurate Band Gap Prediction

The open-access material databases allowed us to approach scientific questions from a completely new perspective with machine learning methods. Here, on the basis of open-access databases, we focus on the classical band gap problem for predicting accurately the band gap of a crystalline compound using a machine learning approach with newly developed tuplewise graph neural networks (TGNN), which is devised to automatically generate input representation of crystal structures in tuple types and to exploit crystal-level properties as one of the input features. Our method brings about a highly accurate prediction of the band gaps at hybrid functionals and GW approximation levels for multiple material data sets without heavy computational cost. Furthermore, to demonstrate the applicability of our prediction model, we provide a data set of GW band gaps for 45835 materials predicted by TGNN posing higher accuracy than standard density functional theory calculations.

Toward Efficient GW Calculations Using Numerical Atomic Orbitals: Benchmarking and Application to Molecular Dynamics Simulations

The use of atomic orbitals in Hedin's GW approximation provides, in principle, an inexpensive alternative to plane-wave basis sets, especially when modeling large molecules. However, benchmarking of the algorithms and basis sets is essential for a careful balance between cost and accuracy. In this paper, we present an implementation of the GW approximation using numerical atomic orbitals and a pseudopotential treatment of core electrons. The combination of a contour deformation technique with a one-shot extraction of quasiparticle energies provides an efficient scheme for many applications. The performance of the implementation with respect to the basis set convergence and the effect of the use of pseudopotentials has been tested for the 117 closed-shell molecules from the G2/97 test set and 24 larger acceptor molecules from another recently proposed test set. Moreover, to demonstrate the potential of our method, we compute the thermally averaged GW density of states of a large photochromic compound by sampling ab initio molecular dynamics trajectories at different temperatures. The computed thermal line widths indicate approximately twice as large electron-phonon couplings with GW than with standard DFT-GGA calculations. This is further confirmed using frozen-phonon calculations.

The GW Compendium: A Practical Guide to Theoretical Photoemission Spectroscopy

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 measured 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, physicists 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.

Accurate description of the electronic structure of organic semiconductors by GW methods

Electronic properties associated with charged excitations, such as the ionization potential (IP), the electron affinity (EA), and the energy level alignment at interfaces, are critical parameters for the performance of organic electronic devices. To computationally design organic semiconductors and functional interfaces with tailored properties for target applications it is necessary to accurately predict these properties from first principles. Many-body perturbation theory is often used for this purpose within the GW approximation, where G is the one particle Green's function and W is the dynamically screened Coulomb interaction. Here, the formalism of GW methods at different levels of self-consistency is briefly introduced and some recent applications to organic semiconductors and interfaces are reviewed.

Vertex Corrections for Positive-Definite Spectral Functions of Simple Metals

We present a systematic study of vertex corrections in a homogeneous electron gas at metallic densities. The vertex diagrams are built using a recently proposed positive-definite diagrammatic expansion for the spectral function. The vertex function not only provides corrections to the well known plasmon and particle-hole scatterings, but also gives rise to new physical processes such as the generation of two plasmon excitations or the decay of the one-particle state into a two-particle-one-hole state. By an efficient Monte Carlo momentum integration we are able to show that the additional scattering channels are responsible for a reduction of the bandwidth, the appearance of a secondary plasmon satellite below the Fermi level, and a substantial redistribution of spectral weights. The feasibility of the approach for first-principles band-structure calculations is also discussed.

Insights and challenges of applying the GW method to transition metal oxides

The ab initio GW method is considered as the most accurate approach for calculating the band gaps of semiconductors and insulators. Yet its application to transition metal oxides (TMOs) has been hindered by the failure of traditional approximations developed for conventional semiconductors. In this work, we examine the effects of these approximations on the values of band gaps for ZnO, Cu2O, and TiO2. In particular, we explore the origin of the differences between the two widely used plasmon-pole models. Based on the comparison of our results with the experimental data and previously published calculations, we discuss which approximations are suitable for TMOs and why.

Valence band x-ray emission spectra of compressed germanium

We report measurements of the valence band width in compressed Ge determined from x-ray emission spectra below the Ge K edge. The width of the valence band does not show any pressure dependence in the semiconducting diamond-type structure of Ge below 10 GPa. On the other hand, in the metallic beta-Sn phase above 10 GPa the valence band width increases under compression. Density-functional calculations show an increasing valence band width under compression both in the semiconducting phase (contrary to experiment) and in the metallic beta-Sn phase of Ge (in agreement with observed pressure-induced broadening). The pressure-independent valence band width in the semiconducting phase of Ge appears to require theoretical advances beyond the density-functional theory or the GW approximation.

Many-body approximation scheme beyond GW

We explore the combination of the extended dynamical mean field theory (EDMFT) with the GW approximation (GWA); the former sums the local contributions to the self-energies to infinite order in closed form and the latter handles the nonlocal ones to lowest order. We investigate the different levels of self-consistency that can be implemented within this method by comparing to the exact quantum Monte Carlo solution of a finite-size model Hamiltonian. We find that using the EDMFT solution for the local self-energies as input to the GWA for the nonlocal self-energies gives the best result.

Inclusion of vertex corrections in the self-consistent calculation of quasiparticles in metals

Based on the Bethe-Salpeter equation and the Ward identity derived from it, we provide a scheme for constructing the vertex function in the self-consistent iteration loop to determine the electron self-energy. The scheme is implemented in the homogeneous electron gas at the sodium density.