Sum-rule Conserving Spectral Functions from the Numerical   Renormalization Group
  We show how spectral functions for quantum impurity models can be calculated very accurately using a complete set of ``discarded'' numerical renormalization group eigenstates, recently introduced by Anders and Schiller. The only approximation is to judiciously exploit energy scale separation. Our derivation avoids both the overcounting ambiguities and the single-shell approximation for the equilibrium density matrix prevalent in current methods, ensuring that relevant sum rules hold rigorously and spectral features at energies below the temperature can be described accurately. 
Andreas Weichselbaum and Jan von Delft