Monte Carlo simulations for two-dimensional Ising system far from equilibrium

With a dynamic Monte Carlo simulation, which is free of critical slowing down, the critical behavior of two-dimensional Ising model at non-equilibrium states is investigated. We focus on the two-time autocorrelation function A(t, t') to identify the dynamic exponent z in two different evolution stages, quenched from a high temperature state and a completely ordered state to the critical state at T-C, respectively. By using the heat-bath algorithm, the dynamic critical exponent is estimated as z approximate to 2.16. As a result, the universality of the scaling behavior is verified numerically, and the exponent lambda(C) is determined complete by lambda(C) = beta/v for the evolutions from an initial completely ordered state.