Temperature evolution of the spin correlation and excitation spectrum is investigated for the Kitaev model defined on a decorated honeycomb lattice by using the quantum Monte Carlo simulation in the Majorana fermion representation. The ground state of this quantum spin model is given by two kinds of chiral spin liquids: one is topologically trivial with Abelian anyon excitations, and the other is topologically nontrivial with non-Abelian anyon excitations. While lowering temperature, the model exhibits several crossovers in the paramagnetic state, which originate from the fractionalization of quantum spins into Majorana fermions, in addition to a phase transition associated with time reversal symmetry breaking. We show that the spin correlation develops around the crossover temperatures, whereas it shows a slight change at the critical temperature, as in other Kitaev-type models. We also calculate the excitation spectrum in terms of Majorana fermions, and find that the excitation gap in the non-Abelian phase is fragile against thermal fluctuations of the Z 2 fluxes, while that in the Abelian phase remains open.
Thermodynamic properties of chiral spin liquids are investigated for a variant of the Kitaev model defined on a decorated honeycomb lattice. Using the quantum Monte Carlo simulation, we find that the model exhibits a finite-temperature phase transition associated with the time reversal symmetry breaking, in both topologically trivial and nontrivial regions. Numerical results for the Chern number and the thermal Hall conductivity indicate that the phase transition changes from a continuous to a discontinuous transition as we vary the coupling constants to reach the non-Abelian phase coming from the Abelian phase of the model. In addition, we find as a diagnostic of the chiral spin liquids, successive crossovers with multistage entropy release above the critical temperature, which indicates that the hierarchical fractionalization of a quantum spin occurs differently between the two regions.
A three-dimensional Kitaev model on a hyperhoneycomb lattice is investigated numerically at finite temperature. The Kitaev model is one of the solvable quantum spin models, where the ground state is given by gapped and gapless spin liquids, depending on the anisotropy of the interactions. This model can be rewritten as a free Majorana fermion system coupled with Z 2 variables. The density of states of Majorana fermions shows an excitation gap in the gapped region, while it is semimetallic in the gapless region reflecting the Dirac node. Performing the Monte Carlo simulation, we calculate the temperature dependence of the Majorana spectra. We find that the semimetallic dip is filled as temperature increases in the gapless region, but surprisingly, the spectrum develops an excitation gap in the region near the gapless-gapped boundary. Such changes of the low-energy spectrum appear sharply at the transition temperature from the spin liquid to the paramagnetic state. The results indicate that thermal fluctuations of the Z 2 fields significantly influence the low-energy state of Majorana fermions, especially in the spin liquid formation.
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.
Finite-temperature (
T) properties of a Kitaev model defined on a honeycomb lattice are investigated by a quantum Monte Carlo simulation, from the viewpoint of fractionalization of quantum
S=1/2 spins into two types of Majorana fermions, itinerant and localized. In this system, the entropy is released successively at two well-separated
T scales, as a clear indication of the thermal fractionalization. We show that the high-
T crossover, which is driven by itinerant Majorana fermions, is closely related with the development of nearest-neighbor spin correlations. On the other hand, the low-
T crossover originates in thermal fluctuations of fluxes composed of localized Majorana fermions, by which the spectrum of itinerant Majorana fermions is significantly disturbed. As a consequence, in the intermediate-
T range, the system exhibits coherent transport of Majorana fermions and
T-linear behavior in the specific heat, despite the Dirac semimetallic spectrum in the low-
T limit. We also show that the flux fluctuations affect the gapless-gapped phase boundary. Our results indicate that the fractionalization is experimentally observable in spin correlations, specific heat, and transport properties.
An effective model in the isolated dimer limit of the Kitaev-type model on a decorated honeycomb lattice is investigated at finite temperature. The ground state of this model is exactly shown to be a chiral spin liquid with spontaneous breaking of time reversal symmetry. We elab- orate the finite-temperature phase diagram by using the mean-field approximation and Monte Carlo simulation. We find that the phase transition between the high-temperature paramagnetic phase and the low-temperature chiral spin liquid phase is always of second order in the Monte Carlo results, although a tricritical point appears in the mean-field phase diagram. The finite-size scaling analysis of the Monte Carlo data suggests that the phase transition belongs to the two-dimensional Ising universality class.