INTERNATIONAL JOURNAL OF MODERN PHYSICS C, cilt.14, sa.10, ss.1305-1320, 2003 (SCI-Expanded)
The two-dimensional antiferromagnetic spin-1 Ising model with positive biquadratic interaction is simulated on a cellular automaton which based on the Creutz cellular automaton for square lattice. Phase diagrams characterizing phase transition of the model are presented for a comparison with those obtained from other calculations. We confirm the existence of the intermediate phase observed in previous works for some values of J/K and D/K. The values of the static critical exponents ( beta, gamma and upsilon) are estimated within the framework of the finite-size scaling theory for D/K < 2J/K. Although the results are compatible with the universal Ising critical behavior in the region of D/K < 2J/K - 4, the model does not exhibit any universal behavior in the interval 2J/K - 4 < D/K < 2J/K.