Monte Carlo study of low/high spin transitions

We study the finite temperature property of a model on two dimensional square lattices with two Ising spins at each lattice site by Monte Carlo simulations. When those Ising spins at a lattice site are parallel the site is said to be in the high-spin state (HS), while when they are antiparallel the site is said to be in the low-spin state (LS). Throughout the study, the energy of HS is presumed to be higher than that of LS. Two Ising spins at each site are added to form a total spin, which interacts with its nearest neighbour total spins via spin-spin couplings. The spin-phonon coupling also is introduced via harmonic springs between nearest neighbour sites with spring constants and equilibrium distances depending on the spin states of the sites involved. In this system, we investigate the feature of transitions between LS and HS (to be called low/high spin transition (LHST)) by varying the temperature. As for the ferromagnetic interaction between total spins, the second order phase transition: pure HSmixed state of HS and LS is possible to occur in a pure spin system, as is expected from mean field calculations. The role of lattice distortions by the change of lattice spacings is shown to be essential for LHST: pure LS(pure)HS. In the model investigated, there appears an indication of the strong first order phase transition which reveals a conspicuous hysteresis.