*What is a critical behavior?*

In layman's term:Examples:

- A drastic change following small variation of one or more parameters.
Phase transition of water:

- liquid <-> ice transition at freezing temperature.
Percolation:Qualitative behaviors of the percolation probability P and average cluster size S:

An interesting thing happens at the percolation threshold, pc; this point is also referred as critical point, and the phenomena associated with it are called critical phenomena.The modern theory of critical phenomena is a well-developed theory and a great success of statistical mechanics. K. Wilson won his Nobel Prize for his renormalization group theory of critical phenomena.

Critical phenomena will be studied again in molecular dynamics and Monte Carlo methods later in the course.

*Critical Exponents:*

Near the critical point pc the two quantities can be described empirically as:

0 , p < p_{c};

P ~ {

(p-p_{c}), p > p^{beta}_{c}, p -> p_{c}.

S ~ | p - p

_{c}|, p -> p^{gamma}_{c}.Where

beta and gamma are some constants called critical exponents.The proportionality sign ~ takes the usual meaning:

The ratio of the left-hand side and the right-hand side approaches a constant as p approaches p_{c}.

- Critical behavior only occurs near a critical point:

The equation or laws that the quantity P or S obeys is called asymptotic law since it takes that form only when p is close to the critical value pc.

- Critical exponents like beta and gamma are universal.

More precisely, the value beta is the same for a whole class of models and experimental samples independent of the details of the models or samples. The value depends only on the dimensionality of the problem.Two dimensional percolation is different from three dimensional percolation. It depends whether the percolation is short ranged or long ranged. But it does not depend whether the lattice is a square lattice or triangular lattice. The critical exponents are the same on both lattices.

The universality of exponents justifies the study of the simplest possible model, and the result also applies to more general situations. One of the tasks of this study is to calculate the critical exponents by computer simulation.

*Scaling assumption of the cluster number*

The percolation probability P and average cluster size S have a simple power-law behavior near the critical point pc. What about the cluster number n_{s}?This quantity is more complicated because it depends on two variables, the cluster size s and occupation probability p. So we should write more precisely as n

_{s}(p).Except in one dimension and a special lattice (Bethe lattice), there are no exact solutions for the percolation problems. Based on many examples and experience in critical phenomena, the following functional form has been assumed:

n_{s}(p) = sf[(p-p^{-tau}_{c})s], p -> p^{sigma}_{c}, s -> infinity

Wheretau and sigma are also critical exponents. The above equation is called thescaling assumption for the cluster number.The exact values for the exponents tau, sigma and the exact form for the function f(z) are not known in general. They must be determined by (computer) experiment or by other methods. In fact, the actual assumption must be tested. People believe that the assumption indeed holds. There has been no rigorous proof of the assumption. Do you want to take on the challenge?

We have introduced four different exponents beta, gamma, tau, and sigma. Since P and S can be derived from n_{s}(p), the exponents are not independent. They are related. This is one of the consequence of the scaling assumption.

A test of the exponent equality is also a test for the scaling assumption. From the equation:

P = p - sumtogether with the scaling assumption, we have:(n_{s = 1, 2, 3, ...}_{s}(p) s)

P = p - sumWe'll use a trick to write the equation in the variable p-ps_{s=1 to infinity}f((p-p^{1-tau}_{c}) s).^{sigma}_{c}. Since P=0 at p=pc, setting p=pc and subtracting it from the above equation, we have:P = (p-p

_{c}) - sums_{s=1 to infinity}[f((p-p^{1-tau}_{c}) s) - f(0)].^{sigma}

We'll convert the summation to integral. The error introduced is small near the percolation threshold.

P = (p-p

_{c}) - ints_{s=s0 to infinity}[f((p-p^{1-tau}_{c}) s) - f(0)]ds^{sigma}Only the upper limit matters. The lower limit does not influence the finite result very much. So we set it at some arbitrary value s0. By changing the integration variable

s=(x/(p-pone finds:_{c}))^{1/sigma}

P = (p-pwhere C is a constant._{c}) - ( p -p_{c})C ~ (p-p^{(tau - 2)/sigma}_{c})^{beta}Comparing the exponents, using the fact that beta < 1, (neglecting (p-p

_{c}) term), we have thefirst relation:

beta = tau - 2/sigmaNow we calculate the average cluster size. This quantity diverges at the percolation threshold. The divergence comes from the numerator: sum (s2 ns ), since the denominator sum (s n_{s}) simply approaches a constant p_{c}.To leading order of the divergence term, we can approximate it by p

_{c}and thus:

S = sum s^{2}n_{s}(p) / sum s n_{s}(p) ~ sum s^{2}n_{s}(p) / p_{c}~ int s

