We study the two-sublattice model in the mean field theory by expanding the Gibbs free energy in terms of the magnetizations M1 (Mup) and M2 (Mdown) with the quadratic coupling M12M22 (quadrupolar interactions) for the ...