The Kitaev model exhibits a quantum spin liquid hosting emergent fractionalized excitations. We study the Kitaev model on the honeycomb lattice coupled to a magnetic field along the [111] axis. Utilizing large-scale matrix product-based numerical models, we confirm three phases with transitions at different field strengths depending on the sign of the Kitaev exchange: a non-Abelian topological phase at low fields, an enigmatic intermediate regime only present for antiferromagnetic Kitaev exchange, and a field-polarized phase. For the topological phase, we numerically observe the expected cubic scaling of the gap and extract the quantum dimension of the non-Abelian anyons. Furthermore, we investigate dynamical signatures of the topological and the field-polarized phase using a matrix product operator-based time evolution method.
Quantum spin liquids (QSLs) host a variety of fractionalized particles. In Kitaev's paradigmatic honeycomb model a spin-fractionalizes intoflux due to emergentgauge field and matter Majorana fermions. Although these excitations have well-defined dynamics in the integrable limit, their direct experimental identification is notoriously challenging: realistic materials inevitably host additional symmetry-allowed interactions that break integrability and hybridize gauge and matter sectors, while magnetic fields, which are often required to suppress competing order and stabilize a putative QSL regime, further entangle the responses of different fractionalized quasiparticles and may even drive the system into field-induced spin-liquid phases that are not adiabatically connected to the integrable limit. A prominent example is the quantum Majorana metal, in which the distinct dynamics of fractionalized Majorana fermions can become directly visible in scattering. This review highlights recent progress on these related questions: in which field-stabilized QSL regimes and nearby emergent phases, and under what conditions, can the response of a specific fractionalized quasiparticle be isolated and positively understood, thereby clarifying the existence and the experimental scope of putative spin liquids? We review the progress on these questions across Abelian, non-Abelian, and emergent quantum phases under magnetic field that are not perturbatively connected to the integrable limit. We connect these field-induced dynamical phenomena to concrete experimental observables, relevant for neutron scattering, resonant inelastic x-ray scattering, and pump-probe spectroscopy, techniques that are capable of identifying the nature of different magnetic-field-induced phases and resolving specific types of fractionalized particles, including Majoranas andfluxes.
The Kitaev model on a honeycomb lattice with bond-dependent Ising interactions offers an exactly solvable model of a quantum spin liquid (QSL) with fractionalized excitations: gapped
Z2 fluxes and gapless linearly dispersing Majorana fermions in the isotropic limit (
Kx=Ky=Kz). We explore the phase diagram along two axes, an external magnetic field
h, applied out of plane of the honeycomb, and anisotropic interactions
Kz larger than the other two (${K}_{x}={K}_{y}\ensuremath{\equiv}K$). For ${K}_{z}/K\ensuremath{\gg}2$ and
h=0, the matter Majorana fermions have the largest gap, and the system is described by a gapped
Z2 toric code. One of the central questions we address is whether the fractionalized excitations in the different phases have sharp signatures that can be detected in experiments. We show that while the response to single-spin excitations is broad, the spectral function corresponding to two-spin excitations across a bond has sharp signatures that can be attributed to specific anyons. In the toric code regime, the $\ensuremath{\epsilon}=e\ifmmode\times\else\texttimes\fi{}m$ fermion, formed from the bosonic Ising electric (
e) and magnetic (
m) charges, disperses along a specific one-dimensional direction that provides a fingerprint of fractionalization. At lower
Kz in the center of the Abelian phase, in a regime we dub the primordial fractionalized regime, the field generates a hybridization between the $\ensuremath{\epsilon}$ fermion and the Majorana matter fermion, resulting in a $\ensuremath{\psi}$ fermion which too has a distinct quasi-one-dimensional dispersion. All the other phases in the field-anisotropy plane are naturally obtained from this primordial soup. These highly constrained fractonlike dispersions can be observable by inelastic light and neutron scattering, thereby providing ``smoking gun'' signatures of fractionalization in the QSL phase. Our analysis is based on calculations of susceptibilities, topological entanglement entropy, and excitation dynamics, obtained using exact diagonalization and density matrix renormalization group, and supported by perturbation theory.
The Kitaev model with an applied magnetic field in the [Formula: see text] direction shows two transitions: from a nonabelian gapped quantum spin liquid (QSL) to a gapless QSL at [Formula: see text] and a second transition at a higher field [Formula: see text] to a gapped partially polarized phase, where K is the strength of the Kitaev exchange interaction. We identify the intermediate phase to be a gapless U(1) QSL and determine the spin structure function [Formula: see text] and the Fermi surface [Formula: see text] of the gapless spinons using the density matrix renormalization group (DMRG) method for large honeycomb clusters. Further calculations of static spin-spin correlations, magnetization, spin susceptibility, and finite temperature-specific heat and entropy corroborate the gapped and gapless nature of the different field-dependent phases. In the intermediate phase, the spin-spin correlations decay as a power law with distance, indicative of a gapless phase.
It is now well established that the Kitaev honeycomb model in a magnetic field along the [111] direction harbors an intermediate gapless quantum spin liquid (QSL) phase sandwiched between a gapped non-Abelian QSL at low fields
H<Hc1 and a partially polarized phase at high fields
H>Hc2. Here, we analyze the low-field and high-field phases and phase transitions in terms of single- and two-magnon excitations using exact diagonalization and density matrix renormalization group methods. We find that the energy to create a bound state of two-magnon ${\mathrm{\ensuremath{\Delta}}}_{p}$ becomes lower than the energy to create a single spin flip ${\mathrm{\ensuremath{\Delta}}}_{s}$ near
Hc2. In the entire Kitaev spin liquid, ${\mathrm{\ensuremath{\Delta}}}_{p}<{\mathrm{\ensuremath{\Delta}}}_{s}$ and both gaps vanish at
Hc2. We make testable predictions for magnon pairing that could be observable in Raman scattering measurements on Kitaev QSL candidate materials.
The authors investigate dynamical signatures of fractionalization in the Kitaev honeycomb model under an out-of-plane magnetic field using advanced tensor network methods. By analyzing spin-spin and dimer-dimer correlations through the infinite Projected Entangled Pair States (iPEPS) simulations, they demonstrate distinct spectral features separating Majorana fermions from Z₂ flux excitations in the chiral spin liquid phase. Crucially, they identify the controversial intermediate phase as a gapless Majorana metal driven by strong flux fluctuations, characterized by logarithmic divergence in the low-energy density of states. This resolves a longstanding debate about the regime in question.