We identify fingerprints of a proximate quantum spin liquid (QSL), observable by finite-temperature dynamical thermal transport within a minimal version of the idealized Kitaev model on a two-leg ladder if subjected to inevitably present Heisenberg couplings. Using exact diagonalization and quantum typicality, we uncover (i) an insulator-conductor crossover induced by recombination of fractionalized excitations at small Heisenberg couplings, (ii) low- and high-energy signatures of fractionalized excitations, which survive far off the pure QSL point, and (iii) a nonmonotonous current lifetime versus Heisenberg couplings. Guided by perturbation theory, we find (iv) a Kitaev-exchange-induced ``one-magnon'' contribution to the dynamical heat transport in the strong Heisenberg rung limit.
We identify fingerprints of a proximate quantum spin liquid (QSL), observable by finite-temperature dynamical thermal transport within a minimal version of the idealized Kitaev model on a two-leg ladder if subjected to inevitably present Heisenberg couplings. Using exact diagonalization and quantum typicality, we uncover (i) an insulator-conductor crossover induced by recombination of fractionalized excitations at small Heisenberg couplings, (ii) low- and high-energy signatures of fractionalized excitations, which survive far off the pure QSL point, and (iii) a nonmonotonous current lifetime versus Heisenberg couplings. Guided by perturbation theory, we find (iv) a Kitaev-exchange-induced ``one-magnon'' contribution to the dynamical heat transport in the strong Heisenberg rung limit.
We study the spin transport through the quantum spin liquid (QSL) by investigating the real-time and real-space dynamics of the Kitaev spin system with zigzag edges using the time-dependent Majorana mean-field theory. After the magnetic-field pulse is introduced to one of the edges, spin moments are excited in the opposite edge region although spin moments are never induced in the Kitaev QSL region. This unusual spin transport originates from the fact that the S=1/2 spins are fractionalized into the itinerant and localized Majorana fermions in the Kitaev system. Although both Majorana fermions are excited by the magnetic pulse, only the itinerant ones flow through the bulk regime without spin excitations, resulting in the spin transport in the Kitaev system despite the presence of a nonzero spin gap. We also demonstrate that this phenomenon can be observed in the system with small Heisenberg interactions using the exact diagonalization.
The Kitaev model offers a platform for quantum spin liquids (QSLs) with fractional excitations, itinerant Majorana fermions and localized fluxes. Since these fractional excitations could be utilized for quantum computing, how to create, observe, and control them through the spin degree of freedom is a central issue. Here, we study dynamical spin transport in a wide range of frequencies for the Kitaev-Heisenberg model, by applying an AC magnetic field to an edge of the system. We find that, in the Kitaev QSL phase, spin polarizations at the other edge are resonantly induced in a specific spin component, even though the static spin correlations are vanishingly small. This interedge spin resonance appears around the input frequency over the broad frequency range. Comparing with the dynamical spin correlations, we clarify that the resonance is governed by the itinerant Majorana fermions with a broad continuum excitation spectrum, which can propagate over long distances, although it vanishes for the pure Kitaev model because of accidental degeneracy and requires weak Heisenberg interactions. We also find that the spin polarizations in the other spin components are weakly induced at an almost constant frequency close to the excitation gap of the localized fluxes, irrespective of the input frequency. These results demonstrate that the dynamical spin transport is a powerful probe of the fractional excitations in the Kitaev QSL. Possible experimental realization of the interedge spin resonance is discussed.
Quantum spin liquids are long-range entangled states of matter with emergent gauge fields and fractionalized excitations. While candidate materials, such as the Kitaev honeycomb ruthenate $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$, show magnetic order at low temperatures T, here we demonstrate numerically a dynamical crossover from magnonlike behavior at low T and frequencies $\ensuremath{\omega}$ to long-lived fractionalized fermionic quasiparticles at higher T and $\ensuremath{\omega}$. This crossover is akin to the presence of spinon continua in quasi-1D spin chains. It is further shown to go hand in hand with persistent typicality down to very low T. This aspect, which has also been observed in the spin-1/2 kagome Heisenberg antiferromagnet, is a signature of proximate spin liquidity and emergent gauge degrees of freedom more generally, and can be the basis for the numerical study of many finite-T properties of putative spin liquids.
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.
We study the Heisenberg-Kitaev spin chain in order to uncover the interplay between two qualitatively different integrable points in the physics of heat transport in one dimension. Focusing on high temperatures and using analytical as well as numerical approaches within linear response theory, we explore several directions in parameter space including exchange-coupling ratios, anisotropies, and external magnetic fields. We show the emergence of purely ballistic energy transport at all integrable points, manifest in pronounced Drude weights and low-frequency suppression of regular-conductivity contributions. Moreover, off integrability, we find extended quantum chaotic regions with vanishing Drude weights and well-defined dc conductivities. In the vicinity of the Kitaev point, we observe clear signatures of the topological gap in the response function. This gap coexists with a nonzero Drude weight in the Kitaev chain.
We study the ground state and low-lying excited states of the Kitaev-Heisenberg model on a two-leg ladder geometry using the density-matrix renormalization-group and Lanczos exact diagonalization methods. The Kitaev and Heisenberg interactions are parametrized as $K=sin\ensuremath{\phi}$ and $J=cos\ensuremath{\phi}$ with an angle parameter $\ensuremath{\phi}$. Based on the results for several types of order parameters, excitation gaps, string order parameter, and entanglement spectra, the $\ensuremath{\phi}$-dependent ground state phase diagram is determined. Remarkably, the phase diagram is quite similar to that of the Kitaev-Heisenberg model on a honeycomb lattice, exhibiting the same long-range-ordered states, namely rung singlet (analog to N\'eel in three dimensions), zigzag, ferromagnetic, and stripy, and the presence of gapped spin liquids around the exactly solvable Kitaev points $\ensuremath{\phi}=\ifmmode\pm\else\textpm\fi{}\ensuremath{\pi}/2$. We also calculate the expectation value of a plaquette operator corresponding to a $\ensuremath{\pi}$-flux state in order to establish how the gapped Kitaev spin liquid extends away from the $\ensuremath{\phi}=\ifmmode\pm\else\textpm\fi{}\ensuremath{\pi}/2$. Furthermore, we determine the dynamical spin structure factor and discuss the effect of the Kitaev interaction on the spin-triplet dispersion.