Generalized numerical renormalization group for dynamical quantities

In this paper we introduce a new approach for calculating dynamical properties within the numerical renormalization group. It is demonstrated that the method previously used fails for the Anderson impurity in a magnetic field due to the absence of energy scale separation. The problem is solved by evaluating the Green function with respect to the reduced density matrix of the full system, leading to accurate spectra in agreement with the static magnetization. The new procedure provides a unifying framework for calculating dynamics at any temperature and represents the correct extension of Wilson's original thermodynamic calculation.