In a many-body perturbation treatment of electronic excitations, one has to solve the so-called Bethe–Salpeter equation (BSE) for the kernel (vertex) describing the interaction between electrons. In general, the BSE exhibits a non-local and frequency-dependent screened interaction and consequently it is extremely difficult to solve. We have developed a scheme that solves the BSE equation iteratively, including dynamically screening, by considering the change in the electron self-energy upon applying a time-dependent field. The BSE is obtained from the self-energy in the GW approximation (GWA) using the Schwinger functional technique. We apply our formalism to a simple model system and discuss briefly changes in the charge response, when dynamical screening is taken into account.
In a many-body perturbation treatment of electronic excitations, one has to solve the so-called Bethe–Salpeter equation (BSE) for the kernel (vertex) describing the interaction between electrons. In general, the BSE exhibits a non-local and frequency-dependent screened interaction and consequently it is extremely difficult to solve. We have developed a scheme that solves the BSE equation iteratively, including dynamically screening, by considering the change in the electron self-energy upon applying a time-dependent field. The BSE is obtained from the self-energy in the GW approximation (GWA) using the Schwinger functional technique. We apply our formalism to a simple model system and discuss briefly changes in the charge response, when dynamical screening is taken into account.
The helium atom is the simplest many-body electronic system provided by nature. The exact solution to the Schrödinger equation is known for helium ground and excited states, and it represents a benchmark for any many-body methodology. Here, we check the ab initio many-body GW approximation and the Bethe-Salpeter equation (BSE) against the exact solution for helium. Starting from the Hartree-Fock method, we show that the GW and the BSE yield impressively accurate results on excitation energies and oscillator strength, systematically improving the time-dependent Hartree-Fock method. These findings suggest that the accuracy of the BSE and GW approximations is not significantly limited by self-interaction and self-screening problems even in this few electron limit. We further discuss our results in comparison to those obtained by time-dependent density-functional theory.
We present a real-time linear-response approach to solving the dynamical Bethe-Salpeter equation (BSE). The polarization part of the screened interaction is obtained from time-dependent Hartree propagation of the one-particle density matrix within an orbital basis-set representation. The frequency convolution defining the dynamical kernel is evaluated as a product in time, avoiding numerical integration and storage of the full screened Coulomb operator. The resulting nonlinear eigenvalue problem is then solved directly, going beyond the static screening approximation with full-frequency dependence in the screened interaction. We validate the method by calculating the low-lying singlet and triplet excitation energies of water, formaldehyde, and acetylene. While the deterministic scaling remains O(N6), the time propagation formulation is readily compatible with grid-based stochastic sampling methods, which will open the possibility for dynamical BSE calculations of very large systems.
In this thesis, an approximate method introduced to solve time-dependent many-body problems known as time-dependent many-body perturbation theory is studied. Many-body perturbation theory for interacting electrons and phonons is reviewed. In particular, the electron propagator G and an unconventional two-component phonon propagator, which satisfy coupled integral Dyson equations, are introduced. In practice, the associated integral kernels known as the electron Σ and phonon self-energies need to be approximated. The conserving approximations known as the Hartree (-Fock) and the first and second Born approximations, which respect the continuity equation between the electron density and current, are considered in this work. Time-dependent systems of interest are studied in this thesis by using the integro-differential forms of the Dyson equations referred to as the Kadanoff-Baym Equations (KBE). The Kadanoff-Baym equations are introduced for the electron and phonon propagators unconventionally as coupled firstorder integro-differential equations. It is reviewed how these equations are solved numerically by describing an integration rule, a class of practical methods and a parallel implementation of the numerical method. In addition, documentation of how the Kadanoff-Baym equations allow to solve the Bethe-Salpeter Equation (BSE) with the kernel δΣ/δG for the density response function, is provided. In two of the enclosed publications, we benchmarked observables obtained by using the Hartree and partially and fully self-consistent Born approximations against numerically exact results for the two- site, two-electron Holstein model. In this model, the two electrons which are constrained to move between two lattice sites interact indirectly via the electron-phonon coupling. It was observed that only the fully self-consistent Born approximation could cope qualitatively correctly with the competition between the delocalizing and localizing effects of the kinetic and interaction energies. For the other two approximations, spurious bifurcative symmetry breaking with an associated unbounded density response was observed. Nevertheless, also the self-consistent Born approximation was concluded to fail in describing the bipolaronic behavior of the true system. In the third publication, we benchmarked the Hartree-Fock and second Born approximations against an exact method for the few-site Hubbard and Pariser-Parr-Pople models in which the underlying lattice is inert and the electrons interact amongst themselves directly. It was found that the second Born approximation is capable of describing the so-called correlation induced doubly-excited states. This is not possible for time-local approximations such as Hartree-Fock. In addition to the qualitative results, which highlight successes of the applied simple self-energy approximations, the approximate and exact results were also compared on a more quantitative level. It is the quantitative and qualitative results combined which are used in this thesis to assess the quality of the many-body approximations used, with the aim to better understand when these methods are predictive.
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 .
The correlation energy of an inhomogeneous electron gas is ruled by the density correlation function L. It obeys a Bethe-Salpeter equation (BSE), whose kernel is the functional derivative of the particle self-energy. It is solved on equal footing for two-particle and single-particle Green's functions within Hedin's GW approximation. According to the Hartree and exchange-correlation contributions to the kernel, the correlation expression decays into random-phase-approximation (RPA) -like and screened-exchange-like parts. The actual treatment of the screened potential W in the BSE for L leads to different second-order screened exchange corrections (SOSEX). Including correlation in the kernel even yields approximations beyond SOSEX. The progress beyond RPA is illustrated for spinless and spin-polarized homogeneous electron gases, indicating that the screening treatment is substantial.
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.