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

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 $\beta$ and spectral width ${\omega }_{\max }$ from $\mathcal{O}((\beta {\omega }_{\max }{)}^{2M-1})$ for a direct quadrature to $\mathcal{O}(M(\log (\beta {\omega }_{\max }){)}^{M+1})$, with controllable high-order accuracy. We benchmark our algorithm using third-order expansions for multi-band 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 Ca${}_2$RuO${}_4$, demonstrating the promise of the method for modeling realistic strongly correlated multi-band 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 black-box 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.

References in corpus (20)

• Mott Insulators in the Strong Spin-Orbit Coupling Limit: From Heisenberg to a Quantum Compass and Kitaev Models
• Continuous-time Monte Carlo methods for quantum impurity models
• Quantum Monte Carlo Impurity Solver for Cluster DMFT and Electronic Structure Calculations in Adjustable Base
• Hybridization expansion impurity solver: General formulation and application to Kondo lattice and two-orbital models
• The GW compendium: A practical guide to theoretical photoemission spectroscopy
• Spin freezing transition and non-Fermi-liquid self-energy in a 3-orbital model
• The coherence-incoherence crossover and the mass-renormalization puzzles in Sr2RuO4
• Strongly Correlated Superconductivity: a plaquette Dynamical mean field theory study
• Bold Diagrammatic Monte Carlo: When Sign Problem is Welcome
• Magnetic couplings, optical spectra, and spin-orbit exciton in 5d electron Mott insulator Sr2IrO4
• Multitier self-consistent $GW$+EDMFT
• The Fermi surface of Sr2RuO4: spin-orbit and anisotropic Coulomb interaction effects
• Bold Line Diagrammatic Monte Carlo Method: General formulation and application to expansion around the Non-Crossing Approximation
• Subband filling and Mott transition in Ca_{2-x}Sr_xRuO_4
• Currents and Green's functions of impurities out of equilibrium -- results from inchworm Quantum Monte Carlo
• Worm Improved Estimators in Continuous-time Quantum Monte Carlo
• Magnetic Excitations in Spin-Orbit Coupled $d^4$ Mott Insulator on Square Lattice
• Spectral functions of Sr${}_2$IrO${}_4$: theory versus experiment
• Low-Scaling Algorithm for the Random Phase Approximation using Tensor Hypercontraction with k-point Sampling
• Quantitative comparison of Anderson impurity solvers applied to transport in quantum dots

Cited by in corpus (9)

• Minimal pole representation for spectral functions
• Stabilizing the calculation of the self-energy in dynamical mean-field theory using constrained residual minimization
• Low-rank quantics tensor train representations of Feynman diagrams for multiorbital electron-phonon models
• Inchworm quasi Monte Carlo for quantum impurities
• Discrete Lehmann representation of three-point functions
• Third-order strong-coupling impurity solver for real-frequency DMFT: Accurate spectral functions for antiferromagnetic and photo-doped states
• Computation Kernel for Feynman Diagrams
• Inchworm tensor train hybridization expansion quantum impurity solver
• Automated evaluation of imaginary time strong coupling diagrams by sum-of-exponentials hybridization fitting