Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1605.04315 · Physical review. B./Physical review. B

Kitaev honeycomb tensor networks: Exact unitary circuits and applications

Philipp Schmoll, Román Orús

Published Jan 10, 2017• 24 citations
arXivHTMLPDF

Abstract

The Kitaev honeycomb model is a paradigm of exactly solvable models, showing nontrivial physical properties such as topological quantum order, Abelian and non-Abelian anyons, and chirality. Its solution is one of the most beautiful examples of the interplay of different mathematical techniques in condensed matter physics. In this paper, we show how to derive a tensor network (TN) description of the eigenstates of this spin-1/2 model in the thermodynamic limit, and in particular for its ground state. In our setting, eigenstates are naturally encoded by an exact 3d TN structure made of fermionic unitary operators, corresponding to the unitary quantum circuit building up the many-body quantum state. In our derivation we review how the different ``solution ingredients'' of the Kitaev honeycomb model can be accounted for in the TN language, namely, Jordan-Wigner transformation, braidings of Majorana modes, fermionic Fourier transformation, and Bogoliubov transformation. The TN built in this way allows for a clear understanding of several properties of the model. In particular, we show how the fidelity diagram is straightforward both at zero temperature and at finite temperature in the vortex-free sector. We also show how the properties of two-point correlation functions follow easily. Finally, we also discuss the pros and cons of contracting of our 3d TN down to a 2d projected entangled pair state (PEPS) with finite bond dimension. The results in this paper can be extended to generalizations of the Kitaev model, e.g., to other lattices, spins, and dimensions.

DOI
10.1103/physrevb.95.045112

Related papers

Complete Entanglement Structure of the Kitaev Honeycomb Spin Liquid from Exact Tensor Networks

Hidehiro Saito, Chisa Hotta

While most of the features of the Kitaev honeycomb model are already clarified from the combined ${\mathbb Z}_2$ plaquette and Majorana fermion picture, we show that retrieving the exact quantum many-body form gives a clue to fix the remaining issues on its spin liquid property. We provide an exact two-dimensional tensor network description of a full set of solutions from a random product state, by sequential quantum-number projection and Gutzwiller projection, both expressed as local matrix-product operators of dimension $D=2$. The ground and excited states share the same bond dimension yet exhibit distinct Schmidt spectra, distinguishing their area-law and volume-law entanglement. We demonstrate that the topological entanglement entropy emerges {\it for all eigenstates only after} the Gutzwiller projection onto the physical spin Hilbert space, where the gauge constraint removes exactly half of the Schmidt states. This establishes that excited states are likewise topologically ordered. We further show that there exists a projection-induced collapse of the wave function into a single gauge sector in the Parton mean-field states that is known to give a good description of the spin liquids. The demonstrations on the $S=1/2$ kagome and triangular Heisenberg models suggest that the present perspective is applied to a wider class of models hosting ${\mathbb Z}_2$ spin liquids. Our finding visualizes a bulk-boundary correspondence with the Wilson-loop characterization of topological order.

arXiv:2610.05395•Oct 4, 2026•0 citations

Thermodynamics of three-dimensional Kitaev quantum spin liquids via tensor networks

Saeed S. Jahromi, Hadi Yarloo, Román Orús

We study the 3D Kitaev and Kitaev-Heisenberg models, respectively, on the hyperhoneycomb and hyperoctagon lattices, both at zero and finite-temperature, in the thermodynamic limit. Our analysis relies on advanced tensor network (TN) simulations based on graph projected entangled-pair states (gPEPS). We map out the TN phase diagrams of the models and characterize their underlying gapped and gapless phases both at zero and finite temperature. In particular, we demonstrate how cooling down the hyperhoneycomb system from high temperature leads to fractionalization of spins to Majorana fermions and gauge fields that occurs in two separate temperature regimes, leaving their fingerprint on specific heat as a double-peak feature as well as on other quantities such as the thermal entropy, spin-spin correlations, and bond entropy. Using the Majorana representation of the Kitaev model, we further show that the low-temperature thermal transition to the Kitaev quantum spin liquid (QSL) phase is associated with the nontrivial Majorana band topology and the presence of Weyl nodes, which manifests itself via nonvanishing Chern number and finite thermal Hall conductivity. Beyond the pure Kitaev limit, we study the 3D Kitaev-Heisenberg (KH) model on the hyperoctagon lattice and extract the full phase diagram for different Heisenberg couplings. We further explore the thermodynamic properties of the magnetically-ordered regions in the KH model and show that, in contrast to the QSL phase, here the thermal phase transition follows the standard Landau symmetry-breaking theory.

arXiv:2011.11577•Physical Review Research•Sep 2, 2021•13 citations

