Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2307.08566 · Physical Review X

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

Jason Kaye, Zhen Huang, Hugo U. R. Strand, Denis Golež

Published Aug 26, 2024• 12 citations
arXivHTMLPDF

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 ) 2 M − 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 Ca 2 RuO 4 , 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. Published by the American Physical Society 2024

DOI
10.1103/physrevx.14.031034

Related papers

Automated evaluation of imaginary time strong coupling diagrams by sum-of-exponentials hybridization fitting

Zhen Huang, Denis Golež, Hugo U. R. Strand +1 more

We present an efficient separation of variables algorithm for the evaluation of imaginary time Feynman diagrams appearing in the bold pseudo-particle strong coupling expansion of the Anderson impurity model. The algorithm uses a fitting method based on AAA rational approximation and numerical optimization to obtain a sum-of-exponentials expansion of the hybridization function, which is then used to decompose the diagrams. A diagrammatic formulation of the algorithm leads to an automated procedure for diagrams of arbitrary order and topology. We also present methods of stabilizing the self-consistent solution of the pseudo-particle Dyson equation. The result is a low-cost and high-order accurate impurity solver for quantum embedding methods using general multi-orbital hybridization functions at low temperatures, appropriate for low-to-intermediate expansion orders. In addition to other benchmark examples, we use our solver to perform a dynamical mean-field theory study of a minimal model of the strongly correlated compound Ca _2 2 RuO _4 4 , describing the anti-ferromagnetic transition and the in- and out-of-plane anisotropy induced by spin-orbit coupling.

arXiv:2503.19727•SciPost Physics•Nov 7, 2025•0 citations

Strong Coupling Quantum Impurity Solver on the Real and Imaginary Axes

Yueyi Wang, Kristjan Haule

The diagrammatic Monte Carlo method has so far been used mainly for weak-coupling expansions. Here we show that the strong-coupling expansion offers a key advantage: it can be implemented efficiently on both the real and imaginary frequency axes at finite temperature. Using a quantum-impurity solver for dynamical mean-field theory (DMFT) as an example, we find rapid convergence with expansion order. We derive closed-form real-axis Feynman rules for diagrams of arbitrary order, and implement them in a bold hybridization-expansion quantum Monte Carlo (BHQMC) impurity solver. Benchmarking against state-of-the-art numerical renormalization group (NRG) results for the DMFT Mott transition of the Hubbard model, we obtain a highly accurate frequency-dependent scattering rate at finite temperature. This enables reliable spectroscopy and provides benchmark transport results within DMFT and cluster-DMFT.

arXiv:2311.09412•Nov 15, 2023•1 citation

Diagrammatic Methods for Strongly-Correlated Electrons in Realistic Systems

Jia Li

This thesis is devoted to the development of numerical methods for correlated electrons in realistic systems based on perturbation theories using Feynman diagrams. Electron correlations in realistic systems give rise to a wide range of interesting physical properties. Diagrammatic expansions and its low-order approximations are key building blocks of theoretical and numerical methods for studying correlated electrons. Recent developments of the diagrammatic quantum Monte Carlo (DiagMC) methods enable efficient stochastic evaluations of diagrams and does not suffer from numerical sign problem which worsens with increasing system size. It is therefore a promising direction to apply DiagMC techniques to realistic electron systems in the pursuit of ab initio calculations of correlated materials. This thesis starts by reviewing basic concepts and definitions of realistic electron systems, the perturbation theory based on Feynman diagrams, Monte Carlo evaluations of the series, as well as self-consistent diagrammatic methods. Low-order diagrammatic methods are applied to realistic calculations and benchmarked against large collections of state-of-the-art many-body methods. The thesis then proposes a diagrammatic Monte Carlo solver for realistic molecules and impurities based on a recursive evaluation of diagrams, which is tested on realistic molecules and impurities with varying temperature, system size, and interaction strength. We observe that this method is ideal for problems with many orbitals with moderate correlations, but may suffer from series divergence at strong correlations. Next, we propose the inchworm Monte Carlo algorithm for interaction expansion, which overcomes limitations of divergent diagrammatic series by replacing the standard perturbation series with a sequence of incremental series expansions. A prototype implementation is presented for system with on-site Hubbard interactions. Finally, we discuss technical aspects of numerical representations of the Green's function, which has significant impact on the efficiency of diagrammatic methods. We review several common numerical representations, and introduce our development of the sparse sampling method, which provides compact representations of the Green's function in both imaginary time and frequency as well as fast transformation formula between the two.

•Jan 1, 2021•0 citations

Discrete Lehmann representation of imaginary time Green's functions

