Yu-Shiba-Rusinov States in s-Wave Kagome Superconductors: Self-Consistent Bogoliubov-de Gennes Calculations

Significant research has recently been conducted into the Yu-Shiba-Rusinov (YSR) states in kagome superconductors through theoretical modeling and experimental investigations. However, additional efforts are still needed to further understand the local superconductivity near magnetic impurities in the kagome lattice and clarify how relevant quantities depend on the interaction strength, J, between such impurities and electrons. In this study, we explore a self-consistent numerical solution of the Bogoliubov-de Gennes equations for an s-wave superconducting Kagome model with a single classical magnetic impurity. Our study reveals that with increasing J, the local pair potential is systematically depressed in the vicinity of the impurity, similar to previous results obtained for the square and triangular lattices. Moreover, when J is further increased, the system undergoes a first-order phase transition with the appearance of stable and metastable states, reflecting the presence of the hysteresis loop in the pertinent quantities. As a consequence of this transition, the minimal energy of the stable YSR state is nonzero at any J, contrary to the expectations based on the assumption of a constant pair potential. A distinctive feature of the kagome lattice is that characteristics of the first-order transition are very sensitive to the position of the chemical potential within the kagome energy spectrum.