Spectral decomposition and high-accuracy Greens functions: Overcoming the Nyquist-Shannon limit via complex-time Krylov expansion

Sebastian Paeckel

2024-11-15

Abstract:The accurate computation of low-energy spectra of strongly correlated quantum many-body systems, typically accessed via Greens-functions, is a long-standing problem posing enormous challenges to numerical methods. When the spectral decomposition is obtained from Fourier transforming a time series, the Nyquist-Shannon theorem limits the frequency resolution $\Delta\omega$ according to the numerically accessible time domain size $T$ via $\Delta\omega = 2/T$. In tensor network methods, increasing the domain size is exponentially hard due to the ubiquitous spread of correlations, limiting the frequency resolution and thereby restricting this ansatz class mostly to one-dimensional systems with small quasi-particle velocities. Here, we show how this fundamental limitation can be overcome using complex-time Krylov spaces. At the example of the critical $S-1/2$ Heisenberg model and light bipolarons in the two-dimensional Su-Schrieffer-Heeger model, we demonstrate the enormous improvements in accuracy, which can be achieved using this method.

Strongly Correlated Electrons,Quantum Physics

What problem does this paper attempt to address?

This paper aims to solve the problem of accurate calculation of the low - energy spectrum of strongly correlated quantum many - body systems, which is usually accessed through the Green's function. When obtaining the spectral decomposition by Fourier transform from the time series, the Nyquist - Shannon sampling theorem limits the frequency resolution \(\Delta\omega\), and its relationship with the numerically accessible time - domain size \(T\) is \(\Delta\omega = 2/T\). In the tensor network method, due to the ubiquitous correlation diffusion, increasing the domain size is computationally exponentially difficult, which limits the frequency resolution and thus mainly restricts this method to one - dimensional systems with a small quasiparticle velocity.

This paper proposes a method using the complex - time Krylov space to overcome this fundamental limitation. Through examples of the critical \(S - 1/2\) Heisenberg model and the formation of optical bipolarons in the two - dimensional Su - Schrieffer - Heeger model, the great improvement in accuracy that can be achieved using this method is demonstrated. Specifically, this method combines real - time and complex - time evolutions to significantly improve the frequency resolution of the Green's function, even in two - dimensional systems.