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.
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.
There has been great interest in identifying a Kitaev quantum spin liquid state in frustrated magnets with bond-dependent interactions. In particular, the experimental report of a half-quantized thermal Hall conductivity in $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$ in the presence of a magnetic field has generated excitement as it could be strong evidence for a field-induced chiral spin liquid. More recent experiments, however, provide a conflicting interpretation advocating for topological magnons in the field-polarized state as the origin of the nonquantized thermal Hall conductivity observed in their experiments. An inherent difficulty in distinguishing between the two scenarios is that the phase transition between a putative two-dimensional spin liquid and the field-polarized state exists only at zero temperature, while the behavior at finite temperature is mostly crossover phenomena. In this work, we provide insights into the finite-temperature crossover behavior between the spin excitation continuum in a quantum spin liquid and topological magnons in the field-polarized state in three different theoretical models with large Kitaev interactions. These models allow for a field-induced phase transition from a spin liquid (or an intermediate field-induced spin liquid) to the field-polarized state in the quantum model. We obtain the dynamical spin structure factor as a function of magnetic field using molecular dynamics simulations and compute thermal Hall conductivity in the field-polarized regime. We demonstrate the gradual evolution of the dynamical spin structure factor exhibiting crossover behavior near magnetic fields where zero-temperature phase transitions occur in the quantum model. We also examine nonlinear effects on topological magnons and the validity of thermal Hall conductivity computed using linear spin wave theory. We discuss the implications of our results for existing and future experiments.
A Majorana fermion is a fermionic particle that is its own antiparticle. Since the theoretical discovery in 1937, the exotic particle has long been searched in particle physics. In the last few decades, however, it has attracted renewed interest in condensed matter physics, where it can be realized as an elementary excitation (quasiparticle) in quantum states of matter. In this review, we discuss another platform for Majorana fermions, the quantum spin liquid, in which interacting magnetic moments remain disordered down to the lowest temperature under strong quantum fluctuations. They are characterized by topological entanglement and fractional excitations, whose possible application to topological quantum computation is recently discussed intensively. As a prime candidate for such exotic states, we here focus on the Kitaev magnets, a subgroup of the spin-orbit Mott insulators. After a brief overview of the Kitaev model and the fractionalization of spins in the exact ground state, we review recent explosive development in this rapidly growing field, with a focus on numerical solutions of the Kitaev model at finite temperatures and the comparison with experiments. The key concept is thermal fractionalization --- two types of fractional excitations manifest themselves at largely different temperatures. This leads to distinct thermodynamics and spin dynamics in a variety of experimentally measurable quantities. We discuss such peculiar behaviors as the signatures of fractional quasiparticles, in careful comparison with the available experimental data for the candidate materials of the Kitaev magnets. Our review gives an overview of the current status of the identification of Majorana fermions in the Kitaev magnets, which would serve as a basis for further experimental and theoretical studies toward the manipulation of the exotic particles for topological quantum computation.
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.
We determine the phase diagrams of anisotropic Kitaev-Heisenberg models on the honeycomb lattice using parton mean-field theories based on different Majorana fermion representations of the S=1/2 spin operators. First, we use a two-dimensional Jordan-Wigner transformation (JWT) involving a semi-infinite snake string operator. To ensure that the fermionized Hamiltonian remains local, we consider the limit of extreme Ising exchange anisotropy in the Heisenberg sector. Second, we use the conventional Kitaev representation in terms of four Majorana fermions subject to local constraints, which we enforce through Lagrange multipliers. For both representations, we self-consistently decouple the interaction terms in the bond and magnetization channels and determine the phase diagrams as a function of the anisotropy of the Kitaev couplings and the relative strength of the Ising exchange. While both mean-field theories produce identical phase boundaries for the topological phase transition between the gapless and gapped Kitaev quantum spin liquids, the JWT fails to correctly describe the magnetic instability and finite-temperature behavior. Our results show that the magnetic phase transition is first order at low temperatures but becomes continuous above a certain temperature. At this energy scale, we also observe a finite-temperature crossover on the quantum-spin-liquid side, from a fractionalized paramagnet at low temperatures, in which gapped flux excitations are frozen out, to a conventional paramagnet at high temperatures. Published by the American Physical Society 2024