We critique a Pad\'e 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 symbolic computation algorithm. As an example of this method in action, it is applied to the problem of determining the spectral function of a single-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.We critique a Pad\'e 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 symbolic computation algorithm. As an example of this method in action, it is applied to the problem of determining the spectral function of a single-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.

arXiv:cond-mat/9908477Physical review. B, Condensed matterFeb 15, 2000140 citations

Comment on the diagonalization of Green functions

Peter A. Henning

arXiv:hep-ph/9502216Physics Letters BJul 1, 19951 citation

Approximation of Green's functions

F. Beck, Benno Fuchssteiner

Zeitschrift für Physik A Hadrons and NucleiJun 1, 19700 citations

Stochastic pole expansion method for analytic continuation of the Green's function

Li Huang, Shuang Liang

In this paper, we propose an 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 this 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 multipeak features in both low- 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 by-product, 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 method could become a pivotal workhorse for numerically analytic continuation and be broadly useful in many applications.In this paper, we propose an 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 this 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 multipeak features in both low- 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 by-product, 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 method could become a pivotal workhorse for numerically analytic continuation and be broadly useful in many applications.

arXiv:2307.11324Physical review. B./Physical review. BDec 13, 202311 citations

Photoemission spectral functions from the three-body Green's function

Gabriele Riva, Timothée Audinet, Matthieu Vladaj +2 more

We present an original strategy for the calculation of direct and inverse photo-emission spectra from first principles. The main goal is to go beyond the standard Green's function approaches, such as the $GW$ method, in order to find a good description not only of the quasiparticles but also of the satellite structures, which are of particular importance in strongly correlated materials. To this end we use as a key quantity the three-body Green's function, or, more precisely, its hole-hole-electron and electron-electron-hole parts, and we show how the one-body Green's function, and hence the corresponding spectral function, can be retrieved from it. We show that, contrary to the one-body Green's function, information about satellites is already present in the non-interacting three-body Green's function. Therefore, simple approximations to the three-body self-energy, which is defined by the Dyson equation for the three-body Green's function and which contains many-body effects, can still yield accurate spectral functions. In particular, the self-energy can be chosen to be static which could simplify a self-consistent solution of the Dyson equation. We give a proof of principle of our strategy by applying it to the Hubbard dimer, for which the exact self-energy is available.We present an original strategy for the calculation of direct and inverse photo-emission spectra from first principles. The main goal is to go beyond the standard Green's function approaches, such as the $GW$ method, in order to find a good description not only of the quasiparticles but also of the satellite structures, which are of particular importance in strongly correlated materials. To this end we use as a key quantity the three-body Green's function, or, more precisely, its hole-hole-electron and electron-electron-hole parts, and we show how the one-body Green's function, and hence the corresponding spectral function, can be retrieved from it. We show that, contrary to the one-body Green's function, information about satellites is already present in the non-interacting three-body Green's function. Therefore, simple approximations to the three-body self-energy, which is defined by the Dyson equation for the three-body Green's function and which contains many-body effects, can still yield accurate spectral functions. In particular, the self-energy can be chosen to be static which could simplify a self-consistent solution of the Dyson equation. We give a proof of principle of our strategy by applying it to the Hubbard dimer, for which the exact self-energy is available.

arXiv:2110.05623SciPost PhysicsMar 15, 202215 citations

On the concept of spectral singularities

Gusein Sh. Guseinov

PramanaSep 1, 200974 citations

GW

method, in order to find a good description not only of the quasiparticles but also of the satellite structures, which are of particular importance in strongly correlated materials. To this end we use as a key quantity the three-body Green's function, or, more precisely, its hole-hole-electron and electron-electron-hole parts, and we show how the one-body Green's function, and hence the corresponding spectral function, can be retrieved from it. We show that, contrary to the one-body Green's function, information about satellites is already present in the non-interacting three-body Green's function. Therefore, simple approximations to the three-body self-energy, which is defined by the Dyson equation for the three-body Green's function and which contains many-body effects, can still yield accurate spectral functions. In particular, the self-energy can be chosen to be static which could simplify a self-consistent solution of the Dyson equation. We give a proof of principle of our strategy by applying it to the Hubbard dimer, for which the exact self-energy is available.