The phase transitions of the frustrated $J_1$-$J_2$ Ising model on the honeycomb lattice

Abstract

We study the phase transitions of the frustrated $J_1$-$J_2$ Ising model on the honeycomb lattice using the non-perturbative first principle Monte Carlo simulations. Here $J_1<0$ and $J_2>0$ are the nearest and next-to-nearest couplings, respectively. In particular, the values of $J_2/|J_1|=0.20,0.22,0.23,0.24,0.3,0.5,0.8,$ and 1.0 are considered in our study. Based on the numerical outcomes, we find that the phase transitions for $J_2/|J_1|=0.20,0.22,0.23,$ and 0.24 are second order and are governed by the 2D Ising universality class. In addition, we find evidence to support the facts that there are transitions for $J_2/|J_1|=0.5,0.8$ and 1.0 and these phase transitions are second order. Our results also indicate phase transition is unlikely to take place for $g=0.3$. We are not able to obtain results for $J_2/|J_1|$ $\in$ (0.24, 0.3) because the associated integrated autocorrelation times or (and) the equilibrium times are extremely large at the low-temperature region. A comparison between the outcomes presented here and the available results in the literature is briefly conducted as well.

Comments: 20 pages (signle column), 67 figures