Prethermalization phenomena in driven systems are generally understood via a local Floquet Hamiltonian obtained from a high-frequency expansion. Remarkably, recently it has been shown that a driven Kitaev spin liquid with fractionalized excitations can realize a quasistationary state that is not captured by this paradigm. Instead distinct types of fractionalized excitations are characterized by vastly different temperatures—a phenomenon dubbed . In our work, we analyze fractionalized prethermalization in a driven one-dimensional Hubbard model at strong coupling which hosts spin-charge fractionalization. At intermediate frequencies quasisteady states emerge which are characterized by a low spin and high charge temperature with lifetimes set by two competing processes: the lifetime of the quasiparticles determined by Fermi's Golden rule and the exponentially long lifetime of a Floquet prethermal plateau. We classify drives into three categories, each giving rise to distinct (fractional) prethermalization dynamics. Resorting to a time-dependent variant of the Schrieffer-Wolff transformation, we systematically analyze how these drive categories are linked to the underlying driven Hubbard model, thereby providing a general understanding of the emergent thermalization dynamics. We discuss routes toward an experimental realization of this phenomenon in quantum simulation platforms.
The long-sought quantum spin liquid is a quantum-entangled magnetic state leading to the fractionalization of spin degrees of freedom. Quasiparticles emergent from the fractionalization affect not only the ground-state properties but also thermodynamic behavior in a peculiar manner. Here we investigate how the spin dynamics evolves from the high-temperature paramagnet to the quantum spin-liquid ground state, for the Kitaev spin model describing the fractionalization into itinerant matter fermions and localized
Z2 fluxes. Beyond the previous studies [J. Yoshitake, J. Nasu, and Y. Motome, Phys. Rev. Lett. 117, 157203 (2016); J. Yoshitake, J. Nasu, Y. Kato, and Y. Motome, Phys. Rev. B 96, 024438 (2017)], in which the mean-field nature of the cluster dynamical mean-field theory prevented us from studying low-temperature properties, we develop a numerical technique by applying the continuous-time quantum Monte Carlo (CTQMC) method to statistical samples generated by the quantum Monte Carlo (QMC) method in a Majorana fermion representation. This QMC+CTQMC method is fully unbiased and enables us to investigate the low-temperature spin dynamics dominated by thermally excited
Z2 fluxes, including the unconventional phase transition caused by
Z2 flux loops in three dimensions, which was unreachable by the previous methods. We apply this technique to the Kitaev model in both two and three dimensions. Our results clearly distinguish two cases: while the dynamics changes smoothly through the crossover in the two-dimensional honeycomb case, it exhibits singular behaviors at the phase transition in the three-dimensional hyperhoneycomb case. We show that the low-temperature spin dynamics is a sensitive probe for thermally fluctuating
Z2 fluxes that behave very differently between two and three dimensions.
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.
Quantum spin liquid is a disordered but highly entangled magnetic state with fractional spin excitations1. The ground state of an exactly solved Kitaev honeycomb model is perhaps its clearest example2. Under a magnetic field, a spin flip in this model fractionalizes into two types of anyon, a quasiparticle with more complex exchange statistics than standard fermions or bosons: a pair of gauge fluxes and a Majorana fermion2,3. Here, we demonstrate this kind of fractionalization in the Kitaev paramagnetic state of the honeycomb magnet α-RuCl3. The spin excitation gap determined by nuclear magnetic resonance consists of the predicted Majorana fermion contribution following the cube of the applied magnetic field2,4,5, and a finite zero-field contribution matching the predicted size of the gauge flux gap2,6. The observed fractionalization into gapped anyons survives in a broad range of temperatures and magnetic fields, which establishes α-RuCl3 as a unique platform for future investigations of anyons.
A Kitaev quantum spin liquid involves fractionalized quasiparticles called itinerant Majorana fermions and localized gauge fluxes. The topological character with chiral excitations under a magnetic field is currently in an active debate. For a candidate spin liquid $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$, we report low-energy spin excitations through nuclear magnetic and quadrupole resonance measurements down to temperatures three orders of magnitude lower than the exchange energy. In a field-driven disordered state, we found the evolution of two kinds of gapped spin excitations with distinct field dependences. The result is indicative of a field-induced hybridization between fractionalized quasiparticles, as inferred by the perturbation theory for the Kitaev model. The detailed field dependence of spin dynamics across the quantum critical region shows the emergence of a gapless intermediate phase between the gapped zigzag order phase and spin-disordered phases.
A prominent feature of quantum spin liquids is fractionalization of the spin degree of freedom. Fractionalized excitations have their own dynamics in different energy scales, and hence, affect finite-temperature (
T) properties in a peculiar manner even in the paramagnetic state harboring the quantum spin liquid state. We here present a comprehensive theoretical study of the spin dynamics in a wide
T range for the Kitaev model on a honeycomb lattice, whose ground state is a quantum spin liquid. In this model, the fractionalization occurs to break up quantum spins into itinerant matter fermions and localized
Z2 fluxes, which results in two crossovers at very different
T scales. Extending the previous study for the isotropic coupling case [J. Yoshitake, J. Nasu, and Y. Motome, Phys. Rev. Lett. 117, 157203 (2016)], we calculate the dynamical spin structure factor $S(\mathbf{q},\ensuremath{\omega})$, the NMR relaxation rate
1/T1, and the magnetic susceptibility $\ensuremath{\chi}$ while changing the anisotropy in the exchange coupling constants, by using the dynamical mean-field theory based on a Majorana fermion representation. We describe the details of the methodology including the continuous-time quantum Monte Carlo method for computing dynamical spin correlations and the maximum entropy method for analytic continuation. We confirm that the combined method provides accurate results in a wide
T range including the region where the spins are fractionalized. We find that also in the anisotropic cases the system exhibits peculiar behaviors below the high-
T crossover whose temperature is comparable to the average of the exchange constants: $S(\mathbf{q},\ensuremath{\omega})$ shows an inelastic response at the energy scale of the averaged exchange constant,
1/T1 continues to grow even though the equal-time spin correlations are saturated and almost
T independent, and $\ensuremath{\chi}$ deviates from the Curie-Weiss behavior. In particular, when the exchange interaction in one direction is stronger than the other two, the dynamical quantities exhibit qualitatively different
T dependences from the isotropic case at low
T, reflecting the opposite parity between the flux-free ground state and the flux-excited state, and a larger energy cost for flipping a spin in the strong interaction direction. On the other hand, when the exchange anisotropy is in the opposite way, the results are qualitatively similar to those in the isotropic case. All these behaviors manifest the spin fractionalization in the paramagnetic region. Among them, the dichotomy between the static and dynamical spin correlations is unusual behavior hardly seen in conventional magnets. We discuss the relation between the dichotomy and the spatial configuration of the
Z2 gauge fields. Our results could stimulate further experimental and theoretical analyses of candidate materials for the Kitaev quantum spin liquids.