Zeeman-Kitaev Honeycomb Model

• Zeeman-Kitaev Honeycomb Model is a quantum spin system defined on a honeycomb lattice with bond-directional interactions and an external Zeeman field.
• The model integrates a Kitaev Hamiltonian with Zeeman coupling and field-induced terms, enabling transitions from a Z2 spin liquid to non-Abelian and polarized phases.
• Intermediate-field regimes exhibit competing interpretations including gapless U(1) spin liquids, abelian chiral phases, and symmetry-broken antiferromagnetic states.

The Zeeman–Kitaev honeycomb model is the Kitaev honeycomb spin model supplemented by an external magnetic field, typically treated through a Zeeman coupling to spin-$\tfrac12$ moments on the honeycomb lattice. It is a canonical setting for studying how bond-directional exchange, fractionalization, and emergent gauge structure evolve under field tuning. In the pure Kitaev limit the model realizes a $\mathbb Z_2$ quantum spin liquid with Majorana fermions and flux excitations, while field perturbations generate gapped non-Abelian descendants, polarized states, and—particularly for antiferromagnetic couplings—an intermediate-field regime that has been described in distinct ways, including as a gapless $U(1)$ spin liquid, an abelian chiral phase with fermionic bulk excitations, and a symmetry-broken antiferromagnetic phase associated with anyon gap closings (Hickey et al., 2018).

1. Microscopic definition

A standard form of the Zeeman–Kitaev honeycomb Hamiltonian is

$$H \;=\;-K_x\!\sum_{\langle ij\rangle_x}\!\sigma_i^x\sigma_j^x\,-\,K_y\!\sum_{\langle ij\rangle_y}\!\sigma_i^y\sigma_j^y\,-\,K_z\!\sum_{\langle ij\rangle_z}\!\sigma_i^z\sigma_j^z \;-\;\sum_i\bigl(h_x\,\sigma_i^x+h_y\,\sigma_i^y+h_z\,\sigma_i^z\bigr)\,,$$

where each bond $\langle ij\rangle_\gamma$ carries an Ising-type exchange $K_\gamma$ coupling the $\gamma$-components of two $S=\tfrac12$ spins, and ${\bf h}=(h_x,h_y,h_z)$ is an arbitrary tilted magnetic field (Hickey et al., 2018).

Several closely related formulations appear in the literature. A commonly studied isotropic version sets $K_x=K_y=K_z\equiv K$, often with the field along $\mathbb Z_2$0, so that $\mathbb Z_2$1 (Gao et al., 2019). In perturbative treatments the field also induces a three-spin “Haldane” term,

$\mathbb Z_2$2

or equivalently a next-nearest-neighbor Majorana hopping that gaps the Majorana spectrum (Chen et al., 2022). Extensions relevant to candidate materials further add isotropic Heisenberg exchange $\mathbb Z_2$3, off-diagonal symmetric $\mathbb Z_2$4 exchange, or second-neighbor Dzyaloshinskii–Moriya interaction $\mathbb Z_2$5 (Wang et al., 2019).

The model is therefore not a single isolated Hamiltonian but a core bond-directional spin model together with a hierarchy of field-induced and symmetry-allowed perturbations. This broader usage is especially important in the literature on honeycomb spin-orbit-entangled Mott insulators, where Zeeman, $\mathbb Z_2$6, $\mathbb Z_2$7, and DM terms are all treated as part of the effective low-energy description.

2. Fractionalization and emergent gauge structure

The exactly solvable Kitaev limit admits a Majorana-$\mathbb Z_2$8 gauge formulation. One representation is

$\mathbb Z_2$9

with the on-site constraint $U(1)$0, and bond variables

$U(1)$1

on an $U(1)$2-link $U(1)$3. The plaquette flux is

$U(1)$4

In a fixed gauge-flux sector, the Kitaev Hamiltonian reduces to a quadratic Majorana problem, and different $U(1)$5 sectors define distinct band structures (Chen et al., 2022).

