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

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

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

Published Jan 31, 2024• 4 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.109.045152

Related papers

Real-frequency diagrammatic Monte Carlo at finite temperature

J. Vučičević, Michel Ferrero

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.

arXiv:1908.11826•Physical review. B./Physical review. B•Feb 12, 2020•35 citations

Algorithmic implementation of particle-particle ladder diagram approximation to study strongly-correlated metals and semiconductors

Aldiansah Prayogi, Muhammad Aziz Majidi

In condensed-matter physics, strongly-correlated systems refer to materials that exhibit variety of fascinating properties and ordered phases, depending on temperature, doping, and other factors. Such unique properties most notably arise due to strong electron-electron interactions, and in some cases due to interactions involving other quasiparticles as well. Electronic correlation effects are non-trivial that one may need a sufficiently accurate approximation technique with quite heavy computation, such as Quantum Monte-Carlo, in order to capture particular material properties arising from such effects. Meanwhile, less accurate techniques may come with lower numerical cost, but the ability to capture particular properties may highly depend on the choice of approximation. Among the many-body techniques derivable from Feynman diagrams, we aim to formulate algorithmic implementation of the Ladder Diagram approximation to capture the effects of electron-electron interactions. We wish to investigate how these correlation effects influence the temperature-dependent properties of strongly-correlated metals and semiconductors. As we are interested to study the temperature-dependent properties of the system, the Ladder diagram method needs to be applied in Matsubara frequency domain to obtain the self-consistent self-energy. However, at the end we would also need to compute the dynamical properties like density of states (DOS) and optical conductivity that are defined in the real frequency domain. For this purpose, we need to perform the analytic continuation procedure. At the end of this study, we will test the technique by observing the occurrence of metal-insulator transition in strongly-correlated metals, and renormalization of the band gap in strongly-correlated semiconductors.

AIP conference proceedings•Jan 1, 2017•0 citations

Dual Bethe-Salpeter equation for the multiorbital lattice susceptibility within dynamical mean-field theory

Erik G. C. P. van Loon, Hugo U. R. Strand

Dynamical mean-field theory describes the impact of strong local correlation effects in many-electron systems. While the single-particle spectral function is directly obtained within the formalism, two-particle susceptibilities can also be obtained by solving the Bethe-Salpeter equation. The solution requires handling infinite matrices in Matsubara frequency space. This is commonly treated using a finite frequency cutoff, resulting in slow linear convergence. A decomposition of the two-particle response in local and nonlocal contributions enables a reformulation of the Bethe-Salpeter equation inspired by the dual boson formalism. The reformulation has a drastically improved cubic convergence with respect to the frequency cutoff, considerably facilitating the calculation of susceptibilities in multi-orbital systems. This improved convergence arises from the fact that local contributions can be measured in the impurity solver. The dual Bethe-Salpeter equation uses the fully reducible vertex which is free from vertex divergences. We benchmark the approach on several systems including the spin susceptibility of strontium ruthenate Sr2RuO4 , a strongly correlated Hund's metal with three active orbitals. Published by the American Physical Society 2024

arXiv:2306.05157•Physical review. B./Physical review. B•Apr 23, 2024•7 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

Real-Frequency Diagrammatics with Pole Representations

Rayan Farid, Lei Zhang, Aiman Al-Eryani +2 more

Finite-temperature many-body theories are typically formulated in imaginary time, where they are amenable to efficient numerical treatment. However, the extraction of spectral properties from the imaginary axis is difficult and requires numerical analytic continuation, while a direct formulation on the real axis requires prohibitively expensive numerical quadrature. In this paper, we represent the objects of many-body perturbation theory using systematically improvable pole and moment expansions and develop the corresponding framework for evaluating diagrammatic contributions directly on the real axis. We introduce a controlled recompression procedure that prevents the proliferation of poles under algebraic and diagrammatic operations. We benchmark the approach on paradigmatic quantum impurity problems and apply it to self-consistent calculations of the uniform electron gas. The resulting spectral quantities converge systematically towards numerical accuracy. This work provides a practical framework for controlled and systematically improvable calculations of spectral quantities within diagrammatic many-body theories, including self-consistent conserving approaches.

arXiv:2609.26703•Sep 22, 2026•0 citations

Real-Time Approach to the Dynamical Bethe-Salpeter Equation for Finite Systems

Tucker Allen, Barry Yangtao Li, Daniel Neuhauser

