Long-range spin correlations in a honeycomb spin model with a magnetic field

We consider the spin-1/2 model on the honeycomb lattice [A. Kitaev, Ann. Phys. 321, 2 (2006)] in the presence of a weak magnetic field h _α ≪ J. Such a perturbation treated in the lowest nonvanishing order over h _α leads [K.S. Tikhonov, M.V. Feigel'man, and A.Yu. Kitaev, Phys. Rev. Lett. 106, 067203 (2011)] to a powerlaw decay of irreducible spin correlations 《 s ^z ( t, r) s ^z (0, 0)》 ∝ h _z ² f( t, r), where f( t, r) ∝ [max( t, Jr)]^-4. We have studied the effects of the next order of perturbation in h _z and found an additional term of the order h _z ⁴ in the correlation function 《 s ^z ( t, r) s ^z (0, 0)》 which scales as h _z ⁴ cosγ/ r ³ at Jt≪ r, where γ is the polar angle in the 2D plane. We demonstrate that such a contribution can be understood as a result of a perturbation of the effective Majorana Hamiltonian by the weak imaginary vector potential A _x ∝ i h _z ².