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:1908.11826 · Physical review. B./Physical review. B

Real-frequency diagrammatic Monte Carlo at finite temperature

J. Vučičević, Michel Ferrero

Published Feb 12, 2020• 35 citations
arXivHTMLPDF

Abstract

Diagrammatic expansions are a central tool for treating correlated electron systems. At thermal equilibrium, they are most naturally defined within the Matsubara formalism. However, extracting any dynamic response function from a Matsubara calculation ultimately requires the ill-defined analytical continuation from the imaginary- to the real-frequency domain. It was recently proposed [A. Taheridehkordi et al., Phys. Rev. B 99, 035120 (2019)] that the internal Matsubara summations of any interaction-expansion diagram can be performed analytically by using symbolic algebra algorithms. The result of the summations is then an analytical function of the complex frequency rather than Matsubara frequency. Here we apply this principle and develop a diagrammatic Monte Carlo technique which yields results directly on the real-frequency axis. We present results for the self-energy $\mathrm{\ensuremath{\Sigma}}(\ensuremath{\omega})$ of the doped $32\ifmmode\times\else\texttimes\fi{}32$ cyclic square-lattice Hubbard model in a nontrivial parameter regime, where signatures of the pseudogap appear close to the antinode. We discuss the behavior of the perturbation series on the real-frequency axis and in particular show that one must be very careful when using the maximum entropy method on truncated perturbation series. Our approach holds great promise for future application in cases when analytical continuation is difficult and moderate-order perturbation theory may be sufficient to converge the result.

DOI
10.1103/physrevb.101.075113

Related papers

Algorithmic approach to diagrammatic expansions for real-frequency evaluation of susceptibility functions

Amir Taheridehkordi, Stephanie H. Curnoe, J. P. F. LeBlanc

We systematically generate the perturbative expansion for the two-particle spin susceptibility in the Feynman diagrammatic formalism and apply this expansion to a model system---the single-band Hubbard model on a square lattice. We make use of algorithmic Matsubara integration (AMI) [A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlanc, Phys. Rev. B 99, 035120 (2019)] to analytically evaluate Matsubara frequency summations, allowing us to symbolically impose analytic continuation to the real-frequency axis. We minimize our computational expense by applying graph invariant transformations [A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlanc, Phys. Rev. B 101, 125109 (2020)]. We highlight extensions of the random-phase approximation and T-matrix methods that, due to AMI, become tractable. We present results for weak interaction strength where the direct perturbative expansion is convergent, and verify our results on the Matsubara axis by comparison to other numerical methods. By examining the spin susceptibility as a function of real frequency via an order-by-order expansion, we can identify precisely what role higher-order corrections play on spin susceptibility and demonstrate the utility and limitations of our approach.

arXiv:2004.11091•Physical review. B./Physical review. B•Jul 9, 2020•28 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

Solving the Bethe-Salpeter equation in real frequencies at finite temperature

I. S. Tupitsyn, Nikolay V. Prokof'ev

Self-consistent Hartree-Fock approximation combined with solutions of the Bethe-Salpeter equation offers a powerful tool for studies of strong correlation effects arising in condensed matter models, nuclear physics, quantum field theories, and real materials. The standard finite-temperature approach would be to first solve the problem in the Matsubara representation and then apply numerical analytic continuation to the real-frequency axis to link theoretical results with experimental probes, but this ill-conditioned procedure often distorts important spectral features even for very accurate imaginary-frequency data. We demonstrate that the ladder-type finite-temperature Bethe-Salpeter equation in the Hartree-Fock basis for the three-point vertex function and, ultimately, system's polarization can be accurately solved directly on the real frequency axis using the diagrammatic Monte Carlo technique and series resummation. We illustrate the method by applying it to the homogeneous electron gas model and demonstrate how multiple scattering events renormalize Landau damping.

arXiv:2309.16124•Physical review. B./Physical review. B•Jan 31, 2024•4 citations

Dynamic Response of an Electron Gas: Towards the Exact Exchange-Correlation Kernel

J. P. F. LeBlanc, Kun Chen, Kristjan Haule +2 more

Precise calculations of dynamics in the homogeneous electron gas (jellium model) are of fundamental importance for design and characterization of new materials. We introduce a diagrammatic Monte Carlo technique based on algorithmic Matsubara integration that allows us to compute frequency and momentum resolved finite temperature response directly in the real frequency domain using a series of connected Feynman diagrams. The data for charge response at moderate electron density are used to extract the frequency dependence of the exchange-correlation kernel at finite momenta and temperature. These results are as important for development of the time-dependent density functional theory for materials dynamics as ground state energies are for the density functional theory.

