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

Calculation of Feynman Diagrams from their Small Momentum Expansion

ALGORITHMIC CALCULATION OF TWO-LOOP FEYNMAN DIAGRAMS

Near threshold expansion of Feynman diagrams

Photon rockets and gravitational radiation

Simplifying the large-mass expansion of Feynman integrals

3-loop Feynman Integral Extrapolations for the Baseball Diagram

Feynman Diagrams as Computational Graphs

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

An exact solution method for the enumeration of connected Feynman diagrams

The development of computational methods for Feynman diagrams

Glue-and-Cut at Five Loops

Ladder and zig-zag Feynman diagrams, operator formalism and conformal triangles

Analytic two-loop results for selfenergy- and vertex-type diagrams with one non-zero mass

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

Vanishing cycles and analysis of singularities of Feynman diagrams

From equations in coordinate space to Picard-Fuchs and back

The asymptotic Hopf Algebra of Feynman Integrals

Bold Feynman diagrams and the Luttinger-Ward formalism via Gibbs measures. Part II: Non-perturbative analysis

Nonconvergence of the Feynman-Dyson diagrammatic perturbation expansion of propagators

Feynman diagrammatics based on discrete pole representations: A path to renormalized perturbation theories

Combinatorial summation of Feynman diagrams: Equation of state of the 2D SU(N) Hubbard model