M. E. J. Newman, R. M. Ziff
Published 2000-05-16, updated 2000-10-26Version 2
We present a new Monte Carlo algorithm for studying site or bond percolation on any lattice. The algorithm allows us to calculate quantities such as the cluster size distribution or spanning probability over the entire range of site or bond occupation probabilities from zero to one in a single run which takes an amount of time scaling linearly with the number of sites on the lattice. We use our algorithm to determine that the percolation transition occurs at occupation probability 0.59274621(13) for site percolation on the square lattice and to provide clear numerical confirmation of the conjectured 4/3-power stretched-exponential tails in the spanning probability functions.
Comments: 8 pages, including 3 postscript figures, minor corrections in this version, plus updated figures for the position of the percolation transition
Journal: Phys. Rev. Lett. 85, 4104-4107 (2000)
Efficient Monte Carlo algorithm in quasi-one-dimensional Ising spin systems
The boundary between long-range and short-range critical behavior
Models and algorithms for the next generation of glass transition studies
