Quantum spin impurities coupled to superconductors are under intense investigation for their relevance to fundamental research as well as the prospects to engineer novel quantum phases of matter. Here we develop a large-NN mean-field theory of a strongly coupled spin-\tfrac{1}{2}12 quantum impurity in a conventional ss-wave superconductor. The approach is benchmarked against Wilson’s numerical renormalization group (NRG). While the large-NN method is not applicable in the weak-coupling regime where the Kondo temperature T_KTK is smaller than the superconducting gap \DeltaΔ, it performs very well in the strong coupling regime where T_K \gtrsim \DeltaTK≳Δ, thus allowing us to obtain a reasonably accurate description of experimentally relevant quantities. The latter include the energy of the Yu-Shiba-Rusinov subgap states, their spectral weight, as well as the local density of continuum states. The method provides a reliable analytical tool that complements other perturbative and non-perturbative methods, and can be extended to more complex impurity models for which NRG may not be easily applicable.
Read full abstract