Develops the many-body one-particle Green function, and explains its physical interpretation in terms of the spectral function, self-energy, and quasiparticle. lifetime. Its application in angle-resolved photoemission spectroscopy is presented in detail. The time-evolution operator in the interaction picture is derived, and time-ordering and adiabatic switching-on are introduced as precursor tools to construct the Feynman–Dyson many-body perturbation theory. A detailed account of Wick’s theorem, normal ordering, and contractions is outlined. Feynman diagrams are constructed, and the emergence of the infinite Dyson series from irreducible diagrams is outlined. Two-particle Green function and the particle-hole excitation spectra are developed. Diagramatic application of RPA for interacting systems is described. The finite-temperature Matsubara Green function is introduced and developed, together with its Fourier series representation and the evaluation of Matsubara sums.
Develops the many-body one-particle Green function, and explains its physical interpretation in terms of the spectral function, self-energy, and quasiparticle. lifetime. Its application in angle-resolved photoemission spectroscopy is presented in detail. The time-evolution operator in the interaction picture is derived, and time-ordering and adiabatic switching-on are introduced as precursor tools to construct the Feynman–Dyson many-body perturbation theory. A detailed account of Wick’s theorem, normal ordering, and contractions is outlined. Feynman diagrams are constructed, and the emergence of the infinite Dyson series from irreducible diagrams is outlined. Two-particle Green function and the particle-hole excitation spectra are developed. Diagramatic application of RPA for interacting systems is described. The finite-temperature Matsubara Green function is introduced and developed, together with its Fourier series representation and the evaluation of Matsubara sums.
We present results for many-body perturbation theory for the one-body Green's function at finite temperatures using the Matsubara formalism. Our method relies on the accurate representation of the single-particle states in standard Gaussian basis sets, allowing to efficiently compute, among other observables, quasiparticle energies and Dyson orbitals of atoms and molecules. In particular, we challenge the second-order treatment of the Coulomb interaction by benchmarking its accuracy for a well-established test set of small molecules, which includes also systems where the usual Hartree-Fock treatment encounters difficulties. We discuss different schemes how to extract quasiparticle properties and assess their range of applicability. With an accurate solution and compact representation, our method is an ideal starting point to study electron dynamics in time-resolved experiments by the propagation of the Kadanoff-Baym equations.
Abstract Green-function techniques are among the most versatile tools for investigating quantum many-body systems at finite temperature. Yet, standard textbook treatments often fragment the formalism—separating bosons and fermions, systems with massless and massive particles, or restricting proofs to special Hamiltonians—thereby obscuring the structural unity of the theory. This work develops a general, Hamiltonian-independent framework for both real- and imaginary-time Green-functions, valid for any combination of bosonic and fermionic degrees of freedom, including mixed systems such as interacting electron–phonon models. The key conceptual advance is the nature index η , which complements the conventional statistical sign ζ . While ζ encodes particle statistics, η distinguishes canonical quantization of systems with massless particles ( η = 0) from second quantization of systems with massive particles ( η = 1). This classification removes ambiguities in chemical-potential terms and unifies the treatment of photons, phonons, electrons, and their mixtures under a single notation. From this basis, we formulate generalized ensemble constructions, generalized Lehmann-representations, and a fluctuation-dissipation theorem valid for arbitrary Hamiltonians. Several results—such as the monotonicity of particle number in the chemical potential, an explicit formula for the quadratic contribution to the internal energy from the retarded Green-function, and the Hartree-contribution to the self-energy in electron–phonon systems—are either absent from, or only sparsely mentioned in, existing literature. We further introduce Green’s Chart, a visual and algebraic map that systematically connects all Green-function types through analytic continuation, Fourier-transforms, and Wick-rotations, enabling all to be generated from the Matsubara Green-function. Applications to both non-interacting and interacting models demonstrate the framework’s ability to compute densities of states, thermodynamic potentials, and spectral functions with a consistent formalism. The results provide a complete and reproducible toolkit for equilibrium many-body theory, consolidating scattered results into a unified structure suitable for both pedagogical and computational purposes.
Abstract Multiparticle thermodynamic Green’s functions, defined in terms of grand canonical ensemble averages of time-ordered products of creation and annihilation operators, are interpreted as tracing the amplitude for time-developing correlated interacting particle motions taking place in the background of a thermal ensemble. Under equilibrium conditions, time-translational invariance permits the one-particle thermal Green’s function to be represented in terms of a single frequency, leading to a Lehmann spectral representation whose frequency poles describe the energy spectrum. This Green’s function has finite values for both t>t′ and t<t′ (unlike retarded Green’s functions), and the two parts G1> and G1< (respectively) obey a simple proportionality relation that facilitates the introduction of a spectral weight function: It is also interpreted in terms of a periodicity/antiperiodicity property of a modified Green’s function in imaginary time capable of a Fourier series representation with imaginary (Matsubara) frequencies. The analytic continuation from imaginary time to real time is discussed, as are related commutator/anticommutator functions, also retarded/advanced Green’s functions, and the spectral weight sum rule is derived. Statistical thermodynamic information is shown to be embedded in physical features of the one- and two-particle thermodynamic Green’s functions.
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.
Recently we proposed a multichannel Dyson equation, in which we coupled the one- and the three-body Green's functions, to model photoemission spectra. In this work we extend this concept to higher-order Green's functions and to other spectroscopies. We show the general structure of the equations and how one can systematically approximate the corresponding multichannel self-energy. As a particular case, we focus on the coupling of the two-body and the four-body Green's functions in the electron-hole channel to describe neutral excitations. This formulation allows for the description of important many-body effects, such as biexcitons, in a natural way. We illustrate our approach by applying it to a two-level model system, which, in a one-particle picture, exhibits single and double excitations. Our method can correctly describe both kinds of excitation, unlike standard approaches, and in good agreement with the exact results.
We present the multichannel Dyson equation that combines two or more many-body Green's functions to describe the electronic structure of materials. In this work we use it to model photoemission spectra by coupling the one-body Green's function with the three-body Green's function. We demonstrate that, unlike methods using only the one-body Green's function, our approach puts the description of quasiparticles and satellites on an equal footing. We propose a multichannel self-energy that is static and only contains the bare Coulomb interaction, making frequency convolutions and self-consistency unnecessary. Despite its simplicity, we demonstrate with a diagrammatic analysis that the physics it describes is extremely rich. Finally, we present a framework based on an effective Hamiltonian that can be solved for any many-body system using standard numerical tools. We illustrate our approach by applying it to the Hubbard dimer and show that it is exact both at 1/4 and 1/2 filling.