^{2}n_{s}(p) ds ~ int sf((p-p^{2-tau}_{c}) s) ds^{sigma}~ |p-p

_{c}|~ |p-p^{(tau - 3)/sigma}_{c}|^{-gamma}

Again we have changed integration variable, and the value of the integral is just some number independent of p. Thus we have thesecond relation:

Note that in order that both beta and gamma are positive we need 2 < tau < 3.gamma = 3 - tau/sigmaHere is a list of percolation exponents for d=2, 3, 4, 5, 6-e (e is a small value close to zero), and for the Bethe lattice:

Exponent |
d=2 |
d=3 |
d=4 |
d=5 |
d=6-e |
Bethe |

beta | 5/36 | 0.41 | 0.64 | 0.84 | 1-e/7 | 1 |

gamma | 43/18 | 1.80 | 1.44 | 1.18 | 1+e/7 | 1 |

nu | 4/3 | 0.88 | 0.68 | 0.57 | 1/2+5e/84 | 1/2 |

sigma | 36/91 | 0.45 | 0.48 | 0.49 | 1/2+O(e^{2}) |
1/2 |

tau | 187/91 | 2.18 | 2.31 | 2.41 | 5/2-3e/14 | 5/2 |

D | 91/48 | 2.53 | 3.06 | 3.54 | 4-10e/21 | 4 |

The dimension d=6 is called upper critical dimension, above which the exponents become independent of the dimensionality and take the values of the Bethe lattice results.The exponent nu is called correlation length exponent, defined below. And D is the fractal dimension of clusters at the percolation threshold.

*Critical Exponents Related to Cluster Structure*

Cluster perimeter

We'll consider site percolation. Each individual cluster is called a lattice animal. The empty sites immediately adjacent to the cluster (empty neighbors) are called perimeter sites.

Let's use:

then:

tas the number of perimeter sites of a given lattice animalsas the total number of occupied sites of the lattice animal.- g
_{st}to denote the number of lattice animals (cluster configurations) with size s and perimeter tn_{s}(p) = sum_{t}[ g_{st}p^{s}( 1-p)^{t }]

is the average number of s clusters per lattice site.If we have an analytic formula for gst, then we should have solved the percolation problem exactly. Unfortunately, this is not the case except in one dimension. However, we can find ns for the first few small values of s.

For s=1 and dimension d=2:

n1 = p ( 1 - p)For s=2, d=2:^{4}

g= 1, all other g_{1,4}_{st }= 0.n2 = 2 pFor s=3, d=2:^{2}(1-p)^{6}

g= 2, all other g_{2,6}_{st }= 0.n3 = 2 p^{3}( 1 - p)^{8}+ 4 p^{3}(1-p)^{7}

g_{3,8 }= 2, g= 4, all other g_{3,7}_{st }= 0.The meaning of the formula correlates closely with the figures shown above . With the help from computer, people were able to obtain gst for s upto 30.The average perimeter t of a cluster is not proportional to the surface area of the cluster. This is because there are lot of holes in a percolation cluster. If we have one hole for, say, every thirty sites we have a perimeter proportional to the number of sites. Thus:

t ~ s, s -> infinity

This result can be proved more rigorously.

Correlation length

We have introduced the correlation function g(r), which is the probability that a site at distance r from a given occupied site is also occupied and belongs to the same cluster. The function behaves at large distance in the formCluster radius and fractal dimensiong(r) ~ e

, p < p^{-r/xi}_{c}.For p>p

_{c}it approaches a constant P^{2}. The quantityxi is called correlation length.The correlation length can also be defined as:

xiThe correlation length measures the typical size of clusters. Near the percolation threshold, clusters become very large, and the correlation length diverges,:^{2}= sum_{r}r^{2}g(r) / sum_{r}g(r)

xi ~ | p - pwhere_{c}|^{-nu}nu is called correlation length exponent.

The correlation length is a global quantity for the whole system. Cluster radius measures the dimension of clusters having exactly s sites. Since the clusters are rather irregular, we can only give an average quantity which measures the size of the clusters.The clusters at the percolation threshold are fractals, in the sense that the total size s or mass (imagine each occupied site has certain weight) is related to its size as

s ~ R, p = p_{s}^{D}_{c}, s >>1D is called the fractal dimension of the clusters. This fractal behavior for large but finite clusters also holds for the largest cluster in a finite system. If L is the linear size of the system with Ld lattice sites, the mass M (or total number of sites) of the largest cluster varies with the system size as:

M ~ LAccording to scaling theory, D can be related to other exponents as:^{D}

D = d - beta/nu, d < 6