This $U(1)$6 formulation underlies the conventional low-field Kitaev spin liquid. It also furnishes the language for visons, which are flux excitations of the emergent gauge field. In the exactly solvable model they are static, whereas perturbations such as a Zeeman field induce vison hopping and, depending on the sign of $U(1)$7, distinct symmetry implementations on vison motion (Chen et al., 2022).

A complementary low-energy description employs fermionic partons coupled to an emergent $U(1)$8 gauge field. One introduces Abrikosov fermions $U(1)$9 and writes

$$H \;=\;-K_x\!\sum_{\langle ij\rangle_x}\!\sigma_i^x\sigma_j^x\,-\,K_y\!\sum_{\langle ij\rangle_y}\!\sigma_i^y\sigma_j^y\,-\,K_z\!\sum_{\langle ij\rangle_z}\!\sigma_i^z\sigma_j^z \;-\;\sum_i\bigl(h_x\,\sigma_i^x+h_y\,\sigma_i^y+h_z\,\sigma_i^z\bigr)\,,$$0

with $$H \;=\;-K_x\!\sum_{\langle ij\rangle_x}\!\sigma_i^x\sigma_j^x\,-\,K_y\!\sum_{\langle ij\rangle_y}\!\sigma_i^y\sigma_j^y\,-\,K_z\!\sum_{\langle ij\rangle_z}\!\sigma_i^z\sigma_j^z \;-\;\sum_i\bigl(h_x\,\sigma_i^x+h_y\,\sigma_i^y+h_z\,\sigma_i^z\bigr)\,,$$1 and a chiral $$H \;=\;-K_x\!\sum_{\langle ij\rangle_x}\!\sigma_i^x\sigma_j^x\,-\,K_y\!\sum_{\langle ij\rangle_y}\!\sigma_i^y\sigma_j^y\,-\,K_z\!\sum_{\langle ij\rangle_z}\!\sigma_i^z\sigma_j^z \;-\;\sum_i\bigl(h_x\,\sigma_i^x+h_y\,\sigma_i^y+h_z\,\sigma_i^z\bigr)\,,$$2-wave pairing channel $$H \;=\;-K_x\!\sum_{\langle ij\rangle_x}\!\sigma_i^x\sigma_j^x\,-\,K_y\!\sum_{\langle ij\rangle_y}\!\sigma_i^y\sigma_j^y\,-\,K_z\!\sum_{\langle ij\rangle_z}\!\sigma_i^z\sigma_j^z \;-\;\sum_i\bigl(h_x\,\sigma_i^x+h_y\,\sigma_i^y+h_z\,\sigma_i^z\bigr)\,,$$3. In this language the Kitaev spin liquid corresponds to $$H \;=\;-K_x\!\sum_{\langle ij\rangle_x}\!\sigma_i^x\sigma_j^x\,-\,K_y\!\sum_{\langle ij\rangle_y}\!\sigma_i^y\sigma_j^y\,-\,K_z\!\sum_{\langle ij\rangle_z}\!\sigma_i^z\sigma_j^z \;-\;\sum_i\bigl(h_x\,\sigma_i^x+h_y\,\sigma_i^y+h_z\,\sigma_i^z\bigr)\,,$$4, which Higgses $$H \;=\;-K_x\!\sum_{\langle ij\rangle_x}\!\sigma_i^x\sigma_j^x\,-\,K_y\!\sum_{\langle ij\rangle_y}\!\sigma_i^y\sigma_j^y\,-\,K_z\!\sum_{\langle ij\rangle_z}\!\sigma_i^z\sigma_j^z \;-\;\sum_i\bigl(h_x\,\sigma_i^x+h_y\,\sigma_i^y+h_z\,\sigma_i^z\bigr)\,,$$5 to $$H \;=\;-K_x\!\sum_{\langle ij\rangle_x}\!\sigma_i^x\sigma_j^x\,-\,K_y\!\sum_{\langle ij\rangle_y}\!\sigma_i^y\sigma_j^y\,-\,K_z\!\sum_{\langle ij\rangle_z}\!\sigma_i^z\sigma_j^z \;-\;\sum_i\bigl(h_x\,\sigma_i^x+h_y\,\sigma_i^y+h_z\,\sigma_i^z\bigr)\,,$$6, while an intermediate field can restore full $$H \;=\;-K_x\!\sum_{\langle ij\rangle_x}\!\sigma_i^x\sigma_j^x\,-\,K_y\!\sum_{\langle ij\rangle_y}\!\sigma_i^y\sigma_j^y\,-\,K_z\!\sum_{\langle ij\rangle_z}\!\sigma_i^z\sigma_j^z \;-\;\sum_i\bigl(h_x\,\sigma_i^x+h_y\,\sigma_i^y+h_z\,\sigma_i^z\bigr)\,,$$7 gauge invariance and produce a spinon Fermi surface coupled to a massless gauge photon (Hickey et al., 2018).

