text: Great eﬀorts have been invested to better understand 
  the properties of the Kitaev model. For example, sev- 
  eral groups have pointed out that the fractionalized Ma- 
  jorana fermion excitations may be understood from the 
  more familiar Jordan-Wigner transformation of 1D spin 
  systems2,3. The analogy between the non-Abelian Ising 
  vortices and vortices in p + ip superconductors has been 
  raised in serveral works4–7. Exact diagonalization has 
  been used to study the Kitaev model on small lattices8. 
  And perturbative expansion methods have been devel- 
  oped to study the gapped phases of the Kitaev-type 
  models9. 
  Kitaev’s exactly solvable spin-1/2 honeycomb lattice 
  model1 (noted as the Kitaev model hereafter) has in- 
  spired great interest since its debut, due to its exact 
  solvability, fractionalized excitations, and the potential  Many generalizations of the Kitaev model have been 