Percolation critical exponents

In the context of the physical and mathematical theory of percolation, a percolation transition is characterized by a set of universal critical exponents, which describe the fractal properties of the percolating medium at large scales and sufficiently close to the transition. The exponents are universal in the sense that they only depend on the type of percolation model and on the space dimension. They are expected to not depend on microscopic details such as the lattice structure, or whether site or bond percolation is considered. This article deals with the critical exponents of random percolation.

Percolating systems have a parameter $p$ which controls the occupancy of sites or bonds in the system. At a critical value $p_c$, the mean cluster size goes to infinity and the percolation transition takes place. As one approaches $p_c$, various quantities either diverge or go to a constant value by a power law in $|p-p_c|$, and the exponent of that power law is the critical exponent. While the exponent of that power law is generally the same on both sides of the threshold, the coefficient or "amplitude" is generally different, leading to a universal amplitude ratio.

Description

Thermodynamic or configurational systems near a critical point or a continuous phase transition become fractal, and the behavior of many quantities in such circumstances is described by universal critical exponents. Percolation theory is a particularly simple and fundamental model in statistical mechanics which has a critical point, and a great deal of work has been done in finding its critical exponents, both theoretically (limited to two dimensions) and numerically.

Critical exponents exist for a variety of observables, but most of them are linked to each other by exponent (or scaling) relations. Only a few of them are independent, and the choice of the fundamental exponents depends on the focus of the study at hand. One choice is the set $\{\sigma ,\tau \}$ motivated by the cluster size distribution, another choice is $\{d_{\text{f}},\nu \}$ motivated by the structure of the infinite cluster. So-called correction exponents extend these sets, they refer to higher orders of the asymptotic expansion around the critical point.

Definitions of exponents

Self-similarity at the percolation threshold

Percolation clusters become self-similar precisely at the threshold density $p_c$ for sufficiently large length scales, entailing the following asymptotic power laws:

The fractal dimension $d_{\text{f}}$ relates how the mass of the incipient infinite cluster depends on the radius or another length measure, $M(L)\sim L^{d_{\text{f}}}$ at $p=p_c$ and for large probe sizes, $L\to \mathrm{\infty }$. Other notation: magnetic exponent $y_h=D-d_f$ and co-dimension ${\mathrm{\Delta }}_{\sigma }=d-d_f$.

The Fisher exponent $\tau$ characterizes the cluster-size distribution $n_s$, which is often determined in computer simulations. The latter counts the number of clusters with a given size (volume) $s$, normalized by the total volume (number of lattice sites). The distribution obeys a power law at the threshold, $n_s\sim s^{-\tau }$ asymptotically as $s\to \mathrm{\infty }$.

The probability for two sites separated by a distance $\overrightarrow{r}$ to belong to the same cluster decays as $g(\overrightarrow{r})\sim |\overrightarrow{r}{|}^{-2(d-d_{\text{f}})}$ or $g(\overrightarrow{r})\sim |\overrightarrow{r}{|}^{-d+(2-\eta )}$ for large distances, which introduces the anomalous dimension $\eta$. Also, $\delta =(d+2-\eta )/(d-2+\eta )$ and $\eta =2-\gamma /\nu$.

The exponent $\mathrm{\Omega }$ is connected with the leading correction to scaling, which appears, e.g., in the asymptotic expansion of the cluster-size distribution, $n_s\sim s^{-\tau }(1+\text{const}\times s^{-\mathrm{\Omega }})$ for $s\to \mathrm{\infty }$. Also, $\omega =\mathrm{\Omega }/(\sigma \nu )=\mathrm{\Omega }{d}_f$.

For quantities like the mean cluster size $S\sim a_0|p-p_c{|}^{-\gamma }(1+a_1(p-p_c{)}^{{\mathrm{\Delta }}_1}+\dots )$, the corrections are controlled by the exponent ${\mathrm{\Delta }}_1=\mathrm{\Omega }\beta \delta =\omega \nu$.