Spin–Spin Correlations in the Kitaev Model at Finite Temperatures: Approximate and Exact Results via Green’s Function Equation of Motion

Hibiki Takegami, Takao Morinari

The Kitaev model, defined on a honeycomb lattice, features an exactly solvable ground state with fractionalized Majorana fermion excitations, which can potentially form non-Abelian anyons crucial for fault-tolerant topological quantum computing. Although Majorana fermions are essential for obtaining the exact ground state, their physical interpretation in terms of spin operators remains unclear. In this study, we employ a Green's function approach that maintains SU(2) symmetry to address this issue and explore the model's finite temperature properties. Our results demonstrate that the computed temperature dependence of the correlation functions closely approximates the exact values at zero temperature, confirming the accuracy of our method. We also present several exact results concerning the spin Green's function and spin-spin correlation functions that are specific to the Kitaev model.

arXiv:2405.13309•Journal of the Physical Society of Japan•Jul 17, 2024•1 citation

The Kitaev Honeycomb Model: From Bogoliubov-de Gennes Transformations to the Quantum Spin Liquid Phase Diagram

Demetris Demetriades

These lecture notes focus on the Kitaev honeycomb model, a two-dimensional quantumsystem with three distinct Ising-type interactions, depending on the direction of eachbond. This model has proven particularly important for describing quantum spin liquids,as well as for applications in topological quantum computing, due to its exotictopological properties. It constitutes the first exactly solvable theoretical model thatfully captures the behavior of a quantum spin liquid. We begin by presenting the geometryof the model, the Hamiltonian that describes the system, and Alexei Kitaev's effort toobtain an exact analytical solution, based on the decomposition of the Hilbert spaceinto flux sectors. We then present Ettore Majorana's observation that the Dirac equationalso admits real solutions, leading to the theoretical prediction of Majorana fermions-- particles that are identical to their antiparticles. We examine how the spatialseparation of Majorana operators leads to the rise of unbound Majorana fermions, andshowcase the exotic topological properties that manifest in physical systems such assuperconducting wires. Next, we represent the Hamiltonian in the basis of Majoranafermions, by expressing spin operators using four Majorana operators. The transitionfrom the extended space $\tilde{\mathcal{L}}$ to the physical subspace $\mathcal{L}$is achieved through a gauge transformation of the $\mathbb{Z}_2$ group. We thenreformulate the Hamiltonian in the basis of complex bond and matter fermions, transitionto momentum space, and assume that the system lies in the zero-flux sector. Through theBogoliubov--de Gennes transformation, we arrive at the final diagonalized Hamiltonianin the quasiparticle basis. Finally, we represent the possible phases -- the gaplessphase (phase $B$) and the gapped phases (phase $A_{\alpha_i}$) -- compute thecorresponding energy spectrum, and present the complete phase diagram. By employing thecorrelation function, we show that all phases are characterized by strong quantumentanglement and spin fractionalization, classifying the system as a quantum spin liquid.

•Mar 24, 2026•0 citations

Tensor network wave function of S=1 Kitaev spin liquids

Hyun-Yong Lee, Naoki Kawashima, Yong Baek Kim

The spin-1/2 Kitaev model offers the exactly solvable example of quantum spin liquids.Possible material realizations of the spin-1/2 Kitaev systems and the prospect of using the Majorana fermion excitations for quantum computations have revolutionized quantum spin liquids research.Recently it has been suggested that higher-spin, especially spin-1, Kitaev exchange interactions can be realized in a variety of materials.Numerical computations on small clusters indicate that the ground state of the spin-1 Kitaev model may also be a quantum spin liquid.On the other hand, the nature of the ground state remains elusive since the spin-1 model is not exactly solvable in contrast to the spin-1/2 model.In this work, using the quantum-entanglement based tensor network approach, we construct an explicit ground-state wave function for the spin-1 Kitaev model, which is written only in terms of physical spin operators.We establish the existence of distinct topological sectors on a torus by constructing the minimally entangled states in the degenerate ground-state manifold and evaluating topological entanglement entropy.Our results suggest that the ground state of the spin-1 Kitaev model is a gapped quantum spin liquid with Z 2 gauge structure and Abelian quasiparticles.We explain the subtle differences between the spin-1/2 and spin-1 Kitaev quantum spin liquids.

arXiv:1911.07714•Physical Review Research•Aug 27, 2020•51 citations

Lattice Defects in the Kitaev Honeycomb Model

John Brennan, Jiří Vala

