Nonlinear response of the Kitaev honeycomb lattice model in a weak magnetic field

We investigate the nonlinear response of the Kitaev honeycomb lattice model in a weak magnetic field using the theory of two-dimensional coherent spectroscopy. We observe that at the isotropic point in the non-Abelian phase of this model, the nonlinear spectrum in the 2D frequency domain consists of sharp signals that originate from the flux excitations and Majorana bound states. Signatures of different flux excitations can be clearly observed in this spectrum, such that one can observe evidences of flux states with 4-adjacent, 2-non-adjacent, and 4-far-separated fluxes, which are not visible in linear response spectroscopy such as neutron scattering experiments. Moreover, in the Abelian phase we perceive that the spectrum in the frequency domain is composed of streak signals. These signals, as in the nonlinear response of the pure Kitaev model, represent a distinct signature of itinerant Majorana fermions. However, deep in the Abelian phase whenever a Kitaev exchange coupling is much stronger than the others, the streak signals are weakened and only single sharp spots are seen in the response, which resembles the dispersionless response of the conventional toric code.