“…Previously, this correlation function has been calculated using a variety of different methods, for example, full exact diagonalization (ED, [40]), Lanczos ED [45], CTQMC, and, very recently, the numerical renormalization group (NRG, [46,47]). All these methods have in common that the evaluation of f is computationally expensive, in fact, at the present stage only CTQMC solvers are applicable in multi-orbital settings [48,49].…”