Abstract:We present a deterministic algorithm for the efficient evaluation of imaginary-time diagrams based on the recently introduced discrete Lehmann representation (DLR) of imaginary-time Green's functions. In addition to the efficient discretization of diagrammatic integrals afforded by its approximation properties, the DLR basis is separable in imaginary-time, allowing us to decompose diagrams into linear combinations of nested sequences of one-dimensional products and convolutions. Focusing on the strong-coupling bold-line expansion of generalized Anderson impurity models, we show that our strategy reduces the computational complexity of evaluating an M th-order diagram at inverse temperature β and spectral width ωmax from O((βωmax)2M−1) for a direct quadrature to O(M(log(βωmax))M+1) , with controllable high-order accuracy. We benchmark our algorithm using third-order expansions for multiband impurity problems with off-diagonal hybridization and spin-orbit coupling, presenting comparisons with exact diagonalization and quantum Monte Carlo approaches. In particular, we perform a self-consistent dynamical mean-field theory calculation for a three-band Hubbard model with strong spin-orbit coupling representing a minimal model of Ca2RuO4 , demonstrating the promise of the method for modeling realistic strongly correlated multiband materials. For both strong and weak coupling expansions of low and intermediate order, in which diagrams can be enumerated, our method provides an efficient, straightforward, and robust blackbox evaluation procedure. In this sense, it fills a gap between diagrammatic approximations of the lowest order, which are simple and inexpensive but inaccurate, and those based on Monte Carlo sampling of high-order diagrams. https://doi.org/10.1103/PhysRevX.14.031034 Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI. Published by the American Physical Society

A Davidson-Lanczos iteration method for computation of continued-fraction expansion of the Green's function at very low temperatures: Applications to the dynamical mean field theory

Quantum Computed Green's Functions using a Cumulant Expansion of the Lanczos Method

Gauss Quadrature Based Finite Temperature Lanczos Method

Low rank compression in the numerical solution of the nonequilibrium Dyson equation

Stabilizing the calculation of the self-energy in dynamical mean-field theory using constrained residual minimization

Advanced Finite-Temperature Lanczos Method for anisotropic spin systems

Matrix-Product-State Algorithm for Finite Fractional Quantum Hall Systems

Solutions to the Landau-Lifshitz-Gilbert equation in the frequency space: Discretization schemes for the dynamic-matrix approach

Iterative subspace algorithms for finite-temperature solution of Dyson equation

Construction of the spectral function from non-commuting spectral moment matrices

Single particle properties of the 2D Hubbard model for real frequencies at weak coupling: Breakdown of the Dyson series for partial self-energy expansions

A Novel Hybrid Approach For Local Fractional Landau–Ginzburg–Higgs Equation Describing Fractal Heat Flow In Superconductors

Lanczos algorithm with matrix product states for dynamical correlation functions

Low-frequency Logarithmic Discretization of the Reservoir Spectrum for Improving the Efficiency of Hierarchical Equations of Motion Approach.

Lanczos, the transfer matrix, and the signal-to-noise problem

Lanczos for lattice QCD matrix elements

Decomposing Imaginary-Time Feynman Diagrams Using Separable Basis Functions: Anderson Impurity Model Strong-Coupling Expansion

Decomposing imaginary time Feynman diagrams using separable basis functions: Anderson impurity model strong coupling expansion

An Efficient Implementation of Analytical Nuclear Gradients for Linear-Response Time-Dependent Density Functional Theory in the Plane Wave Basis

Deflation of eigenvalues for iterative methods in lattice QCD