Using a histogram Monte Carlo simulation method (HMCSM), Hu, Lin, and Chen found that bond and site percolation models on planar lattices have universal finite-size scaling functions for the existence probability
Ep, the percolation probability
P, and the probability
Wn for the appearance of
n percolating clusters in these models. In this paper we extend above study to percolation on three-dimensional lattices with various linear dimensions
L. Using the HMCSM, we calculate the existence probability
Ep and the percolation probability
P for site and bond percolation on a simple-cubic (sc) lattice, and site percolation on body-centered-cubic and face-centered-cubic lattices; each lattice has the same linear dimension in three dimensions. Using the data of
Ep and
P in a percolation renormalization group method, we find that the critical exponents obtained are quite consistent with the universality of critical exponents. Using a small number of nonuniversal metric factors, we find that
Ep and
P have universal finite-size scaling functions. This implies that the critical
Ep is a universal quantity, which is $0.265\ifmmode\pm\else\textpm\fi{}0.005$ for free boundary conditions and $0.924\ifmmode\pm\else\textpm\fi{}0.005$ for periodic boundary conditions. We also find that
Wn for site and bond percolation on sc lattices have universal finite-size scaling functions.
We propose a site random-cluster model by introducing an additional cluster weight in the partition function of the traditional site percolation. To simulate the model on a square lattice, we combine the color-assignation and the Swendsen-Wang methods to design a highly efficient cluster algorithm with a small critical slowing-down phenomenon. To verify whether or not it is consistent with the bond random-cluster model, we measure several quantities, such as the wrapping probability Re, the percolating cluster density P∞, and the magnetic susceptibility per site χp, as well as two exponents, such as the thermal exponent yt and the fractal dimension yh of the percolating cluster. We find that for different exponents of cluster weight q=1.5, 2, 2.5, 3, 3.5, and 4, the numerical estimation of the exponents yt and yh are consistent with the theoretical values. The universalities of the site random-cluster model and the bond random-cluster model are completely identical. For larger values of q, we find obvious signatures of the first-order percolation transition by the histograms and the hysteresis loops of percolating cluster density and the energy per site. Our results are helpful for the understanding of the percolation of traditional statistical models.
We study the bond-percolation model in two and three dimensions by Monte Carlo simulation, and investigate the finite-size scaling behavior of several quantities that account for fluctuations of the total numbers of clusters and occupied bonds,
Nc and
Nb, respectively. These quantities include ${C}_{2c}=(⟨{N}_{c}^{2}⟩\ensuremath{-}{⟨{N}_{c}⟩}^{2})∕{L}^{d}$ and ${C}_{cb}=(⟨{N}_{c}{N}_{b}⟩\ensuremath{-}⟨{N}_{c}⟩⟨{N}_{b}⟩)∕{L}^{d}$, where
L is the linear system size and
d is the spatial dimensionality. In statistical models with thermal fluctuations,
C2c and
Ccb are specific heatlike quantities. Despite the absence of thermal fluctuations in percolation, we find that the leading finite-size scaling of
C2c and
Ccb is described by the thermal critical exponent ${y}_{t}\text{\ensuremath{-}}d$. We also measure quantity ${\ensuremath{\kappa}}_{b}=2⟨{N}_{b}{S}_{2}⟩∕⟨{S}_{2}⟩\ensuremath{-}⟨{N}_{b}{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{S}}_{4}⟩∕⟨{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{S}}_{4}⟩\ensuremath{-}⟨{N}_{b}⟩$ and an analogous quantity ${\ensuremath{\kappa}}_{c}$ for
Nc, where
S2 and ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{S}}_{4}$ are quantities associated with the second and the fourth moments of cluster sizes, respectively. At criticality, we show that ${\ensuremath{\kappa}}_{b}$ and ${\ensuremath{\kappa}}_{c}$ diverge as
Lyt for $L\ensuremath{\rightarrow}\ensuremath{\infty}$. The analysis of the data of ${\ensuremath{\kappa}}_{b}$ and ${\ensuremath{\kappa}}_{c}$ yields
yt=1.145(2) for the three-dimensional percolation, in good agreement with existing results.
In two space dimensions, the percolation point of the pure-site clusters of the Ising model coincides with the critical point
Tc of the thermal transition and the percolation exponents belong to a special universality class. By introducing a bond probability
pB<1, the corresponding site-bond clusters keep on percolating at
Tc and the exponents do not change, until ${p}_{B}{=p}_{\mathrm{CK}}=1\ensuremath{-}\mathrm{exp}(\ensuremath{-}2J/kT):$ for this special expression of the bond weight the critical percolation exponents switch to the two-dimensional Ising universality class. We show here that the result is valid for a wide class of bidimensional models with a continuous magnetization transition: there is a critical bond probability
pc such that, for any
pB> pc, the onset of percolation of the site-bond clusters coincides with the critical point of the thermal transition. The percolation exponents are the same for
pc<pB< 1 but, for
pB=pc, they suddenly change to the thermal exponents, so that the corresponding clusters are critical droplets of the phase transition. Our result is based on Monte Carlo simulations of various systems near criticality.
We obtain the critical parameters for the site-percolation problem on the square lattice to a high degree of accuracy (comparable to that of series expansions) by using a Monte Carlo position-space renormalization-group procedure directly on the site-occupation probability. Our method involves calculating recursion relations using progressively larger lattice rescalings,
b. We find smooth sequences for the value of the critical percolation concentration
pc(b) and for the scaling powers
yp(b) and
yh(b). Extrapolating these sequences to the limit $b\ensuremath{\rightarrow}\ensuremath{\infty}$ leads to quite accurate numerical predictions. Further, by considering other weight functions or "rules" which also embody the essential connectivity feature of percolation, we find that the numerical results in the infinite-cell limit are in fact "rule independent." However, the actual fashion in which this limit is approached does depend upon the rule chosen. A connection between extrapolation of our renormalization-group results and finite-size scaling is made. Furthermore, the usual finite-size scaling arguments lead to independent estimates of
pc and
yp. Combining both the large-cell approach and the finite-size scaling results, we obtain ${y}_{p}=0.7385\ifmmode\pm\else\textpm\fi{}0.0080$ and ${y}_{h}=1.898\ifmmode\pm\else\textpm\fi{}0.003$. Thus we find ${\ensuremath{\alpha}}_{p}=\ensuremath{-}0.708\ifmmode\pm\else\textpm\fi{}0.030$, ${\ensuremath{\beta}}_{p}=0.138(+0.006,\ensuremath{-}0.005)$, ${\ensuremath{\gamma}}_{p}=2.432\ifmmode\pm\else\textpm\fi{}0.035$, ${\ensuremath{\delta}}_{p}=18.6\ifmmode\pm\else\textpm\fi{}0.6$, ${\ensuremath{\nu}}_{p}=1.354\ifmmode\pm\else\textpm\fi{}0.015$, and $2\ensuremath{-}{\ensuremath{\eta}}_{p}=1.796\ifmmode\pm\else\textpm\fi{}0.006$. The site-percolation threshold is found for the square lattice at ${p}_{c}=0.5931\ifmmode\pm\else\textpm\fi{}0.0006$. We note that our calculated value of ${\ensuremath{\nu}}_{p}$ is in considerably better agreement with the proposal of Klein et al. that ${\ensuremath{\nu}}_{p}=\frac{\mathrm{ln}\sqrt{3}}{\mathrm{ln}(\frac{3}{2})}\ensuremath{\cong}1.3548$, than with den Nijs' recent conjecture, which predicts ${\ensuremath{\nu}}_{p}=\frac{4}{3}$. However, our results cannot entirely rule out the latter possibility.