Physical review. B, Condensed matter
Improved treatment of retarded self-energy calculations in perturbation theory
S. Wermbter
Abstract
I present a numerical treatment for a certain class of Feynman diagrams in perturbation theory with respect to the Coulomb correlation U. This method is formally exact for the frequency integrals occurring in the expression for the retarded Green's function. I specify this method to the self-consistent second-order perturbation theory (SOPT) for the single-impurity Anderson model and to the self-consistent fluctuation exchange approximation (FLEX) for the Hubbard model in infinite dimensions d\ensuremath{\rightarrow}\ensuremath{\infty}. The presented treatment is not only restricted to the limit d\ensuremath{\rightarrow}\ensuremath{\infty}. I present the resulting equations for the retarded self-energy, numerical results for the self-energy, density of states and the interaction spectrum in the FLEX for different scalings of the real frequency axis. The numerical solution of the SOPT and the FLEX could be obtained with very little numerical effort. In this form the calculation is nearly independent of the temperature and over a wide range nearly independent of the value of U. The computational time and the memory requirement for the solution of the self-energy equation are reduced to a minimum in this numerical treatment. \textcopyright{} 1996 The American Physical Society.