Anatomy of fermionic entanglement and criticality in Kitaev spin liquids

We analyze in detail the effect of nontrivial band topology on the area-law behavior of the entanglement entropy in Kitaev's honeycomb model. By mapping the translationally invariant 2D spin model onto 1D fermionic subsystems, we identify those subsystems responsible for universal entanglement contributions in the gapped phases and those responsible for critical entanglement scaling in the gapless phases. For the gapped phases, we analytically show how the topological edge states contribute to the entanglement entropy and provide a universal lower bound for it. For the gapless semimetallic phases and topological phase transitions, the identification of the critical subsystems shows that they fall always into the Ising or the XY universality classes. As our study concerns the fermionic degrees of freedom in the honeycomb model, qualitatively similar results are expected to apply also to generic topological insulators and superconductors.