Solution (source code)

= Solution

Let $\eta_i=+1$ and $-1$ on the two sublattices of the square lattice. A unit-vector decomposition that separates staggered and uniform magnetization is
$$
\boxed{
\mathbf n_i
=\eta_i\widetilde{\mathbf n}(\mathbf x_i)
\sqrt{1-|\mathbf m(\mathbf x_i)|^2}
+\mathbf m(\mathbf x_i),}
$$
where
$$
|\widetilde{\mathbf n}|^2=1,
\qquad
\mathbf m\mathbin\cdot\widetilde{\mathbf n}=0,
\qquad
|\mathbf m|\ll1.
$$
The alternating part is the <Néel order parameter>, while $\mathbf m$ is the slowly varying <uniform magnetization> generated by canting the two sublattices.