The minimum or chemical distance or shortest-path exponent $d_{\mathrm{m}\mathrm{i}\mathrm{n}}$ describes how the average minimum distance $⟨\ell ⟩$ relates to the Euclidean distance $r$, namely $⟨\ell ⟩\sim r^{d_{\mathrm{m}\mathrm{i}\mathrm{n}}}$ Note, it is more appropriate and practical to measure average $r$, <$r$> for a given $\ell$. The elastic backbone has the same fractal dimension as the shortest path. A related quantity is the spreading dimension $d_{\ell }$, which describes the scaling of the mass M of a critical cluster within a chemical distance $\ell$ as $M\sim {\ell }^{d_{\ell }}$, and is related to the fractal dimension $d_f$ of the cluster by $d_{\ell }=d_f/d_{\mathrm{m}\mathrm{i}\mathrm{n}}$. The chemical distance can also be thought of as a time in an epidemic growth process, and one also defines ${\nu }_t$ where $d_{\mathrm{m}\mathrm{i}\mathrm{n}}={\nu }_t/\nu =z$, and $z$ is the dynamical exponent. One also writes ${\nu }_{\parallel }={\nu }_t$.

Also related to the minimum dimension is the simultaneous growth of two nearby clusters. The probability that the two clusters coalesce exactly in time $t$ scales as $p(t)\sim t^{-\lambda }$ with $\lambda =1+5/(4d_{\mathrm{m}\mathrm{i}\mathrm{n}})$.

The dimension of the backbone, which is defined as the subset of cluster sites carrying the current when a voltage difference is applied between two sites far apart, is $d_{\text{b}}$ (or $d_{\text{BB}}$). One also defines $\xi =d-d_{\text{b}}$.

The fractal dimension of the random walk on an infinite incipient percolation cluster is given by $d_w$.

The spectral dimension $\tilde{d}$ such that the average number of distinct sites visited in an $N$-step random walk scales as $N^{\tilde{d}}$.

Critical behavior close to the percolation threshold

The approach to the percolation threshold is governed by power laws again, which hold asymptotically close to $p_c$:

The exponent $\nu$ describes the divergence of the correlation length $\xi$ as the percolation transition is approached, $\xi \sim |p-p_c{|}^{-\nu }$. The infinite cluster becomes homogeneous at length scales beyond the correlation length; further, it is a measure for the linear extent of the largest finite cluster. Other notation: Thermal exponent $y_t=1/\nu$ and dimension ${\mathrm{\Delta }}_{ϵ}=d-1/\nu$.

Off criticality, only finite clusters exist up to a largest cluster size $s_{max}$, and the cluster-size distribution is smoothly cut off by a rapidly decaying function, $n_s\sim s^{-\tau }f(s/s_{max})$. The exponent $\sigma$ characterizes the divergence of the cutoff parameter, $s_{max}\sim |p-p_c{|}^{-1/\sigma }$. From the fractal relation we have $s_{max}\sim {\xi }^{d_{\text{f}}}$, yielding $\sigma =1/\nu d_{\text{f}}$.

The density of clusters (number of clusters per site) $n_c$ is continuous at the threshold but its third derivative goes to infinity as determined by the exponent $\alpha$: $n_c\sim A+B(p-p_c)+C(p-p_c{)}^2+D_{\pm }|p-p_c{|}^{2-\alpha }+\cdots$, where $D_{\pm }$ represents the coefficient above and below the transition point.

The strength or weight of the percolating cluster, $P$ or $P_{\mathrm{\infty }}$, is the probability that a site belongs to an infinite cluster. $P$ is zero below the transition and is non-analytic. Just above the transition, $P\sim (p-p_c{)}^{\beta }$, defining the exponent $\beta$. $P$ plays the role of an order parameter.

The divergence of the mean cluster size $S=\sum _s{s}^2{n}_s/p_c\sim |p-p_c{|}^{-\gamma }$ introduces the exponent $\gamma$.

The gap exponent Δ is defined as Δ = 1/(β+γ) = 1/σ and represents the "gap" in critical exponent values from one moment $M_n$ to the next $M_{n+1}$ for $n>2$.

