Date of Original Version
We introduce a novel variance-reducing Monte Carlo algorithm for accurate determination of correlation times. We apply this method to two-dimensional Ising systems with sizes up to 15×15, using single-spin flip dynamics, random site selection, and transition probabilities according to the heat-bath method. From a finite-size scaling analysis of these correlation times, the dynamic critical exponent z is determined as z=2.1665(12).
M. P. Nightingale and H. W. J. Blöte. (1996). "Dynamic Exponent of the Two-Dimensional Ising Model and Monte Carlo Computation of the Subdominant Eigenvalue of the Stochastic Matrix." Physical Review Letters, 76(24), 4548. Available at: http://dx.doi.org/10.1103/PhysRevLett.76.4548