arXiv:2205.13595•Physical Review Letters•Dec 7, 2022•31 citations

Diagrammatic Monte Carlo approach for diagrammatic extensions of dynamical mean-field theory: Convergence analysis of the dual fermion technique

Jan Gukelberger, Evgeny Kozik, Hartmut Hafermann

The dual fermion approach provides a formally exact prescription for calculating properties of a correlated electron system in terms of a diagrammatic expansion around dynamical mean-field theory (DMFT). Most practical implementations, however, neglect higher-order interaction vertices beyond two-particle scattering in the dual effective action and further truncate the diagrammatic expansion in the two-particle scattering vertex to a leading-order or ladder-type approximation. In this work, we compute the dual fermion expansion for the two-dimensional Hubbard model including all diagram topologies with two-particle interactions to high orders by means of a stochastic diagrammatic Monte Carlo algorithm. We benchmark the obtained self-energy against numerically exact diagrammatic determinant Monte Carlo simulations to systematically assess convergence of the dual fermion series and the validity of these approximations. We observe that, from high temperatures down to the vicinity of the DMFT N\'eel transition, the dual fermion series converges very quickly to the exact solution in the whole range of Hubbard interactions considered ($4\ensuremath{\le}U/t\ensuremath{\le}12$), implying that contributions from higher-order vertices are small. As the temperature is lowered further, we observe slower series convergence, convergence to incorrect solutions, and ultimately divergence. This happens in a regime where magnetic correlations become significant. We find, however, that the self-consistent particle-hole ladder approximation yields reasonable and often even highly accurate results in this regime.

arXiv:1611.07523•Physical review. B./Physical review. B•Jul 28, 2017•55 citations

Algorithmic Matsubara integration for Hubbard-like models

Amir Taheridehkordi, Stephanie H. Curnoe, J. P. F. LeBlanc

We present an algorithm to evaluate Matsubara sums for Feynman diagrams composed of bare Green's functions with single-band dispersions and local $U$ Hubbard interaction vertices. The algorithm provides an exact construction of the analytic result for the frequency integrals of a diagram that can then be evaluated for all parameters $U$, temperature $T$, chemical potential $\ensuremath{\mu}$, external frequencies, and internal/external momenta. This method allows for symbolic analytic continuation of results to the real frequency axis, avoiding any ill-posed numerical procedure. This method can also be used to simultaneously evaluate diagrams throughout the entire $T\text{\ensuremath{-}}U\text{\ensuremath{-}}\ensuremath{\mu}$ phase space of Hubbard-like models even in the $T=0$ limit at minimal computational expense.

