Pad\'e approximants for improved finite-temperature spectral functions in the numerical renormalization group

TL;DR

This paper presents an improved method using Padé approximants for obtaining smooth, high-resolution finite-temperature spectral functions in quantum impurity models via NRG, reducing artifacts and resolving fine spectral details.

Contribution

The authors introduce a Padé approximant-based analytic continuation method for NRG spectral functions, enhancing resolution and reliability over traditional broadening techniques.

Findings

Better resolution of spectral features at low frequencies.

Reduced artifacts compared to conventional broadening.

Effective in both Anderson and Hubbard models.

Abstract

We introduce an improved approach for obtaining smooth finite-temperature spectral functions of quantum impurity models using the numerical renormalization group (NRG) technique. It is based on calculating first the Green's function on the imaginary-frequency axis, followed by an analytic continuation to the real-frequency axis using Pad\'e approximants. The arbitrariness in choosing a suitable kernel in the conventional broadening approach is thereby removed and, furthermore, we find that the Pad\'e method is able to resolve fine details in spectral functions with less artifacts on the scale of omega ~ T. We discuss the convergence properties with respect to the NRG calculation parameters (discretization, truncation cutoff) and the number of Matsubara points taken into account in the analytic continuation. We test the technique on the the single-impurity Anderson model and the Hubbard model (within the dynamical mean-field theory). For the Anderson impurity model, we discuss the shape of the Kondo resonance and its temperature dependence. For the Hubbard model, we discuss the inner structure of the Hubbard bands in metallic and insulating solutions at half-filling, as well as in the doped Mott insulator. Based on these test cases we conclude that the Pad\'e approximant approach provides more reliable results for spectral functions at low-frequency scales of omega < T and that it is capable of resolving sharp spectral features also at high frequencies. It outperforms broadening in most respects.