A fast time domain solver for the equilibrium Dyson equation

TL;DR

This paper introduces a fast, FFT-based numerical solver for the real-time equilibrium Dyson equation, enabling efficient simulation of quantum many-body systems over large timescales.

Contribution

It presents a novel quasilinear-scaling algorithm that respects the nonlinear structure of the Dyson equation, improving computational efficiency for quantum dynamics.

Findings

The solver efficiently handles large propagation times.

It is applicable to models like Bethe graph and Sachdev-Ye-Kitaev.

The method reduces computational complexity compared to traditional approaches.

Abstract

We consider the numerical solution of the real time equilibrium Dyson equation, which is used in calculations of the dynamical properties of quantum many-body systems. We show that this equation can be written as a system of coupled, nonlinear, convolutional Volterra integro-differential equations, for which the kernel depends self-consistently on the solution. As is typical in the numerical solution of Volterra-type equations, the computational bottleneck is the quadratic-scaling cost of history integration. However, the structure of the nonlinear Volterra integral operator precludes the use of standard fast algorithms. We propose a quasilinear-scaling FFT-based algorithm which respects the structure of the nonlinear integral operator. The resulting method can reach large propagation times, and is thus well-suited to explore quantum many-body phenomena at low energy scales. We demonstrate the solver with two standard model systems: the Bethe graph, and the Sachdev-Ye-Kitaev model.