Effective Hamiltonian and Exact Solutions for Two Anyons in the Kitaev Toric Code

We study the spectrum of two quasiparticle excitations in the Kitaev toric code model subjected to an external magnetic field. Using the perturbative continuous unitary transformation (PCUT) method, we derive an effective Hamiltonian that captures the low-energy behavior of the system. The spin degrees of freedom are explicitly mapped onto creation and annihilation operators corresponding to each anyon type. We address the subtleties associated with implementing braiding statistics through a careful construction of the basis states. For magnetic fields applied along the x or z directions, the system supports either e or m type excitations, which we describe as hardcore bosons and solve numerically. In the case where both e and m excitations coexist, we obtain an exact solution to the two-anyon problem, fully accounting for the braiding phase and its impact on the excitation spectrum.