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.
In the field of frustrated magnetism, Kitaev models provide prototypical examples of exactly solvable quantum spin liquids, in which the spin degrees of freedom fractionalize into Majorana fermions coupled to an emergent Z2 gauge field. While the ground states of these models are already well understood as exhibiting different kinds of Majorana (semi)metals, it has still been an open task to achieve a detailed understanding of the thermodynamics, and, namely, the low-temperature ordering of the gauge field into different Z2 flux configurations. In this thesis, we investigate the thermodynamics of Kitaev systems in two and three spatial dimensions with sign-problem-free quantum Monte Carlo simulations. In a first study of elementary 3D Kitaev models, we verify that the ground state Z2 flux sectors of these systems are entirely determined by their elementary plaquette length – a result which shows the validity of Lieb’s theorem for lattice geometries which lack the geometric requirements for its proof. We closely investigate the low-temperature phase transition associated with gauge-ordering, which is a particular realization of an inverted Ising phase transition, as it occurs in general lattice Z2 gauge theories. Our results corroborate the understanding of this transition as separating different regimes in terms of vison-loop excitations, by showing a clear correlation between the critical temperature and the vison gap. We also introduce the concept of “gauge frustration”, which, for a particular 3D Kitaev model, leads to the suppression of the phase transition and a more complex physical behavior at low temperatures. A second quantum Monte Carlo study focuses on a generalized 2D Kitaev model on a five-coordinated lattice system. Here, the gauge ordering is accompanied by the spontaneous breaking of time-reversal symmetry, a scenario which not only allows for a variety of topological ground states, but also for the occurrence of a phase transition in two spatial dimensions. We show that this model exhibits such a transition at a particularly high temperature scale. Moreover, the system possesses a number of phases where the Z2 fluxes are only partially ordered.
Kitaev quantum spin liquids host Majorana fermions via the fractionalization of spins. In a magnetic field, the Majorana fermions were predicted to comprise a topological state, which has attracted great attention by the discovery of the half-quantized thermal Hall conductivity. Nevertheless, a reliable theory remains elusive for the field effect, especially at finite temperature. Here we present unbiased large-scale numerical results for the Kitaev model in a wide range of magnetic field and temperature. We find that the unconventional paramagnetic region showing fractional spin dynamics extends at finite temperature, far beyond the field range where the topological state is expected at zero temperature. Our results show the confinement-deconfinement behavior between the fractional Majorana excitations and the conventional magnons.
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.
The quantum spin liquid is an enigmatic quantum state in insulating magnets, in which conventional long-range order is suppressed by strong quantum fluctuations. Recently, an unconventional phase transition was reported between the low-temperature quantum spin liquid and the high-temperature paramagnet in the Kitaev model on a three-dimensional hyperhoneycomb lattice. Here, we show that a similar ``liquid-gas'' transition takes place in another three-dimensional lattice, the hyperoctagon lattice. We investigate the critical phenomena by adopting the Green-function based Monte Carlo technique with the kernel polynomial method, which enables systematic analysis of up to 2048 sites. The critical temperature is lower than that in the hyperhoneycomb case, reflecting the smaller flux gap. We also discuss the transition on the basis of an effective model in the anisotropic limit.