Abstract
Using the Monte Carlo Simulation (MCS) method combined with the Metropolis algorithm, we explore the magnetic properties of a two-dimensional (2D) spin system with spins 1, 2 and 3. The adapted simulation to the Ising model aims to study the impact of the correction term a of the interaction with second neighbors and the crystal field D on the critical temperature Tc of the 2D anisotropic square lattice for different integer values of spins. The obtained results demonstrate that the critical temperature Tc is significantly affected by the competing interactions between the anisotropy parameter, second-neighbor interaction and the crystal field. The considered model exhibits a critical phase transition regardless of the spin S value where the critical temperature increases with the value of S. Note that this phase transition occurs under the condition D≥−2. Moreover, the obtained results show that the complex interactions give rise to a variety of phase diagrams, thereby revealing the richness of possible magnetic behaviors in the studied system.
