Percolation analysis: Difference between revisions
No edit summary |
|||
(14 intermediate revisions by 3 users not shown) | |||
Line 1: | Line 1: | ||
This entry focuses on the application of '''percolation analysis''' to problems in [[statistical mechanics]]. For a general discussion see Refs. <ref name="Stauffer"> Dietrich Stauffer and Ammon Aharony "Introduction to Percolation Theory", CRC Press (1994) ISBN 9780748402533</ref> <ref name="Torquato"> Salvatore Torquato "Random Heterogeneous Materials, Microscopic and Macroscopic Properties", Springer, New York (2002) ISBN 9780387951676</ref> | This entry focuses on the application of '''percolation analysis''' to problems in [[statistical mechanics]]. For a general discussion see Refs. <ref name="Stauffer"> Dietrich Stauffer and Ammon Aharony "Introduction to Percolation Theory", CRC Press (1994) ISBN 9780748402533</ref> <ref name="Torquato">Salvatore Torquato "Random Heterogeneous Materials, Microscopic and Macroscopic Properties", Springer, New York (2002) ISBN 9780387951676</ref> | ||
==Sites, bonds, and clusters == | ==Sites, bonds, and clusters == | ||
Line 9: | Line 9: | ||
The sites of the system can belong to different types (species in the chemistry language). | The sites of the system can belong to different types (species in the chemistry language). | ||
Bonds are usually permitted only between near sites. | Bonds are usually permitted only between near sites. | ||
== Lattice and continuum (off-lattice) models == | |||
Attending to the spatial distribution of the sites, one can classify the models into lattice models and continuum (or off-lattice) models. | |||
Off-lattice models are more difficult to deal with from the numerical point of view, but in many applications they are expected to be more ''realistic'' than lattice models to capture the physics of a number of real systems <ref name=lee > [http://dx.doi.org/10.1063/1.455411 Sang Bub Lee and S. Torquato, "Pair connectedness and mean cluster size for continuum-percolation models: Computer-simulation results", Journal of Chemical Physics '''89''', 6427 (1988)] </ref> <ref name=becker > [http://dx.doi.org/10.1103/PhysRevE.80.041101 Adam M. Becker and Robert M. Ziff, "Percolation thresholds on two-dimensional Voronoi networks and Delaunay triangulations", Physical Review E '''80''', 041101 (2009)]</ref> | |||
== Connectivity rules == | == Connectivity rules == | ||
The connectivity rules (bonding criteria) that permit the creation of bonds can be of different nature: distance between sites, energetic | The connectivity rules (bonding criteria) that permit the creation of bonds can be of different nature: distance between sites, energetic | ||
interaction, types of sites, ''etc''. | interaction, types of sites, ''etc''. | ||
In addition the bonding criteria can be deterministic or probabilistic. | In addition the bonding criteria can be either deterministic or probabilistic. | ||
In | In [[statistical mechanics]] applications one can find different bonding criteria, for example: | ||
* Geometric distance: Two sites, <math>i</math>, <math>j</math>, are bonded if the distance between then satisfies: <math> r_{ij} < R_p </math>. This criterium is used in the so-called site percolation models (see below). | * Geometric distance: Two sites, <math>i</math>, <math>j</math>, are bonded if the distance between then satisfies: <math> r_{ij} < R_p </math>. This criterium is used in the so-called site percolation models (see below). | ||
* Energetic criteria: As an example, in the simulation of Ising | * Probabilistic criteria: Two sites located at a certain distance <math> r </math> are bonded with a given probability <math> b(r) </math>, with | ||
: <math> 0 \le b(r) < 1 </math>. This is the case of the so-called bond-percolation models. | |||
* Energetic criteria (probabilistic): As an example, in the simulation of [[Ising model]]s using [[Cluster algorithms|cluster algorithms]]; two sites <math>i</math>, <math>j</math>, have a bonding probability given by <math> b(r_{ij}) = 1- \min \left\{ 1, \exp \left[ 2 \beta u_{ij}(r_{ij}) \right] \right\} </math>; where <math> u_{ij} </math> is the interaction energy between sites, and <math> \beta = 1/ k_BT</math> (See [[Cluster algorithms|cluster algorithms]] for details). | |||
==Percolation threshold== | ==Percolation threshold== | ||
Line 38: | Line 45: | ||
== Percolation and finite-size scaling analysis == | == Percolation and finite-size scaling analysis == | ||
=== Example: Site-percolation on a square lattice === | === Example: Site-percolation on a square lattice === | ||
Let us consider a standard example of percolation theory, <ref name=deng > [http://dx.doi.org/10.1103/PhysRevE.72.016126 Youjin Deng and Henk W. J. Blöte, "Monte Carlo study of the site-percolation model in two and three dimensions", Physical Review E '''72''' 016126 (2005)]</ref> | Let us consider a standard example of percolation theory, <ref name=deng> [http://dx.doi.org/10.1103/PhysRevE.72.016126 Youjin Deng and Henk W. J. Blöte, "Monte Carlo study of the site-percolation model in two and three dimensions", Physical Review E '''72''' 016126 (2005)]</ref> | ||
a two-dimensional [[building up a square lattice|square lattice]] in which: | a two-dimensional [[building up a square lattice|square lattice]] in which: | ||
* Each site of the lattice can be occupied (by one ''particle'') or empty. | * Each site of the lattice can be occupied (by one ''particle'') or empty. | ||
Line 69: | Line 76: | ||
== Computation of the percolation threshold == | == Computation of the percolation threshold == | ||
A couple of simple procedures to estimate the percolation threshold (<math> x_c </math> in the example introduced above) are described here. | A couple of simple procedures to estimate the percolation threshold (<math> x_c </math> in the example introduced above) are described here. | ||
These procedures are similar to those used in the analysis of critical thermodynamic transitions<ref>[http://dx.doi.org/10.2277/0521842387 David P. Landau and Kurt Binder "A Guide to Monte Carlo Simulations in Statistical Physics", Cambridge University Press (2005)] </ref>. More sophisticated methods can be found in the literature (See Refs. <ref name= | These procedures are similar to those used in the analysis of critical thermodynamic transitions<ref>[http://dx.doi.org/10.2277/0521842387 David P. Landau and Kurt Binder "A Guide to Monte Carlo Simulations in Statistical Physics", Cambridge University Press (2005)] </ref>. More sophisticated methods can be found in the literature (See Refs. <ref name=deng/> <ref name=lin> [http://dx.doi.org/10.1103/PhysRevE.58.1521 Chai-Yu Lin and Chin-Kun Hu, "Universal finite-size scaling functions for percolation on three-dimensional lattices", Physical Review E '''58''', 1521 - 1527 (1998)] </ref> | ||
<ref name= | <ref name=Newman> [http://dx.doi.org/10.1103/PhysRevE.64.016706 M. E. J. Newman and R. M. Ziff, "Fast Monte Carlo algorithm for site or bond percolation", Physical Review E '''64''', 016706 (2001)] </ref> for details). | ||
=== Crossing of the <math> X_{\rm per}(x,L) </math> for different system sizes === | === Crossing of the <math> X_{\rm per}(x,L) </math> for different system sizes === | ||
Line 92: | Line 99: | ||
:<math> x_c \left( L \right) = x_c \left( \infty \right) + a L^{- b} </math> | :<math> x_c \left( L \right) = x_c \left( \infty \right) + a L^{- b} </math> | ||
where <math> b </math> is a [[Critical exponents |critical exponent]] (See Refs. <ref name="Stauffer" | where <math> b </math> is a [[Critical exponents |critical exponent]] (See Refs. <ref name="Stauffer"/> <ref name=Torquato/>for details). Therefore, by fitting the results of <math> x_c(L) </math> it is | ||
possible to estimate the percolation transition location: <math> x_c = x_c ( \infty ) </math>. | possible to estimate the percolation transition location: <math> x_c = x_c ( \infty ) </math>. | ||
Line 98: | Line 105: | ||
In some systems, with an appropriate definition of bonding criteria, the percolation transition occurs at the same value of the control parameter (density, temperature, [[chemical potential]]) as the thermodynamic transition <ref name=fortunato > [http://dx.doi.org/10.1103/PhysRevB.67.014102 Santo Fortunato, "Critical droplets and phase transitions in two dimensions", Physical Review B ''' 67''' 014102 (2003)] </ref> | In some systems, with an appropriate definition of bonding criteria, the percolation transition occurs at the same value of the control parameter (density, temperature, [[chemical potential]]) as the thermodynamic transition <ref name=fortunato > [http://dx.doi.org/10.1103/PhysRevB.67.014102 Santo Fortunato, "Critical droplets and phase transitions in two dimensions", Physical Review B ''' 67''' 014102 (2003)] </ref> | ||
<ref name=fortunato_2> [http://dx.doi.org/10.1088/0305-4470/36/15/304 Santo Fortunato, "Cluster percolation and critical behaviour in spin models and SU(N) gauge theories",Journal of Physics A: Mathematical and Theoretical '''36''' pp. 4269-4281 (2002)] </ref> . In these case [[cluster algorithms|cluster algorithms]] become very efficient, and moreover, the percolation analysis can be useful to develop algorithms to locate the transition (see the [[cluster algorithms|cluster algorithms]] page for more details). | <ref name=fortunato_2> [http://dx.doi.org/10.1088/0305-4470/36/15/304 Santo Fortunato, "Cluster percolation and critical behaviour in spin models and SU(N) gauge theories", Journal of Physics A: Mathematical and Theoretical '''36''' pp. 4269-4281 (2002)] </ref> | ||
<ref name=hu >[http://dx.doi.org//10.1103/PhysRevB.40.5007 Chin-Kun Hu and Kit-Sing Ma, "Monte Carlo study of the Potts model on the square and the simple cubic lattices", Physical Review B '''40''', 5007-5014 (1989)] </ref> | |||
. In these case [[cluster algorithms|cluster algorithms]] become very efficient, and moreover, the percolation analysis can be useful to develop algorithms to locate the transition (see the [[cluster algorithms|cluster algorithms]] page for more details). | |||
==References== | ==References== | ||
<references/> | <references/> | ||
;Related reading | |||
*Dietrich Stauffer and Ammon Aharony "Introduction To Percolation Theory", 2nd Edition, CRC Press (1994) ISBN 9780748402533 | |||
[[Category: Confined systems]] | [[Category: Confined systems]] |
Latest revision as of 19:12, 15 October 2023
This entry focuses on the application of percolation analysis to problems in statistical mechanics. For a general discussion see Refs. [1] [2]
Sites, bonds, and clusters[edit]
This topic concerns the analysis of connectivity of elements (sites) distributed in different positions of a given large system. Using some connectivity rules it is possible to define bonds between pairs of sites. These bonds can be used to build up clusters of sites. Two sites in a cluster can be connected directly by a bond between them or indirectly by one or more sequences of bonds between pairs of sites. The sites of the system can belong to different types (species in the chemistry language). Bonds are usually permitted only between near sites.
Lattice and continuum (off-lattice) models[edit]
Attending to the spatial distribution of the sites, one can classify the models into lattice models and continuum (or off-lattice) models. Off-lattice models are more difficult to deal with from the numerical point of view, but in many applications they are expected to be more realistic than lattice models to capture the physics of a number of real systems [3] [4]
Connectivity rules[edit]
The connectivity rules (bonding criteria) that permit the creation of bonds can be of different nature: distance between sites, energetic interaction, types of sites, etc. In addition the bonding criteria can be either deterministic or probabilistic. In statistical mechanics applications one can find different bonding criteria, for example:
- Geometric distance: Two sites, , , are bonded if the distance between then satisfies: . This criterium is used in the so-called site percolation models (see below).
- Probabilistic criteria: Two sites located at a certain distance are bonded with a given probability , with
- . This is the case of the so-called bond-percolation models.
- Energetic criteria (probabilistic): As an example, in the simulation of Ising models using cluster algorithms; two sites , , have a bonding probability given by ; where is the interaction energy between sites, and (See cluster algorithms for details).
Percolation threshold[edit]
The sizes of the clusters of a given system depend on different control parameters: density and distribution of sites, bonding criteria (which could include the effect of temperature and energy interactions), etc.
Let us consider some initial conditions of the control parameters, in which, the bonding criteria leads to the formation of small clusters, i.e. all the cluster contain a small number of particles and the cluster size is much smaller than the linear dimension of the system. Now, if one varies gradually some control parameter(s) to increase the number of bonds in the system, then the number of clusters is expected to decrease, the number of sites per cluster and the cluster size will increase; and, eventually, the largest cluster size(s) (in one or several directions) will be similar to the overall system size (the system reaches the percolation threshold of the percolation transition).
Percolation and boundary conditions[edit]
There are different possible criteria to consider that a cluster has percolated. The choice of percolation criteria usually depends on computational convenience, type of boundary conditions, dimensionality of the space, etc. In the particular case of considering periodic boundary conditions, a cluster realization is usually considered as percolating when, at least, one of the clusters becomes of infinite size (length) in, at least, one direction. This infinite size occurs, obviously, via the replication of the system that appears due to the periodic boundary conditions.
Percolation and finite-size scaling analysis[edit]
Example: Site-percolation on a square lattice[edit]
Let us consider a standard example of percolation theory, [5] a two-dimensional square lattice in which:
- Each site of the lattice can be occupied (by one particle) or empty.
- The probability of occupancy of each site is , with .
- Two sites are considered to be bonded if and only if:
- They are nearest neighbours and
- Both sites are occupied.
Fraction of percolating realizations[edit]
On such a system, it is possible to perform simulations considering different system sizes (with sites), using periodic boundary conditions. In such simulations one can generate different system realizations for given values of , and compute the fraction, , of realizations with percolating clusters. For low values of one will have , whereas when , then . Considering the behavior of as a function of , for different values of the transition between and occurs more abruptly as increases. In addition, it is possible to compute the value of the occupancy probability at which the transition would take place for an infinite system (that is to say, in the thermodynamic limit).
Finite-size scaling[edit]
Considering the functions the percolation theory predicts for large system sizes:
In addition, at , it is expected that the fraction of percolating realizations do not
depend on the system size:
- ; for large values of .
Computation of the percolation threshold[edit]
A couple of simple procedures to estimate the percolation threshold ( in the example introduced above) are described here. These procedures are similar to those used in the analysis of critical thermodynamic transitions[6]. More sophisticated methods can be found in the literature (See Refs. [5] [7] [8] for details).
Crossing of the for different system sizes[edit]
In practice, one has to compute the fraction of percolating realizations for different values of the control parameter and different system sizes . The critical value is then estimated by plotting as a function of for several values of . The crossing of the curves with different values of provide estimates of both and .
Computation of pseudo-critical parameters and extrapolation[edit]
Given the results of for a given system size , a pseudo-critical size dependent variable is computed by matching .
If the universal value value is unknown for the type of transition considered, an alternative definition for can be taken, for instance:
- .
The percolation theory predicts that the pseudo-critical values will scale as:
where is a critical exponent (See Refs. [1] [2]for details). Therefore, by fitting the results of it is possible to estimate the percolation transition location: .
Percolation threshold and critical thermodynamic transitions[edit]
In some systems, with an appropriate definition of bonding criteria, the percolation transition occurs at the same value of the control parameter (density, temperature, chemical potential) as the thermodynamic transition [9] [10] [11] . In these case cluster algorithms become very efficient, and moreover, the percolation analysis can be useful to develop algorithms to locate the transition (see the cluster algorithms page for more details).
References[edit]
- ↑ 1.0 1.1 Dietrich Stauffer and Ammon Aharony "Introduction to Percolation Theory", CRC Press (1994) ISBN 9780748402533
- ↑ 2.0 2.1 Salvatore Torquato "Random Heterogeneous Materials, Microscopic and Macroscopic Properties", Springer, New York (2002) ISBN 9780387951676
- ↑ Sang Bub Lee and S. Torquato, "Pair connectedness and mean cluster size for continuum-percolation models: Computer-simulation results", Journal of Chemical Physics 89, 6427 (1988)
- ↑ Adam M. Becker and Robert M. Ziff, "Percolation thresholds on two-dimensional Voronoi networks and Delaunay triangulations", Physical Review E 80, 041101 (2009)
- ↑ 5.0 5.1 Youjin Deng and Henk W. J. Blöte, "Monte Carlo study of the site-percolation model in two and three dimensions", Physical Review E 72 016126 (2005)
- ↑ David P. Landau and Kurt Binder "A Guide to Monte Carlo Simulations in Statistical Physics", Cambridge University Press (2005)
- ↑ Chai-Yu Lin and Chin-Kun Hu, "Universal finite-size scaling functions for percolation on three-dimensional lattices", Physical Review E 58, 1521 - 1527 (1998)
- ↑ M. E. J. Newman and R. M. Ziff, "Fast Monte Carlo algorithm for site or bond percolation", Physical Review E 64, 016706 (2001)
- ↑ Santo Fortunato, "Critical droplets and phase transitions in two dimensions", Physical Review B 67 014102 (2003)
- ↑ Santo Fortunato, "Cluster percolation and critical behaviour in spin models and SU(N) gauge theories", Journal of Physics A: Mathematical and Theoretical 36 pp. 4269-4281 (2002)
- ↑ Chin-Kun Hu and Kit-Sing Ma, "Monte Carlo study of the Potts model on the square and the simple cubic lattices", Physical Review B 40, 5007-5014 (1989)
- Related reading
- Dietrich Stauffer and Ammon Aharony "Introduction To Percolation Theory", 2nd Edition, CRC Press (1994) ISBN 9780748402533