arXiv:1808.05188•Physical review. B./Physical review. B•Jan 10, 2019•49 citations
Diagrammatic expansions are a central tool for treating correlated electron systems. At thermal equilibrium, they are most naturally defined within the Matsubara formalism. However, extracting any dynamic response function from a Matsubara calculation ultimately requires the ill-defined analytical continuation from the imaginary- to the real-frequency domain. It was recently proposed [A. Taheridehkordi et al., Phys. Rev. B 99, 035120 (2019)] that the internal Matsubara summations of any interaction-expansion diagram can be performed analytically by using symbolic algebra algorithms. The result of the summations is then an analytical function of the complex frequency rather than Matsubara frequency. Here we apply this principle and develop a diagrammatic Monte Carlo technique which yields results directly on the real-frequency axis. We present results for the self-energy \ensuremathΣ(\ensuremathω)\mathrm{\ensuremath{\Sigma}}(\ensuremath{\omega})\ensuremathΣ(\ensuremathω) of the doped 32\ifmmode×\else\texttimes\fi3232\ifmmode\times\else\texttimes\fi{}3232\ifmmode×\else\texttimes\fi32 cyclic square-lattice Hubbard model in a nontrivial parameter regime, where signatures of the pseudogap appear close to the antinode. We discuss the behavior of the perturbation series on the real-frequency axis and in particular show that one must be very careful when using the maximum entropy method on truncated perturbation series. Our approach holds great promise for future application in cases when analytical continuation is difficult and moderate-order perturbation theory may be sufficient to converge the result.
We systematically generate the perturbative expansion for the two-particle spin susceptibility in the Feynman diagrammatic formalism and apply this expansion to a model system---the single-band Hubbard model on a square lattice. We make use of algorithmic Matsubara integration (AMI) [A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlanc, Phys. Rev. B 99, 035120 (2019)] to analytically evaluate Matsubara frequency summations, allowing us to symbolically impose analytic continuation to the real-frequency axis. We minimize our computational expense by applying graph invariant transformations [A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlanc, Phys. Rev. B 101, 125109 (2020)]. We highlight extensions of the random-phase approximation and T-matrix methods that, due to AMI, become tractable. We present results for weak interaction strength where the direct perturbative expansion is convergent, and verify our results on the Matsubara axis by comparison to other numerical methods. By examining the spin susceptibility as a function of real frequency via an order-by-order expansion, we can identify precisely what role higher-order corrections play on spin susceptibility and demonstrate the utility and limitations of our approach.
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.
Self-consistent Hartree-Fock approximation combined with solutions of the Bethe-Salpeter equation offers a powerful tool for studies of strong correlation effects arising in condensed matter models, nuclear physics, quantum field theories, and real materials. The standard finite-temperature approach would be to first solve the problem in the Matsubara representation and then apply numerical analytic continuation to the real-frequency axis to link theoretical results with experimental probes, but this ill-conditioned procedure often distorts important spectral features even for very accurate imaginary-frequency data. We demonstrate that the ladder-type finite-temperature Bethe-Salpeter equation in the Hartree-Fock basis for the three-point vertex function and, ultimately, system's polarization can be accurately solved directly on the real frequency axis using the diagrammatic Monte Carlo technique and series resummation. We illustrate the method by applying it to the homogeneous electron gas model and demonstrate how multiple scattering events renormalize Landau damping.
Precise calculations of dynamics in the homogeneous electron gas (jellium model) are of fundamental importance for design and characterization of new materials. We introduce a diagrammatic Monte Carlo technique based on algorithmic Matsubara integration that allows us to compute frequency and momentum resolved finite temperature response directly in the real frequency domain using a series of connected Feynman diagrams. The data for charge response at moderate electron density are used to extract the frequency dependence of the exchange-correlation kernel at finite momenta and temperature. These results are as important for development of the time-dependent density functional theory for materials dynamics as ground state energies are for the density functional theory.
The dual fermion approach provides a formally exact prescription for calculating properties of a correlated electron system in terms of a diagrammatic expansion around dynamical mean-field theory (DMFT). Most practical implementations, however, neglect higher-order interaction vertices beyond two-particle scattering in the dual effective action and further truncate the diagrammatic expansion in the two-particle scattering vertex to a leading-order or ladder-type approximation. In this work, we compute the dual fermion expansion for the two-dimensional Hubbard model including all diagram topologies with two-particle interactions to high orders by means of a stochastic diagrammatic Monte Carlo algorithm. We benchmark the obtained self-energy against numerically exact diagrammatic determinant Monte Carlo simulations to systematically assess convergence of the dual fermion series and the validity of these approximations. We observe that, from high temperatures down to the vicinity of the DMFT N\'eel transition, the dual fermion series converges very quickly to the exact solution in the whole range of Hubbard interactions considered ($4\ensuremath{\le}U/t\ensuremath{\le}12$), implying that contributions from higher-order vertices are small. As the temperature is lowered further, we observe slower series convergence, convergence to incorrect solutions, and ultimately divergence. This happens in a regime where magnetic correlations become significant. We find, however, that the self-consistent particle-hole ladder approximation yields reasonable and often even highly accurate results in this regime.
We present an algorithm to evaluate Matsubara sums for Feynman diagrams composed of bare Green's functions with single-band dispersions and local UUU Hubbard interaction vertices. The algorithm provides an exact construction of the analytic result for the frequency integrals of a diagram that can then be evaluated for all parameters UUU, temperature TTT, chemical potential $\ensuremath{\mu}$, external frequencies, and internal/external momenta. This method allows for symbolic analytic continuation of results to the real frequency axis, avoiding any ill-posed numerical procedure. This method can also be used to simultaneously evaluate diagrams throughout the entire $T\text{\ensuremath{-}}U\text{\ensuremath{-}}\ensuremath{\mu}$ phase space of Hubbard-like models even in the T=0T=0T=0 limit at minimal computational expense.