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.

arXiv comment19 pages, 1 table, 15 figures

References in corpus (49)

• The GW method
• Continuous-time Monte Carlo methods for quantum impurity models
• Functional renormalization group approach to correlated fermion systems
• Diagrammatic routes to nonlocal correlations beyond dynamical mean field theory
• Non-perturbative many-body approach to the Hubbard model and single-particle pseudogap
• The Bethe-Salpeter Equation Formalism: From Physics to Chemistry
• Local Electronic Correlation at the Two-Particle Level
• Tracking the Footprints of Spin Fluctuations: A MultiMethod, MultiMessenger Study of the Two-Dimensional Hubbard Model
• Compressing Green's function using intermediate representation between imaginary-time and real-frequency domains
• From infinite to two dimensions through the functional renormalization group
• Divergent Precursors of the Mott-Hubbard Transition at the Two-Particle Level
• Sparse sampling approach to efficient ab initio calculations at finite temperature
• Sparse modeling approach to analytical continuation of imaginary-time quantum Monte Carlo data
• QCD modeling of hadron physics
• High-frequency asymptotics of the vertex function: diagrammatic parametrization and algorithmic implementation
• Superconductivity, antiferromagnetism and phase separation in the two-dimensional Hubbard model: A dual-fermion approach
• Efficient implementation of the parquet equations -- role of the reducible vertex function and its kernel approximation
• Non-perturbative landscape of the Mott-Hubbard transition: Multiple divergence lines around the critical endpoint
• Ab initio dynamical vertex approximation
• Multiloop functional renormalization group that sums up all parquet diagrams
• Formation mechanism of helical Q structure in Gd-based skyrmion materials
• Efficient treatment of two-particle vertices in dynamical mean-field theory
• Divergences of the irreducible vertex functions in correlated metallic systems: Insights from the Anderson Impurity Model
• Fingerprints of the local moment formation and its Kondo screening in the generalized susceptibilities of many-electron problems
• Collective Charge Excitations of Strongly Correlated Electrons, Vertex Corrections and Gauge Invariance
• Truncated Unity Parquet Solver
• Generic optical excitations of correlated systems: $\pi$-tons
• Single-boson exchange decomposition of the vertex function
• Analytical investigation of singularities in two-particle irreducible vertex functions of the Hubbard atom
• Quantitative functional renormalization-group description of the two-dimensional Hubbard model
• Ab-Initio self-energy embedding for the photoemission spectra of NiO and MnO
• Sparse Modeling in Quantum Many-Body Problems
• irbasis: Open-source database and software for intermediate-representation basis functions of imaginary-time Green's function
• Two-particle response in Cluster Dynamical Mean-Field Theory: Formalism and application to the Raman Response of High-temperature Superconductors
• Overcomplete compact representation of two-particle Green's functions
• Performance analysis of a physically constructed orthogonal representation of imaginary-time Green's function
• Magnetic exchange coupling in cuprate-analog $d^9$ nickelates
• Effect of non-local correlations on the electronic structure of LiFeAs
• Symmetric Improved Estimators for Continuous-time Quantum Monte Carlo
• Efficient evaluation of the polarization function in dynamical mean-field theory
• Sparse sampling and tensor network representation of two-particle Green's functions
• Dynamically screened vertex correction to $GW$
• Efficient ab initio Migdal-Eliashberg calculation considering the retardation effect in phonon-mediated superconductors
• 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
• Bethe--Salpeter equation in QCD
• Truncated-Unity Parquet Equations: Application to the Repulsive Hubbard Model
• Local force method for the ab initio tight-binding model with spin-dependent hopping
• Functional renormalization group study of an eight-band model for the iron arsenides
• 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