We study the spin S=1 antiferromagnetic Heisenberg model on the square lattice with, in addition to the nearest-neighbor interaction, a three-site interaction of the form (S_i S_j)*(S_j S_k) + h.c.. This interaction appears naturally in a strong coupling exansion of the two-orbital, half-filled Hubbard model. For spin 1/2, this model... Show more