The Kitaev honeycomb lattice system is an important model of topological materials whose phase diagram exhibits both abelian and non-abelian topological phases. The latter, a so-called Ising phase, is related to topological superconductors. Its quasiparticle excitations, which are formed by Majorana fermions attached to vortices, show non-abelian fractional statistics and are known as Ising anyons. We investigate dislocation defects in the Ising phase of the Kitaev honeycomb model. After introducing them to the system, we accordingly generalize our solution of this model to the situation with the defects. The important part of this effort is developing an appropriate Jordan-Wigner fermionization procedure. It is expected that the presence of defects manifests itself by the formation of fermionic zero-energy modes around the defect end points. We numerically confirm this expectation and further investigate properties of these modes. The computational potential of our technique is demonstrated for both diagonalization and dynamical simulations. The latter focuses on the process of fusion of the vortex zero-energy modes with the Majorana fermions attached to the defect. This process simulates fusion of non-abelian Ising anyons.

The Journal of Physical Chemistry A•Feb 17, 2016•13 citations
The Kitaev honeycomb model is a paradigm of exactly solvable models, showing nontrivial physical properties such as topological quantum order, Abelian and non-Abelian anyons, and chirality. Its solution is one of the most beautiful examples of the interplay of different mathematical techniques in condensed matter physics. In this paper, we show how to derive a tensor network (TN) description of the eigenstates of this spin-1/2 model in the thermodynamic limit, and in particular for its ground state. In our setting, eigenstates are naturally encoded by an exact 3d TN structure made of fermionic unitary operators, corresponding to the unitary quantum circuit building up the many-body quantum state. In our derivation we review how the different ``solution ingredients'' of the Kitaev honeycomb model can be accounted for in the TN language, namely, Jordan-Wigner transformation, braidings of Majorana modes, fermionic Fourier transformation, and Bogoliubov transformation. The TN built in this way allows for a clear understanding of several properties of the model. In particular, we show how the fidelity diagram is straightforward both at zero temperature and at finite temperature in the vortex-free sector. We also show how the properties of two-point correlation functions follow easily. Finally, we also discuss the pros and cons of contracting of our 3d TN down to a 2d projected entangled pair state (PEPS) with finite bond dimension. The results in this paper can be extended to generalizations of the Kitaev model, e.g., to other lattices, spins, and dimensions.
While most of the features of the Kitaev honeycomb model are already clarified from the combined Z2{\mathbb Z}_2Z2​ plaquette and Majorana fermion picture, we show that retrieving the exact quantum many-body form gives a clue to fix the remaining issues on its spin liquid property. We provide an exact two-dimensional tensor network description of a full set of solutions from a random product state, by sequential quantum-number projection and Gutzwiller projection, both expressed as local matrix-product operators of dimension D=2D=2D=2. The ground and excited states share the same bond dimension yet exhibit distinct Schmidt spectra, distinguishing their area-law and volume-law entanglement. We demonstrate that the topological entanglement entropy emerges {\it for all eigenstates only after} the Gutzwiller projection onto the physical spin Hilbert space, where the gauge constraint removes exactly half of the Schmidt states. This establishes that excited states are likewise topologically ordered. We further show that there exists a projection-induced collapse of the wave function into a single gauge sector in the Parton mean-field states that is known to give a good description of the spin liquids. The demonstrations on the S=1/2S=1/2S=1/2 kagome and triangular Heisenberg models suggest that the present perspective is applied to a wider class of models hosting Z2{\mathbb Z}_2Z2​ spin liquids. Our finding visualizes a bulk-boundary correspondence with the Wilson-loop characterization of topological order.
We study the 3D Kitaev and Kitaev-Heisenberg models, respectively, on the hyperhoneycomb and hyperoctagon lattices, both at zero and finite-temperature, in the thermodynamic limit. Our analysis relies on advanced tensor network (TN) simulations based on graph projected entangled-pair states (gPEPS). We map out the TN phase diagrams of the models and characterize their underlying gapped and gapless phases both at zero and finite temperature. In particular, we demonstrate how cooling down the hyperhoneycomb system from high temperature leads to fractionalization of spins to Majorana fermions and gauge fields that occurs in two separate temperature regimes, leaving their fingerprint on specific heat as a double-peak feature as well as on other quantities such as the thermal entropy, spin-spin correlations, and bond entropy. Using the Majorana representation of the Kitaev model, we further show that the low-temperature thermal transition to the Kitaev quantum spin liquid (QSL) phase is associated with the nontrivial Majorana band topology and the presence of Weyl nodes, which manifests itself via nonvanishing Chern number and finite thermal Hall conductivity. Beyond the pure Kitaev limit, we study the 3D Kitaev-Heisenberg (KH) model on the hyperoctagon lattice and extract the full phase diagram for different Heisenberg couplings. We further explore the thermodynamic properties of the magnetically-ordered regions in the KH model and show that, in contrast to the QSL phase, here the thermal phase transition follows the standard Landau symmetry-breaking theory.
The Kitaev model, defined on a honeycomb lattice, features an exactly solvable ground state with fractionalized Majorana fermion excitations, which can potentially form non-Abelian anyons crucial for fault-tolerant topological quantum computing. Although Majorana fermions are essential for obtaining the exact ground state, their physical interpretation in terms of spin operators remains unclear. In this study, we employ a Green's function approach that maintains SU(2) symmetry to address this issue and explore the model's finite temperature properties. Our results demonstrate that the computed temperature dependence of the correlation functions closely approximates the exact values at zero temperature, confirming the accuracy of our method. We also present several exact results concerning the spin Green's function and spin-spin correlation functions that are specific to the Kitaev model.
These lecture notes focus on the Kitaev honeycomb model, a two-dimensional quantumsystem with three distinct Ising-type interactions, depending on the direction of eachbond. This model has proven particularly important for describing quantum spin liquids,as well as for applications in topological quantum computing, due to its exotictopological properties. It constitutes the first exactly solvable theoretical model thatfully captures the behavior of a quantum spin liquid. We begin by presenting the geometryof the model, the Hamiltonian that describes the system, and Alexei Kitaev's effort toobtain an exact analytical solution, based on the decomposition of the Hilbert spaceinto flux sectors. We then present Ettore Majorana's observation that the Dirac equationalso admits real solutions, leading to the theoretical prediction of Majorana fermions-- particles that are identical to their antiparticles. We examine how the spatialseparation of Majorana operators leads to the rise of unbound Majorana fermions, andshowcase the exotic topological properties that manifest in physical systems such assuperconducting wires. Next, we represent the Hamiltonian in the basis of Majoranafermions, by expressing spin operators using four Majorana operators. The transitionfrom the extended space L~\tilde{\mathcal{L}}L~ to the physical subspace L\mathcal{L}Lis achieved through a gauge transformation of the Z2\mathbb{Z}_2Z2​ group. We thenreformulate the Hamiltonian in the basis of complex bond and matter fermions, transitionto momentum space, and assume that the system lies in the zero-flux sector. Through theBogoliubov--de Gennes transformation, we arrive at the final diagonalized Hamiltonianin the quasiparticle basis. Finally, we represent the possible phases -- the gaplessphase (phase BBB) and the gapped phases (phase AαiA_{\alpha_i}Aαi​​) -- compute thecorresponding energy spectrum, and present the complete phase diagram. By employing thecorrelation function, we show that all phases are characterized by strong quantumentanglement and spin fractionalization, classifying the system as a quantum spin liquid.
The spin-1/2 Kitaev model offers the exactly solvable example of quantum spin liquids.Possible material realizations of the spin-1/2 Kitaev systems and the prospect of using the Majorana fermion excitations for quantum computations have revolutionized quantum spin liquids research.Recently it has been suggested that higher-spin, especially spin-1, Kitaev exchange interactions can be realized in a variety of materials.Numerical computations on small clusters indicate that the ground state of the spin-1 Kitaev model may also be a quantum spin liquid.On the other hand, the nature of the ground state remains elusive since the spin-1 model is not exactly solvable in contrast to the spin-1/2 model.In this work, using the quantum-entanglement based tensor network approach, we construct an explicit ground-state wave function for the spin-1 Kitaev model, which is written only in terms of physical spin operators.We establish the existence of distinct topological sectors on a torus by constructing the minimally entangled states in the degenerate ground-state manifold and evaluating topological entanglement entropy.Our results suggest that the ground state of the spin-1 Kitaev model is a gapped quantum spin liquid with Z 2 gauge structure and Abelian quasiparticles.We explain the subtle differences between the spin-1/2 and spin-1 Kitaev quantum spin liquids.
The Kitaev honeycomb lattice system is an important model of topological materials whose phase diagram exhibits both abelian and non-abelian topological phases. The latter, a so-called Ising phase, is related to topological superconductors. Its quasiparticle excitations, which are formed by Majorana fermions attached to vortices, show non-abelian fractional statistics and are known as Ising anyons. We investigate dislocation defects in the Ising phase of the Kitaev honeycomb model. After introducing them to the system, we accordingly generalize our solution of this model to the situation with the defects. The important part of this effort is developing an appropriate Jordan-Wigner fermionization procedure. It is expected that the presence of defects manifests itself by the formation of fermionic zero-energy modes around the defect end points. We numerically confirm this expectation and further investigate properties of these modes. The computational potential of our technique is demonstrated for both diagonalization and dynamical simulations. The latter focuses on the process of fusion of the vortex zero-energy modes with the Majorana fermions attached to the defect. This process simulates fusion of non-abelian Ising anyons.