Superstatistical two-temperature Ising model

The previous approach of the nonequilibrium Ising model was based on the local temperature in which each site or part of the system has its own specific temperature. We introduce an approach of the two-temperature Ising model as a prototype of the superstatistic critical phenomena. The model is described by two temperatures (T_{1},T_{2}) in a zero magnetic field. To predict the phase diagram and numerically estimate the exponents, we develop the Metropolis and Swendsen-Wang Monte Carlo method. We observe that there is a nontrivial critical line, separating ordered and disordered phases. We propose an analytic equation for the critical line in the phase diagram. Our numerical estimation of the critical exponents illustrates that all points on the critical line belong to the ordinary Ising universality class.