Monte Carlo study of the Spin-1 Baxter-Wu model in a crystal field

The two-dimensional Baxter-Wu model with spin-1, in the presence of a crystal field, is studied by using Monte Carlo simulations. The standard single-spin-flip Metropolis algorithm is used to generate the configurations from which the order parameter, specific heat and magnetic susceptibility are measured. Single histogram techniques are used to get the results close to the phase transitions. The finite-size scaling procedure is employed in order to obtain the critical behavior. The simulations have shown that the critical exponents are different from the spin-1/2 case and are crystal field dependent.