Jason Kaye, Kun Chen, Olivier Parcollet

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 $\ensuremath{\delta}$ functions. The basis is determined only by an upper bound on the product $\ensuremath{\beta}{\ensuremath{\omega}}_{max}$, with $\ensuremath{\beta}$ the inverse temperature and ${\ensuremath{\omega}}_{max}$ an energy cutoff, and a user-defined error tolerance $\ensuremath{\epsilon}$. The number $r$ of basis functions scales as $O(log(\ensuremath{\beta}{\ensuremath{\omega}}_{max})log(1/\ensuremath{\epsilon}))$. 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.

arXiv:2107.13094•Physical review. B./Physical review. B•Jun 13, 2022•61 citations

Renormalized diagrammatic perturbation theory using algorithmic Matsubara integration

Daria Gazizova

This thesis develops and applies numerical methods for computing correlated electrons at finite temperatures within perturbation theory, formulated in terms of Feynman diagrams. First we compute the fully renormalized effective interaction W in the 2D Hubbard model beyond the conventional low-order random phase approximation (RPA) and additionally resolve W into same-spin and opposite-spin channels. This analysis shows that starting from a purely local repulsive bare interaction U we get attractive effective nearest-neighbor interaction. We identify regimes of inverse temperature t and interaction U=t, where t denotes the nearest-neighbor hopping amplitude, in which this attraction appears. Such an emergent of attractive interaction is known to be able to promote superconductivity in the repulsive Hubbard model. Next, we introduce a method for incorporating higher-order diagrams within selfconsistent perturbative schemes. We merge algorithmic Matsubara integration (AMI) with discrete pole representations, including the discrete Lehmann representation and Prony (minimal pole) representation. AMI automates nested internal Matsubara frequency summations for arbitrary diagrams topologies through symbolic contour integration, which requires explicit access to the analytic pole structure of the input Green’s functions. By compressing the Green’s function into a small set of poles and corresponding weights, we regain access to its analytic pole structure, which allows us to iterate the calculation within self-consistent schemes. Benchmarks on the Bethe lattice demonstrate accurate single-shot and self-consistent self-energies on both Matsubara and real axes. For self-consistent calculations, we use two standard schemes, where AMI combined with the discrete pole representation works as a selfenergy and impurity solver. We further extend calculations to third and fourth order of interaction in self-energy diagrams. This approach yields the ability to perform self-consistent calculations at a fraction of the cost compared to other Matsubara-axis methods without moving to the real axis at every step of the iteration. Finally we demonstrate the universality of the AMI+DLR framework, by decomposing any diagram evaluation into a computational kernel that depends only on the diagram topology and DLR parameters, and a set of coefficients determined by the physical system. This separation enables precomputation and reuse of the computational kernel for a certain diagram topology and any physical problem, both in single-shot and iterative calculations. We illustrate the universality by reusing the same computational kernel for two simple systems: the 2D Hubbard square lattice and the H₂ molecule.

•May 1, 2026•0 citations

Multipoint Correlation Functions: Spectral Representation and Numerical Evaluation

Fabian B. Kugler, Seung‐Sup B. Lee, Jan von Delft

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.

