Phonon dynamics in the site-disordered Kitaev spin liquid

The Kitaev honeycomb model provides a paradigmatic example of an exactly solvable quantum spin liquid (QSL), where spins fractionalize into itinerant Majorana fermions coupled to a static background of ${\mathbb{Z}}_{2}$ gauge fluxes. This model has attracted significant interest due to its potential experimental realization in spin-orbit Mott insulators like $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$. Among various experimental techniques, ultrasound measurements of sound attenuation have emerged as a promising approach to detect spin fractionalization in these materials. However, deviations from the ideal Kitaev model, often due to disorder, introduce localized modes that dominate the low-energy physics. To investigate these effects, we calculate the sound attenuation coefficient in the site-disordered Kitaev honeycomb model under an applied magnetic field that breaks time-reversal symmetry. Our analysis reveals that quasilocalized modes from quasivacancies impact the phonon self-energy, yet the sound attenuation coefficient maintains its sixfold symmetry and linear temperature dependence at low temperatures, even with disorder. This robustness underscores sound attenuation's reliability as a probe for spin fractionalization. Furthermore, we show that this linear behavior persists under an external field and in the random-flux sector, highlighting the persistent influence of fractionalized excitations, despite disorder.