Abstract Motivated by experiments on spin–orbit coupled magnets with Kitaev exchange in magnetic fields, we present an analysis of the classical Kitaev honeycomb model in the presence of a magnetic field. We show that there is a spin liquid regime that exists within a finite window of fields from zero up to a finite threshold before transitioning into the polarized paramagnet. We uncover constraints that spins need to satisfy in the ground state and show that they determine the exact limiting zero temperature behavior of the heat capacity and magnetic susceptibility within the spin liquid as a function of field. When the field is finite, both the two-point spin and the quadrupolar correlations are short-ranged, in contrast to the zero-field case. We rationalize an effective mass for the quadrupolar correlations in terms of a coarse-grained theory with fluctuating effective charge degrees of freedom. Finally, we show that weak site-dilution does not change the magnetization within the spin liquid – a kind of “perfect” compensation of the site dilution.
McClarty et al. (2026) studied this question.