The conductivity exponent $t=\nu t^{\prime }$ describes how the electrical conductivity $C$ goes to zero in a conductor-insulator mixture, $C\sim (p-p_c{)}^t$. Also, $t^{\prime }=\zeta$.

Surface critical exponents

The probability a point at a surface belongs to the percolating or infinite cluster for $p\ge p_c$ is $P_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}\sim (p-p_c{)}^{{\beta }_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}}$.

The surface fractal dimension is given by $d_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}=d-1-{\beta }_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}/\nu$.

Correlations parallel and perpendicular to the surface decay as $g_{\parallel }(\overrightarrow{r})\sim |\overrightarrow{r}{|}^{2-d-{\eta }_{\parallel }}$ and $g_{\perp }(\overrightarrow{r})\sim |\overrightarrow{r}{|}^{2-d-{\eta }_{\perp }}$.

The mean size of finite clusters connected to a site in the surface is ${\chi }_1\sim |p-p_c{|}^{-{\gamma }_1}$.

The mean number of surface sites connected to a site in the surface is ${\chi }_{1,1}\sim |p-p_c{|}^{-{\gamma }_{1,1}}$.

Scaling relations

Hyperscaling relations

$\tau =\frac{d}{d_{\text{f}}}+1$

$d_{\text{f}}=d-\frac{\beta }{\nu }$

$\eta =2+d-2d_{\text{f}}$

Relations based on Template:Math

$\alpha =2-\frac{\tau -1}{\sigma }$

$\beta =\frac{\tau -2}{\sigma }$

$\gamma =\frac{3-\tau }{\sigma }$

$\beta +\gamma =\frac{1}{\sigma }=\nu d_f$

$\nu =\frac{\tau -1}{\sigma d}=\frac{2\beta +\gamma }{d}$

$\delta =\frac{1}{\tau -2}$

Relations based on Template:Math

$\alpha =2-\nu d$

$\beta =\nu (d-d_{\text{f}})$

$\gamma =\nu (2d_{\text{f}}-d)$

$\sigma =\frac{1}{\nu d_{\text{f}}}$

Conductivity scaling relations

$d_w=d+\frac{t-\beta }{\nu }$

$t^{\prime }=d_w-d_f$

$\tilde{d}=2d_f/d_w$

Surface scaling relations

${\eta }_{\parallel }=2-d+2{\beta }_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}/\nu$

$d_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}=d-1-{\beta }_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}/\nu$

${\eta }_{\parallel }=2{\eta }_{\perp }-\eta$

${\gamma }_1=\nu (2-{\eta }_{\perp })$

${\gamma }_{1,1}=\nu (d-1-2{\beta }_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}/\nu )=\nu (1-{\eta }_{\parallel })$

$\gamma +\nu =2{\gamma }_1-{\gamma }_{1,1}$

$x_1={\beta }_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}/\nu$

Exponents for standard percolation

