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).
Nightingale, M. P., & Blöte, H. W.J. (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-4551. doi: 10.1103/PhysRevLett.76.4548
Available at: http://dx.doi.org/10.1103/PhysRevLett.76.4548