We present a real-time linear-response approach to solving the dynamical Bethe-Salpeter equation (BSE). The polarization part of the screened interaction is obtained from time-dependent Hartree propagation of the one-particle density matrix within an orbital basis-set representation. The frequency convolution defining the dynamical kernel is evaluated as a product in time, avoiding numerical integration and storage of the full screened Coulomb operator. The resulting nonlinear eigenvalue problem is then solved directly, going beyond the static screening approximation with full-frequency dependence in the screened interaction. We validate the method by calculating the low-lying singlet and triplet excitation energies of water, formaldehyde, and acetylene. While the deterministic scaling remains $\mathcal{O}(N^6)$, the time propagation formulation is readily compatible with grid-based stochastic sampling methods, which will open the possibility for dynamical BSE calculations of very large systems.

arXiv:2608.08948•Aug 9, 2026•0 citations
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.
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.
In condensed-matter physics, strongly-correlated systems refer to materials that exhibit variety of fascinating properties and ordered phases, depending on temperature, doping, and other factors. Such unique properties most notably arise due to strong electron-electron interactions, and in some cases due to interactions involving other quasiparticles as well. Electronic correlation effects are non-trivial that one may need a sufficiently accurate approximation technique with quite heavy computation, such as Quantum Monte-Carlo, in order to capture particular material properties arising from such effects. Meanwhile, less accurate techniques may come with lower numerical cost, but the ability to capture particular properties may highly depend on the choice of approximation. Among the many-body techniques derivable from Feynman diagrams, we aim to formulate algorithmic implementation of the Ladder Diagram approximation to capture the effects of electron-electron interactions. We wish to investigate how these correlation effects influence the temperature-dependent properties of strongly-correlated metals and semiconductors. As we are interested to study the temperature-dependent properties of the system, the Ladder diagram method needs to be applied in Matsubara frequency domain to obtain the self-consistent self-energy. However, at the end we would also need to compute the dynamical properties like density of states (DOS) and optical conductivity that are defined in the real frequency domain. For this purpose, we need to perform the analytic continuation procedure. At the end of this study, we will test the technique by observing the occurrence of metal-insulator transition in strongly-correlated metals, and renormalization of the band gap in strongly-correlated semiconductors.
Dynamical mean-field theory describes the impact of strong local correlation effects in many-electron systems. While the single-particle spectral function is directly obtained within the formalism, two-particle susceptibilities can also be obtained by solving the Bethe-Salpeter equation. The solution requires handling infinite matrices in Matsubara frequency space. This is commonly treated using a finite frequency cutoff, resulting in slow linear convergence. A decomposition of the two-particle response in local and nonlocal contributions enables a reformulation of the Bethe-Salpeter equation inspired by the dual boson formalism. The reformulation has a drastically improved cubic convergence with respect to the frequency cutoff, considerably facilitating the calculation of susceptibilities in multi-orbital systems. This improved convergence arises from the fact that local contributions can be measured in the impurity solver. The dual Bethe-Salpeter equation uses the fully reducible vertex which is free from vertex divergences. We benchmark the approach on several systems including the spin susceptibility of strontium ruthenate Sr2RuO4 , a strongly correlated Hund's metal with three active orbitals. Published by the American Physical Society 2024
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.
Finite-temperature many-body theories are typically formulated in imaginary time, where they are amenable to efficient numerical treatment. However, the extraction of spectral properties from the imaginary axis is difficult and requires numerical analytic continuation, while a direct formulation on the real axis requires prohibitively expensive numerical quadrature. In this paper, we represent the objects of many-body perturbation theory using systematically improvable pole and moment expansions and develop the corresponding framework for evaluating diagrammatic contributions directly on the real axis. We introduce a controlled recompression procedure that prevents the proliferation of poles under algebraic and diagrammatic operations. We benchmark the approach on paradigmatic quantum impurity problems and apply it to self-consistent calculations of the uniform electron gas. The resulting spectral quantities converge systematically towards numerical accuracy. This work provides a practical framework for controlled and systematically improvable calculations of spectral quantities within diagrammatic many-body theories, including self-consistent conserving approaches.
We present a real-time linear-response approach to solving the dynamical Bethe-Salpeter equation (BSE). The polarization part of the screened interaction is obtained from time-dependent Hartree propagation of the one-particle density matrix within an orbital basis-set representation. The frequency convolution defining the dynamical kernel is evaluated as a product in time, avoiding numerical integration and storage of the full screened Coulomb operator. The resulting nonlinear eigenvalue problem is then solved directly, going beyond the static screening approximation with full-frequency dependence in the screened interaction. We validate the method by calculating the low-lying singlet and triplet excitation energies of water, formaldehyde, and acetylene. While the deterministic scaling remains O(N6)\mathcal{O}(N^6)O(N6), the time propagation formulation is readily compatible with grid-based stochastic sampling methods, which will open the possibility for dynamical BSE calculations of very large systems.