Monte Carlo simulations of the spin-2 Blume-Emery-Griffiths model with four-spin interactions

The magnetic properties of a spin S = 2 Ising system with bilinear exchange interaction J1, the biquadratic exchange interaction K, four-spin exchange interactions J4 and crystal field Δ are discussed using the Monte Carlo simulation. The lattice is divided into two sublattices: A and B, for which we compute the magnetizations mA and mB. The phase obtained diagrams of this system are deduced in the planes: (T, Δ/J1), (K/J1, Δ/J1), (Δ/J1, J4/J1) and (J4/J1, K/J1). In addition to the usual phases, we found a new phase called nonmagnetic quadratic, for which the magnetizations are mA ≠ mB and the quadrupolar moments are so that are qA = qB. Furthermore, the behavior of the magnetizations as a function of temperature, crystal field, four-spin exchange interactions and biquadratic exchange interaction are deduced.