arXiv:2101.00707•Physical Review X•Oct 11, 2021•48 citations
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 ) 2 M − 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 Ca 2 RuO 4 , 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. Published by the American Physical Society 2024
We present an efficient separation of variables algorithm for the evaluation of imaginary time Feynman diagrams appearing in the bold pseudo-particle strong coupling expansion of the Anderson impurity model. The algorithm uses a fitting method based on AAA rational approximation and numerical optimization to obtain a sum-of-exponentials expansion of the hybridization function, which is then used to decompose the diagrams. A diagrammatic formulation of the algorithm leads to an automated procedure for diagrams of arbitrary order and topology. We also present methods of stabilizing the self-consistent solution of the pseudo-particle Dyson equation. The result is a low-cost and high-order accurate impurity solver for quantum embedding methods using general multi-orbital hybridization functions at low temperatures, appropriate for low-to-intermediate expansion orders. In addition to other benchmark examples, we use our solver to perform a dynamical mean-field theory study of a minimal model of the strongly correlated compound Ca _2 2 RuO _4 4 , describing the anti-ferromagnetic transition and the in- and out-of-plane anisotropy induced by spin-orbit coupling.
The diagrammatic Monte Carlo method has so far been used mainly for weak-coupling expansions. Here we show that the strong-coupling expansion offers a key advantage: it can be implemented efficiently on both the real and imaginary frequency axes at finite temperature. Using a quantum-impurity solver for dynamical mean-field theory (DMFT) as an example, we find rapid convergence with expansion order. We derive closed-form real-axis Feynman rules for diagrams of arbitrary order, and implement them in a bold hybridization-expansion quantum Monte Carlo (BHQMC) impurity solver. Benchmarking against state-of-the-art numerical renormalization group (NRG) results for the DMFT Mott transition of the Hubbard model, we obtain a highly accurate frequency-dependent scattering rate at finite temperature. This enables reliable spectroscopy and provides benchmark transport results within DMFT and cluster-DMFT.
This thesis is devoted to the development of numerical methods for correlated electrons in realistic systems based on perturbation theories using Feynman diagrams. Electron correlations in realistic systems give rise to a wide range of interesting physical properties. Diagrammatic expansions and its low-order approximations are key building blocks of theoretical and numerical methods for studying correlated electrons. Recent developments of the diagrammatic quantum Monte Carlo (DiagMC) methods enable efficient stochastic evaluations of diagrams and does not suffer from numerical sign problem which worsens with increasing system size. It is therefore a promising direction to apply DiagMC techniques to realistic electron systems in the pursuit of ab initio calculations of correlated materials. This thesis starts by reviewing basic concepts and definitions of realistic electron systems, the perturbation theory based on Feynman diagrams, Monte Carlo evaluations of the series, as well as self-consistent diagrammatic methods. Low-order diagrammatic methods are applied to realistic calculations and benchmarked against large collections of state-of-the-art many-body methods. The thesis then proposes a diagrammatic Monte Carlo solver for realistic molecules and impurities based on a recursive evaluation of diagrams, which is tested on realistic molecules and impurities with varying temperature, system size, and interaction strength. We observe that this method is ideal for problems with many orbitals with moderate correlations, but may suffer from series divergence at strong correlations. Next, we propose the inchworm Monte Carlo algorithm for interaction expansion, which overcomes limitations of divergent diagrammatic series by replacing the standard perturbation series with a sequence of incremental series expansions. A prototype implementation is presented for system with on-site Hubbard interactions. Finally, we discuss technical aspects of numerical representations of the Green's function, which has significant impact on the efficiency of diagrammatic methods. We review several common numerical representations, and introduce our development of the sparse sampling method, which provides compact representations of the Green's function in both imaginary time and frequency as well as fast transformation formula between the two.
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 $\ensuremath{\delta}$ functions. The basis is determined only by an upper bound on the product $\ensuremath{\beta}{\ensuremath{\omega}}_{max}$, with $\ensuremath{\beta}$ the inverse temperature and ${\ensuremath{\omega}}_{max}$ an energy cutoff, and a user-defined error tolerance $\ensuremath{\epsilon}$. The number rrr of basis functions scales as $O(log(\ensuremath{\beta}{\ensuremath{\omega}}_{max})log(1/\ensuremath{\epsilon}))$. The discrete Lehmann representation of a particular imaginary time Green's function can be recovered by interpolation at a set of rrr 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.
This thesis develops and applies numerical methods for computing correlated electrons at finite temperatures within perturbation theory, formulated in terms of Feynman diagrams. First we compute the fully renormalized effective interaction W in the 2D Hubbard model beyond the conventional low-order random phase approximation (RPA) and additionally resolve W into same-spin and opposite-spin channels. This analysis shows that starting from a purely local repulsive bare interaction U we get attractive effective nearest-neighbor interaction. We identify regimes of inverse temperature t and interaction U=t, where t denotes the nearest-neighbor hopping amplitude, in which this attraction appears. Such an emergent of attractive interaction is known to be able to promote superconductivity in the repulsive Hubbard model. Next, we introduce a method for incorporating higher-order diagrams within selfconsistent perturbative schemes. We merge algorithmic Matsubara integration (AMI) with discrete pole representations, including the discrete Lehmann representation and Prony (minimal pole) representation. AMI automates nested internal Matsubara frequency summations for arbitrary diagrams topologies through symbolic contour integration, which requires explicit access to the analytic pole structure of the input Green’s functions. By compressing the Green’s function into a small set of poles and corresponding weights, we regain access to its analytic pole structure, which allows us to iterate the calculation within self-consistent schemes. Benchmarks on the Bethe lattice demonstrate accurate single-shot and self-consistent self-energies on both Matsubara and real axes. For self-consistent calculations, we use two standard schemes, where AMI combined with the discrete pole representation works as a selfenergy and impurity solver. We further extend calculations to third and fourth order of interaction in self-energy diagrams. This approach yields the ability to perform self-consistent calculations at a fraction of the cost compared to other Matsubara-axis methods without moving to the real axis at every step of the iteration. Finally we demonstrate the universality of the AMI+DLR framework, by decomposing any diagram evaluation into a computational kernel that depends only on the diagram topology and DLR parameters, and a set of coefficients determined by the physical system. This separation enables precomputation and reuse of the computational kernel for a certain diagram topology and any physical problem, both in single-shot and iterative calculations. We illustrate the universality by reusing the same computational kernel for two simple systems: the 2D Hubbard square lattice and the H₂ molecule.
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.