“…Based on the Heisenberg model and Monte Carlo method, the code had been developed by Zhang et al and has been successfully used in estimating the critical temperature of ferromagnetic monolayers. 34,43–45 For the ferromagnetic MnN 4 C 6 , FeN 4 C 6 , and CoN 4 C 6 monolayers, the couplings J m ( m = 1, 2, 3, and 4) are derived from the energy differences among FM, AFM-I, AFM-II, AFM-III, and Fi states and the values are listed in Table 4. Meanwhile, we estimate their Curie temperatures, and the magnetization and susceptibility variations with temperature are shown in Fig.…”