Abstract:In this paper, we propose a new analytic continuation method to extract real frequency spectral functions from imaginary frequency Green's functions of quantum many-body systems. This method is based on the pole representation of Matsubara Green's function and a stochastic sampling procedure is utilized to optimize the amplitudes and locations of poles. In order to capture narrow peaks and sharp band edges in the spectral functions, a constrained sampling algorithm and a self-adaptive sampling algorithm are developed. To demonstrate the usefulness and performance of the new method, we at first apply it to study the spectral functions of representative fermionic and bosonic correlators. Then we employ this method to tackle the analytic continuation problems of matrix-valued Green's functions. The synthetic Green's functions, as well as realistic correlation functions from finite temperature quantum many-body calculations, are used as input. The benchmark results demonstrate that this method is capable of reproducing most of the key characteristics in the spectral functions. The sharp, smooth, and multi-peak features in both low-frequency and high-frequency regions of spectral functions could be accurately resolved, which overcomes one of the main limitations of the traditional maximum entropy method. More importantly, it exhibits excellent robustness with respect to noisy and incomplete input data. The causality of spectral function is always satisfied even in the presence of sizable noises. As a byproduct, this method could derive a fitting formula for the Matsubara data, which provides a compact approximation to the many-body Green's functions. Hence, we expect that this new method could become a pivotal workhorse for numerically analytic continuation and be broadly useful in many applications.

Minimal Pole Representation and Controlled Analytic Continuation of Matsubara Response Functions

Minimal pole representation and analytic continuation of matrix-valued correlation functions

Stochastic pole expansion method

Analytic continuation of multipoint correlation functions

A reliable Pade analytical continuation method based on a high accuracy symbolic computation algorithm

Feynman diagrammatics based on discrete pole representations: A path to renormalized perturbation theories

Using Analytic Continuation for the Hadronic Vacuum Polarization Computation

Hadronic Structure, Conformal Maps, and Analytic Continuation

Analytic Continuation in the Coupling Constant Method for the Dirac Equation

Reconstructing lattice QCD spectral functions with stochastic pole expansion and Nevanlinna analytic continuation

Nevanlinna Analytic Continuation for Migdal-Eliashberg Theory

Computing the hadronic vacuum polarization function by analytic continuation

ana_cont: Python package for analytic continuation

Bosonic Nevanlinna Analytic Continuation

From Experimental Data to Pole Parameters in a Model Independent Way (Fixed-t Single-Energy Partial Wave Analysis (SE PWA) and Laurent + Pietarinen Model (L+P) for Pole Extraction)

A convergence framework for Airy$_β$ line ensemble via pole evolution

Leading-logarithmic approximation by one-loop renormalization group within Matsubara formalism

Analytic continuation of the finite-volume three-particle amplitudes

The asymptotic approach to the continuum of lattice QCD spectral observables