After fifty years of lattice gauge theories (LGTs), the nature of the transition between their topological phases (confinement or deconfinement) remains challenging due to the absence of a local order parameter. In this work, we conduct a percolation analysis of Wegner's three-dimensional Z_{2} lattice gauge model using intensive Monte Carlo simulations and finite-size scaling, offering fresh insights into the topological phase transitions of gauge-invariant systems. We demonstrate that, regardless of the connection rules, geometrical loops, constructed by piercing excited plaquettes percolate precisely at the thermal critical point T_{c}, with critical exponents coinciding with those of the loop representation of the dual 3D Ising model. Further, we construct Fortuin-Kasteleyn (FK) clusters in a random-cluster representation, showing that they also percolate at T_{c}, enabling access to all thermal critical exponents. Strikingly, the Binder cumulants of the percolation order parameters for both loops and FK clusters reveal a pseudo-first-order transition. This work sheds new light on the critical behavior of pure LGTs, with potential implications for condensed matter systems and quantum error correction.
After fifty years of lattice gauge theories (LGTs), the nature of the transition between their topological phases (confinement or deconfinement) remains challenging due to the absence of a local order parameter. In this work, we conduct a percolation analysis of Wegner's three-dimensional Z_{2} lattice gauge model using intensive Monte Carlo simulations and finite-size scaling, offering fresh insights into the topological phase transitions of gauge-invariant systems. We demonstrate that, regardless of the connection rules, geometrical loops, constructed by piercing excited plaquettes percolate precisely at the thermal critical point T_{c}, with critical exponents coinciding with those of the loop representation of the dual 3D Ising model. Further, we construct Fortuin-Kasteleyn (FK) clusters in a random-cluster representation, showing that they also percolate at T_{c}, enabling access to all thermal critical exponents. Strikingly, the Binder cumulants of the percolation order parameters for both loops and FK clusters reveal a pseudo-first-order transition. This work sheds new light on the critical behavior of pure LGTs, with potential implications for condensed matter systems and quantum error correction.
Abstract Wegner’s Z 2 gauge model is the earliest formulation of pure lattice gauge theory and predicts the topological nature of the confinement–deconfinement transition. In three dimensions ( D = 3 ), the equilibrium critical behavior of the model is understood in terms of geometrically defined objects, namely loop excitations and Fortuin–Kasteleyn clusters. This work investigates the critical quench dynamics of this model from a geometric perspective, following quenches from both a high-temperature percolation phase and the zero-temperature ground state. Using time-dependent finite-size scaling analysis, we find that the critical non-equilibrium relaxation of the percolation order parameter is governed by a dynamical exponent z p ∼ 2.6 , consistent with that associated with the energy density, z c . Importantly, the value of z p is robust with respect to the initial quench condition and the choice of geometrical objects. Furthermore, we provide a detailed characterization of the kinetics of different geometrical objects during the evolution from the percolation phase. Notably, we observe that the quench dynamics obeys dynamic scaling in terms of a growing lengthscale, ξ p ( t ) ∼ t 1 / z p , despite the absence of a local order parameter.
Lattice gauge theories (LGTs) were introduced in 1974 by Wilson to study quark confinement. These models have been shown to exhibit (de)confined phases, yet it remains challenging to define experimentally accessible order parameters. Here we propose percolation-inspired order parameters (POPs) to probe confinement of dynamical matter in Z2 LGTs using electric field basis snapshots accessible to quantum simulators. We apply the POPs to study a classical Z2 LGT and find a confining phase up to temperature $T=\ensuremath{\infty}$ in two dimensions (critical Tc, i.e., finite-T phase transition, in three dimensions) for any nonzero density of Z2 charges. Further, using quantum Monte Carlo we demonstrate that the POPs reproduce the square lattice Fradkin-Shenker phase diagram at T=0 and explore the phase diagram at T>0. The correlation length exponent coincides with the one of the three-dimensional Ising universality class and we determine the POP critical exponent characterizing percolation. Our proposed POPs provide a geometric perspective of confinement and are directly accessible to snapshots obtained in quantum simulators, making them suitable as a probe for quantum spin liquids.
We present a Monte Carlo study of the geometric properties of Fortuin-Kasteleyn (FK) clusters of the Ising model on square [two-dimensional (2D)] and simple-cubic [three-dimensional (3D)] lattices. The wrapping probability, a dimensionless quantity characterizing the topology of the FK clusters on a torus, is found to suffer from smaller finite-size corrections than the well-known Binder ratio and yields a high-precision critical coupling as Kc(3D)=0.221654631(8). We then study other geometric properties of FK clusters at criticality. It is demonstrated that the distribution of the critical largest-cluster size C1 follows a single-variable function as $P({C}_{1},L)\phantom{\rule{0.16em}{0ex}}d{C}_{1}=\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{P}(x)\phantom{\rule{0.16em}{0ex}}dx$ with $x\ensuremath{\equiv}{C}_{1}/{L}^{{d}_{\mathrm{F}}}$ (L is the linear size), where the fractal dimension dF is identical to the magnetic exponent. An interesting bimodal feature is observed in distribution $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{P}(x)$ in three dimensions, and attributed to the different approaching behaviors for $K\ensuremath{\rightarrow}{K}_{c}+{0}^{\ifmmode\pm\else\textpm\fi{}}$. To characterize the compactness of the FK clusters, we measure their graph distances and determine the shortest-path exponents as dmin(3D)=1.25936(12) and dmin(2D)=1.0940(2). Further, by excluding all the bridges from the occupied bonds, we obtain bridge-free configurations and determine the backbone exponents as dB(3D)=2.1673(15) and dB(2D)=1.7321(4). The estimates of the universal wrapping probabilities for the 3D Ising model and of the geometric critical exponents dmin and dB either improve over the existing results or have not been reported yet. We believe that these numerical results would provide a testing ground in the development of further theoretical treatments of the 3D Ising model.
Self-dual criticality in three-dimensional Z2 gauge theory with matter
The simplest topologically ordered phase in 2+1D is the deconfined phase of Z2 gauge theory, realized for example in the toric code. This phase permits a duality that exchanges electric and magnetic excitations (``e'' and ``m'' particles). Condensing either particle while the other remains gapped yields a phase transition with 3D Ising exponents. More mysterious, however, is the transition out of the deconfined phase when self-duality symmetry is preserved. If this transition is continuous, which has so far been unclear, then it may be the simplest critical point for which we still lack any useful continuum Lagrangian description. This transition also has a soft matter interpretation, as a multicritical point for classical membranes in 3D. We study the self-dual transition with Monte Carlo simulations of the Z2 gauge-Higgs model on cubic lattices of linear size L≤96. Our results indicate a continuous transition: for example, cumulants show a striking parameter-free scaling collapse. We estimate scaling dimensions by using duality symmetry to distinguish the leading duality-odd/duality-even scaling operators A and S. All local operators have large scaling dimensions, making standard techniques for locating the critical point ineffective. We develop an alternative using ``renormalization group trajectories'' of cumulants. We check that two- and three-point functions, and temporal correlators in the Monte-Carlo dynamics, are scale-invariant, with scaling dimensions xA and xS and dynamical exponent z. We also give a picture for emergence of 1-form symmetries, in some parts of the phase diagram, in terms of ``patching'' of membranes/worldsurfaces. We relate this to the percolation of anyon worldlines in spacetime. Analyzing percolation yields a fourth exponent for the self-dual transition. We propose variations of the model for further investigation.
Abstract In the ordered phase of the 3D Ising model, minority spin clusters are surrounded by a boundary of dual plaquettes. As the temperature is raised, these spin clusters become more numerous, and it is found that eventually their boundaries undergo a percolation transition when about 13% of spins are minority. Boundary percolation differs from the more commonly studied site and link percolation, although it is related to an unusual type of site percolation that includes next to nearest neighbor relationships. Because the Ising model can be reformulated in terms of the domain boundaries alone, there is reason to believe boundary percolation should be relevant here. A symmetry-breaking order parameter is found in the dual theory, the 3D gauge Ising model. It is seen to undergo a phase transition at a coupling close to that predicted by duality from the boundary percolation. This transition lies in the disordered phase of the gauge theory and has the nature of a spin-glass transition. Its critical exponent ν ∼ 1.3 is seen to match the finite-size shift exponent of the percolation transition further cementing their connection. This predicts a very weak specific heat singularity with exponent α ∼ − 1.9 . The third energy cumulant fits well to the expected non-infinite critical behavior in a manner consistent with both the predicted exponent and critical point, indicating a true thermal phase transition. Unlike random boundary percolation, the Ising boundary percolation has two different ν exponents, one associated with largest-cluster scaling and the other with finite-size transition-point shift. This suggests there may be two different correlation lengths present.
The fractal structure and scaling properties of a 2d slice of the 3d Ising model is studied using Monte Carlo techniques. The percolation transition of geometric spin (GS) clusters is found to occur at the Curie point, reflecting the critical behavior of the 3d model. The fractal dimension and the winding angle statistics of the perimeter and external perimeter of the geometric spin clusters at the critical point suggest that, if conformally invariant in the scaling limit, they can be described by the theory of Schramm-Löwner evolution (SLE κ ) with diffusivity of κ=5 and 16/5, respectively, putting them in the same universality class as the interfaces in 2d tricritical Ising model. It is also found that the Fortuin-Kasteleyn (FK) clusters associated with the cross-sections undergo a nontrivial percolation transition, in the same universality class as the ordinary 2d critical percolation.