Emergent gauge theories take a prominent role in the description of quantum matter, supporting deconfined phases with topological order and fractionalized excitations. A common construction of \mathbb{Z}_2 ℤ 2 lattice gauge theories, first introduced by Wegner, involves Ising gauge spins placed on links and subject to a discrete \mathbb{Z}_2 ℤ 2 Gauss law constraint. As shown by Kitaev, \mathbb{Z}_2 ℤ 2 lattice gauge theories also emerge in the exact solution of certain spin systems with bond-dependent interactions. In this context, the \mathbb{Z}_2 ℤ 2 gauge field is constructed from Majorana fermions, with gauge constraints given by the parity of Majorana fermions on each site. In this work, we provide an explicit Jordan-Wigner transformation that maps between these two formulations on the square lattice, where the Kitaev-type gauge theory emerges as the exact solution of a spin-orbital (Kugel-Khomskii) Hamiltonian. We then apply our mapping to study local perturbations to the spin-orbital Hamiltonian, which correspond to anisotropic interactions between electric-field variables in the \mathbb{Z}_2 ℤ 2 gauge theory. These are shown to induce anisotropic confinement that is characterized by emergence of weakly-coupled one-dimensional spin chains. We study the nature of these phases and corresponding confinement transitions in both absence and presence of itinerant fermionic matter degrees of freedom. Finally, we discuss how our mapping can be applied to the Kitaev spin-1/2 model on the honeycomb lattice.
While phase transitions between magnetic analogs of three states of matter -- a long-range ordered state, paramagnet, and spin liquid -- have been extensively studied, the possibility of liquid-liquid transitions, namely, between different spin liquids, remains elusive. By introducing the Ising coupling into the honeycomb Kitaev model with bond asymmetry, we discover that the Kitaev spin liquid becomes a bond-nematic quantum paramagnet before magnetically ordered. The phase transition between the two liquid-like states with different topological nature is of first order, driven by delocalization of the
Z2 gauge fluxes, and persists to a critical point at finite temperature located inside the regime where quantum spins are fractionalized. It is suggested that similar transitions may occur in other perturbed Kitaev magnets with bond asymmetry.
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.
We investigate the phase diagram of a bilayer Kitaev honeycomb model with Ising interlayer interactions, deriving effective models via perturbation theory and performing Majorana mean-field theory calculations. We show that a diverse array of magnetic and topological phase transitions occur, depending on the direction of the interlayer Ising interaction and the relative sign of Kitaev interactions. When two layers have the same sign of the Kitaev interaction, a first-order transition from a Kitaev spin liquid to a magnetically ordered state takes place. The magnetic order points along the Ising axis and it is (anti)ferromagnetic for (anti)ferromagnetic Kitaev interactions. However, when two layers have opposite signs of the Kitaev interaction, we observe a notable weakening of magnetic ordering tendencies and the Kitaev spin liquid survives up to a remarkably larger interlayer exchange. Our mean-field analysis suggests the emergence of an intermediate gapped
Z2 spin-liquid state, which eventually becomes unstable upon vison condensation. The confined phase is described by a highly frustrated ${120}^{\ensuremath{\circ}}$ compass model. We furthermore use perturbation theory to study the model with the Ising axis pointing along the $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{z}$ axis or lying in the
xy plane. In both cases, our analysis reveals the formation of one-dimensional Ising chains, which remain decoupled in perturbation theory, resulting in a subextensive ground-state degeneracy. Our results highlight the interplay between topological order and magnetic ordering tendencies in bilayer quantum spin liquids.
Kitaev's honeycomb-lattice spin-
1/2 model has become a paradigmatic example for
Z2 quantum spin liquids, both gapped and gapless. Here we study the fate of these spin-liquid phases in differently stacked bilayer versions of the Kitaev model. Increasing the ratio between the interlayer Heisenberg coupling ${J}_{\ensuremath{\perp}}$ and the intralayer Kitaev couplings
Kx,y,z destroys the topological spin liquid in favor of a paramagnetic dimer phase. We study phase diagrams as a function of ${J}_{\ensuremath{\perp}}/K$ and Kitaev coupling anisotropies using Majorana-fermion mean-field theory, and we employ different expansion techniques in the limits of small and large ${J}_{\ensuremath{\perp}}/K$. For strongly anisotropic Kitaev couplings, we derive effective models for the different layer stackings that we use to discuss the quantum phase transition out of the Kitaev phase. We find that the phase diagrams depend sensitively on the nature of the stacking and anisotropy strength. While in some stackings and at strong anisotropies we find a single transition between the Kitaev and dimer phases, other stackings are more involved. Most importantly, we prove the existence of two novel macrospin phases, which can be understood in terms of Ising chains that can be either coupled ferromagnetically or remain degenerate, thus realizing a classical spin liquid. In addition, our results suggest the existence of a flux phase with spontaneous interlayer coherence. We discuss prospects for experimental realizations.
Kitaev quantum spin liquids (KQSLs) host fractionalized excitations described by itinerant Majorana quasiparticles and gapped
Z2 fluxes (visons), providing a platform for emergent topological matter. Whether such a state survives under strong magnetic fields, however, remains an open question. The layered honeycomb magnet
α-RuCl
3 is a leading candidate material: an in-plane field of
∼ 7 T suppresses antiferromagnetic order and induces a quantum-disordered phase exhibiting signatures consistent with Majorana excitations, including an anomalous thermal Hall effect and field-angle-dependent specific heat. At higher fields, the magnetization approaches saturation, suggesting a transition to a spin-polarized state, yet the microscopic evolution between these limits remains unresolved. Here we report high-field specific heat measurements up to 24 T that reveal a distinct crossover at $μ_0H^*\approx$15 T, beyond which the perturbative Kitaev description breaks down. Above
H∗, the characteristic six-fold angular modulation of the specific heat collapses and the excitation gap deviates from the predicted
H3 scaling. Meanwhile, the gap decreases with increasing field and the in-plane magnetization anisotropy persists up to
∼ 24 T, both in sharp contrast to a trivial spin-polarized state, indicating that KQSL signatures are preserved even at
∼ 90 % of magnetization saturation. These results reveal that the KQSL in
α-RuCl
3 extends well beyond the perturbative window, persisting as a nonperturbative regime in which the Majorana and vison energy scales merge, before eventually giving way to spin polarization. This thermodynamic roadmap provides a basis for understanding how fractionalized phases evolve under strong magnetic fields.