$d$	$1$	$2$	$3$	$4$	$5$	$6-\epsilon$Template:Refn	$6+$
$\alpha$	1	–2/3	-0.625(3) -0.64(4)	-0.756(40) -0.75(2)	-0.870(1)	$-1+\frac{\epsilon }{7}-\frac{443}{2^2{3}^2{7}^3}{\epsilon }^2$	-1
$\beta$	0	0.14(3) 5/36	0.39(2) 0.4181(8) 0.41(1) 0.405(25), 0.4273 0.4053(5) 0.429(4)	0.52(3) 0.639(20) 0.657(9) 0.6590 0.658(1)	0.66(5) 0.835(5) 0.830(10) 0.8457 0.8454(2)	$1-\frac{\epsilon }{7}-\frac{61}{2^2{3}^2{7}^3}{\epsilon }^2$	1
$\gamma$	1	43/18	1.6 1.80(5) 1.66(7) 1.793(3) 1.805(20) 1.8357 1.819(3) 1.78(3)	1.6(1) 1.48(8) 1.422(16) 1.4500 1.435(15) 1.430(6)	1.3(1) 1.18(7) 1.185(5) 1.1817 1.1792(7)	$1+\frac{\epsilon }{7}+\frac{565}{2^2{3}^2{7}^3}{\epsilon }^2$	1
$\delta$	$\mathrm{\infty }$	91/5, 18	5.29(6) Template:Ref 5.3 5.16(4)	3.9 3.198(6) 3.175(8)	3.0 2.3952(12)	$2+\frac{2}{7}\epsilon +\frac{565}{2\cdot 3^2{7}^3}{\epsilon }^2$	2
$\eta$	1	5/24	-0.046(8) -0.059(9) -0.07(5) -0.0470 −0.03(1) −0.0431⁢(14)	-0.12(4) -0.0944(28) -0.0929(9) -0.0954 -0.084(4)	-0.075(20) -0.0565 −0.0547(10)	$-\frac{\epsilon }{21}-\frac{206}{3^3{7}^3}{\epsilon }^2$	0
$\nu$	1	1.33(5) 4/3	0.8(1), 0.80(5), 0.872(7) 0.875(1) 0.8765(18) 0.8960 0.8764(12) 0.8751(11) 0.8762(12) 0.8774(13) 0.88(2) 0.8762(7)	0.6782(50) 0.689(10) 0.6920 0.693 0.6852(28) 0.6845(23) 0.6845(6) 0.686(2) 0.6842(16)	0.51(5) 0.569(5) cited in 0.571(3) 0.5746 0.5723(18) 0.5737(33) 0.5757(7) 0.5739(1) 0.5720(43)	$\frac{1}{2}+\frac{5}{84}\epsilon +\frac{589}{2^2{3}^3{7}^3}{\epsilon }^2$	1/2
$\sigma$	1	36/91	0.42(6) 0.445(10) 0.4522(8) 0.4524(6) 0.4419 0.452(7)	0.476(5) 0.4742 0.4789(14)	0.496(4) 0.4933 0.49396(13)	$\frac{1}{2}-\frac{1}{98}{\epsilon }^2$	1/2
$\tau$	2	187/91	2.186(2) 2.1888 2.189(2) 2.190(2) 2.189(1) 2.18906(8) 2.18909(5) 2.1892(1) 2.1938(12)	2.26 2.313(3) 2.3127(6) 2.313(2) 2.3124 2.3142(5) 2.3150(8)	2.33 2.412(4) 2.4171 2.419(1) 2.4175(2)	$\frac{5}{2}-\frac{\epsilon }{14}+\frac{313}{2^3{3}^2{7}^3}{\epsilon }^2$	5/2
$d_{\text{f}}$	1	91/48	2.523(4) Template:Ref 2.530(4) Template:Ref 2.5230(1) 2.5226(1) 2.52293(10)	3.12(2), 3.05(5), 3.003 3.0472(14) 3.046(7) 3.046(5) 3.0479 3.0437(11) 3.0446(7)	3.54(4) 3.69(2) 3.528 3.524(2) 3.5260(14)	$4-\frac{10}{21}\epsilon +\frac{103}{3^3{7}^3}{\epsilon }^2$	4
$\mathrm{\Omega }$		0.70(2) 0.77(4) 0.77(2) 72/91 0.44(9)	0.50(9) 0.64(2) 0.73(8) 0.65(2) 0.60(8) 0.77(3) 0.64(5)	0.31(5) 0.5(1) 0.37(4) 0.4008	0.27(7) 0.2034 0.210(2)	$\frac{\epsilon }{4}-\frac{283}{2^2{3}^2{7}^2}{\epsilon }^2$
$\omega$		3/2	1.26(23) 1.6334 1.62(13) 1.61(5)	0.94(15) 1.2198 1.13(10) 1.0(2)	0.96(26) 0.7178	$\epsilon -\frac{671}{2\cdot 3^2{7}^2}{\epsilon }^2$	0
$t^{\prime }=t/\nu$		0.9479 0.995(1) 0.977(8)) 0.9825(8)	2.276(12) 2.26(4) 2.305(15) 2.283(3)				3
$d_w$		2.8784(8)
$d_s$		4/3 1.327(1) 1.3100(11)	1.32(6)
$d_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}$		2/3	1.04(5) 1.030(6) 1.0246(4)	1.32(7)	1.65(3)	$2-\frac{5}{21}\epsilon$	2
${\beta }_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}/\nu =x_1$		1/3	0.98(2) 0.970(6) 0.975(4) 0.9754(4) 0.974(2)	1.64(2)	2.408(5)	3
${\eta }_{\parallel }$ (surf)		2/3	1.02(12) 1.08(10)	1.37(13)	1.7(6)
$d_{\text{b}}$		1.60(5) 1.64(1) 1.647(4) 1.6432(8) 1.6434(2) 1.64336(10) 1.64333316328711...*	1.8, 1.77(7) 1.855(15)	1.95(5) 1.9844(11)	2.00(5) 2.0226(27)		2
$d_{\text{min}}=z$		1.132(2) 1.130(3) 1.1307(4) 1.1303(8) 1.1306(3) 1.130 77(2)	1.35(5) 1.34(1) 1.374(6) 1.3756(6) 1.3756(3) 1.3755(3)	1.607(5) 1.6042(5)	1.812(6) 1.8137(16)		2
$\lambda =\mu +1$		2.1055(10) 2.1056(3) 2.1045(10) 2.105

