Phase diagrams of the FCC Blume-Emery-Griffiths model on a cellular automaton

The spin-1 Ising (BEG) model has been simulated on a cellular automaton improved from the Creutz cellular automaton (CCA) for a face-centered cubic lattice. The simulations have been made in the 0 <= d = D/J <= 7 and -1.25 < k = KIJ <= 0 parameter region. In this region, the ground state diagram (k, d) has ferromagnetic and perfect zero ordering regions. The ferromagnetic ordering region separates into four regions which exhibit different phase transition types as the first order, the second order, the reentrant, the double-reentrant and the successive phase transitions. The simulation results show that the model has the tricritical points, the critical end points and the bicritical points on the (kT(C)/zJ, d) and (kT(C)/zJ, k) planes as indicated by the Mean Field approximation results.