The study of the Ising model from a percolation perspective has played a significant role in the modern theory of critical phenomena. We consider the celebrated square-lattice Ising model and construct percolation clusters by placing bonds, with probability p, between any pair of parallel spins within an extended range beyond nearest neighbors. At the Ising criticality, we observe two percolation transitions as p increases: starting from a disordered phase with only small clusters, the percolation system enters into a stable critical phase that persists over a wide range p_{c_{1}}<p<p_{c_{2}}, and then develops a long-ranged percolation order with giant clusters for both up and down spins. At p_{c1} and for the stable critical phase, the critical behaviors agree well with those for the Fortuin-Kasteleyn random clusters and the spin domains of the Ising model, respectively. At p_{c2}, the fractal dimension of clusters and the scaling exponent along p direction are estimated as y_{h2}=1.9580(6) and y_{p2}=0.552(9), of which the exact values remain unknown. These findings reveal interesting geometric properties of the two-dimensional Ising model that has been studied for more than 100 years.
The study of the Ising model from a percolation perspective has played a significant role in the modern theory of critical phenomena. We consider the celebrated square-lattice Ising model and construct percolation clusters by placing bonds, with probability p, between any pair of parallel spins within an extended range beyond nearest neighbors. At the Ising criticality, we observe two percolation transitions as p increases: starting from a disordered phase with only small clusters, the percolation system enters into a stable critical phase that persists over a wide range p_{c_{1}}<p<p_{c_{2}}, and then develops a long-ranged percolation order with giant clusters for both up and down spins. At p_{c1} and for the stable critical phase, the critical behaviors agree well with those for the Fortuin-Kasteleyn random clusters and the spin domains of the Ising model, respectively. At p_{c2}, the fractal dimension of clusters and the scaling exponent along p direction are estimated as y_{h2}=1.9580(6) and y_{p2}=0.552(9), of which the exact values remain unknown. These findings reveal interesting geometric properties of the two-dimensional Ising model that has been studied for more than 100 years.
Geometric representations provide a useful perspective on critical phenomena in the Ising model. In a recent study [Phys. Rev. E 112, 034118 (2025)], we found that the two-dimensional critical Ising model exhibits two consecutive percolation transitions for geometric spin clusters as the bond-occupation probability p between parallel spins increases. Here, through extensive Monte Carlo simulations, we show that this phenomenon does not persist in three dimensions, where we observe only a single percolation transition on critical Ising configurations. Further theoretical analysis of the Ising model on the complete graph also yields the same scenario. In addition, we study percolation on a two-dimensional layer embedded in the three-dimensional critical Ising model. For this layer system, we estimate the red-bond exponent yp=0.426(6) and the fractal dimensions of the largest cluster, hull, and shortest path as df=1.8926(20), dhull=1.663(4), and dmin=1.080(10), respectively. These values indicate a distinct universality class induced by coupling to out-of-plane critical correlations.
In this paper we consider the Ising model (called the dynamical Ising model with dynamical temperature T) on the Ising-correlated percolation lattice. To construct this lattice, another Ising model (which is called the quenched Ising model with quench artificial temperature Tq) is employed to model the correlations between imperfections. We argue that for each quenched temperature, there is a particular dynamical temperature Tc(Tq) for which the dynamical Ising model becomes critical. We present some finite-size arguments (based on the moments of the magnetism and the energy as well as the spanning cluster probability) to extract the critical points and also show that they are compatible with the finite-size scaling of the singular point of the magnetic susceptibility. The model is thoroughly characterized in and out of these points. We find that the critical behaviors of the model change significantly with respect to the regular Ising model as well as the Ising model on the uncorrelated percolation lattice. It is shown that for the critical lattice, ${T}_{q}={T}_{c}^{\text{square}\phantom{\rule{4.pt}{0ex}}\text{Ising}\phantom{\rule{4.pt}{0ex}}\text{model}}\ensuremath{\approx}2.269$, the critical temperature for the dynamical Ising model in the thermodynamic limit is ${T}_{c}({T}_{q}={T}_{c}^{\text{square}\phantom{\rule{4.pt}{0ex}}\text{Ising}\phantom{\rule{4.pt}{0ex}}\text{model}})=1.94\ifmmode\pm\else\textpm\fi{}0.005$, and the fractal dimension of the exterior perimeter of geometrical spin clusters is ${D}_{f}^{T={T}_{c}({T}_{q}=2.269)}=1.408\ifmmode\pm\else\textpm\fi{}0.002$. Many quantities, such as the dynamical critical temperatures, all local and global critical exponents, and the fractal dimension of loops Df, scale with the quench temperature in a power-law fashion, with some critical exponents that are reported. Significantly we see that ${D}_{f}^{T={T}_{c}({T}_{q})}\ensuremath{-}{D}_{f}^{T={T}_{c}({T}_{q}=2.269)}\ensuremath{\sim}\frac{1}{\sqrt{\ensuremath{\zeta}({T}_{q})}}$ in which $\ensuremath{\zeta}({T}_{q})$ is the correlation length of the quenched Ising model at temperature Tq.
Large-deviation tails of critical probability distributions provide a sensitive probe of universality beyond standard finite-size scaling. We study these tails for critical percolation and Fortuin--Kasteleyn Ising models on two-dimensional lattices, three-dimensional lattices, and complete graphs. We consider two rescaled order parameters: the magnetization-like variable xm=∣M∣/⟨∣M∣⟩, including a signed cluster-mass analogue for percolation, and the largest-cluster variable xC=C1/⟨C1⟩. For xm, we test the expected stretched-exponential large-deviation tail and show that the same form applies to the percolation analogue. For xC, guided by the exact complete-graph result and scaling arguments, we propose universal scaling forms for both tails of the cumulative distribution and test them by extensive Monte Carlo simulations. In the complete-graph FK-Ising model, the left tail is governed by rare configurations with percolation-like scaling rather than by the typical Ising scaling. Our results show that the tails of order-parameter distributions reveal universal features of critical fluctuations that are not captured by averaged observables alone.
We simulate the two-dimensional XY model in the flow representation by a worm-type algorithm, up to linear system size L=4096, and study the geometric properties of the flow configurations. As the coupling strength K increases, we observe that the system undergoes a percolation transition K_{perc} from a disordered phase consisting of small clusters into an ordered phase containing a giant percolating cluster. Namely, in the low-temperature phase, there exhibits a long-ranged order regarding the flow connectivity, in contrast to the quasi-long-range order associated with spin properties. Near K_{perc}, the scaling behavior of geometric observables is well described by the standard finite-size scaling ansatz for a second-order phase transition. The estimated percolation threshold K_{perc}=1.1053(4) is close to but obviously smaller than the Berezinskii-Kosterlitz-Thouless (BKT) transition point K_{BKT}=1.1193(10), which is determined from the magnetic susceptibility and the superfluid density. Various interesting questions arise from these unconventional observations, and their solutions would shed light on a variety of classical and quantum systems of BKT phase transitions.
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.
To demonstrate the application of computational physics in modern statistical physics, this paper investigates the critical behavior of a finite-size ferromagnetic system based on the two-dimensional square-lattice Ising model.A nearest-neighbor Ising Hamiltonian is constructed, and periodic boundary conditions are adopted to reduce finite-boundary effects.The Metropolis Monte Carlo algorithm is used to sample spin configurations at different temperatures, while energy, magnetization, specific heat, magnetic susceptibility and Binder cumulant are calculated to describe the transition from an ordered state to a disordered state.The simulation results show that the system exhibits clear long-range order at low temperatures, where the order parameter remains high.As temperature increases, thermal fluctuations become stronger, magnetization decreases rapidly, and the system gradually evolves into a disordered state.The specific heat curve shows a pronounced peak at approximately (T = 2.269), while the susceptibility peak shifts slightly toward the high-temperature side due to finite-size effects and the discrete temperature grid.Overall, the simulated critical region is consistent with the theoretical critical temperature of the two-dimensional square-lattice Ising model.Spin-configuration analysis further reveals enhanced fluctuations and domain reconstruction near the critical region.In addition, a thermalization criterion and an extended sampling test for L = 48 near the critical region are added to examine the influence of finite sampling and equilibration.These results indicate that Monte Carlo simulation of the two-dimensional Ising model can effectively reproduce continuous phase transition, critical fluctuation and finite-size effects, providing a basic computational case for understanding critical phenomena in complex systems and data-driven physical analysis.