Discrete Lehmann representation of imaginary time Green's functions

Abstract

We present an efficient basis for imaginary time Green's functions based on a low rank decomposition of the spectral Lehmann representation. The basis functions are simply a set of well-chosen exponentials, so the corresponding expansion may be thought of as a discrete form of the Lehmann representation using an effective spectral density which is a sum of $\delta$ functions. The basis is determined only by an upper bound on the product $\beta {\omega }_{\max }$, with $\beta$ the inverse temperature and ${\omega }_{\max }$ an energy cutoff, and a user-defined error tolerance $\epsilon$. The number $r$ of basis functions scales as $\mathcal{O}\left(\log (\beta {\omega }_{\max })\log (1/\epsilon )\right)$. The discrete Lehmann representation of a particular imaginary time Green's function can be recovered by interpolation at a set of $r$ imaginary time nodes. Both the basis functions and the interpolation nodes can be obtained rapidly using standard numerical linear algebra routines. Due to the simple form of the basis, the discrete Lehmann representation of a Green's function can be explicitly transformed to the Matsubara frequency domain, or obtained directly by interpolation on a Matsubara frequency grid. We benchmark the efficiency of the representation on simple cases, and with a high precision solution of the Sachdev-Ye-Kitaev equation at low temperature. We compare our approach with the related intermediate representation method, and introduce an improved algorithm to build the intermediate representation basis and a corresponding sampling grid.

References in corpus (3)

• Sparse Modeling in Quantum Many-Body Problems
• libdlr: Efficient imaginary time calculations using the discrete Lehmann representation
• Efficient ab initio many-body calculations based on sparse modeling of Matsubara Green's function

Cited by in corpus (47)

• sparse-ir: optimal compression and sparse sampling of many-body propagators
• Multiscale space-time ansatz for correlation functions of quantum systems based on quantics tensor trains
• Fully Self-Consistent Finite-Temperature $GW$ in Gaussian Bloch Orbitals for Solids
• libdlr: Efficient imaginary time calculations using the discrete Lehmann representation
• Efficient ab initio many-body calculations based on sparse modeling of Matsubara Green's function
• Low-Scaling Algorithm for the Random Phase Approximation using Tensor Hypercontraction with k-point Sampling
• Deep Learning the Functional Renormalization Group
• A fast time domain solver for the equilibrium Dyson equation
• Minimal Pole Representation and Controlled Analytic Continuation of Matsubara Response Functions
• Excitations and spectra from equilibrium real-time Green's functions
• Superconductivity in the Uniform Electron Gas: Irrelevance of Kohn-Luttinger Mechanism
• Minimal pole representation and analytic continuation of matrix-valued correlation functions
• Compactness of quantics tensor train representations of local imaginary-time propagators
• Effective Modeling of Open Quantum Systems by Low-rank Discretization of Structured Environments
• Origin of the Coulomb Pseudopotential
• Adaptive Time Stepping for the Two-Time Integro-Differential Kadanoff-Baym Equations
• Decomposing imaginary time Feynman diagrams using separable basis functions: Anderson impurity model strong coupling expansion
• Equivariant neural network for Green's functions of molecules and materials
• Low rank Green's function representations applied to dynamical mean-field theory
• Minimal pole representation for spectral functions
• Stabilizing the calculation of the self-energy in dynamical mean-field theory using constrained residual minimization
• Quasiparticle and fully self-consistent GW methods: an unbiased analysis using Gaussian orbitals
• Low-rank quantics tensor train representations of Feynman diagrams for multiorbital electron-phonon models
• Chirped amplitude mode in photo-excited superconductors
• Inchworm quasi Monte Carlo for quantum impurities
• Analytic continuation of multipoint correlation functions
• Discrete Lehmann representation of three-point functions
• Compressing the two-particle Green's function using wavelets: Theory and application to the Hubbard atom
• Efficient Pseudomode Representation and Complexity of Quantum Impurity Models
• cppdlr: Imaginary time calculations using the discrete Lehmann representation
• Overcomplete intermediate representation of two-particle Green's functions and its relation to partial spectral functions
• Subleading logarithmic behavior in the parquet formalism
• Intertwined fluctuations and isotope effects in the Hubbard-Holstein model on the square lattice from functional renormalization
• Precise Low-Temperature Expansions for the Sachdev-Ye-Kitaev model
• Heating and cooling in self-consistent many-body simulations
• Compact representation and long-time extrapolation of real-time data for quantum systems using the ESPRIT algorithm
• Precursory Cooper Flow in Ultralow-Temperature Superconductors
• High-temperature series expansion of the dynamic Matsubara spin correlator
• From strong to weak correlations in breathing-mode kagome van der Waals materials: Nb${}_3$(F,Cl,Br,I)${}_8$ as a robust and versatile platform for many-body engineering
• Strange Metal to Insulator Transitions in the Lowest Landau Level
• Finite Temperature Minimal Entangled Typical Thermal States Impurity Solver
• Orbital selective order and ${\mathbb{Z}}_3$ Potts nematicity from a non-Fermi liquid
• Computation Kernel for Feynman Diagrams
• Inchworm tensor train hybridization expansion quantum impurity solver
• Diagrammatic bosonization, aspects of criticality, and the Hohenberg-Mermin-Wagner theorem in parquet approaches
• Onset of spin-valley order and Stoner boundaries in twisted WSe${}_2$
• Automated evaluation of imaginary time strong coupling diagrams by sum-of-exponentials hybridization fitting