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 modeling for two-particle objects on the Matsubara axis, we develop an algorithm that solves the Bethe-Salpeter equation in time with memory, where 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 multiorbital weak-coupling limit, where we observe the expected exponential convergence to the analytical results. We then showcase the method for a realistic impurity problem.
8 More- Received 10 December 2020
- Accepted 23 July 2021
DOI:https://doi.org/10.1103/PhysRevResearch.3.033168
Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.
Published by the American Physical Society