The many-body problem is usually approached from one of two perspectives: the first originates from an action and is based on Feynman diagrams, the second is centered around a Hamiltonian and deals with quantum states and operators.The connection between results obtained in either way is made through spectral (or Lehmann) representations, well known for two-point correlation functions.Here, we complete this picture by deriving generalized spectral representations for multipoint correlation functions that apply in all of the commonly used many-body frameworks: the imaginary-frequency Matsubara and the realfrequency zero-temperature and Keldysh formalisms.Our approach separates spectral from time-ordering properties and thereby elucidates the relation between the three formalisms.The spectral representations of multipoint correlation functions consist of partial spectral functions and convolution kernels.The former are formalism independent but system specific; the latter are system independent but formalism specific.Using a numerical renormalization group method described in the accompanying paper, we present numerical results for selected quantum impurity models.We focus on the four-point vertex (effective interaction) obtained for the single-impurity Anderson model and for the dynamical mean-field theory solution of the one-band Hubbard model.In the Matsubara formalism, we analyze the evolution of the vertex down to very low temperatures and describe the crossover from strongly interacting particles to weakly interacting quasiparticles.In the Keldysh formalism, we first benchmark our results at weak and infinitely strong interaction and then reveal the rich real-frequency structure of the dynamical mean-field theory vertex in the coexistence regime of a metallic and insulating solution.
The many-body problem is usually approached from one of two perspectives: the first originates from an action and is based on Feynman diagrams, the second is centered around a Hamiltonian and deals with quantum states and operators.The connection between results obtained in either way is made through spectral (or Lehmann) representations, well known for two-point correlation functions.Here, we complete this picture by deriving generalized spectral representations for multipoint correlation functions that apply in all of the commonly used many-body frameworks: the imaginary-frequency Matsubara and the realfrequency zero-temperature and Keldysh formalisms.Our approach separates spectral from time-ordering properties and thereby elucidates the relation between the three formalisms.The spectral representations of multipoint correlation functions consist of partial spectral functions and convolution kernels.The former are formalism independent but system specific; the latter are system independent but formalism specific.Using a numerical renormalization group method described in the accompanying paper, we present numerical results for selected quantum impurity models.We focus on the four-point vertex (effective interaction) obtained for the single-impurity Anderson model and for the dynamical mean-field theory solution of the one-band Hubbard model.In the Matsubara formalism, we analyze the evolution of the vertex down to very low temperatures and describe the crossover from strongly interacting particles to weakly interacting quasiparticles.In the Keldysh formalism, we first benchmark our results at weak and infinitely strong interaction and then reveal the rich real-frequency structure of the dynamical mean-field theory vertex in the coexistence regime of a metallic and insulating solution.
Abstract Conceptually, the Matsubara formalism (MF), using imaginary frequencies, and the Keldysh formalism (KF), formulated in real frequencies, give equivalent results for systems in thermal equilibrium. The MF has less complexity and is thus more convenient than the KF. However, computing dynamical observables in the MF requires the analytic continuation from imaginary to real frequencies. The analytic continuation is well‐known for two‐point correlation functions (having one frequency argument), but, for multipoint correlators, a straightforward recipe for deducing all Keldysh components from the MF correlator had not been formulated yet. Recently, a representation of MF and KF correlators in terms of formalism‐independent partial spectral functions and formalism‐specific kernels was introduced by Kugler, Lee, and von Delft [Phys. Rev. X 11, 041006 (2021)]. This representation is used to formally elucidate the connection between both formalisms. How a multipoint MF correlator can be analytically continued to recover all partial spectral functions and yield all Keldysh components of its KF counterpart is shown. The procedure is illustrated for various correlators of the Hubbard atom.
A major challenge in the field of correlated electrons is the computation of dynamical correlation functions. For comparisons with experiment, one is interested in their real-frequency dependence. This is difficult to compute because imaginary-frequency data from the Matsubara formalism require analytic continuation, a numerically ill-posed problem. Here, we apply quantum field theory to the single-impurity Anderson model using the Keldysh instead of the Matsubara formalism with direct access to the self-energy and dynamical susceptibilities on the real-frequency axis. We present results from the functional renormalization group (fRG) at the one-loop level and from solving the self-consistent parquet equations in the parquet approximation. In contrast with previous Keldysh fRG works, we employ a parametrization of the four-point vertex which captures its full dependence on three real-frequency arguments. We compare our results to benchmark data obtained with the numerical renormalization group and to second-order perturbation theory. We find that capturing the full frequency dependence of the four-point vertex significantly improves the fRG results compared with previous implementations, and that solving the parquet equations yields the best agreement with the numerical renormalization group benchmark data but is only feasible up to moderate interaction strengths. Our methodical advances pave the way for treating more complicated models in the future.
We develop a self-consistent spectral quadrature (sc-SQ) framework for the calculation of many-body Green functions. The method approximates the Källén--Lehmann spectral measure by Gauss--Christoffel (GC) quadrature, yielding a rational Green function representation with guaranteed spectral positivity that exactly reproduces the first 2N spectral moments at pole order N. A key component is an SVD-based rank-selection criterion on the Hankel matrix, which identifies the numerically resolvable pole rank N∗ from the singular-value gap and acts as a precision-guided diagnostic of correlation complexity. The scheme is made self-consistent by requiring that the spectral function used to evaluate expectation values coincides with the spectral function generated by the quadrature reconstruction. This defines a fixed-point hierarchy that connects systematically to established approximations, including Hartree--Fock and Hubbard-I, and incorporates non-perturbative features such as multi-peak spectral structure. We benchmark the approach for the Anderson impurity model against numerical renormalization group (NRG) results and apply it within dynamical mean-field theory for the Hubbard model on the Bethe lattice. The method captures the three-peak Anderson impurity spectrum and the suppression of quasiparticle weight in the half-filled Hubbard model on the Bethe lattice, including Mott-gap formation on the insulating branch for N⩾5, in qualitative agreement with NRG references.
We introduce an alternative approach for obtaining smooth finite-temperature spectral functions of quantum impurity models using the numerical renormalization group (NRG) technique. It is based on calculating first the Green's function on the imaginary-frequency axis, followed by an analytic continuation to the real-frequency axis using Pad\'e approximants. The arbitrariness in choosing a suitable kernel in the conventional broadening approach is thereby removed and, furthermore, we find that the Pad\'e method is able to resolve fine details in spectral functions with less artifacts on the scale of $\ensuremath{\omega}\ensuremath{\sim}T$. We discuss the convergence properties with respect to the NRG calculation parameters (discretization $\ensuremath{\Lambda}$, z-averaging, truncation cutoff) and the number of Matsubara points taken into account in the analytic continuation. We test the technique on the single-impurity Anderson model and the Hubbard model (within the dynamical mean-field theory). For the Anderson impurity model, we discuss the shape of the Kondo resonance and its temperature dependence. For the Hubbard model, we discuss the inner structure of the Hubbard bands in metallic and insulating solutions at half-filling, as well as in the doped Mott insulator. Based on these test cases, we conclude that the Pad\'e approximant approach provides improved results for spectral functions at low-frequency scales of $\ensuremath{\omega}\ensuremath{\lesssim}T$ and that it is capable of resolving sharp spectral features also at high frequencies.
An essential role for the description of correlated quantum many-particle systems is played by impurity problems. These consist of a small number of localized orbitals, where the electrons are subject to a Coulomb-interaction, coupled to noninteracting bath degrees of freedom. Impurity models are used for the description of magnetic impurities in metals, of quantum points in nanostructures, and in the context of the dynamical mean field theory (DMFT). For their solution a large number of numerical methods exists. We mention here quantum Monte Carlo methods, the numerical renormalization group, or the exact diagonalization method. In this thesis we introduce a new approach to impurity problems, which is based on the functional renormalization group (fRG). In contrast to the methods mentioned above this approach is approximative, but it can be applied with comparatively lower numerical effort. The aim of this thesis is to examine the advantages and limits of this method. An important role in this connection is played by the single impurity Anderson model (SIAM), which is already understood very well with other methods and is seen as a prototype for models with strong local correlations. The spin fluctuations of this model are governed by the so-called Kondo energy scale and a central question of this thesis is whether this scale can be reproduced with our approach. We begin with a short introduction to the Green's function formalism for correlated quantum many particle systems. Then the fRG flow equations for the one-particle irreducible vertex functions are deduced, and several approximation schemes are discussed. In chapter 4 we introduce the fRG scheme used in this thesis. At first the SIAM is mapped to a semi-infinite chain, in which the interacting orbital is given by the first site of this chain. The system is then subdivided into two parts. The first part, which is called 'core', consists of the interacting orbital and the first L bath sites. The remaining bath sites form the second part. At the beginning both parts are decoupled from each other such that the core can be solved exactly. Starting with this exact solution the coupling between the core and the remaining bath sites is switched on slowly in a renormalization group flow. This way it is possible to calculate the one-particle and two-particle correlation functions on the correlated site. We call this flow scheme 'hybridization flow'. The flow equations are formulated in an effective theory on the first bath site outside the core, which turns out to be advantageous compared to other implementations of the hybridization flow. In chapter 5 this hybridization flow method is applied to the SIAM with semi-elliptic bath density of states. We discuss the differences which arise for different core-sizes (L = 0,1,2,3) and for different truncations of the fRG flow equations. In the cases L=0 and L=2 the local Fermi liquid properties of the SIAM and the associated Kondo scale can not be reproduced in both approximation schemes. For the core sizes L=1 and L=3 the ground state of the core is already a spin singlet state in the beginning of the flow and the Fermi liquid properties are reproduced. The Kondo energy scale, which determines the width of the central quasi-particle peak and the size of the local spin-susceptibility, is however not accurately resolved. This is partly due to the finite temperature used in our implementation such that for larger interaction strength the temperature is above the Kondo energy scale. The application of the hybridization flow method to quantum impurity problems in the context of the DMFT is discussed in chapter 6. In this approach lattice models with a local interaction, like the Hubbard model, are mapped to an effective impurity problem for a single lattice site embedded in a dynamical mean field representing the influence of the other electrons. This scheme is exact in infinitely many dimensions and a nonperturbative approximation method for finite-dimensional systems. In particular the Mott metal-insulator transition in the Hubbard model is successfully described by this method. One class of quantities, which is easily accessible from the hybridization flow scheme, is given by the local two-particle correlation functions. These play an essential role in non-local extensions of the DMFT. Therefore we mainly focus on these quantities. At first we apply our flow scheme to the Hubbard model on the Bethe lattice in infinite dimensions and calculate the local one-particle irreducible vertex functions in the insulating and the metallic phase in good agreement with previous calculations, which use the exact diagonalization method as impurity solver. After this we extend our flow scheme to a cluster DMFT method, which includes short-range antiferromagnetic correlations and calculate again the local and non-local vertex functions for the two-dimensional Hubbard model. We compare these to the results of the single-site DMFT.
In the first part of the Thesis we mostly concentrate on spectral properties of strongly correlated systems and on their equilibrium properties. This is accomplished by the general concept of imaginary-time dynamics which we apply to a number of different problems in which different strengths of this approach emerge. In Chapter 1 we introduce the formalism that allows for a connection between the quantum and the classical worlds. The connection is established by means of the imaginary-time quantum evolution which, under certain circumstances, is shown to be equivalent to a classical stochastic process. It is further shown that exact static and spectral properties of correlated systems can be obtained when this mapping is feasible. The relationship between the imaginary-time dynamics in different frameworks such as the path-integral and the perturbative one is also underlined. In Chapter 2 we present a specific implementation of the general ideas previously presented. In particular we introduced an extension to lattice systems of the Reptation Monte Carlo algorithm [30] which benefits of a sampling scheme based on directed updates. Specific improvements over the existing methodologies consist in the unbiased evaluation of the imaginary-time path integrals for bosons and a systematic scheme to improve over the Fixed-node approximation for fermions. Applications to the Hubbard and the Heisenberg models are presented. In Chapter 3 we demonstrate the application of the imaginary-time dynamics to the exact study of spectral properties. Subject of our attention is a highly anharmonic and correlated quantum crystal such as Helium 4 at zero temperature.[33] Concerning this system, we have obtained the first ab-initio complete phonon dispersion in good agreement with neutron spectroscopy experiments. Moreover, we have also studied the density excitations of solid helium in a region of wave-vectors in between the collective (phonon) and the single-particle regimes, where the presence of residual coherence in the dynamics shows analogies between the highly anharmonic crystal and the superfluid phase. In Chapter 4 we introduce a novel method, based on the imaginary-time dynamics, to obtain unbiased estimates of fermionic properties.[34] By means of this method and of a very accurate variational state, we provide strong evidence for the stability of a saturated ferromagnetic phase in the high-density regime of the two-dimensional infinite-U Hubbard model. By decreasing the electron density, we observe a discontinuous transition to a paramagnetic phase, accompanied by a divergence of the susceptibility on the paramagnetic side. This behavior, resulting from a high degeneracy among different spin sectors, is consistent with an infinite-order phase transition scenario. In Chapter 5 the use of imaginary-time dynamics in the context of finite-temperature response functions is highlighted. As an application, we study an intriguing quantum phase featuring both glassy order and Bose-Einstein condensation. [35] We introduce and validate a model for the role of geometrical frustration in the coexistence of off-diagonal long range order with an amorphous density profile. The exact characterization of the response of the system to an external density perturbation is what allows here to establish the existence of a spin-glass phase. The differences between such a phase and the otherwise insulating "Bose glasses" are further elucidated in the Chapter. In the second part of the Thesis we focus our attention on the dynamics of closed systems out of equilibrium. This is accomplished by both non-stochastic exact methods for the dynamics and the introduction of a novel time-dependent Variational Monte Carlo scheme. In Chapter 6 exact diagonalization schemes and renormalization-based methods for one-dimensional systems are introduced. We identify key phenomenological traits resulting from the many-body correlation in closed systems driven sufficiently away from equilibrium.[31] We provide evidences that the dynamics of interacting lattice bosons away from equilibrium can be trapped into extremely long-lived inhomogeneous metastable states. The slowing down of incoherent density excitations above a threshold energy, much reminiscent of a dynamical arrest on the verge of a glass transition, is identified as the key feature of this phenomenon. In Chapter 7 we present an extension to dynamical properties of the Variational Quantum Monte Carlo method.[32] This is accomplished by introducing a general class of time-dependent variational states which is based on the mapping of the many-body dynamics onto an instantaneous ground-state problem. The application of the method to the experimentally relevant quantum quenches of interacting bosons reveals the accuracy and the reliability of the introduced numerical scheme. We indeed obtain for the first time a consistent variational description of the approach to the equilibrium of local observables and underline the origin of the metastability and glassy behavior previously identified. In the very last part we draw our conclusions and show some possible paths for stimulating future research.