Effects of bond randomness and site dilution are systematically investigated for the Kitaev model describing a quantum spin liquid with fractional excitations of itinerant Majorana fermions and localized fluxes. We find that, in the high-temperature region where the itinerant Majorana fermions release their entropy, both types of disorders suppress the longitudinal thermal conductivity while keeping the specific heat almost unchanged. This suggests that both disorders reduce the mean-free path of the Majorana fermions. On the other hand, in the low-temperature region, the other specific heat peak associated with the entropy release from the localized fluxes is suppressed for both cases, but it is broadened and shifted to the lower-temperature side by the bond randomness, while the position and the width are almost unchanged against the site dilution. Contrasting behavior is also found in the thermal Hall effect under a magnetic field; the half quantization of the thermal Hall conductivity is fragile against the site dilution, while it remains for the bond randomness despite the reduced onset temperature. We discuss the contrasting behavior from the stability of the topological nature by calculating flux condensation and the Majorana excitation gap.
Effects of bond randomness and site dilution are systematically investigated for the Kitaev model describing a quantum spin liquid with fractional excitations of itinerant Majorana fermions and localized fluxes. We find that, in the high-temperature region where the itinerant Majorana fermions release their entropy, both types of disorders suppress the longitudinal thermal conductivity while keeping the specific heat almost unchanged. This suggests that both disorders reduce the mean-free path of the Majorana fermions. On the other hand, in the low-temperature region, the other specific heat peak associated with the entropy release from the localized fluxes is suppressed for both cases, but it is broadened and shifted to the lower-temperature side by the bond randomness, while the position and the width are almost unchanged against the site dilution. Contrasting behavior is also found in the thermal Hall effect under a magnetic field; the half quantization of the thermal Hall conductivity is fragile against the site dilution, while it remains for the bond randomness despite the reduced onset temperature. We discuss the contrasting behavior from the stability of the topological nature by calculating flux condensation and the Majorana excitation gap.
Thermal Fractionalization of Quantum Spins in a Kitaev Model: Coherent Transport of Majorana Fermions and T-linear Specific Heat
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.
We study how stable the Majorana-mediated spin transport in a quantum spin Kitaev model is against thermal fluctuations. Using the time-dependent thermal pure quantum state method, we examine finite-temperature spin dynamics in the Kitaev model. The model exhibits two characteristic temperatures, TL and TH, which correspond to energy scales of the local flux and the itinerant Majorana fermion, respectively. At low temperatures $(T\ensuremath{\ll}{T}_{L})$, an almost flux-free state is realized, and the spin excitation propagates in a similar way to that for the ground state. Namely, after the magnetic pulse is introduced at one of the edges, the itinerant Majorana fermions propagate the spin excitations even through the quantum spin liquid state region, and oscillations in the spin moment appear in the other edge with a tiny magnetic field. When $T\ensuremath{\sim}{T}_{L}$, larger oscillations in the spin moments are induced in the other edge, compared with the results at the ground state. At higher temperatures, excited Z2 fluxes disturb the coherent motion of the itinerant Majorana fermions, which suppresses the spin propagation. Our results demonstrate a crucial role of thermal fluctuations in the Majorana-mediated spin transport.
We study the 3D Kitaev and Kitaev-Heisenberg models, respectively, on the hyperhoneycomb and hyperoctagon lattices, both at zero and finite-temperature, in the thermodynamic limit. Our analysis relies on advanced tensor network (TN) simulations based on graph projected entangled-pair states (gPEPS). We map out the TN phase diagrams of the models and characterize their underlying gapped and gapless phases both at zero and finite temperature. In particular, we demonstrate how cooling down the hyperhoneycomb system from high temperature leads to fractionalization of spins to Majorana fermions and gauge fields that occurs in two separate temperature regimes, leaving their fingerprint on specific heat as a double-peak feature as well as on other quantities such as the thermal entropy, spin-spin correlations, and bond entropy. Using the Majorana representation of the Kitaev model, we further show that the low-temperature thermal transition to the Kitaev quantum spin liquid (QSL) phase is associated with the nontrivial Majorana band topology and the presence of Weyl nodes, which manifests itself via nonvanishing Chern number and finite thermal Hall conductivity. Beyond the pure Kitaev limit, we study the 3D Kitaev-Heisenberg (KH) model on the hyperoctagon lattice and extract the full phase diagram for different Heisenberg couplings. We further explore the thermodynamic properties of the magnetically-ordered regions in the KH model and show that, in contrast to the QSL phase, here the thermal phase transition follows the standard Landau symmetry-breaking theory.
The Kitaev model exhibits a canonical quantum spin liquid as a ground state and hosts two fractional quasiparticles, i.e., itinerant Majorana fermion and localized flux excitation. The former can carry heat and spin modulations in the quantum spin liquid, but the role of the latter remains unknown for the transport phenomena. Here, we focus on spin transport in the presence of excited fluxes and report that they yield strong interference in the propagation of the Majorana fermions, which have a gaugelike potential emergent around the fluxes. We examine the transient spin dynamics triggered by a pulsed magnetic field at an edge. In the absence of excited fluxes, the magnetic-field pulse creates the plane wave of the Majorana fermions, which flows in the quantum spin liquid. Although this wave does not accompany the change of local spin moments in bulk, it induces local moments at the side opposite to the edge under the magnetic-field pulse. We observe the spatial modulation of induced spin moments when fluxes are excited in the bulk region. This behavior is more striking than the case of lattice defects. Moreover, we find that although the amplitude of the spatial change is almost independent of the distance between lattice defects, it is strongly enhanced by increasing the distance for the case of excited fluxes. The difference is understood from the influence on the itinerant Majorana fermions; the lattice defects change the system locally, but flux excitations alter all the transfer integrals on the string connecting them. The present results will provide another route to observing intrinsic flux excitations distinguished from extrinsic effects such as lattice defects.
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.
We study the effects of the flux configurations on the emergent Majorana fermions in the S=21 Kitaev model on a honeycomb lattice, where quantum spins are fractionalized into itinerant Majorana fermions and localized fluxes. A quantum spin liquid appears as the ground state of the Kitaev model in the flux-free sector, which has intensively been investigated so far. In this flux sector, the Majorana fermion system has linear dispersions and shows power-law behavior in the Majorana correlations. On the other hand, periodically arranged flux configurations yield low-energy excitations in the Majorana fermion system, which are distinctly different from those in the flux-free state. We find that one of the periodically arranged flux states results in the gapped Majorana dispersion and the exponential decay in the Majorana correlations. The Kitaev system with another flux configuration exhibits a semi-Dirac-type dispersion, leading to the power-law decay with a smaller power than that in the flux-free sector along symmetry axes. We also examine the effect of the randomness in the flux configurations and clarify that the Majorana density of states is filled by increasing the flux density, and power-law decay in the Majorana correlations remains. The present results could be important to control the motion of Majorana fermions, which carries the spin excitations, in the Kitaev candidate materials.