Within the framework of many-particle perturbation theory, we develop an analytical approach that allows us to determine the small distance behavior of Green's functions and related quantities in electronic structure theory. As a case study, we consider the one-particle Green's function up to 2nd order in the perturbation approach. We derive explicit expressions for the leading order terms of the asymptotic small distance behavior. In particular, we demonstrate the appearance of a logarithmic term in the corresponding 2nd order Feynman diagrams. Our asymptotic analysis leads to an improved classification scheme for the diagrams, which takes into account not only the perturbation order, but also the asymptotic smoothness properties near their diagonals. Such a classification may be useful in the design of numerical algorithms and helps to improve their efficiency.
Within the framework of many-particle perturbation theory, we develop an analytical approach that allows us to determine the small distance behavior of Green's functions and related quantities in electronic structure theory. As a case study, we consider the one-particle Green's function up to 2nd order in the perturbation approach. We derive explicit expressions for the leading order terms of the asymptotic small distance behavior. In particular, we demonstrate the appearance of a logarithmic term in the corresponding 2nd order Feynman diagrams. Our asymptotic analysis leads to an improved classification scheme for the diagrams, which takes into account not only the perturbation order, but also the asymptotic smoothness properties near their diagonals. Such a classification may be useful in the design of numerical algorithms and helps to improve their efficiency.
Abstract Within the framework of many-particle perturbation theory, we develop an analytical approach that allows us to determine the small distance behaviour of Green’s functions and related quantities in electronic structure theory. As a case study, we consider the one-particle Green’s function up to second-order in the perturbation approach. We derive explicit expressions for the leading order terms of the asymptotic small distance behaviour. In particular, we demonstrate the appearance of logarithmic terms in the corresponding second-order Feynman diagrams. Our asymptotic analysis leads to an improved classification scheme for the diagrams, which takes into account not only the perturbation order, but also the asymptotic smoothness properties near their diagonals. Such a classification may be useful in the design of numerical algorithms and helps to improve their efficiency.
A thorough analytical and numerical characterization of the whole perturbation series of one-particle many-body Green's function (MBGF) theory is presented in a pedagogical manner. Three distinct but equivalent algebraic (first-quantized) recursive definitions of the perturbation series of the Green's function are derived, which can be combined with the well-known recursion for the self-energy. Six general-order algorithms of MBGF are developed, each implementing one of the three recursions, the ΔMPn method (where n is the perturbation order) [S. Hirata et al., J. Chem. Theory Comput. 11, 1595 (2015)], the automatic generation and interpretation of diagrams, or the numerical differentiation of the exact Green's function with a perturbation-scaled Hamiltonian. They all display the identical, nondivergent perturbation series except ΔMPn, which agrees with MBGF in the diagonal and frequency-independent approximations at 1≤n≤3 but converges at the full-configuration-interaction (FCI) limit at n=∞ (unless it diverges). Numerical data of the perturbation series are presented for Koopmans and non-Koopmans states to quantify the rate of convergence towards the FCI limit and the impact of the diagonal, frequency-independent, or ΔMPn approximation. The diagrammatic linkedness and thus size-consistency of the one-particle Green's function and self-energy are demonstrated at any perturbation order on the basis of the algebraic recursions in an entirely time-independent (frequency-domain) framework. The trimming of external lines in a one-particle Green's function to expose a self-energy diagram and the removal of reducible diagrams are also justified mathematically using the factorization theorem of Frantz and Mills. Equivalence of ΔMPn and MBGF in the diagonal and frequency-independent approximations at 1≤n≤3 is algebraically proven, also ascribing the differences at n = 4 to the so-called semi-reducible and linked-disconnected diagrams.
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.
Using a general-order ab initio many-body Green's function method, we numerically illustrate several pathological behaviors of the Feynman-Dyson diagrammatic perturbation expansion of one-particle many-body Green's functions as electron Feynman propagators. (i) The perturbation expansion of the frequency-dependent self-energy is not convergent at the exact self-energy in many frequency domains. (ii) An odd-perturbation-order self-energy has a qualitatively wrong shape and, as a result, many roots of the corresponding Dyson equation are nonphysical in that the poles may be complex or residues can exceed unity or be negative. (iii) A higher even-order self-energy consists of vertical lines at many frequencies, predicting numerous phantom poles with zero residues. (iv) Infinite partial resummations of diagrams by vertex or edge renormalization tend to exacerbate these pathologies. (v) The nonconvergence is caused by the nonanalyticity of the rational-function form of the exact Green's function at many frequencies, where the radius of convergence of its Taylor expansion is zero. This is consistent with the fact that (vi) Pad\'e approximants (power-series expansions of a rational function) can largely restore the correct shape and poles of the Green's function. Nevertheless, not only does the nonconvergence render higher-order Feynman-Dyson diagrammatic perturbation theory useless for many lower-lying ionization or higher-lying electron-attachment states but it also calls into question the validity of its combined use with the Ans\"atze requiring the knowledge of all poles and residues. Such Ans\"atze include the Galitskii-Migdal identity, the self-consistent Green's function methods, and some models of the algebraic diagrammatic construction.
This is the first of a series of papers dealing with many-particle systems from a unified, nonperturbative point of view. It contains derivations and discussions of various field-theoretical techniques which will be applied in subsequent papers. In a short introduction the general method of approach is summarized, and its relationship to other field-theoretic problems indicated. In the second section the macroscopic properties of the spectra of many-particle systems are described. Asymptotic evaluations are performed which characterize these macroscopic features in terms of intensive parameters, and the relationship of these parameters to thermodynamics is discussed. The special characteristics of the ground state are shown to follow as a limiting case of the asymptotic evaluations. The third section is devoted to the time-dependent field correlation functions, or Green's functions, which describe the microscopic behavior of a multiparticle system. These functions are defined, and related to intensive macroscopic variables when the energy and number of particles are large. Spectral representations and other properties of various one-particle Green's functions are derived. In the fourth section the treatment of non-equilibrium processes is considered. As a particular example, the electromagnetic properties of a system are expressed in terms of the special two-particle Green's function which describes current correlation. The discussion yields specifically a fluctuation-dissipation theorem, a sum rule for conductivity, and certain dispersion relations. The fifth section deals with the differential equations which determine the Green's functions. The boundary conditions that characterize the Green's function equations are exhibited without reference to adiabatic decoupling. A method for solving the equations approximately, by treating the correlations among successively larger numbers of particles, is considered. The first approximation in this sequence is shown to yield a generalized Hartree-like equation. A related, but rigorous, identity for the single-particle Green's function is then derived. A second approximation, which takes certain two-particle correlations into account, is shown to produce various additional effects: The interaction between particles is altered in a manner characterized by the intensive macroscopic parameters, and the modification and spread of the energy-momentum relation come into play. In the final section compact formal expressions for the Green's functions and other physical quantities are derived. Alternative equations and systematic approximations for the Green's functions are obtained.
This dissertation concerns the quantum many-body problem, which is the problem of predicting the properties of systems of several quantum particles from the first principles of quantum mechanics. Included under this umbrella are various problems of fundamental importance in quantum chemistry, condensed matter physics, and materials science. Of particular interest is the electronic structure problem, the problem of determining the state of the electrons in a system with fixed atomic nuclei. Since direct numerical solution of the many-body Schrödinger equation is intractable even for systems of moderate size, a diverse array of approximate methods has been developed. The broad goals of this dissertation are to improve the mathematical understanding of certain widely-used approximations, as well as to propose new methods. Roughly speaking, we consider three (overlapping) categories of methods: Green's function methods, embedding methods, and variational methods.One can understand Green's function methods in terms of many-body perturbation theory, which computes series expansions of physical quantities about a non-interacting reference system. These expansions can be expressed graphically in terms of Feynman diagrams, which can in turn be reorganized, in some cases, into an expansion in terms of so-called bold diagrams. Green's function methods can be specified by choosing a subset of bold diagrams to approximate the sum. At the same time, such methods can be understood in terms of an object known as the Luttinger-Ward (LW) functional, which admits a representation in terms of the bold diagrams. Many aspects of these constructions are purely formal, and indeed the existence of the fermionic LW functional as a single-valued functional has recently been called into question. To contribute to the understanding of these issues, we provide rigorous proofs of the combinatorial construction and analytic interpretation of the bold diagrams in the simplified setting of a classical field theory. In this setting we also provide a rigorous non-perturbative construction of the LW functional via convex duality and prove several key properties, including continuity up to the boundary of its domain and asymptotics in the limit of large interaction.Quantum embedding methods, meanwhile, view a large system as being composed of smaller fragments that are treated with high accuracy and embedded in the larger system in a mutually consistent way. Inspired by a connection between the boundary analysis of the LW functional and embedding, we perform similar analysis for the 1-RDM theory for fermionic systems, which is also developed via convex duality, illustrating a relation to fermionic embedding methods such as the density matrix embedding theory (DMET).Another embedding method of note is the dynamical mean-field theory (DMFT), which is at the same time a Green's function method that can be understood in terms of the LW functional. DMFT relies on the solution of impurity problems, which specify the embedding of an interacting system into a non-interacting bath. Underlying DMFT is a result about the sparsity pattern of the self-energy matrix for impurity problems, which to our knowledge has not been proved in the literature. We provide a rigorous proof of this result in various classical and quantum settings. We go on to investigate the fermionic DMFT in depth, identifying the key mathematical structures that appear in the algorithmic loop for solving it and using these to prove the well-posedness of this loop, in a certain sense.Finally, we introduce a suite of new approaches to the quantum many-body problem that provide variational lower bounds to the ground-state energy. These methods, which combine the themes of convexity and embedding, are based on novel convex relaxations of the variational principles for the ground-state energies of many-body systems. To begin, we recover a second-quantized version of the formalism of strictly correlated electrons (SCE), which yields an exact expression for the exchange-correlation functional in Kohn-Sham density functional theory in the limit of infinite Coulomb repulsion in terms of the solution of a multi-marginal optimal transport problem. We introduce a semidefinite relaxation method for approximately solving this problem and obtaining a lower bound for the ground-state energy. The ideas underlying this relaxation are generalized considerably, outside the context of SCE, to yield much tighter lower bounds, which we validate numerically for both quantum spin systems and fermionic systems. We also describe how these relaxation methods can be interpreted as embedding methods via convex duality.