Candidate materials for the Kitaev spin liquid generically have residual interactions beyond the Kitaev coupling. It therefore becomes necessary to understand how signatures of the quantum spin liquid, e.g., fractionalization of the spin excitations, are affected by the presence of these interactions. Recently it was shown that phonon dynamics is an indirect but effective probe to study fractionalized excitations in the Kitaev spin liquid. Ultrasound experiments can measure sound attenuation, which should show characteristic temperature and angular dependence of the sound-attenuation coefficient if the scattering of phonons happens predominantly on Majorana fermions. So far the computation of the sound attenuation was only done in the pure spin-phonon coupled Kitaev model, without taking into account residual interactions. In order to understand experimental signatures, here we present a mean-field study of the sound attenuation in the generalized $J\text{\ensuremath{-}}K\text{\ensuremath{-}}\mathrm{\ensuremath{\Gamma}}$ model, which is relevant to candidate materials. Our findings show that as long as the system is in the spin liquid phase, characteristic features of the sound attenuation remain observable even in the presence of residual interactions.
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
Z2 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.
Motivated by the appearance of Majorana fermions in a broad range of correlated and topological electronic systems, we develop a general method to compute the dynamical response of interacting Majorana fermions in the random-phase approximation (RPA). This can be applied self-consistently on top of Majorana mean-field theory backgrounds, thereby in particular providing a powerful tool to analyze behavior in the vicinity of (various heavily studied) exactly soluble models. Prime examples are quantum spin liquids (QSL) with emergent Majorana excitations, with the celebrated exact solution of Kitaev. We employ the RPA to study in considerable detail phase structure and dynamics of the extended Kitaev honeycomb K J Γ model, with and without an applied field. First, we benchmark our method with Kitaev's exactly soluble model, finding a remarkable agreement. The interactions between Majorana fermions even turn out to mimic the effect of local Z 2 flux excitations, which we explain analytically. Second, we show how small non-Kitaev couplings J and Γ induce Majorana bound states, resulting in sharp features in the dynamical structure factor in the presence of fractionalization: such “spinon excitons” naturally appear, and can coexist and interact with the broad Majorana continuum. Third, for increasing couplings or field, our theory predicts instabilities of the Kitaev QSL (KQSL) triggered by the condensation of the sharp modes. From the high-symmetry momenta of the condensation we can deduce which magnetically ordered phases surround the KQSL, in good agreement with previous finite-size numerics. We discuss implications for experiments and the broad range of applicability of our method to other QSL and Majorana systems.
The possible existence of a quantum spin liquid (QSL) phase, an exotic state of matter with long-range quantum entanglement and fractionalized excitations, in
α-RuCl
3 has sparked widespread interests in exploring QSLs in various Kitaev models under magnetic fields. Recently, a
K-
J-
Γ-
Γ′ model has been proposed to accurately describe the compound by fitting the measured thermodynamic data, where
J is Heisenberg interaction, and
Γ,
Γ′ are off-diagonal exchanges on top of the dominant Kitaev coupling
K. Based on this effective model, an intermediate QSL phase in presence of an out-of-plane fields along the
[111] direction, between the low-field zigzag order and high-field polarized phase, has been predicted. By combining density matrix renormalization group (DMRG), exponential tensor renormalization group (XTRG), and variational Monte Carlo (VMC) calculations, we address the nature of this QSL phase in the honeycomb
K-
J-
Γ-
Γ′ model under the
[111]-direction field. Our DMRG calculations find the algebraic-like decay of spin correlation function, the finite spin scalar chiral order, and lattice nematic order. Together with the XTRG results of power-law specific heat at low temperature, our findings naturally suggest a gapless nematic chiral spin liquid. On the other hand, our VMC study finds a gapped nematic chiral spin liquid with the variational energy very close to that obtained in DMRG. As VMC finds a very small gap that is beyond the current resolution of both ground-state DMRG and finite-temperature XTRG calculations on finite-width cylinders, we resort the full clarification for the nature of the intermediate QSL to future studies. Lastly, we discuss the implications of our results to the recent experiment on
α-RuCl
3 and the QSL-like phase in a similar
K-
Γ-
Γ′ model.
Quantum spin liquids (QSLs) host a variety of fractionalized particles. In Kitaev's paradigmatic honeycomb model a spin-fractionalizes intoflux due to emergentgauge field and matter Majorana fermions. Although these excitations have well-defined dynamics in the integrable limit, their direct experimental identification is notoriously challenging: realistic materials inevitably host additional symmetry-allowed interactions that break integrability and hybridize gauge and matter sectors, while magnetic fields, which are often required to suppress competing order and stabilize a putative QSL regime, further entangle the responses of different fractionalized quasiparticles and may even drive the system into field-induced spin-liquid phases that are not adiabatically connected to the integrable limit. A prominent example is the quantum Majorana metal, in which the distinct dynamics of fractionalized Majorana fermions can become directly visible in scattering. This review highlights recent progress on these related questions: in which field-stabilized QSL regimes and nearby emergent phases, and under what conditions, can the response of a specific fractionalized quasiparticle be isolated and positively understood, thereby clarifying the existence and the experimental scope of putative spin liquids? We review the progress on these questions across Abelian, non-Abelian, and emergent quantum phases under magnetic field that are not perturbatively connected to the integrable limit. We connect these field-induced dynamical phenomena to concrete experimental observables, relevant for neutron scattering, resonant inelastic x-ray scattering, and pump-probe spectroscopy, techniques that are capable of identifying the nature of different magnetic-field-induced phases and resolving specific types of fractionalized particles, including Majoranas andfluxes.
Experimental identification of quantum spin liquids remains a challenge, as the pristine nature is to be seen in asymptotically low temperatures. We here theoretically show that the precursor of quantum spin liquids appears in the spin dynamics in the paramagnetic state over a wide temperature range. Using the cluster dynamical mean-field theory and the continuous-time quantum Monte Carlo method, which are newly developed in the Majorana fermion representation, we calculate the dynamical spin structure factor, relaxation rate in nuclear magnetic resonance, and magnetic susceptibility for the honeycomb Kitaev model whose ground state is a canonical example of the quantum spin liquid. We find that dynamical spin correlations show peculiar temperature and frequency dependence even below the temperature where static correlations saturate. The results provide the experimentally accessible symptoms of the fluctuating fractionalized spins evincing the quantum spin liquids.