A comprehensive study of the phase transitions of the frustrated J1-J2 Ising model on the square lattice

We study the phase transitions of the two-dimensional frustrated J1-J2 Ising model on the square lattice with Monte Carlo simulations. Here J1 < 0 and J2 > 0 are the couplings for the nearest and next nearest neighbor Ising spins, respectively. In addition, we consider J1 = −1 and J2/|J1| ≡ g = 0.56, 0.58, 0.60, 0.62, 0.64 and 0.7 in our investigation. Several physical quantities including the Binder ratio U, the specific heat Cv, and the energy density E are obtained. The analyses of these calculated observables indicate that the considered phase transitions for g ≥ 0.58 are likely second order, and the nature of the phase transition associated with g = 0.56 needs further investigation. A comparison between the histograms of E for the phase transitions considered here and the one associated with the 5-state ferromagnetic Potts model on the square lattice is conducted. The outcomes imply that if the studied phase transitions are weakly first order, then they have much larger correlation lengths than that of the 5-state Potts model.