The Andreev states of a superconducting quantum dot: mean field vs exact   numerical results
  We analyze the spectral density of a single level quantum dot coupled to superconducting leads focusing on the Andreev states appearing within the superconducting gap. We use two complementary approaches: the numerical renormalization group and the Hartree-Fock approximation. Our results show the existence of up to four bound states within the gap when the ground state is a spin doublet (\pi\ phase). Furthermore the results demonstrate the reliability of the mean field description within this phase. This is understood from a complete correspondence that can be established between the exact and the mean field quasiparticle excitation spectrum 
A. Mart\'in-Rodero and A. Levy Yeyati
cond-mat.mes-hall