---
title: Kitaev-Heisenberg Model Overview
url: https://www.emergentmind.com/topics/kitaev-heisenberg-model
type: topic
---

# Kitaev-Heisenberg Model Overview

Searching arXiv for the cited Kitaev–Heisenberg papers to ground the article in published work.
The Kitaev-Heisenberg model is a class of spin Hamiltonians that combines bond-directional Kitaev exchange with isotropic Heisenberg exchange, originally formulated on the honeycomb lattice as an effective description of the low-energy magnetism of layered iridates \(A_2\)IrO\(_3\) with \(A=\)Li, Na. In that setting, strong spin-orbit coupling acting on \(\mathrm{Ir}^{4+}\) ions in edge-sharing octahedra produces effective \(j_{\mathrm{eff}}=1/2\) moments and highly anisotropic nearest-neighbor exchange. The model interpolates between the antiferromagnetic Heisenberg limit and the exactly solvable Kitaev limit, while also exhibiting an intermediate exact stripy antiferromagnet at a special point of the interpolation [1004.2964].

## 1. Microscopic origin and canonical formulation

In the iridate construction, each \(\mathrm{Ir}^{4+}\) ion sits in an octahedral environment, so the \(t_{2g}\) manifold carries an effective orbital angular momentum \(l=1\). Strong spin-orbit coupling then produces a Kramers doublet with total angular momentum \(j_{\mathrm{eff}}=1/2\), treated as an effective spin-\(\tfrac12\) degree of freedom. Because the octahedra share edges, neighboring ions exchange through bond-dependent \(90^\circ\) Ir–O–Ir hopping paths and also direct \(dd\) overlap, generating highly anisotropic exchange interactions [1004.2964].

