Realistic Spin Model for Multiferroic NiI_{2}

A realistic first-principle-based spin Hamiltonian is constructed for the type-II multiferroic NiI_{2}, using a symmetry-adapted cluster expansion method. Besides single ion anisotropy and isotropic Heisenberg terms, this model further includes the Kitaev interaction and a biquadratic term, and can well reproduce striking features of the experimental helical ground state, that are, e.g., a proper screw state, canting of rotation plane, propagation direction, and period. Using this model to build a phase diagram, it is demonstrated that, (i) the in-plane propagation direction of ⟨11[over ¯]0⟩ is determined by the Kitaev interaction, instead of the long-believed exchange frustrations and (ii) the canting of rotation plane is also dominantly determined by Kitaev interaction, rather than interlayer couplings. Furthermore, additional Monte Carlo simulations reveal three equivalent domains and different topological defects. Since the ferroelectricity is induced by spins in type-II multiferroics, our work also implies that Kitaev interaction is closely related to the multiferroicity of NiI_{2}.