Coupled Majorana modes in a dual vortex of the Kitaev honeycomb model
Phys. Rev. B 114, 125137 – Published 31 August, 2026
DOI: https://doi.org/10.1103/qhsx-q63d
Abstract
The Kitaev model is exactly solvable in terms of Majorana fermions hopping on a honeycomb lattice and coupled to a static gauge field, giving the possibility of vortices in hexagonal plaquettes. In the vortex-full sector and in the presence of a time-reversal–breaking three-spin term of strength , the energy spectrum is gapped, and the ground state possesses an even Chern number. An isolated vortex-free plaquette acts as a dual vortex and binds a fermionic mode at finite energy in the bulk gap. This mode is equivalent to two coupled Majorana zero modes located on the same dual vortex. In a continuum approximation, we analytically compute the Majorana wave functions and their coupling in the two limits of small or large . The analytical approach is confirmed by numerical perturbation theory directly on the lattice. The latter is in excellent agreement with the full numerics on a finite-size system. We contrast our results with states bound to an isolated vortex in a topological superconductor with even Chern number.