Alternative method to characterize continuous and discontinuous phase transitions in surface reaction models

In this paper we revisited the Ziff-Gulari-Barshad model to study its phase transitions and critical exponents through time-dependent Monte Carlo simulations. We use a method proposed recently to locate the nonequilibrium second-order phase transitions and that has been successfully used in systems with defined Hamiltonians and with absorbing states. This method, which is based on optimization of the coefficient of determination of the order parameter, was able to characterize the continuous phase transition of the model, as well as its upper spinodal point, a pseudocritical point located near the discontinuous phase transition. The static critical exponents β, ν_{∥}, and ν_{⊥}, as well as the dynamic critical exponents θ and z for the continuous transition point, were also estimated and are in excellent agreement with results found in literature.