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

K. S. D. Beach; R. J. Gooding; F. Marsiglio

doi:10.1103/PhysRevB.61.5147

← Back to search

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

K. S. D. Beach, R. J. Gooding, F. Marsiglio

arXiv:cond-mat/9908477 · DOI 10.1103/PhysRevB.61.5147 · 1999-08-31 · Physical Review B

T3112 citations

Active bibliographic source — not scientific approval

Bibliographic access preserves source history; it does not approve extracted materials or validate reported claims. Review warnings on each occurrence separately.

Abstract

We critique a Pade analytic continuation method whereby a rational polynomial function is fit to a set of input points by means of a single matrix inversion. This procedure is accomplished to an extremely high accuracy using a novel symbolic computation algorithm. As an example of this method in action we apply it to the problem of determining the spectral function of a one-particle thermal Green's function known only at a finite number of Matsubara frequencies with two example self energies drawn from the T-matrix theory of the Hubbard model. We present a systematic analysis of the effects of error in the input points on the analytic continuation, and this leads us to propose a procedure to test quantitatively the reliability of the resulting continuation, thus eliminating the black magic label frequently attached to this procedure.