The same emergent $$H \;=\;-K_x\!\sum_{\langle ij\rangle_x}\!\sigma_i^x\sigma_j^x\,-\,K_y\!\sum_{\langle ij\rangle_y}\!\sigma_i^y\sigma_j^y\,-\,K_z\!\sum_{\langle ij\rangle_z}\!\sigma_i^z\sigma_j^z \;-\;\sum_i\bigl(h_x\,\sigma_i^x+h_y\,\sigma_i^y+h_z\,\sigma_i^z\bigr)\,,$$8 framework is used to describe field-induced internal gauge flux in the presence of second-neighbor DM interactions. The scalar spin chirality

$$H \;=\;-K_x\!\sum_{\langle ij\rangle_x}\!\sigma_i^x\sigma_j^x\,-\,K_y\!\sum_{\langle ij\rangle_y}\!\sigma_i^y\sigma_j^y\,-\,K_z\!\sum_{\langle ij\rangle_z}\!\sigma_i^z\sigma_j^z \;-\;\sum_i\bigl(h_x\,\sigma_i^x+h_y\,\sigma_i^y+h_z\,\sigma_i^z\bigr)\,,$$9

generates an internal flux $\langle ij\rangle_\gamma$0 through second-neighbor triangles according to

$\langle ij\rangle_\gamma$1

with $\langle ij\rangle_\gamma$2 for small flux (Gao et al., 2019).

3. Field-driven phases and their diagnostics

Large-scale exact diagonalization on clusters up to $\langle ij\rangle_\gamma$3 sites has been used to map the phase diagram in tilted magnetic fields. The standard diagnostics include the low-lying spectrum $\langle ij\rangle_\gamma$4, the static structure factor

$\langle ij\rangle_\gamma$5

the dynamical structure factor

$\langle ij\rangle_\gamma$6

the plaquette flux operator $\langle ij\rangle_\gamma$7, and the finite-temperature specific heat

$\langle ij\rangle_\gamma$8

Phase boundaries are identified from peaks in $\langle ij\rangle_\gamma$9 and minima in the ground-state fidelity $K_\gamma$0 (Hickey et al., 2018).

Regime	Principal characterization	Representative hallmarks
Small field	Kitaev spin liquid (KSL)	three-fold quasi-degenerate torus manifold; finite vison gap; two peaks in $K_\gamma$1
Intermediate field	gapless phase / intermediate phase	enhanced low-energy levels; flux-gap collapse; no sharp Bragg peaks
Large field	polarized (PL) phase	product state; magnon excitations; one Schottky-like peak in $K_\gamma$2

For antiferromagnetic $K_\gamma$3, three regimes were reported as $K_\gamma$4 increases. At small field, the Kitaev spin liquid survives and becomes gapped once the field has a component along $K_\gamma$5, acquiring non-Abelian Ising topological order. At large field, the system enters a trivial polarized phase. Between them, for a wide range of tilt angles, exact diagonalization identifies an intermediate gapless phase with a strong enhancement of low-energy states, a featureless $K_\gamma$6, strong $K_\gamma$7 fluctuations, collapse of the zero-field flux gap in $K_\gamma$8, and a downward shift of the low-$K_\gamma$9 specific-heat peak while the high-$\gamma$0 peak remains almost unchanged (Hickey et al., 2018).

