Yu-Shiba-Rusinov states of single magnetic molecule in an $s-wave$ superconductor

Abstract: We use the numerical renormalization group theory to investigate the Yu-Shiba-Rusinov (YSR) bound state properties of single magnetic molecules placed in an s-wave superconducting substrate. The molecule consist of a large core spin and a single orbital, coupled via an exchange interaction. The critical Coulomb interaction for the singlet/doublet transition decreases in the presence of this exchange interaction for both ferro and anti-ferromagnetic couplings. The number of YSR states also increase to two pairs, however, in the singlet phase, one of the pairs have zero spectral weight. We explore the evolution of the in-gap states using the Anderson model. Away from the particle-hole symmetry point, the results suggest a doublet-singlet-doublet transition as the on-site energy is lowered while keeping the Coulomb interaction fixed. To understand these results, we write down an effective model for the molecule in the limit of large superconducting order parameter. Qualitatively, it explains the various phase transitions and spectral nature of the in-gap states. Finally, we analyze the effects of magnetic anisotropic fields of the core spin on in-gap states. Due to internal degrees of freedom of the spin excited states, a multitude of new states emerges within the gap. Depending on the sign and strength of the uniaxial anisotropic field, the results indicate up to three pairs of YSR states.