Z2 flux binding to higher-spin impurities in the Kitaev spin liquid

Stabilizing Z2 fluxes in Kitaev spin liquids (KSLs) is crucial for both characterizing candidate materials and identifying Ising anyons. In this study, we investigate the effects of spin-S magnetic impurities embedded in the spin-1/2 KSL. Utilizing exact diagonalization and density matrix renormalization group methods, we examine the impurity magnetization and ground-state flux sector with varying impurity coupling and spin size. Our findings reveal that impurity magnetization exhibits an integer/half-integer spin dependence, which aligns with analytical predictions, and a flux-sector transition from bound-flux to zero-flux occurs at low coupling strengths, independent of the impurity spin. Notably, for spin-3/2 impurities, we observe a reentrant bound-flux sector, which remains stable under magnetic fields. By considering fermionic representations of our spin Hamiltonian, we provide phenomenological explanations for the transitions. Our results suggest a novel way of binding a flux in KSLs, beyond the proposals of vacancies or Kondo impurities.