A New Numerical Algorithm for the Analytic Continuation of Green's Functions

The need to calculate the spectral properties of a Hermitian operatorHfrequently arises in the technical sciences. A common approach to its solution involves the construction of the Green's function operatorG(z) = [z−H]−1in the complexzplane. For example, the energy spectrum and other physical...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational physics Vol. 126; no. 1; pp. 99 - 108
Main Authors: Natoli, Vincent D., Cohen, Morrel H., Fornberg, Bengt
Format: Journal Article
Language:English
Published: Elsevier Inc 01.06.1996
ISSN:0021-9991, 1090-2716
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:The need to calculate the spectral properties of a Hermitian operatorHfrequently arises in the technical sciences. A common approach to its solution involves the construction of the Green's function operatorG(z) = [z−H]−1in the complexzplane. For example, the energy spectrum and other physical properties of condensed matter systems can often be elegantly and naturally expressed in terms of the Kohn–Sham Green's functions. However, the nonanalyticity of resolvents on the real axis makes them difficult to compute and manipulate. The Herglotz property of a Green's function allows one to calculate it along an arc with a small but finite imaginary part, i.e.,G(x+iy), and then to continue it to the real axis to determine quantities of physical interest. In the past, finite-difference techniques have been used for this continuation. We present here a fundamentally new algorithm based on the fast Fourier transform which is both simpler and more effective.
Bibliography:ObjectType-Article-2
SourceType-Scholarly Journals-1
ObjectType-Feature-1
content type line 23
ISSN:0021-9991
1090-2716
DOI:10.1006/jcph.1996.0123