Competition between two- and three-sublattice ordering for S=1 spins on the square lattice

We provide strong evidence that the S = 1 bilinear-biquadratic Heisenberg model with nearest-neighbor interactions on the square lattice possesses an extended three-sublattice phase induced by quantum fluctuations for sufficiently large biquadratic interactions, in spite of the bipartite nature of the lattice. The argumentation relies on exact diagonalizations of finite clusters and on a semiclassical treatment of quantum fluctuations within linear flavor-wave theory. In zero field, this three-sublattice phase is purely quadrupolar, and upon increasing the field it replaces most of the plateau at 1/2 that is predicted by the classical theory.