try with 2D Heisenberg
Hamiltonian checks out,

i.e. 1/2 \sum_k \lambda_k v^2_k = H

so the AFMPO's are correct