Two-particle Green's functions and the vertex functions play a critical role in theoretical frameworks for describing strongly correlated electron systems. However, numerical calculations at the two-particle level often suffer from large computation time and massive memory consumption. We derive a general expansion formula for the two-particle Green's functions in terms of an overcomplete representation based on the recently proposed ``intermediate representation'' basis. The expansion formula is obtained by decomposing the spectral representation of the two-particle Green's function. We demonstrate that the expansion coefficients decay exponentially, while all high-frequency and long-tail structures in the Matsubara-frequency domain are retained. This representation therefore enables efficient treatment of two-particle quantities and opens a route to the application of modern many-body theories to realistic strongly correlated electron systems.
Two-particle Green's functions and the vertex functions play a critical role in theoretical frameworks for describing strongly correlated electron systems. However, numerical calculations at the two-particle level often suffer from large computation time and massive memory consumption. We derive a general expansion formula for the two-particle Green's functions in terms of an overcomplete representation based on the recently proposed ``intermediate representation'' basis. The expansion formula is obtained by decomposing the spectral representation of the two-particle Green's function. We demonstrate that the expansion coefficients decay exponentially, while all high-frequency and long-tail structures in the Matsubara-frequency domain are retained. This representation therefore enables efficient treatment of two-particle quantities and opens a route to the application of modern many-body theories to realistic strongly correlated electron systems.
Two-particle response functions are a centerpiece of both experimental and theoretical quantum many-body physics. Yet, due to their size and discontinuity structure, they are challenging to handle numerically. Recently, two advances were made to tackle this problem: first, the overcomplete intermediate representation (OIR), which provides a highly efficient compression of Green's functions in imaginary frequency, and second, partial spectral functions (PSFs), which allow for an efficient evaluation in real frequency. We show that there is a two-to-one correspondence between PSFs and OIR coefficients and exploit this fact to construct the OIR for three-or-more-particle propagators. We then use OIR to fit and compress imaginary-frequency data obtained from the numerical renormalization group (NRG), reaching a compression ratio of more than 400. Finally, we attempt to match the OIR data to partial Green's functions from NRG. Due to the overcompleteness, we achieve only qualitative agreement. Published by the American Physical Society 2024
Many-body calculations at the two-particle level require a compact representation of two-particle Green’s functions. In this paper, we introduce a sparse sampling scheme in the Matsubara frequency domain as well as a tensor network representation for two-particle Green’s functions. The sparse sampling is based on the intermediate representation basis and allows an accurate extraction of the generalized susceptibility from a reduced set of Matsubara frequencies. The tensor network representation provides a system independent way to compress the information carried by two-particle Green’s functions. We demonstrate efficiency of the present scheme for calculations of static and dynamic susceptibilities in single- and two-band Hubbard models in the framework of dynamical mean-field theory.
The imaginary-time Green's function is a building block of various numerical methods for correlated electron systems. Recently, it was shown that a model-independent compact orthogonal representation of the Green's function can be constructed by decomposing its spectral representation. We investigate the performance of this so-called intermediate representation (IR) from several points of view. First, we develop an efficient algorithm for computing the IR basis functions of arbitrary high degree. Second, for two simple models, we study the number of coefficients required to represent the Green's function within a given tolerance. We show that the number of coefficients grows only as $O(ln\ensuremath{\beta})$ for fermions and converges to a constant for bosons as temperature $T=1/\ensuremath{\beta}$ decreases. Third, we show that this remarkable feature is ascribed to the properties of the physically constructed basis functions. The fermionic basis functions on the real-frequency axis have features whose width is scaled as O(T), which is consistent with the low-T properties of quasiparticles in a Fermi liquid state. On the other hand, the properties of the bosonic basis functions are consistent with those of spin/orbital susceptibilities at low T. These results demonstrate the potential wide applications of the IR to calculations of correlated systems.
We study the expansion of single-particle and two-particle imaginary-time Matsubara Green's functions of quantum impurity models in the basis of Legendre orthogonal polynomials. We discuss various applications within the dynamical mean-field theory (DMFT) framework. The method provides a more compact representation of the Green's functions than standard Matsubara frequencies, and therefore significantly reduces the memory-storage size of these quantities. Moreover, it can be used as an efficient noise filter for various physical quantities within the continuous-time quantum Monte Carlo impurity solvers recently developed for DMFT and its extensions. In particular, we show how to use it for the computation of energies in the context of realistic DMFT calculations in combination with the local density approximation to the density functional theory (LDA+DMFT) and for the calculation of lattice susceptibilities from the local irreducible vertex function.
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.
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.