• For $d=2$, $d_{\mathrm{b}}=2-(z^2-1)/12$ where $z$ satisfies $\sqrt{3}z/4+\sin (2\pi z/3)=0$ near $z=2.3$.

Exponents for protected percolation

In protected percolation, bonds are removed one at a time only from the percolating cluster. Isolated clusters are no longer modified. Scaling relations: ${\beta }^{\prime }=\beta /(1+\beta )$, ${\gamma }^{\prime }=\gamma /(1+\beta )$, ${\nu }^{\prime }=\nu /(1+\beta )$, ${\tau }^{\prime }=\tau$ where the primed quantities indicated protected percolation

$d$	$1$	$2$	$3$	$4$	$5$	$6-\epsilon$	$6+$
${\beta }^{\prime }$		5/41	0.288 71(15)
${\gamma }^{\prime }$		86/41	1.3066(19)
${\tau }^{\prime }$		187/91	2.1659(21)

Exponents for directed percolation

Directed percolation (DP) refers to percolation in which the fluid can flow only in one direction along bonds—such as only in the downward direction on a square lattice rotated by 45 degrees. This system is referred to as "1 + 1 dimensional DP" where the two dimensions are thought of as space and time.

${\nu }_{\perp }$ and ${\nu }_{\parallel }$ are the transverse (perpendicular) and longitudinal (parallel) correlation length exponents, respectively. Also $\zeta =1/z={\nu }_{\perp }/{\nu }_{\parallel }$. It satisfies the hyperscaling relation $d/z=\eta +2\delta$.

Another convention has been used for the exponent $z$, which here we call $z^{\prime }$, is defined through the relation $⟨R^2⟩\sim t^{z^{\prime }}$, so that $z^{\prime }={\nu }_{\parallel }/{\nu }_{\perp }=2/z$. It satisfies the hyperscaling relation $d{z}^{\prime }=2\eta +4\delta$.

$\delta$ is the exponent corresponding to the behavior of the survival probability as a function of time: $P(t)\sim t^{-\delta }$.

$\eta$ (sometimes called $\mu$) is the exponent corresponding to the behavior of the average number of visited sites at time $t$ (averaged over all samples including ones that have stopped spreading): $N(t)\sim t^{-\eta }$.

The d(space)+1(time) dimensional exponents are given below.

