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 introduce a framework of constructing topological subsystem codes based on the class of anticommuting quantum spin liquids described in [Phys. Rev. B 113, 064402 (2026)]. A canonical model from this class can be considered as a spatial modification of the toric code that voids its stabilizer code property. Rather, these models contain an extensive set of anticommuting local conserved operators that lead to an extensive ground state degeneracy. This degeneracy forms the basis of the subsystem degrees of freedom in the associated quantum error correcting code. The code inherits the many-body topological order of the quantum spin liquid, making it a topological subsystem code. We present two concrete and detailed examples for constructing these codes on a square lattice and a kagome lattice geometry, requiring weight-4 and weight-3 local check operator measurements respectively. In contrast to other subsystem codes, a unique property of these codes is the presence of an extensive number of local gauge qubits that are left undisturbed by the check operators apart from the logical qubits. Our construction provides a template for generating this new category of topological subsystem codes on different lattice or graph geometries, suitable for implementation on various quantum hardware platforms.
In this chapter, we consider Kitaev's honeycomb lattice model (Kitaev, 2006). This is an analytically tractable spin model that gives rise to quasiparticles with Abelian as well as non-Abelian statistics. Some of its properties are similar to the fractional quantum Hall effect, which has been studied experimentally in great detail even though it evades exact analytical treatment (Moore and Read, 1991). Due to its simplicity, the honeycomb lattice model is likely to be the first topological spin model to be realised in the laboratory, e.g., with optical lattice technology (Micheli et al. , 2006). Understanding its properties can facilitate its physical realisation and can provide a useful insight into the mechanisms underlining topological insulators and the fractional quantum Hall effect. The honeycomb lattice model comprises interacting spin-½ particles arranged on the sites of a honeycomb lattice. It is remarkable that such a simple model can support a rich variety of topological behaviours. For certain values of its couplings, Abelian anyons emerge that behave like the toric code anyons. For another coupling regime, non-Abelian anyons emerge that correspond to the Ising anyons. The latter are manifested as vortex-like configurations of the original spin model that can effectively be described by Majorana fermions. These are fermionic fields that are antiparticles of themselves. They were first introduced in the context of high-energy physics (Majorana, 1937) and become increasingly important in the analysis of solid state phenomena (Wilczek, 2009).
In the field of quantum magnetism, the exactly solvable Kitaev honeycomb model serves as a paradigm for the fractionalization of spin degrees of freedom and the formation of [Formula: see text] quantum spin liquids. An intense experimental search has led to the discovery of a number of spin-orbit entangled Mott insulators that realize its characteristic bond-directional interactions and, in the presence of magnetic fields, exhibit no indications of long-range order. Here, we map out the complete phase diagram of the Kitaev model in tilted magnetic fields and report the emergence of a distinct gapless quantum spin liquid at intermediate field strengths. Analyzing a number of static, dynamical, and finite temperature quantities using numerical exact diagonalization techniques, we find strong evidence that this phase exhibits gapless fermions coupled to a massless U(1) gauge field. We discuss its stability in the presence of perturbations that naturally arise in spin-orbit entangled candidate materials.
Abstract The Kitaev honeycomb model is an approximate topological quantum error correcting code in the same phase as the toric code, but requiring only a 2-body Hamiltonian. As a frustrated spin model, it is well outside the commuting models of topological quantum codes that are typically studied, but its exact solubility makes it more amenable to analysis of effects arising in this noncommutative setting than a generic topologically ordered Hamiltonian. Here we study quantum error correction in the honeycomb model using both analytic and numerical techniques. We first prove explicit exponential bounds on the approximate degeneracy, local indistinguishability, and correctability of the code space. These bounds are tighter than can be achieved using known general properties of topological phases. Our proofs are specialized to the honeycomb model, but some of the methods may nonetheless be of broader interest. Following this, we numerically study noise caused by thermalization processes in the perturbative regime close to the toric code renormalization group fixed point. The appearance of non-topological excitations in this setting has no significant effect on the error correction properties of the honeycomb model in the regimes we study. Although the behavior of this model is found to be qualitatively similar to that of the standard toric code in most regimes, we find numerical evidence of an interesting effect in the low-temperature, finite-size regime where a preferred lattice direction emerges and anyon diffusion is geometrically constrained. We expect this effect to yield an improvement in the scaling of the lifetime with system size as compared to the standard toric code.
We develop an exactly solvable model with Kitaev-type interactions and study its phase diagram on the dual lattice of the quasicrystalline Ammann-Beenker lattice. Our construction is based on the $\mathrm{\ensuremath{\Gamma}}$-matrix generalization of the Kitaev model and utilizes the cut-and-project correspondence between the four-dimensional simple cubic lattice and the Ammann-Beenker lattice to designate four types of bonds. We obtain a rich phase diagram with gapped (chiral and abelian) and gapless spin liquid phases via Monte Carlo simulations and variational analysis. We show that the ground state can be further tuned by the inclusion of an onsite term that selects 21 different vison configurations while maintaining the integrability of the model. Our results highlight the rich physics at the intersection of quasicrystals and quantum magnetism.