Re-entrant phase transitions of the Blume-Emery-Griffiths model

We have investigated the spin-1 Ising model on the simple cubic lattice with bilinear, biquadratic interaction and anisotropic energy (BEG model). We have been specially interested in the case of antiferro biquadratic interaction, because the interaction will cause the competition with bilinear interaction and anisotropy. A two-sublattice ordering, so called the staggered quadrupole (SQ) phase, occurs as long as biquadratic interaction is negative large enough. We have obtained a full phase diagram in the whole interaction parameter space (for the positive bilinear interaction) by the Bethe approximation, and found several kinds of phase transitions, such as successive, re-entrant and double re-entrant transitions. These transitions are also confirmed by Monte Carlo simulations on simple cubic lattices.