$d+1$	$1+1$	$2+1$	$3+1$	$4-\epsilon$	$\text{Mean Field}$
$\beta$	0.276486(8) 0.276 7(3)	0.5834(30) 0.580(4)	0.813(9) 0.818(4) 0.82205	$1-\frac{\epsilon }{6}$	1
$\delta ,\alpha$	0.159464(6) 0.15944(2)	0.4505(1) 0.451(3) 0.4509(5) 0.4510(4) 0.460(6)	0.732(4) 0.7398(10) 0.73717	$1-\frac{\epsilon }{4}$	1
$\eta ,\theta$	0.313686(8) 0.31370(5)	0.2303(4) 0.2307(2) 0.2295(10) 0.229(3) 0.214(8)	0.1057(3) 0.114(4) 0.12084	$\frac{\epsilon }{12}$	0
${\nu }_{\parallel }$	1.733847(6) 1.733825(25) 1.7355(15) 1.73(2)	1.16(5) 1.287(2) 1.295(6)	1.106(3) 1.11(1) 1.10571	$1+\frac{\epsilon }{12}$	1
${\nu }_{\perp }$	1.096854(4) 1.096844(14) 1.0979(10)	0.7333(75) 0.729(1)	0.584(5) 0.582(2) 0.58360	$\frac{1}{2}+\frac{\epsilon }{16}$	$\frac{1}{2}$
$z=2/z^{\prime }$	1.580745(10) 1.5807(2)	1.7660(16) 1.765(3) 1.766(2) 1.7665(2) 1.7666(10)	1.88746 1.8990(4) 1.901(5)	$2-\frac{\epsilon }{12}$	2
$\gamma$	2.277730(5) = 41/18?, 2.278(2)	1.595(18)	1.237(23)	$1+\frac{\epsilon }{6}$	1
$\tau$	2.112(5), 2.1077(13), 2.10825(8)

Scaling relations for directed percolation

$\beta =\frac{\tau -2}{\sigma }$

$\gamma =\frac{3-\tau }{\sigma }$

$\tau =2+\frac{2}{1+\gamma /\beta }$

$\tilde{\tau }={\nu }_{\parallel }-\beta$

$\eta =\gamma /{\nu }_{\parallel }-1$

$d_{\mathrm{D}\mathrm{P}}=2-\beta /{\nu }_{\parallel }$

$d_{b,\mathrm{D}\mathrm{P}}=2-2\beta /{\nu }_{\parallel }$

$\mathrm{\Delta }=\beta +\gamma$

$d{z}^{\prime }=2\eta +4\delta$

$d/z=\eta +2\delta$

Exponents for dynamic percolation

For dynamic percolation (epidemic growth of ordinary percolation clusters), we have

$P(t)\sim L^{-\beta /\nu }\sim (t^{1/d_{\mathrm{m}\mathrm{i}\mathrm{n}}}{)}^{-\beta /\nu }=t^{-\delta }$, implying

$\delta =\frac{\beta }{\nu d_{\mathrm{m}\mathrm{i}\mathrm{n}}}=\frac{d-d_f}{d_{\mathrm{m}\mathrm{i}\mathrm{n}}}$

For $N(t)\sim t^{\eta }$, consider $N(\le s)\sim s^{3-\tau }\sim R^{d_f(3-\tau )}\sim t^{d_f(3-\tau )/d_{\mathrm{m}\mathrm{i}\mathrm{n}}}$, and taking the derivative with respect to $t$ yields $N(t)\sim t^{d_f(3-\tau )/d_{\mathrm{m}\mathrm{i}\mathrm{n}}-1}$, implying

$\eta =\frac{d_f(3-\tau )}{d_{\mathrm{m}\mathrm{i}\mathrm{n}}}-1=\frac{2d_f-d}{d_{\mathrm{m}\mathrm{i}\mathrm{n}}}-1$

Also, $z=d_{\mathrm{m}\mathrm{i}\mathrm{n}}$

Using exponents above, we find

$d:$	$2$	$3$	$4$	$5$	$6-\epsilon$	$\text{Mean Field}$
$\delta$	0.09212	0.34681	0.59556	0.8127		1
$\eta ,\mu$	0.584466	0.48725	0.30233	0.1314		0

See also

• Critical exponent
• Universality class
• Percolation theory
• Percolation threshold
• Percolation surface critical behavior
• Conductivity near the percolation threshold

References

Further reading

• Template:Citation