{'fraction_non_alphanumeric': 0.024950007691124442, 'fraction_numerical': 0.0072604214736194436, 'mean_word_length': 4.455857670359181, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Random percolation as a gauge theory

Three-dimensional bond or site percolation theory on a lattice can be interpreted as a gauge theory in which the Wilson loops are viewed as counters of topological linking with random clusters. Beyond the percolation threshold large Wilson loops decay with an area law and show the universal shape effects due to flux tube quantum fluctuations like in ordinary confining gauge theories. Wilson loop correlators define a non-trivial spectrum of physical states of increasing mass and spin, like the glueballs of ordinary gauge theory. The crumbling of the percolating cluster when the length of one periodic direction decreases below a critical threshold accounts for the finite temperature deconfinement, which belongs to 2-D percolation universality class.


Introduction
Percolation is a purely geometrical phenomenon which in many respects resembles a continuous thermal phase transition. The theoretical description of the percolation processes is conventionally given in terms of the cluster sizes [1], and most of the universal scalings deal with size distribution of clusters. The point of view which is taken in this paper is different. We focus on topological entanglement of random clusters and use it to describe how percolation theory in three dimensions can be viewed as a full-fledged gauge theory.
Though the gauge group of the theory in question is trivial (it is the q → 1 limit of the symmetric group S q ), the occurrence of a confining phase yields some new hints on the mechanisms of quark confinement of more general theories.
Transcribing percolation in terms of gauge theory has also some important consequences for percolation itself. In three-dimensional systems there are almost no exact results whatsoever, but in gauge theories we have a number of well-verified conjectures that can be translated into the language of percolation. In this way we shall, for instance, relate certain linking properties of the closed paths within a percolating cluster to the universal quantum fluctuations of the chromoelectric flux tube joining a quark pair in the confining phase of whatever gauge theory.
An unsuspected property of random percolation which emerges from this new viewpoint is that its spectrum should be composed by a (possibly infinite) tower of physical states of increasing mass and spin -the glue-balls of the corresponding gauge theory. Their mass ratios near the percolation point define a totally new set of critical amplitude ratios belonging to the universality class of 3-D random percolation.
Another piece of useful information comes from considering percolation in a slab which is infinite in two dimensions, but of finite length ℓ and periodic in the remaining direction. The associated gauge model describes a system at finite temperature T = 1/ℓ. This transition is accurately described by the universality class of two-dimensional random percolation, but the corresponding deconfining temperature T c may be used to define a new critical amplitude of the three-dimensional system.
We test the universality of this new set of critical amplitudes by performing large scale numerical experiments in three different kinds of lattices: the simple cubic (SC) with bond or site percolation and the body centered cubic lattice (BCC) with bond percolation. The numerical implementation of these systems is straightforward in comparison with simulations of ordinary gauge models: no Markovian process is needed and there are neither thermalization problems nor critical slowing down. Preliminary results have been presented in [2].
The organisation of the paper is as follows. In the next Section we define a new class of observables of the percolation theory to be identified with the Wilson loops of the corresponding gauge theory. In Section 3 we describe a method to study the linking properties of the random clusters which is at the heart of our analysis. Section 4 is devoted to the comparison of confinement mechanisms and in Section 5 we extract the string tension and show the relevance of the universal terms generated by quantum fluctuations of an underlying effective string. In Section 6 we study the plaquette correlators in order to find the low-lying states of the spectrum. Section 7 is devoted to the transcription of the percolation on a slab into a gauge system at finite temperature and Section 8 discusses the universality class of the finite temperature deconfinement transition. Finally in Section 9 we draw some concluding remarks.

Observables
The most basic observables of any gauge theory are the Wilson loops. These are operators which assign to each pair (C, γ) formed by an arbitrary gauge configuration C and any closed path γ of the space a suitably defined complex number W γ (C). Their importance stems from the fact that they serve as order parameters for confinement: in the confining phase the vacuum expectation value of large Wilson loops exhibits area law.
One can define similar observables in the framework of random percolation. In this context the configurations are generated simply by occupying each site or bond on an initially empty 3-D lattice Λ with independent probability p. Two sites are considered neighbours if they share one bond. The resulting configuration is a graph G drawn on the lattice, composed by the occupied bonds (bond percolation problem) or by the bonds joining occupied neighbour sites (site percolation problem). The connected components of G are the clusters of the configuration. We choose as γ's the closed paths of the dual latticeΛ. The value W γ (G) measures the topological entanglement between γ and G. More precisely we apply the following rule This definition did not come out of the blue. Starting from the Fortuin-Kasteleyn random cluster representation of the 3-D Ising model [3] combined with the Kramers-Wannier duality [4], it is easy to express the Wilson loop W γ belonging to a Z 2 gauge system in terms of the winding numbers modulo 2 of the Fortuin-Kasteleyn clusters across the γ loop [5]. Generalising this result to q-state Potts model one is led to the above rule in the case of non-integer q. A recently developed algorithm could even allow to evaluate explicitly these quantities for real q > 0 [6]. We extend the above recipe to any percolating system, owing to the fact that bond percolation can be viewed as the q → 1 limit of the q-state Potts model.
The removal of an occupied bond b from a graph G can lead to two different issues. If the number c(G) of clusters is kept invariant, then b is necessarily a step of a closed path (or loop) of G, whereas if c(G) increases by one, b cannot lie in a loop and is called a bridge (see Fig. 1). Clearly only the former bonds contribute to W γ (G). If two graphs G and G ′ have the same loops and differ only in the bridges, then they yield the same value of W γ for any γ. In more technical terms we may write the double implication hence the transformation G → G ′ has some resemblance to a gauge transformation. Note that the connected correlator among occupied bonds is exactly zero by construction, but this is not an invariant quantity under the above G → G ′ transformation. Cutting all bridges of G yields the maximal invariant subset B G of G, made by bonds belonging to some loop. Of course we have W γ (G) = W γ (B G ) for any γ. The connected correlator among bonds belonging to B G is by no means trivial and is directly related to the connected correlator of the

Cutting all bridges
In our approach the only bonds which play a role in the evaluation of the observables defining the gauge theory are those belonging to loops. Thus, once a new configuration G is generated, we first get rid of all bridges. One way of achieving this goal is the following.
1. Eliminate all the dangling ends (see Fig. 2a and b). At this stage the remaining graph is formed by loops and lines of bridges connecting them.
2. Build a reduced graph in which the only vertices are the lattice sites with more than two incident occupied bonds. The edges of the graph are formed by lines of occupied bonds (Fig. 2b).
3. Erase one edge at a time and apply each time a cluster reconstruction algorithm (for instance Hoshen-Kopelman [7]) in order to pick off the remaining bridges (Fig. 2c).
In order to check whether a planar loop γ ⊂Λ is linked to a configuration G, we first project out all bridges as discussed. Then we switch off the layer of occupied bonds which pierce the planar surface Σ encircled by γ and rebuild the cluster structure of the cut graph. For a non-trivial linking there must be a c b at least one cluster which reaches the layer on either side (see the graph on the left of Fig. 3). In such a case we build an auxiliary graph in which the vertices represent the clusters of sites on either side of the layer; a cluster connecting sites of opposite sides is represented by two vertices, one for each side (see the graph on the right of Fig. 3). We draw an edge between two clusters in the opposite sides of the layer if they are connected by switched off bonds. The configuration is truly linked to γ if there is a path (at least) joining two vertices lying on opposite margins of the layer, but belonging to the same cluster of the cut graph (like the vertices a and a ′ of the figure). The whole procedure of cutting all bridges of a configuration and then evaluating its linking property with a set of Wilson loops is time demanding, thus a good implementation is mandatory. However the variance of the measured quantities turns out to be small, hence one can reach more precise results than the corresponding estimates in ordinary gauge theories.

Confinement mechanisms
According to the recipe given in the previous Section, the vacuum expectation value of the Wilson operator W γ is defined as What is the functional form of this quantity for large loops? It depends on the value of the occupation probability p. If p is below the percolation threshold p c , then there are only finite-size clusters. If the loop γ has much a larger size than the clusters, then the configurations where W γ (G i ) = 0 necesarily have some cluster located near the loop. The number of these clusters grows linearly with the perimeter |γ| of γ and produces the exponential decay W γ ∝ e −α|γ| . We say that the theory is deconfined. On the contrary if p > p c the theory is confined. Confinement is expected to show up in an area law for the expectation value of large Wilson loops. Indeed in such a case there is an infinite cluster and the number of closed paths linked to γ grows with the area A of the minimal surface Σ encircled by γ. These paths will pierce Σ at points; let N = αA be their mean number. Assuming these points to be randomly distributed on the surface, the probability of finding k such points inside Σ is binomial, Note that only the k = 0 term contributes to the numerator of (4.1), so the expectation value of the Wilson loop becomes One thus apparently obtains an area law decay with string tension σ for any Wilson loop, including those of small size. There is however a flaw in the argument; even if the configurations are obtained by populating each bond (or site) of the lattice independently with a probability p, when all the bridges are erased it is no longer true that the remaining bonds are randomly distributed, as anticipated previously. In fact, since the interaction among the intersection points of Σ is rather weak, we expect an area law only for large enough Wilson loops.
In gauge theories two different confining mechanisms were proposed. One is based on the condensation of center vortices [8]. These objects are stringlike structures which are created by gauge transformations with a non-trivial homotopy associated to the center of the gauge group C(G). In a 3-D lattice the center vortices are represented by a skein of loops forming an infinite network in the confining phase [9]. Each center vortex linked with a Wilson loop contributes to it with a factor of ζ ∈ C(G), thus each configuration contributes with a factor ζ n , where n is the number of linked loops. Although at the end the overall effect is again a decay of W γ with an area law, such a mechanism is slightly different from the one we have described in pure percolation, where a single linked loop suffices to gives a weight zero to the configuration. Another difference is that center vortices carry some conserved charge. For instance, if C(G) = Z N the vortex flux is conserved modulo N , while paths of random percolation do not carry any conserved charge and can intersect freely.
The other confining mechanism is based on the old conjecture [10] that the vacuum behaves like a dual superconductor. The key element of this picture is the monopole condensate which squeezes the gauge field generated by a pair of sources (quarks) into a thin flux tube (the dual version of the Abrikosov vortex). This causes the Wilson loop to decay with an area law. One is led to conjecture that such a thin flux tube should vibrate as a free string [11]. As a consequence, the expectation value of a rectangular Wilson loop of size R × T is expected to have the following asymptotic functional form in the continuum limit [12] W where C, p and σ are functions of the coupling constant and η is the Dedekind function The factor under the square root accounts for the universal quantum contributions of the supposedly string-like flux tube describing the interaction between far-away sources. Random percolation, lacking any non-trivial conserved charge, can hardly account for effects which play the role of magnetic monopole condensation. Notwithstanding this difficulty, we get indirect evidence of the formation of a vibrating string-like flux tube by measuring the universal shape effects it produces, as discussed in the next Section.  We estimated the string tension σ by fitting the mean values of the Wilson loops associated to squares of side R to Eq. (4.4), that in such a case becomes Typically, in a lattice of size L 3 we considered all the squares with R ≤ L/2. The fits for not too small R are very good (see Fig. 4), nevertheless the parameters slightly depend on the value R min of the smallest square included in the fit. Since these formulae are expected to be valid only asymptotically for large values of R, we progressively eliminated the data of lower R until stable parameters were obtained. In order to check for presence of universal shape effects ascribed to the quantum fluctuations of the effective string, we considered, as in [13], the quantity which asymptotically (i.e. large R and R − n) should be only a function of the ratio t = n R , namely, .

(5.3)
Note that it does not contain any adjustable parameters. This function is plotted in Figure 5 and compared with the numerical data for three different values of p. The presence of the expected universal shape effects seems uncontroversial.
The string tension σ is a physical quantity with the dimensions of an inverse square length, hence it is expected to exhibit the following power law behaviour sufficiently close to the percolation threshold where ν is the critical exponent of the correlation length in 3-D percolation. We used the value ν = 0.8765 (16) Table 1.
The percolation threshold p c depends on the lattice and on the kind of percolation (site or bond) which is studied. We checked Eq. (5.4) in the percolating region of three different lattices: simple cubic (SC) with site or bond percolation and body-centered cubic (BCC) with bond percolation. Precise estimates of p c are known and are reported in Tab. 1. In the SC cases we used a lattice of size L 3 = 64 3 , while in the BCC case we had L 3 = 55 3 . In all cases the one-parameter fit to Eq. (5.4) is very good and the scaling window seems rather wide (see Fig. 6). The resulting amplitudes S for the three lattices are reported in Tab. 1.
Noteworthy, in order to extract the string tension σ we have thus far assumed that the square Wilson loops obey the asymptotic form (5.1), where the factor R 1 4 accounts for the contribution of the string fluctuations. If one neglected this factor and only took the area term into account, the χ 2 test of the critical power law (5.4) would grow worse by one or two orders of magnitude, depending on the kind of lattice. We consider this fact another strong evidence of a vibrating, confining string.

Spectrum
As for standard gauge theories, we expect that the confining phase of percolation possesses a rich spectrum of physical states with increasing mass and spin, that we still call glueballs. . These data have been displaced downwards by two orders of magnitude for clarity. The other line is a two-exponential fit to the plaquette data (+). The semi-logarithmic plot makes it evident that in the latter case a single exponential does not suffice.
The basic method that goes into the computation of such mass spectrum is very simple. One first constructs a linear combination of Wilson loops on a fixed time slice of the three-dimensional lattice carrying the quantum numbers of the state one wishes to investigate. One then builds zero momentum operators by summing such a linear combination over the entire spatial lattice. The simplest example of a zero momentum operator coupling to the spin 0 + states is given by where W j ( x) denotes an elementary plaquette variable with base at x = (x, y, t) and orthogonal to the j coordinate axis. According to Section 2, in random percolation the plaquette variable is replaced by the dual link variable ℓ j ( x), defined as equal to 1 if the corresponding bond belongs to the subset B G of loops and null otherwise. The low-lying mass spectrum can be extracted by studying the exponential decay of the connected correlator which is expected to have the following asymptotic expansion where c n denote positive constants and m n are the glueball masses. An example of such a correlator in a SC bond lattice is reported in Fig. 7, where it is evident that at least two different scalar states contribute to C(t). The estimates of the lowest mass in the range 0.258 ≤ p ≤ 0.270 fit well with the scaling form (see Fig. 8) Also the first excited state seems to follow the same power law as expected, though the errors are rather big. Its mass M ′ 0 is about twice M 0 . Such a behaviour is very different from the one observed in the standard twopoint correlation function of the percolation problem, defined as the probability G( x, x ′ ) that the sites x and x ′ are in the same cluster 2 . The corresponding connected, zero-momentum, projection exhibits a single exponential behaviour with a mass term which coincides, within the numerical accuracy, with that of the lowest energy state coupled to the plaquette operator (see Fig. 7). This indicates that the standard cluster correlator couples only to the lowest glueball, while the new observables suggested by the gauge theory interpretation of the percolation disclose a totally unexpected spectrum of physical states. The lightest spinning glueball is the 2 + state. It can be observed in the exponential decay of the correlation function of the operator (6.6) A difficulty encountered in this case is that the signal is drowned within the statistical noise for values of t beyond three or four lattice spacings. In spite of this accuracy problem, one can still verify that in a SC bond lattice the scaling form (6.3) is approximately obeyed with an amplitude In order to have more accurate results on the mass spectrum the basis of the operators should be enlarged to Wilson loops of different shapes, trying to enhance their overlap with the glueball states.

Deconfinement at finite temperature
In quantum field theory the concept of temperature is introduced by simply compactifying the Euclidean time direction and identifying the inverse temperature with the temporal extension of the space-time manifold. Lattice field theories at a temperature T are formulated in a slab which is infinite in the spatial dimensions, but of finite length ℓ = 1/T (in lattice spacing units) and periodic in the remaining temporal direction. In any confining gauge theory there is a critical temperature T c above which the system is deconfined in the sense that for T ≥ T c the string tension σ vanishes. In this Section we demonstrate that the same phenomenon also occurs in random percolation. In the latter case the deconfinement mechanism is particularly transparent, showing its purely geometric origin: Non-vanishing string tension requires an infinite, percolating, cluster. Shrinking the width of the slab reduces the number of possible percolating paths along the spatial directions until the infinite cluster crumbles away, yielding a vanishing σ.
To put it in different terms, note that as the temperature varies from zero to infinity a three-dimensional system is gradually dominated by two-dimensional behaviour; in particular the percolation threshold is a decreasing function of the space dimensions. For instance in the SC bond lattice at T = 0 the percolation threshold is at p c (D = 3) = 0.2488... (see Table 1). At T = ∞ the system reduces to a square lattice, where p c (D = 2) = 1 2 , hence heating a system which at zero temperature lies in the percolating phase with p < 1 2 inevitably undergoes a deconfinement transition at a finite temperature T c . In order to estimate T c in various site or bond percolation lattices we considered a slab of size L x × L y × ℓ with L x = L y = L ≫ ℓ and periodic boundary conditions in all directions (see Fig. 12). We calculated the probability R L (p) for a cluster to wrap around one of the large dimensions. Wrapping probabilities can be defined in different ways: wrapping around the x direction, around the y direction, around either direction, around both directions, etc. To be definite, we consider wrapping around the x direction. For large L it coincides with the probability that the system percolates along x.
We evaluated R L (p) using a very efficient algorithm described by Newman and Ziff [17]. For a bond percolation problem it consists of repeatedly adding a random bond to an initially empty lattice, identifying the clusters joined by the bond and merging them if they are different. At each step one checks whether the touched cluster wraps around x using a clever method described in Ref. [18].
The process stops as soon as a wrapped cluster is detected. In this way one can evaluate the probability Q L (n) that a random configuration with n occupied bonds is wrapping around L x for any n ≤ N , where N is the total number of bonds of the lattice. This method may be adapted to site percolation in a straightforward way. Then one simply finds the required quantity R L (p) for any value of p by convolution with the binomial distribution (7.1) Figure 9 shows some examples of R L (p) in the case of SC site percolating lattice. In the L → ∞ limit this wrapping probability becomes a step function. We have where the threshold value p = p ℓ depends on the type of lattice and on its width ℓ. When L = ∞ the slab system is equivalent to a 2-D torus, where the wrapping probability at criticality has been calculated exactly by Pinson [19] and is R ∞ (p ℓ ) = 0.521058290... . We can use his result to measure the value of p ℓ in a slab of width ℓ by finding the value of p for which This method was first applied in the case of the Ising model in Ref. [20] (with the appropriate value of R ∞ , of course). These estimates in the case of 2-D percolation turn out to scale particularly well with the system size: Newman and Ziff argued that in planar lattices the leading order finite size correction goes like p = p ℓ + c L −2−1/ν2 = p ℓ + c L −11/4 , (7.4) where ν 2 = 4 3 is the thermal exponent of 2-D random percolation. We checked it in the slab geometry finding good agreement for L large enough, as Figure 10 shows in the case of a slab of width ℓ = 6.
We estimated in this way the threshold p ℓ for seven different lattices. The results are reported in Tab. 2. The slab widths (or equivalently the inverse temperatures) were chosen in such a way that the values of p ℓ lie in the scaling region of the string tension σ, as determined in Section 5. In this manner we were able to evaluate also the ratio T c / √ σ using the amplitudes of Tab. 1. It turns out that these ratios for different lattices and widths are clearly compatible with a common value, as required by universality (see last column of Tab. 2).
In (d + 1)-dimensional field theory at finite temperature there is a characteristic interplay between d + 1 and d critical behaviours. This is particularly evident in the present instance: the p ℓ values are extracted through a 2-D percolation power law (7.4) in order to take into account the finite size scaling tied to L. However the tower of p ℓ values as a function of the width ℓ obeys a typical power law of 3-D percolation:  where ν indicates, as in all the other formulae of this paper, the thermal exponent of 3-D percolation, p c is the critical threshold as listed in Tab. 1 and the amplitude T c depends on the kind of lattice. Figure 11 shows the four p ℓ values of the bond SC lattice as a function of ℓ −1/ν . A one-parameter fit to Eq. (7.5) yields T c = 4.45 ± 0.02 .
It should be noted parenthetically that, as the p ℓ 's are essentially two-dimensional quantities, they can be evaluated with high precision. Comparison between Table 1 and Table 2 shows that in some cases the level of precision of those overcomes that of the best estimates of p c . Moreover, with a modest additional computational effort, it would be possible to further improve their precision. Thus one could perhaps envisage to apply systematically Eq. (7.5) to improve the estimates of p c and/or ν.

Universality class of deconfinement
The deconfined phase of any (d + 1)-dimensional gauge theory at finite temperature is characterised by the vanishing of the string tension σ. The interaction between static sources (quarks) is described in terms of Polyakov operators. These are straight Wilson loops wrapped around the short periodic direction ℓ (see Fig. 12). The Polyakov-Polyakov correlator P (0)P ( R) of two parallel Polyakov operators only depends on their relative positions in the d-dimensional sub-lattice. At the deconfining point it is expected to obey a power law dictated by the universality class of the transition. The critical behaviour of gauge theories at the deconfining temperature is well described by the Svetitsky-Yaffe (SY) conjecture [21] which can be formulated as follows. Suppose a d + 1 dimensional gauge theory with gauge group G has a second-order deconfinement transition at a certain temperature T c ; then its universality class is the same of the order-disorder transition of a d-dimensional spin system with a global symmetry group coinciding with the center of the gauge group. In particular, the Polyakov-Polyakov correlator corresponds to the spin-spin correlator of such a d-dimensional system. Therefore, at the critical point, it should decay as where η is the magnetic exponent of the spin model. The validity of this conjecture has been well established in a large number of numerical studies. In the present case the above conjecture requires some generalisation, owing to the fact that not only the center, but the whole gauge group is trivial. Actually the SY conjecture is somehow related to the dimensional crossover in a layered lattice system [22]. The universality class of such a system, as it approaches a critical point, depends on the number of spatial directions which are going to infinity in the thermodynamic limit. This simple observation provides the basis for arguing that the critical behaviour of percolation in a slab of finite thickness is well described by the 2-D percolation universality class. This is also supported by the fact that finite size scaling of threshold probability p ℓ is driven by the thermal exponent ν 2 , as Eq. (7.4) and Figure 10 show. Such a conclusion is much less obvious, and more interesting, when considering the Polyakov-Polyakov correlator. This quantity can be evaluated using exactly the same method described in Section 3 for the Wilson loops. One simply has to take into account that in the Wilson loop the surface Σ (see Fig. 3) used to evaluate cluster wrapping is a rectangle, while in the latter case is a cylinder bounded by the the two Polyakov lines. Periodic boundary conditions put into play another difference: there are now four topologically different surfaces bounded by the two Polyakov lines (see Figure 13), hence the correlator must be written as the sum of these four contributions. When the distance R between Polyakov lines is much smaller than lattice size L the main contribution comes from the top two surfaces of Figure 13.
One instance is reported in Figure 14, where we plot the estimated Polyakov-Polyakov correlator extracted from 10 4 configurations in a 200 × 200 × 6 lattice at the critical value of p 6 as determined in Tab. 2. Plotting these data versus R −η with η = 5 24 shows linear behaviour, as expected for a critical system which lies in the universality class of 2-D percolation.

Conclusions
We have studied some consequences of a new point of view in three-dimensional random percolation which allows us to reinterpret it as a full-fledged gauge theory.
A key difference between the conventional and the present approach is that instead of studying the universal scalings through the size distribution of random clusters we only consider their topological entanglement with suitable closed paths. This suggests a reformulation in the language of gauge theory. In this context a more detailed description of 3-D percolation universality class emerges as a major issue: the gauge-percolation dictionary leads to define new physical quantities that can be used to extract new universal amplitude ratios.
A typical example is the ratio T 2 c /σ between the (square of) deconfining temperature T c and the string tension σ. We evaluated it on seven different lattices (see Tab. 2) finding excellent agreement with universality. Notice that neither T c nor σ are truly foreign concepts of percolation theory: slab percolation has always been a subject of intensive study. Strictly speaking, string tension does not seem to have been considered previously in percolation studies, being a typical notion of gauge theory. There is however an intimately related quantity, the surface tension, which can also be defined in percolation [23]. The novelty introduced by the gauge theory interpretation is a previously unsuspected relationship between slab percolation and surface tension. Another class of universal amplitude ratios came to us as a surprise. It turns out that the Wilson loop correlators receive contribution from a tower of physical states of increasing mass and spin, like in ordinary gauge theories. The ratios of their masses define universal quantities which further characterise the universality class of 3-D percolation.
Transcribing percolation into gauge theory language has also some interesting consequences in the study of quark confinement mechanisms. The confinement generated by a random percolating cluster is similar, but not identical, to that produced by an infinite network of center vortices. In particular the latter carry some conserved charge which regulates their mutual intersections, while the former do not carry any conserved charge and intersect freely. This indicates that the intersection rules of center vortices do not play an important role in producing confinement.
One of the most surprising findings of the present approach is the observation of some universal shape effects in the vacuum expectation value of Wilson loops which have been always ascribed to a different picture of confinement. This picture says that the flux of the gauge field generated by pair of quarks is squeezed by the magnetic monopole condensate into a string-like structure which can vibrate freely. It turns out that these very vibrations generate the universal effects mentioned above. Apparently, the two different pictures of confinement are different descriptions of the same physical phenomenon.
It would be interesting to extend our percolation approach to a 4-D gauge theory. This would require the study of plaquette percolation and the Wilson loops should measure their entanglement with closed surfaces.