Solving the Bethe-Salpeter equation with exponential convergence

Abstract

The Bethe-Salpeter equation plays a crucial role in understanding the physics of correlated fermions, relating to optical excitations in solids as well as resonances in high-energy physics. Yet, it is notoriously difficult to control numerically, typically requiring an effort that scales polynomially with energy scales and accuracy. This puts many interesting systems out of computational reach. Using the intermediate representation and sparse modelling for two-particle objects on the Matsubara axis, we develop an algorithm that solves the Bethe-Salpeter equation in $O(L^8)$ time with $O(L^4)$ memory, where $L$ grows only logarithmically with inverse temperature, bandwidth, and desired accuracy, This opens the door for computations in hitherto inaccessible regimes. We benchmark the method on the Hubbard atom and on the multi-orbital weak-coupling limit, where we observe the expected exponential convergence to the analytical results. We then showcase the method for a realistic impurity problem.

19 pages, 1 table, 15 figures

References in corpus (11)

• Continuous-time Monte Carlo methods for quantum impurity models
• The Bethe-Salpeter Equation Formalism: From Physics to Chemistry
• Superconductivity, antiferromagnetism and phase separation in the two-dimensional Hubbard model: A dual-fermion approach
• Ab initio dynamical vertex approximation
• Collective Charge Excitations of Strongly Correlated Electrons, Vertex Corrections and Gauge Invariance
• Sparse Modeling in Quantum Many-Body Problems
• Two-particle response in Cluster Dynamical Mean-Field Theory: Formalism and application to the Raman Response of High-temperature Superconductors
• Efficient evaluation of the polarization function in dynamical mean-field theory
• Dynamically screened vertex correction to $GW$
• Efficient fluctuation exchange approach to low-temperature spin fluctuations and superconductivity: from the Hubbard model to Na${}_x$CoO${}_2\cdot y$H${}_2$O
• Summation of Parquet diagrams as an ab initio method in nuclear structure calculations

Cited by in corpus (25)

• Discrete Lehmann representation of imaginary time Green's functions
• sparse-ir: optimal compression and sparse sampling of many-body propagators
• Multiscale space-time ansatz for correlation functions of quantum systems based on quantics tensor trains
• Efficient ab initio many-body calculations based on sparse modeling of Matsubara Green's function
• Parametrizations of local vertex corrections from weak to strong coupling: importance of the Hedin three-leg vertex
• Strongly correlated superconductivity with long-range spatial fluctuations
• Robust analytic continuation combining the advantages of the sparse modeling approach and Padé approximation
• Two-particle calculations with quantics tensor trains: Solving the parquet equations
• Real-frequency quantum field theory applied to the single-impurity Anderson model
• Data-driven compression of electron-phonon interactions
• Compactness of quantics tensor train representations of local imaginary-time propagators
• Low rank Green's function representations applied to dynamical mean-field theory
• Multipolar ordering from dynamical mean field theory with application to CeB6
• Dual Bethe-Salpeter equation for the multi-orbital lattice susceptibility within dynamical mean-field theory
• Analytic continuation of multipoint correlation functions
• Discrete Lehmann representation of three-point functions
• Compressing the two-particle Green's function using wavelets: Theory and application to the Hubbard atom
• Overcomplete intermediate representation of two-particle Green's functions and its relation to partial spectral functions
• Subleading logarithmic behavior in the parquet formalism
• Parquet theory for molecular systems: Formalism and static kernel parquet approximation
• Suppression of the charge fluctuations by nonlocal correlations close to the Mott transition
• Entanglement across scales: Quantics tensor trains as a natural framework for renormalization
• Compressing local vertex functions from the multipoint numerical renormalization group using quantics tensor cross interpolation
• Antiferromagnetic magnons and local anisotropy: dynamical mean-field study
• Automated evaluation of imaginary time strong coupling diagrams by sum-of-exponentials hybridization fitting