{ "id": "ep91y9edbv", "slug": "kitaev-honeycomb-model", "title": "Kitaev honeycomb model", "updated": "2026-10-15", "topic_path": [ { "id": "physical", "label": "Physical world and mathematics", "api_url": "https://www.edgepedia.org/api/v1/topics/physical" }, { "id": "physical.physics", "label": "Physics", "api_url": "https://www.edgepedia.org/api/v1/topics/physical.physics" }, { "id": "physical.physics.quantum", "label": "Quantum physics", "api_url": "https://www.edgepedia.org/api/v1/topics/physical.physics.quantum" } ], "geo": [ { "id": "geo.nongeo.t2001.physical.physics", "label": "Non-geographic · 2001 to 2020: Physics", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t2001.physical.physics", "path": [ { "id": "geo.nongeo", "label": "Non-geographic", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo" }, { "id": "geo.nongeo.t2001", "label": "Non-geographic · 2001 to 2020", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t2001" }, { "id": "geo.nongeo.t2001.physical", "label": "Physical world and mathematics", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t2001.physical" }, { "id": "geo.nongeo.t2001.physical.physics", "label": "Physics", "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t2001.physical.physics" } ] } ], "excerpt": "The Kitaev honeycomb model is an exactly solvable spin model on a honeycomb lattice, introduced by Alexei Kitaev in 2006, whose ground state is a topological quantum spin liquid of emergent Majorana fermions.", "snippet": "The Kitaev honeycomb model is an exactly solvable spin model on a honeycomb lattice, introduced by Alexei Kitaev in 2006, whose ground state is a topological quantum spin liquid of emergent Majorana fermions.", "node": "physical.physics.quantum", "markdown": "# Kitaev honeycomb model\n\nThe Kitaev honeycomb model is an exactly solvable spin model on a honeycomb lattice in which every nearest-neighbor bond carries an Ising-type interaction along a bond-dependent axis, and whose ground state is a topological quantum spin liquid built from emergent Majorana fermions. It is solved exactly by a reduction to free fermions in a static \\( \\mathbb{Z}_2 \\) gauge field, a feature that makes it a rare microscopic realization of a quantum spin liquid.[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) • [2](https://journals.aps.org/rmp/abstract/10.1103/3m4m-3v59) Under a magnetic field its gapless phase becomes a [chiral spin liquid](https://www.edgepedia.org/chiral-spin-liquid) with non-Abelian anyons, connecting the model to topological quantum computation.[2](https://journals.aps.org/rmp/abstract/10.1103/3m4m-3v59)\n\n| Key fact | Value |\n|---|---|\n| Introduced by | Alexei Kitaev, \"Anyons in an exactly solved model and beyond\", Annals of Physics 321, 2–111 (2006)[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) |\n| Hamiltonian | \\( H = -J_x \\sum_{x\\text{-links}} \\sigma_j^x \\cdot \\sigma_k^x - J_y \\sum_{y\\text{-links}} \\sigma_j^y \\cdot \\sigma_k^y - J_z \\sum_{z\\text{-links}} \\sigma_j^z \\cdot \\sigma_k^z \\)[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) |\n| Conserved quantity | Plaquette operators \\( W_p \\), a static \\( \\mathbb{Z}_2 \\) gauge field[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) |\n| Phases at zero field | Gapped A phases (toric-code anyons) and gapless B phase under triangle inequalities[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) • [3](https://www.thp.uni-koeln.de/trebst/Lectures/Seminar14/Handout9.pdf) |\n| Field-gapped phase | Non-Abelian chiral spin liquid, Chern number \\( \\nu = \\pm 1 \\)[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) |\n| Minimal material model | \\( J\\text{-}K\\text{-}\\Gamma\\text{-}\\Gamma' \\) exchange for α-RuCl₃ and iridates[4](https://arxiv.org/pdf/2308.01943) |\n| Candidate materials | α-RuCl₃, Na₂IrO₃, α/β/γ-Li₂IrO₃, H₃LiIr₂O₆[5](https://jeffrau.ca/assets/lectures/kitaev-2021.pdf) |\n\n## How it works\n\nEach bond of the honeycomb lattice carries an Ising interaction along one axis: x-links couple \\( \\sigma^x \\sigma^x \\), y-links \\( \\sigma^y \\sigma^y \\), and z-links \\( \\sigma^z \\sigma^z \\), with independently settable strengths \\( J_x, J_y, J_z \\).[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) The model is frustrated, yet it remains exactly solvable.[5](https://jeffrau.ca/assets/lectures/kitaev-2021.pdf) The solution maps each spin to four Majorana operators \\( (b^x, b^y, b^z, c) \\), enlarging the [Hilbert space](https://www.edgepedia.org/hilbert-space) from \\( 2^N \\) to \\( 4^N \\) and requiring a projector back to the physical subspace.[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) • [6](https://thomashodson.com/thesis/2_Background/2.2_HKM_Model.html) The spin operator on a link becomes a Majorana hopping term \\( i \\cdot c_i \\cdot c_j \\) multiplied by a link variable \\( \\hat{u}_{jk} = i \\cdot b^\\alpha_j \\cdot b^\\alpha_\\kappa \\).[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) • [7](https://www.nature.com/articles/s41586-025-09475-0)\n\nThe link variables commute with the Hamiltonian and with each other, so the Hilbert space splits into sectors labeled by a static \\( \\mathbb{Z}_2 \\) gauge field. Within each sector the Hamiltonian is a quadratic form \\( H = \\frac{i}{4} \\sum_{jk} A_{jk} c_j \\cdot c_k \\) in Majorana operators, which is diagonalizable exactly.[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) The gauge-invariant plaquette operators \\( W_p = \\sigma^x_1 \\cdot \\sigma^y_2 \\cdot \\sigma^z_3 \\cdot \\sigma^x_4 \\cdot \\sigma^y_5 \\cdot \\sigma^z_6 \\) commute with the Hamiltonian and measure the flux through each hexagon.[3](https://www.thp.uni-koeln.de/trebst/Lectures/Seminar14/Handout9.pdf) A theorem fixes the ground state to the vortex-free sector, \\( w_p = 1 \\) for all plaquettes.[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) • [5](https://jeffrau.ca/assets/lectures/kitaev-2021.pdf) In that sector the fermion dispersion \\( \\varepsilon(\\mathbf{q}) = \\pm |f(\\mathbf{q})| \\), with \\( f(\\mathbf{q}) = 2(J_x e^{i \\mathbf{q} \\cdot \\mathbf{n}_1} + J_y e^{i \\mathbf{q} \\cdot \\mathbf{n}_2} + J_z) \\), has Dirac cones near the Brillouin-zone corners, the same spectrum as graphene.[3](https://www.thp.uni-koeln.de/trebst/Lectures/Seminar14/Handout9.pdf)\n\n## How it is done\n\nSolving the model by hand or by computer follows the same sequence. First write each spin in the Majorana representation. Second, evaluate the plaquette operators \\( W_p \\) to fix a flux sector; each sector is solvable in \\( O(N^3) \\) by diagonalizing the quadratic Majorana Hamiltonian.[5](https://jeffrau.ca/assets/lectures/kitaev-2021.pdf) Third, use Lieb's theorem to select the flux-free ground sector.[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) Fourth, diagonalize \\( A_{jk} \\) to obtain the fermion spectrum and thermodynamics. For perturbed or non-integrable variants, the same fractionalized structure guides numerical methods: exact diagonalization, variational [Monte Carlo](https://www.edgepedia.org/monte-carlo) with symmetry-guided ansätze, density-matrix renormalization group on cylinders, tensor networks including iPEPS, and quantum Monte Carlo that samples flux sectors at finite temperature.[5](https://jeffrau.ca/assets/lectures/kitaev-2021.pdf) • [8](https://arxiv.org/pdf/2601.14496)\n\n## Origin\n\nThe model was introduced by [Alexei Kitaev](https://www.edgepedia.org/alexei-kitaev) in \"Anyons in an exactly solved model and beyond\", [Annals of Physics](https://www.edgepedia.org/annals-of-physics), volume 321, issue 1 (January 2006), pages 2–111.[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) It built on the earlier toric code paper, published in Annals of Physics 303, 2–30; topologically ordered states can serve as physical analogs of error-correcting codes.[9](https://www.annualreviews.org/content/journals/10.1146/annurev-conmatphys-033117-053934) • [1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) Kitaev's paper also places the work in a longer lineage: the study of anyons was initiated in the early 1980s, the resonating-valence-bond idea was put forward by Philip Anderson, and the weak-pairing p-wave BCS phase with non-Abelian vortices was identified by [Nicholas Read](https://www.edgepedia.org/nicholas-read) and Dmitry Green.[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381)\n\n## Variants\n\nThe model is exactly solvable on a wide range of tricoordinated lattices beyond the honeycomb, including the 3D hyperhoneycomb and stripy-honeycomb geometries relevant to β- and γ-Li₂IrO₃.[4](https://arxiv.org/pdf/2308.01943) For real materials, the Kitaev-Heisenberg model of Jiří Chaloupka, George Jackeli, and Giniyat Khaliullin ([Physical Review Letters](https://www.edgepedia.org/physical-review-letters) 105, 027204, 2010) interpolates between the Heisenberg and exactly solvable Kitaev limits for layered iridates \\( A_2\\text{IrO}_3 \\) (A = Li, Na); exact diagonalization and spin-wave analysis found an extended spin-liquid phase near the Kitaev limit, a Néel state near the [Heisenberg limit](https://www.edgepedia.org/heisenberg-limit), and a stripy antiferromagnetic state that is the exact ground state at the midpoint between the two limits.[10](https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.105.027204) The \\( J\\text{-}K\\text{-}\\Gamma\\text{-}\\Gamma' \\) model, adding symmetric off-diagonal \\( \\Gamma \\) exchange, is broadly considered the minimal model for many Kitaev materials.[4](https://arxiv.org/pdf/2308.01943)\n\n## Applications\n\nThe model's phases carry quantitative topological signatures. The gapped A phase has the same anyonic structure as the toric code, with Abelian anyons.[3](https://www.thp.uni-koeln.de/trebst/Lectures/Seminar14/Handout9.pdf) • [6](https://thomashodson.com/thesis/2_Background/2.2_HKM_Model.html) The gapless B phase occupies the central region of the phase diagram where \\( |J_x|, |J_y|, |J_z| \\) satisfy the triangle inequalities, with exactly two fermion zeros at \\( \\mathbf{q} = \\pm \\mathbf{q}^* \\) when the inequalities are strict.[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) • [3](https://www.thp.uni-koeln.de/trebst/Lectures/Seminar14/Handout9.pdf) A magnetic field generates, at third order in perturbation theory, a three-spin term \\( \\sigma^x \\cdot \\sigma^z \\cdot \\sigma^y \\) that breaks time-reversal symmetry and acts as a Haldane-type second-neighbor hopping, gapping the Dirac cones and producing chiral Majorana edge modes and a half-quantized thermal [Hall effect](https://www.edgepedia.org/hall-effect).[5](https://jeffrau.ca/assets/lectures/kitaev-2021.pdf) Kitaev classified 16 types of chiral spin liquid by [Chern number](https://www.edgepedia.org/chern-number) \\( \\nu \\): \\( \\mathbb{Z}_2 \\) vortices are Abelian anyons for even \\( \\nu \\) and non-Abelian for odd \\( \\nu \\), and the thermal Hall conductance is \\( \\kappa_{xy} = \\Lambda \\cdot T \\cdot c_- \\) with \\( \\Lambda = \\pi k_B^2 / 6h \\) and \\( c_- = \\nu/2 \\).[11](https://www.nature.com/articles/s41535-024-00704-9) The Abelian and non-Abelian phases of the original model correspond to \\( \\nu = 0 \\) and \\( \\nu = \\pm 1 \\).[1](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381) In topological quantum computing, an odd Chern number guarantees that a flux (vortex) binds an unpaired [Majorana zero mode](https://www.edgepedia.org/majorana-zero-mode) with non-Abelian Ising-anyon statistics.[7](https://www.nature.com/articles/s41586-025-09475-0)\n\nThe Jackeli-Khaliullin mechanism, derived by George Jackeli and Giniyat Khaliullin in 2009 (Physical Review Letters 102, 017205), shows that bond-directional Kitaev exchange is the leading interaction in spin-orbit Mott insulators with edge-sharing octahedra.[17](https://www.semanticscholar.org/reader/8df21f6b54ac5aa28240e7f704cf595e8dcad4aa) • [12](https://doi.org/10.1103/physrevlett.102.017205) Candidate Kitaev materials include α-RuCl₃, Na₂IrO₃, α-, β-, and γ-Li₂IrO₃, H₃LiIr₂O₆, and Cu₂IrO₃.[5](https://jeffrau.ca/assets/lectures/kitaev-2021.pdf) In Na₂IrO₃, α-Li₂IrO₃, and α-RuCl₃ the Kitaev coupling is dominant and ferromagnetic, but \\( \\Gamma \\) is comparable in size to, or can exceed, \\( |K| \\) except in Na₂IrO₃.[4](https://arxiv.org/pdf/2308.01943) Most candidate materials magnetically order at low temperature: Na₂IrO₃ shows a zigzag-ordered transition around \\( T_N \\approx 15 \\) K, and α-RuCl₃ orders at \\( T_N \\approx 7 \\) K yet shows Raman fermionic excitations, a diffuse neutron-scattering continuum, and a reported half-quantized thermal Hall effect in a field-induced state.[13](https://cond-mat.de/events/correl23/manuscripts/trebst.pdf) [Experiments](https://www.edgepedia.org/experiments) indicate that H₃LiIr₂O₆ and α-RuCl₃ in an applied magnetic field show signatures of the QSL state.[14](https://www.osti.gov/biblio/1509555) In 2025, a reconfigurable atom-array processor digitally simulated the Kitaev honeycomb model, preparing the non-Abelian spin-liquid phase and verifying it by measuring a change from \\( C = 0 \\) to \\( C = 1 \\) in the Chern number of the lowest energy band, using a Floquet circuit of depth 6 at the simulated system size.[7](https://www.nature.com/articles/s41586-025-09475-0)\n\n## Limitations and alternatives\n\nThe exact solution holds only at the pure Kitaev point; competing interactions of realistic size destabilize the spin liquid in most materials.[4](https://arxiv.org/pdf/2308.01943) • [15](https://google.iopscience.iop.org/article/10.1088/1361-648X/ad6827/meta) The gapped A phase is physically equivalent to Kitaev's own toric code.[3](https://www.thp.uni-koeln.de/trebst/Lectures/Seminar14/Handout9.pdf) Magnetic fields can tune between phases and have been investigated for potentially counteracting destabilizing interactions and reviving the topological phase.[15](https://google.iopscience.iop.org/article/10.1088/1361-648X/ad6827/meta) The field-induced phases of α-RuCl₃ remain debated: thermal Hall measurements have yielded contradictory results, attributed to sample quality issues such as stacking disorder.[11](https://www.nature.com/articles/s41535-024-00704-9) In the pure ferromagnetic Kitaev model, variational Monte Carlo finds the non-Abelian chiral spin liquid terminating at \\( g \\mu_B B / |K| = 0.19 \\), a critical field smaller than the vison gap \\( \\Delta_v = 0.07K \\), while in the antiferromagnetic model the non-Abelian phase survives up to \\( h_c \\sim 0.3K \\).[11](https://www.nature.com/articles/s41535-024-00704-9) • [16](https://www.thp.uni-koeln.de/trebst/pubs/KitaevMaterialsReview.pdf) For the antiferromagnetic model under field, iPEPS calculations locate two critical fields, \\( h_{c1} \\simeq 0.45 \\) and \\( h_{c2} \\simeq 0.70 \\), separating the chiral spin liquid from an intermediate gapless phase with approximate power-law correlations \\( C(R) \\propto R^{-m} \\), \\( m \\approx 1.5 \\), and then a polarized phase.[8](https://arxiv.org/pdf/2601.14496) The interpretation of the field-induced state in α-RuCl₃ remains actively debated, and neural quantum states appear in the published literature only as an emerging tool for the model.[2](https://journals.aps.org/rmp/abstract/10.1103/3m4m-3v59) • [8](https://arxiv.org/pdf/2601.14496)\n\n## References\n\n1. [Anyons in an exactly solved model and beyond (Kitaev, Annals of Physics 321, Issue 1, January 2006, Pages 2-111)](https://www.sciencedirect.com/science/article/abs/pii/S0003491605002381)\n2. [Kitaev quantum spin liquids (Reviews of Modern Physics 97, 045003, published 3 December 2025)](https://journals.aps.org/rmp/abstract/10.1103/3m4m-3v59)\n3. [Kitaev Honeycomb Model (University of Cologne seminar handout)](https://www.thp.uni-koeln.de/trebst/Lectures/Seminar14/Handout9.pdf)\n4. [Kitaev materials at finite fields and temperatures: J-K-Γ-Γ′ model review (arXiv:2308.01943)](https://arxiv.org/pdf/2308.01943)\n5. [Tutorial: Physics of the Kitaev Model and its Realization in Kitaev Materials](https://jeffrau.ca/assets/lectures/kitaev-2021.pdf)\n6. [Background - The Kitaev Honeycomb Model (PhD thesis chapter)](https://thomashodson.com/thesis/2_Background/2.2_HKM_Model.html)\n7. [Digital quantum simulation of Kitaev's honeycomb model on a reconfigurable atom array (Nature, 2025)](https://www.nature.com/articles/s41586-025-09475-0)\n8. [Report on Progress: magnetic-field-induced phenomena in Z2 quantum spin liquids (arXiv preprint, 2026)](https://arxiv.org/pdf/2601.14496)\n9. [Physics of the Kitaev Model: Fractionalization, Dynamic Correlations, and Material Connections (Annual Review of Condensed Matter Physics)](https://www.annualreviews.org/content/journals/10.1146/annurev-conmatphys-033117-053934)\n10. [Kitaev-Heisenberg Model on a Honeycomb Lattice: Possible Exotic Phases in Iridium Oxides A2IrO3 (Chaloupka, Jackeli, Khaliullin, PRL 105, 027204, 2010)](https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.105.027204)\n11. [Multinode quantum spin liquids in extended Kitaev honeycomb models (npj Quantum Materials, 2024)](https://www.nature.com/articles/s41535-024-00704-9)\n12. [G. Jackeli, G. Khaliullin (2009). Mott Insulators in the Strong Spin-Orbit Coupling Limit: From Heisenberg to a Quantum Compass and Kitaev Models. Physical Review Letters.](https://doi.org/10.1103/physrevlett.102.017205)\n13. [Kitaev Magnets (Trebst lecture notes)](https://cond-mat.de/events/correl23/manuscripts/trebst.pdf)\n14. [Concept and realization of Kitaev quantum spin liquids (OSTI.GOV record)](https://www.osti.gov/biblio/1509555)\n15. [Field tuning Kitaev systems for spin fractionalization and topological order (J. Phys.: Condens. Matter, 2024)](https://google.iopscience.iop.org/article/10.1088/1361-648X/ad6827/meta)\n16. [Kitaev materials (Trebst & Catuneanu, Physics Reports 2022, author-hosted PDF)](https://www.thp.uni-koeln.de/trebst/pubs/KitaevMaterialsReview.pdf)\n17. [8df21f6b54ac5aa28240e7f704cf595e8dcad4aa (semanticscholar.org)](https://www.semanticscholar.org/reader/8df21f6b54ac5aa28240e7f704cf595e8dcad4aa)\n\n---\n*Topic: Encyclopedia › Physical world and mathematics › Physics › Quantum physics*\n\n*Initially written Sep 29, 2026 · Reviewed: Sep 30, 2026 · Edited: Sep 30, 2026; Oct 15, 2026 · Last review: Sep 30, 2026*\n\n*Copyright 2026 EdgeChat AI, a subsidiary of Biostate AI.*\n\nLicense: Edgepedia Community License 1.0, https://www.edgepedia.org/edgepedia/license\n", "same_as": [], "url": "https://www.edgepedia.org/kitaev-honeycomb-model", "markdown_url": "https://www.edgepedia.org/kitaev-honeycomb-model.md", "license": { "name": "Edgepedia Community License 1.0", "url": "https://www.edgepedia.org/edgepedia/license", "summary": "Free with credit, commercial use included. AI training is open to everyone. For other uses, organizations over USD 100M in revenue or 100M monthly users license separately.", "spdx": "LicenseRef-Edgepedia-Community-1.0" }, "credit": "\"Kitaev honeycomb model\", Edgepedia (EdgeChat), https://www.edgepedia.org/kitaev-honeycomb-model. Edgepedia Community License 1.0.", "credit_md": "\"[Kitaev honeycomb model](https://www.edgepedia.org/kitaev-honeycomb-model)\", Edgepedia (EdgeChat), [https://www.edgepedia.org/kitaev-honeycomb-model](https://www.edgepedia.org/kitaev-honeycomb-model). [Edgepedia Community License 1.0](https://www.edgepedia.org/edgepedia/license).", "credit_html": "\"Kitaev honeycomb model\", Edgepedia (EdgeChat), https://www.edgepedia.org/kitaev-honeycomb-model. Edgepedia Community License 1.0.", "speakable": "The Kitaev honeycomb model is an exactly solvable spin model on a honeycomb lattice, introduced by Alexei Kitaev in 2006, whose ground state is a topological quantum spin liquid of emergent Majorana fermions." }