Spin model for the Honeycomb $\rm NiPS_3$

TL;DR

This paper develops a microscopic spin Hamiltonian for NiPS3, revealing key exchange interactions and magnetic ground states, which helps explain experimental observations in bulk and monolayer forms.

Contribution

It introduces a minimal J1-J3-K1 spin model derived from a Hubbard framework, elucidating the magnetic interactions in NiPS3.

Findings

Ferromagnetic J1 and antiferromagnetic J3 in bulk NiPS3

Ground state is zig-zag antiferromagnetic order

Potential for spin spiral states in monolayer NiPS3

Abstract

In the Van der Waal material $\rm NiPS_3$, Ni atoms have spin S=1 and realize a honeycomb lattice. Six sulfur atoms surround each Ni and split their d manifold into three filled and two unfilled bands. Aimed to determine the spin Hamiltonian of $\rm NiPS_3$, we study its exchange mechanisms using a two-band half-filled Hubbard model. Hopping between d orbitals is mediated by p orbitals of sulfur and gives rise to bilinear and biquadratic spin couplings in the limit of strong electronic correlations. The microscopic model exposed a ferromagnetic biquadratic spin interaction $\rm K_1$ allowing the completion of a minimal $\rm J_1-J_3-K_1$ spin Hamiltonian for $\rm NiPS_3$. In bulk, a ferromagnetic first nearest neighbor $\rm J_1$ and a more significant antiferromagnetic third nearest neighbor spin coupling $\rm J_3$ agreed with the literature, while in monolayer $\rm J_1$ is positive and very small in comparison. Using a variational scheme we found that a zig-zag antiferromagnetic order is the ground state of bulk samples. The zig-zag pattern is adjacent to commensurate and incommensurate spin spirals, which could hint at the puzzling results reported in $\rm NiPS_3$ monolayers.