In the Zeeman-plus-DM model, exact-diagonalization phase diagrams place an intermediate $\gamma$1 spin liquid in the approximate field window

$\gamma$2

for antiferromagnetic $\gamma$3 when the field is close to $\gamma$4. In that setting the gapless $\gamma$5 phase survives a finite second-neighbor DM interaction up to $\gamma$6 (Gao et al., 2019).

4. Intermediate-field regime and competing descriptions

The intermediate-field regime in the antiferromagnetic model is the main conceptual focus of the Zeeman–Kitaev literature. One interpretation is a field-driven gapless $\gamma$7 spin liquid with fermionic spinons and a massless gauge field. In the parton description, the low-field Kitaev spin liquid is a paired state with $\gamma$8, while increasing field suppresses the pairing amplitude at a critical $\gamma$9, restores $S=\tfrac12$0 gauge invariance, and produces a spinon Fermi surface. The transition into the polarized state is then described as confinement via monopole proliferation once the spinons are gapped (Hickey et al., 2018).

A distinct microscopic gauge-theory description has been formulated in terms of a lattice mutual Chern–Simons theory. In that construction, both ferro- and antiferromagnetic models exhibit low-field nonabelian Ising topological order with $S=\tfrac12$1. For the antiferromagnetic case, an intermediate-field phase appears with effective $S=\tfrac12$2-matrix

$S=\tfrac12$3

no ground-state degeneracy, $S=\tfrac12$4, and fermionic bulk excitations; despite the change in topological order, the chiral central charge remains $S=\tfrac12$5, so the half-quantized thermal Hall response is unchanged across the transition (Das et al., 2024).

A third line of work analyzes field-induced anyon polarons. In that approach, visons, fermionic polarons, and bosonic bound states are tracked perturbatively in small Zeeman field. For the antiferromagnetic model, single visons become gapless at $S=\tfrac12$6, the first fermion closing occurs at $S=\tfrac12$7, and boson condensation occurs at $S=\tfrac12$8. The bosonic excitation carries the quantum numbers of an in-plane Néel order orthogonal to $S=\tfrac12$9, and its soft mode transforms as an ${\bf h}=(h_x,h_y,h_z)$0 irreducible representation of ${\bf h}=(h_x,h_y,h_z)$1; this was used to argue that the intermediate phase is a genuine antiferromagnetic phase rather than a spin liquid (Chen et al., 2024).

These descriptions are not identical. The literature therefore contains a genuine disagreement over whether the antiferromagnetic intermediate-field regime is best understood as a gapless ${\bf h}=(h_x,h_y,h_z)$2 spin liquid, an abelian chiral phase with trivial two-dimensional topological order, or a symmetry-broken ordered phase. What is common across these analyses is that the intermediate regime is tied to vison-gap collapse, strong field sensitivity of gauge-sector excitations, and a nontrivial restructuring of the low-energy spectrum.

5. Thermal Hall response, visons, and band topology

The Zeeman–Kitaev model is closely connected to thermal Hall transport because magnetic field breaks time-reversal symmetry and can induce topological band structure for emergent quasiparticles. In a ${\bf h}=(h_x,h_y,h_z)$3 spin liquid with a spinon Fermi surface, second-neighbor DM interactions and Zeeman coupling generate an internal ${\bf h}=(h_x,h_y,h_z)$4 gauge flux that twists spinon motion. At mean-field level the spinons occupy flux-twisted bands ${\bf h}=(h_x,h_y,h_z)$5 with Berry curvature ${\bf h}=(h_x,h_y,h_z)$6, and the thermal Hall conductivity takes the form

${\bf h}=(h_x,h_y,h_z)$7

with

${\bf h}=(h_x,h_y,h_z)$8