On the honeycomb lattice, the nearest-neighbor Hamiltonian on a bond of type \(\gamma=x,y,z\) is
\[
{\cal H}_{ij}^{(\gamma)}=-J_{1}\,S_{i}^{\gamma}S_{j}^{\gamma} +J_{2}\,\mathbf{S}_{i}\cdot\mathbf{S}_{j}.
\]
The first term is a bond-dependent Ising interaction of Kitaev type, while the second is the ordinary Heisenberg exchange. The three bond types correspond to the three cubic components \(\gamma=x,y,z\) in the local axes of the IrO\(_6\) octahedra. In the microscopic derivation,
\[
J_1=\eta_1+2\eta_2,\qquad J_2=\eta_2+\eta_3,
\]
with \(\eta_1\) from Hund’s coupling on the Ir ion, \(\eta_2\) from oxygen-mediated processes involving \(U_p\), and \(\eta_3=(t'/t)^2\) from direct Ir–Ir \(dd\) hopping. The energy unit is \(4t^2/9U_d\), where \(t=t_{pd\pi}^2/\Delta_{pd}\) [1004.2964].

A commonly used one-parameter form sets
\[
J_1=2\alpha,\qquad J_2=1-\alpha,\qquad 0\le \alpha\le 1,
\]
so that \(\alpha=0\) is the antiferromagnetic Heisenberg limit and \(\alpha=1\) is the pure Kitaev limit. Later studies often use angular parametrizations such as \(K=\sin\varphi,\ J=\cos\varphi\) or \(K=2\sin\alpha,\ J=\cos\alpha\), which sweep the full nearest-neighbor coupling space including ferromagnetic and antiferromagnetic sectors [1608.05333], [1911.12854].

## 2. Solvable limits, hidden symmetry, and duality

The defining structural feature of the model is the coexistence of exactly solvable or symmetry-enhanced limits with intervening frustrated regimes. At \(\alpha=0\), the honeycomb model is the antiferromagnetic Heisenberg model; at \(\alpha=1\), it is the exactly solvable Kitaev model. The original honeycomb study also identified a special solvable point at \(\alpha=\tfrac12\), where a four-sublattice spin rotation maps the Hamiltonian to an isotropic ferromagnet in rotated variables,
\[
{\cal H}_{ij}^{(\gamma)}=-\frac12\,\tilde{\mathbf S}_i\cdot \tilde{\mathbf S}_j,
\]
so the exact ground state is a fully polarized ferromagnet in the rotated basis and a stripy antiferromagnetic pattern in the original spins [1004.2964].

This hidden ferromagnetic structure became the basis for a broader duality framework. On the honeycomb and on other lattices built from edge-sharing \(\mathrm{IrO}_6\) octahedra, site-dependent \(\pi\)-rotations about spin axes define a Klein four-group structure \(\mathbb Z_2\times\mathbb Z_2\). In the conventional exchange language, the duality acts as
\[
J_H \to -J_H,\qquad J_K \to J_K + 2J_H,
\]
and it predicts fluctuation-free ordered states that are analogs of honeycomb stripy order on triangular, kagome, hyperkagome, fcc, and pyrochlore lattices [1303.3290]. On the three-dimensional hyperhoneycomb lattice, a hidden four-sublattice symmetry similarly generates hidden SU(2) points at \(\varphi/\pi=3/4\) and \(7/4\), in addition to the ordinary Heisenberg points at \(\varphi/\pi=0\) and \(1\) [2303.09156].

A recurrent implication is that apparently complicated collinear orders may become simple ferromagnetic or antiferromagnetic states in rotated variables. This suggests that the phase structure of the Kitaev-Heisenberg family is organized not only by competing exchanges, but also by nontrivial transformations of the spin basis that expose hidden isotropic points [1004.2964], [1303.3290].

## 3. Honeycomb phase diagram and critical behavior

For the spin-\(\tfrac12\) honeycomb model in the original \(\alpha\in[0,1]\) interpolation, exact diagonalization on a 24-site cluster and complementary spin-wave analysis found three main zero-temperature regions: a conventional Néel phase near \(\alpha=0\), a stripy antiferromagnetic phase centered around \(\alpha=\tfrac12\), and an extended spin-liquid phase near the Kitaev limit. The Néel-to-stripy transition was described as first-order, while the stripy-to-spin-liquid transition was described as second-order or weakly first-order based on the numerical data. The most likely transition points were identified near \(\alpha\simeq 0.4\) and \(\alpha\simeq 0.8\) [1004.2964].

Semiclassically, the ordered phases are strongly affected by accidental degeneracies and their lifting by fluctuations. In the original honeycomb study, the classical Néel-stripy boundary sits at \(\alpha=\tfrac13\), where the linear spin-wave spectra develop zero-energy lines; quantum fluctuations shift the boundary upward and open a spin-wave gap, with
\[
\Delta \simeq \frac{2}{\alpha}\left(\alpha-\frac12\right)^2
\]
near \(\alpha\sim\tfrac12\). The stripy state at the midpoint is fluctuation-free and has a saturated order parameter despite being antiferromagnetic in the original variables [1004.2964].

When the full nearest-neighbor coupling circle is considered, the honeycomb phase diagram contains four magnetically ordered phases—Néel, zigzag, ferromagnetic, and stripy—separated by two Kitaev spin-liquid phases around the antiferromagnetic and ferromagnetic Kitaev points. Exact diagonalization, cluster mean-field theory, linear spin-wave theory, and second-order perturbation theory all support the sequence
\[
\text{Néel} \;\to\; \text{AF KSL} \;\to\; \text{zigzag} \;\to\; \text{FM} \;\to\; \text{FM KSL} \;\to\; \text{stripy} \;\to\; \text{Néel},
\]
with the ferromagnetic-side KSL substantially broader than the antiferromagnetic-side KSL because the neighboring ferromagnetic and stripy phases have very weak quantum fluctuations [1608.05333].

At finite temperature, the classical honeycomb Kitaev-Heisenberg model exhibits a three-phase structure: a low-temperature magnetically ordered phase with spontaneously broken \(Z_6\) symmetry, an intermediate critical Kosterlitz-Thouless phase with emergent \(U(1)\) symmetry and algebraic correlations, and a high-temperature disordered phase. Thermal fluctuations select collinear order along cubic axes by order-by-disorder, and finite-size scaling gives exponents near \(\eta=1/9\) and \(\eta=1/4\) at the lower and upper boundaries of the critical phase, as expected for six-state clock criticality [1205.3967].

The nature of the quantum transition out of the spin liquid remains nontrivial. A slave-particle mean-field treatment of the transition between the gapless \(Z_2\) spin liquid and stripy antiferromagnet found a discontinuous transition at mean-field level, but also argued that spinon confinement effects associated with the instability of a gapped \(U(1)\) spin liquid in two spatial dimensions may be important, leaving open the possibility of a more exotic continuous transition beyond mean field [1206.5814].

## 4. Generalizations in spin, dimension, and lattice geometry

The honeycomb model has been generalized in several distinct directions. For spin-1 local moments on the two-dimensional honeycomb lattice, iDMRG on infinite cylinders finds two spin-liquid phases and four symmetry-broken phases. The spin liquids occur near the pure antiferromagnetic and ferromagnetic Kitaev points, are gapless according to finite-entanglement scaling with central charge \(c=1\), and show approximate \(Z_2\) local conservations: the plaquette Wilson-loop expectation value stays near \(1\), and the static spin-spin correlations remain short-range in the entire spin-liquid phases [1911.12854].

In three dimensions, the spin-\(\tfrac12\) model on the hyperhoneycomb lattice exhibits a phase diagram closely paralleling the two-dimensional honeycomb case. PFFRG for the isotropic model \(J_x=J_y=J_z=1\) finds four magnetically ordered phases—Néel AFM, zigzag AFM, FM, and stripy AFM—and QSL regions around both pristine Kitaev points, with the ferromagnetic-side QSL wider than the antiferromagnetic-side QSL. Introducing anisotropy through \(J_x=J_y=(3-J_z)/2\) narrows the QSL region and replaces part of it by stripy order, while strong anisotropy can lead to a dimer phase [2303.09156].

An earlier hyperhoneycomb study, motivated by \(\beta\)-Li\(_2\)IrO\(_3\), combined semiclassical analysis, the exact solution at the Kitaev point, and slave-fermion mean-field theory. It found four collinear magnetic phases—Néel, polarized ferromagnet, skew-stripy, and skew-zig-zag—together with an extended three-dimensional \(\mathbb Z_2\) spin liquid around the Kitaev point. In that spin liquid, the Majorana spectrum has a gapless line node described as a deformed Fermi-circle of co-dimension \(d_c=2\) [1308.05333].

The bond-directional mechanism is not restricted to tricoordinated honeycomb-derived systems. Analogous Kitaev-Heisenberg interactions were argued to arise on triangular, kagome, hyperkagome, fcc, and pyrochlore lattices built from edge-sharing \(\mathrm{IrO}_6\) octahedra. On the triangular lattice, 2D DMRG found a fully magnetically ordered phase diagram with \(120^\circ\) antiferromagnetic, \(\mathbb Z_2\)-vortex crystal, nematic, dual \(\mathbb Z_2\)-vortex crystal, \(\mathbb Z_6\) ferromagnetic, and dual ferromagnetic phases, with first-order transitions between them [1512.02334]. On a honeycomb-triangular interpolation, classical and quantum studies found that known honeycomb and triangular phases can merge through coexistence regions such as HN-\(\mathbb Z_2\)VC and extended nematic or stripy regimes, while quantum fluctuations were reported to affect the phase diagram only weakly [1804.06080].

## 5. Perturbations and enriched descendants

Doping converts the spin model into a setting for unconventional superconductivity. In an \(SU(2)\) slave-boson mean-field treatment of the doped honeycomb model, light doping near the Kitaev limit stabilizes a triplet \(p\)-wave superconducting state \(p\)SC\(_1\) that breaks time-reversal symmetry and has Chern number \(+1\), irrespective of the sign of the Kitaev interaction. At larger doping, the antiferromagnetic Kitaev side favors a singlet \(d_{x^2-y^2}+id_{xy}\) state and eventually an \(s\)-wave state on the AF Heisenberg side, whereas the ferromagnetic Kitaev side favors a distinct time-reversal-symmetric triplet state \(p\)SC\(_2\), which may be topologically trivial or nontrivial depending on doping and interactions [1212.5218]. In an itinerant formulation around quarter filling, strong spin-orbit coupling was argued to generate six finite-momentum inversion symmetry centers of the Fermi surface, leading to spin-triplet FFLO superconductivity with three separated degenerate ground states and finite-momentum Cooper pairing [1511.03289].

Magnetic fields strongly reorganize the phase structure. A systematic \(1/S\) expansion for fields along \([001]\) and \([111]\) showed that quantum corrections substantially modify the classical phase diagram, reduce the stability of the high-field polarized phase, and strongly suppress exotic large-unit-cell phases in a \([111]\) field [2007.03717]. A 24-site exact-diagonalization study of the spin-\(\tfrac12\) model in field found overall agreement with nonlinear spin-wave theory and identified an intermediate-field spin-disordered phase, especially robust for \(h\parallel[111]\), where it survives well away from the pure Kitaev point and may end near a quantum tricritical point [2111.11474].

Disorder changes the competition between spin liquid and magnetic order in a direction-dependent way. Lanczos exact diagonalization with random-coupling and singular-coupling disorder found that, in the nearest-neighbor honeycomb model, disorder shrinks the Kitaev spin-liquid window, broadens sharp transitions into more crossover-like features, replaces long-range zigzag and stripy order by their three domains with different ordering direction, and eventually drives the system toward a spin-glass-like state. Disorder can also close the flux gap, generate vortices in the plaquette-flux arrangement, and produce commensurate flux patterns. When second- and third-neighbor Heisenberg interactions are included, however, singular-coupling disorder can instead suppress long-range magnetic order and expand the spin-liquid regime [2210.17198].

Additional interactions and geometry further enrich the phase diagram. Adding the bond-dependent off-diagonal \(\Gamma\) term to obtain the extended Kitaev-Heisenberg model produces a global phase diagram with eight distinct quantum phases—Kitaev spin liquid, FM, AFM, stripy, zigzag, rotated FM, rotated AFM, and a valence bond solid in a quadro-critical region—together with a dual mapping
\[
J'=-\Gamma,\qquad K'=-(K-\Gamma-J),\qquad \Gamma'=-J
\]
that preserves the model form [1501.06990]. In the ferromagnetic phase of the extended model, linear spin-wave theory finds topological magnon bands with chiral zig-zag edge states protected by non-zero Chern numbers, and for \([001]\) polarization a topological phase transition can be driven by anisotropy of the Kitaev couplings [1803.01515]. On the bilayer honeycomb lattice with interlayer Heisenberg exchange, large-scale iPEPS and high-order series expansions find the familiar FM, AFM, zigzag, stripy, and two Kitaev QSL phases, but also a new rung-singlet valence bond solid between AFM and stripy order, adiabatically connected to isolated Heisenberg dimers [2412.17495].

## 6. Materials relevance, interpretation, and open issues

The original physical motivation was that layered iridates may realize a broad swath of the Kitaev-Heisenberg phase diagram because strong spin-orbit coupling converts the orbital structure of the \(t_{2g}\) electrons into bond-dependent magnetic anisotropy. The honeycomb model therefore proposed \(A_2\)IrO\(_3\) as a realistic setting in which Kitaev physics, quantum spin-liquid behavior, and stripy magnetic order can occur in the same microscopic model. At the same time, the original study emphasized that the available experimental information was insufficient to pin down the materials’ exact location in the phase diagram, and that Na/Ir site disorder may complicate interpretation because impurities can induce local moments, including in Kitaev-like systems [1004.2964].

Subsequent work sharpened both the usefulness and the limitations of the model as a materials description. The finite-temperature classical study argued that the thermal phase structure of the nearest-neighbor honeycomb model is relevant to Na\(_2\)IrO\(_3\) and likely to Li\(_2\)IrO\(_3\), but also stated that the model alone is not sufficient to fully describe Na\(_2\)IrO\(_3\), because further-neighbor couplings are needed to capture the zigzag ground state and the detailed excitation spectrum [1205.3967]. This suggests that the nearest-neighbor Kitaev-Heisenberg Hamiltonian is best regarded as a minimal or organizing model rather than a complete material-specific Hamiltonian.

In three-dimensional iridates such as \(\beta\)-Li\(_2\)IrO\(_3\), \(\gamma\)-Li\(_2\)IrO\(_3\), and \(\beta\)-ZnIrO\(_3\), the hyperhoneycomb version provides a reference phase diagram with QSL regions near both Kitaev limits and four commensurate ordered phases, but realistic modeling likely requires additional interactions such as the \(\Gamma\) term to account for experimentally observed incommensurate noncoplanar order [2303.09156]. Disorder studies further indicate that bond randomness can either destabilize or extend proximate spin-liquid regimes depending on whether the dominant ordered state originates from nearest-neighbor or longer-range exchanges, a point highlighted as relevant to \(\alpha\)-RuCl\(_3\) and H\(_3\)LiIr\(_2\)O\(_6\) [2210.17198].

A final interpretive issue is that the “Kitaev-Heisenberg model” now names not one unique Hamiltonian but a family of closely related models: the original spin-\(\tfrac12\) honeycomb form, full-circle nearest-neighbor variants, higher-spin and three-dimensional analogs, doped and field-driven descendants, and extended models with \(\Gamma\), further-neighbor, disorder, or interlayer terms. What remains common across this family is the central competition between bond-directional Kitaev exchange and isotropic Heisenberg exchange, together with the hidden symmetry structures and proximate spin-liquid regimes that emerge from that competition [1004.2964], [1608.05333], [1501.06990].

Source: https://www.emergentmind.com/topics/kitaev-heisenberg-model