One way to get around the critical slowing down is to use Cluster algorithm ranther than Metropolis's algorithm[1]. Metropolis's algorithm tries to flip a single spin each time. SW[2] and Wolff's algorithm[3] are so called Cluster algorithm. Cluster algorithm can dramatically reduce the effects of critical slowing down and have been applied for many models. Wolff's method is probably the best sequential cluster algorithm. We are going to use Metropolis's algorithm and Wolff's Cluster algorithm in our project.
[1] N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth,
A.H.Teller, E.Teller, J. Chem. Phys. 21, 1087 (1953).
[2] R. H. Swendsen and J. S. Wang, Phys. Rev. Lett. 58, 86 (1987).
[3] U. Wolff, Phys. Rev. Lett. 62,361 (1989).