For ${\bf h}=(h_x,h_y,h_z)$9, $K_x=K_y=K_z\equiv K$0 approaches a finite, noninteger plateau as $K_x=K_y=K_z\equiv K$1, then decreases monotonically with temperature and vanishes at high $K_x=K_y=K_z\equiv K$2 (Gao et al., 2019).

Vison motion contributes a separate topological channel. In a Zeeman field, the vison hopping problem maps to a tight-binding model on the triangular dual lattice. The key distinction is between ferromagnetic and antiferromagnetic Kitaev exchange. For $K_x=K_y=K_z\equiv K$3, the vison experiences $K_x=K_y=K_z\equiv K$4, translations act non-projectively, the vison band has zero Berry curvature, and there is no intrinsic vison thermal Hall contribution. For $K_x=K_y=K_z\equiv K$5, the vison experiences $K_x=K_y=K_z\equiv K$6, translations act projectively, the unit cell doubles, and two gapped vison bands appear with Chern numbers $K_x=K_y=K_z\equiv K$7 and $K_x=K_y=K_z\equiv K$8, yielding an intrinsic vison contribution to $K_x=K_y=K_z\equiv K$9 (Chen et al., 2022).

The low-field non-Abelian phase generated by a field-induced three-spin term has its own quantized response. In proximate Kitaev spin liquids of the $\mathbb Z_2$00-$\mathbb Z_2$01-$\mathbb Z_2$02 model, a field normal to the plane gaps the state into chiral spin liquids with Chern numbers $\mathbb Z_2$03 and $\mathbb Z_2$04, corresponding to

$\mathbb Z_2$05

before a transition to a trivial polarized phase (Wang et al., 2019). This shows that thermal Hall plateaus in honeycomb Kitaev systems need not be limited to the conventional $\mathbb Z_2$06 response.

6. Perturbations, extensions, and broader formulations

The field-driven phenomenology of the Zeeman–Kitaev model has been tested against perturbations expected in materials. Exact diagonalization with an additional Heisenberg term $\mathbb Z_2$07 and off-diagonal $\mathbb Z_2$08 exchange found that for modest $\mathbb Z_2$09, both the Kitaev spin liquid and the intermediate gapless phase survive, with the main effect being a shift of critical fields. Larger perturbations favor magnetically ordered states including zig-zag, Néel, and “vortex” patterns (Hickey et al., 2018).

Variational Monte Carlo studies of the $\mathbb Z_2$10-$\mathbb Z_2$11-$\mathbb Z_2$12 model further identify a “generic” Kitaev spin liquid and a single proximate Kitaev spin liquid. The proximate phase is a gapless $\mathbb Z_2$13 state with 14 Majorana cones and a gapless spin response; in a field along $\mathbb Z_2$14, it realizes a non-Abelian chiral spin liquid with $\mathbb Z_2$15, then an Abelian chiral spin liquid with $\mathbb Z_2$16, and only at larger field a trivial polarized phase (Wang et al., 2019). This places the Zeeman perturbation within a broader landscape of nearby honeycomb spin liquids.

Recent work has also recast the model in more general duality frameworks. A staggered-Majorana or “dumbbell fermion” construction expresses Kitaev’s exact solution as a special case of a higher-dimensional duality between Pauli spins and Majorana fermions, with the one-particle problem in each $\mathbb Z_2$17 flux sector governed by

$\mathbb Z_2$18

In this formulation the Zeeman term spoils exact solvability, but third-order perturbation theory again generates the effective chiral term proportional to $\mathbb Z_2$19 that gaps the Majorana spectrum (Banks, 13 Feb 2025).

Taken together, these developments place the Zeeman–Kitaev honeycomb model at the intersection of exactly solvable $\mathbb Z_2$20 gauge theory, emergent $\mathbb Z_2$21 gauge dynamics, topological band theory of fractionalized quasiparticles, and materials-motivated extensions. The recurring technical themes are the competition between Majorana and vison sectors, the role of field-induced chirality and gauge flux, and the sensitivity of the intermediate regime to both microscopic perturbations and the theoretical framework used to describe it.