The many-body problem cannot in general be solved exactly, and one of the most prominent approximations is to build perturbation expansions. A huge variety of expansions are possible, which differ by the quantity to be expanded, the expansion variable, the starting point, and ideas of how to re-sum or terminate the series. Although much has been discussed and much has been done, some choices were made for historical reasons, in particular, limited computation or storage capacities. The present work aims to examine the justification for different choices made in different contexts, by comparing the components of functionals based on Green's functions on the one side, and on the charge density on the other side. Of particular interest will be the question of how to build an optimal starting point for the approximation of non-local quantities, making use of near- or far-sightedness, and daring to consider models beyond the homogeneous electron gas. This will include the use of connector approximations. We will also discuss why it is a good idea to build functionals of the density.
The many-body problem cannot in general be solved exactly, and one of the most prominent approximations is to build perturbation expansions. A huge variety of expansions are possible, which differ by the quantity to be expanded, the expansion variable, the starting point, and ideas of how to re-sum or terminate the series. Although much has been discussed and much has been done, some choices were made for historical reasons, in particular, limited computation or storage capacities. The present work aims to examine the justification for different choices made in different contexts, by comparing the components of functionals based on Green's functions on the one side, and on the charge density on the other side. Of particular interest will be the question of how to build an optimal starting point for the approximation of non-local quantities, making use of near- or far-sightedness, and daring to consider models beyond the homogeneous electron gas. This will include the use of connector approximations. We will also discuss why it is a good idea to build functionals of the density.
We especially need imagination in science. aria Mitchell Summary The many-body problem consists of two parts: the first is the non-interacting system in a materials-specific external potential; the second is the Coulomb interaction that makes the problem so hard to solve. The most straightforward idea is to use perturbation theory, with the Coulomb interaction as perturbation. This is conceptually simple, but it turns out to be difficult in practice, since the Coulomb interaction is often not small compared with typical energy differences, it is long-ranged and in the thermodynamic limit there is an infinite number of particles, contributing with an infinite number of mutual interaction processes. The present chapter outlines how one can deal with this problem. It contains an overview of facts that one can also find in many standard textbooks on the many-body problem, but that are useful to keep in mind in order to look at later chapters from a sound and well-established perspective. The many-body problem is a tough one, and it has many facets. Sorting it out is like putting together a huge puzzle. The eight introductory chapters of this book provide pieces of the puzzle, and ideas on what one might do about it. In the present chapter we choose to go in one of the possible directions, in order to arrive at something tangible. The chapter gives the general framework and the main ideas; specific approximations are the topic of Chs. 10–15. The idea is to start from an independent-particle problem and add the Coulomb interaction as a perturbation. This is not easy: first, the interaction is responsible for a rich variety of phenomena that are completely absent otherwise, such as the finite lifetime of quasiparticles, or additional structures in spectra due to the fact that a quasi-particle excitation may transfer its energy to other elementary excitations, for example plasmons. Second, because of the two-body Coulomb interaction, the problem scales badly with the number of electrons, and straightforward perturbation theory for the many-body hamiltonian with the Coulomb interaction as perturbation rapidly becomes intractable or even useless, especially in large systems. To get started, Sec. 9.1 recalls why things are not so easy. The following sections try to solve one problem after the other, starting from Sec. 9.2 where the Green's function is reformulated in a way that is appropriate for a perturbation expansion.
To avoid confusion, it is necessary to be quite explicit about what is assumed and what is to be proved. A. Klein, Phys. Rev. 121, 950 (1961) Summary The topic of this chapter is functionals that provide a concise formulation for thermodynamic quantities and Green's function methods in interacting many-body systems. Building upon the instructive examples of density functional theory and the Hartree– Fock approximation, functionals of the Green's function are developed to provide a framework for practical methods used extensively in many-body perturbation theory (Chs. 9–15), dynamical mean-field theory (Chs. 16–20), and their combination in Ch. 21. Functional concepts give a firm foundation for those methods, for example, they give the conditions for approximations to obey conservation laws. Additional background and specific aspects are in App. H. Many physical quantities can be expressed as functionals , i.e., they depend on an entire function.1 For example, in the quantum variational method to determine a wavefunction, the total energy of a system of many electrons is a functional of a many-body trial wavefunction ( r 1 , r 2 , …): it is the expectation value of the hamiltonian. At the variational minimum, assuming the trial wavefunction is completely flexible, one finds the many-body Schrödinger equation for the ground state. Many of the properties we are interested in are, formally, simple functionals of the ground- or excited-state wavefunctions. In the search for feasible ways to calculate properties such as excitation spectra, we are led to different types of functional. The idea is to deal only with a physical, measurable quantity Q of interest and a variational functional E [ Q ] of this quantity. The variational character of E allows one to derive equations that determine Q 0 , the value of Q for the system in a desired state, typically the ground state or equilibrium state at T ≠ 0. Other properties should be different functionals F m [ Q ] that can then be evaluated at Q = Q 0 . This strategy is outlined in Fig. 8.1. To appreciate the power of such approaches, one can consider the impact of density functional theory, which is formulated as a functional of the density n ( r ) of electrons; it is defined for a range of densities and the variational minimum for the energy functional E [ n ] leads to equations to determine the actual ground-state density n 0 ( r ) and total energy E 0 .
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.
Many-body perturbation theory is often formulated in terms of an expansion in the dressed instead of the bare Green's function, and in the screened instead of the bare Coulomb interaction. However, screening can be calculated on different levels of approximation, and it is important to define what is the most appropriate choice. We explore this question by studying a zero-dimensional model (so called 'one-point model') that retains the structure of the full equations. We study both linear and non-linear response approximations to the screening. We find that an expansion in terms of the screening in the random phase approximation is the most promising way for an application in real systems. Moreover, by making use of the nonperturbative features of the Kadanoff-Baym equation for the one-body Green's function, we obtain an approximate solution in our model that is very promising, although its applicability to real systems has still to be explored.
Several widely used methods for the calculation of band structures and photo emission spectra, such as the GW approximation, rely on many-body perturbation theory. They can be obtained by iterating a set of functional differential equations (DEs) relating the one-particle Green's function (GF) to its functional derivative with respect to an external perturbing potential. In this work, we apply a linear response expansion in order to obtain insights into various approximations for GF calculations. The expansion leads to an effective screening while keeping the effects of the interaction to all orders. In order to study various aspects of the resulting equations, we discretize them and retain only one point in space, spin and time for all variables. Within this one-point model we obtain an explicit solution for the GF, which allows us to explore the structure of the general family of solutions and to determine the specific solution that corresponds to the physical one. Moreover, we analyze the performances of established approaches like GW over the whole range of interaction strength, and we explore alternative approximations. Finally, we link certain approximations for the exact solution to the corresponding manipulations of the DE which produces them. This link is crucial in view of a generalization of our findings to the real (multidimensional functional) case where only the DE is known.
In this paper, four distinct ideas are combined, which under a wide range of circumstances can give very rapidly converging series expansions for functional integrals. (1) Expansion of the functional being integrated in functional Taylor series. In the familiar case arising in quantum statistical mechanics, that of the Wiener integral of exp[−∫V(x(t))dt],V(x) being the perturbing potential, this is equivalent to expanding the characteristic functional of the probability functional of V(x(t)) in central moments of V(x(t)). The lowest-order term of the series is the approximation obtained by Feynman and Hibbs through a variational method. (2) Transfer of the harmonic term of the potential, when the functional integral is the quantum-statistical density matrix (Green's function of the Bloch equation), to the weighting function. This transforms the functional integral from a Wiener to an Uhlenbeck-Ornstein integral. The formal expressions for the terms of the expansion are somewhat more complicated, but they can be worked out, and the result is a great improvement in the speed of convergence of the series with decreasing temperature and/or decreasing relative magnitude of the anharmonic part of the potential. (3) ``Reservation of variables'' in the integration. This amounts to breaking the averaging process down into an average over subsets of the distributed random function (conditional average), followed by an averaging of these averages. Any step of this kind (it may be repeated within the subsets, etc.) gives an improvement of accuracy. (4) When the quantity being evaluated through functional integration is the partition function, the device introduced by Feynman and Hibbs, of interchanging the functional integration with the integration of the Green's function over the equated initial and final configuration-space points, may be combined with the above techniques. This eliminates one integration in the terms of the expansion and seems to improve accuracy at the same time. The general series obtained is correlated with more conventional operator techniques of quantum-mechanical perturbation theory, in order to answer the perennial question, does the path-integral method bring with it anything that could not be derived by other methods? It is in some sense a Feynman-Dyson expansion of the Green's function, but one that is further modified mathematically in a way characteristic only of the path-integral point of view, and which, moreover, improves its accuracy. It thus appears unlikely that the result is merely one of standard type disguised as a functional integral result. Sample numerical calculations are given to assay the accuracy of the methods, which are shown to compare very favorably with the traditional approximation of finite subdivision of the time interval.