We study the percolation properties of geometrical clusters defined in the overlap space of two statistically independent replicas of a square-lattice Ising model that are simulated at the same temperature. In particular, we consider two distinct types of clusters in the overlap, which we dub soft- and hard-constraint clusters, and which are subsets of the regions of constant spin overlap. By means of Monte Carlo simulations and a finite-size scaling analysis we estimate the transition temperature as well as the set of critical exponents characterizing the percolation transitions undergone by these two cluster types. The results suggest that both soft- and hard-constraint clusters percolate at the critical temperature of the Ising model and their critical behavior is governed by the correlation-length exponent ν=1 found by Onsager. At the same time, they exhibit nonstandard and distinct sets of exponents for the average cluster size and percolation strength.
We study the percolation properties of geometrical clusters defined in the overlap space of two statistically independent replicas of a square-lattice Ising model that are simulated at the same temperature. In particular, we consider two distinct types of clusters in the overlap, which we dub soft- and hard-constraint clusters, and which are subsets of the regions of constant spin overlap. By means of Monte Carlo simulations and a finite-size scaling analysis we estimate the transition temperature as well as the set of critical exponents characterizing the percolation transitions undergone by these two cluster types. The results suggest that both soft- and hard-constraint clusters percolate at the critical temperature of the Ising model and their critical behavior is governed by the correlation-length exponent ν=1 found by Onsager. At the same time, they exhibit nonstandard and distinct sets of exponents for the average cluster size and percolation strength.
We consider geometrical or Ising clusters (i.e., domains of parallel spins) in the square lattice random-field Ising model by varying the strength of the Gaussian random field $\ensuremath{\Delta}$. In agreement with the conclusion of a previous investigation [Phys. Rev. E 63, 066109 (2001)], the geometrical correlation length, i.e., the average size of the clusters $\ensuremath{\xi}$ is finite for $\ensuremath{\Delta}>{\ensuremath{\Delta}}_{c}\ensuremath{\approx}1.65$ and divergent for $\ensuremath{\Delta}\ensuremath{\leqslant}{\ensuremath{\Delta}}_{c}$. The scaling function of the distribution of the mass of the clusters as well as the geometrical correlation function are found to involve the scaling exponents of critical percolation. On the other hand, the divergence of the correlation length, $\ensuremath{\xi}(\ensuremath{\Delta})\ensuremath{\sim}{(\ensuremath{\Delta}\ensuremath{-}{\ensuremath{\Delta}}_{c})}^{\ensuremath{-}\ensuremath{\nu}}$, with $\ensuremath{\nu}\ensuremath{\approx}2$, is related to that of tricritical percolation. It is verified numerically that critical geometrical correlations transform conformally.
The fractal dimensions and the percolation exponents of geometrical spin clusters of like sign at criticality are obtained numerically for an Ising model with temperature dependent annealed bond dilution, also known as the thermalized bond Ising model (TBIM), in two dimensions. For this purpose, a modified Wolff single-cluster Monte Carlo simulation is used to generate equilibrium spin configurations on square lattices in the critical region. A tie-breaking rule is employed to identify nonintersecting spin cluster boundaries along the edges of the dual lattice. The values obtained for the fractal dimensions of the spanning geometrical clusters D c , and their interfaces D I , are in perfect agreement with those reported for the standard two-dimensional ferromagnetic Ising model. Furthermore, the variance of the winding angles results in a diffusivity κ = 3 for the two-dimensional thermalized bond Ising model, thus placing it in the universality class of the regular Ising model. A finite-size scaling analysis of the largest geometrical clusters results in a reliable estimation of the critical percolation exponents for the geometrical clusters in the limit of an infinite lattice size. The percolation exponents thus obtained are also found to be consistent with those reported for the regular Ising model. These consistencies are explained in terms of the Fisher renormalization relations, which express the thermodynamic critical exponents of systems with annealed bond dilution in terms of those of the regular model system.
An extensive Monte Carlo study of the two-dimensional Ising model is made to investigate the statistical behavior of spin clusters and interfaces as a function of temperature, T . We use a tie-breaking rule to define interfaces of spin clusters on a square lattice with strip geometry and show that such definitions are consistent with conformally invariant properties of interfaces at the critical temperature, T c . The effective fractal dimensions of spin clusters and interfaces ( d c and d I , respectively) are obtained as a function of temperature. We find that the effective fractal dimension of the spin clusters varies almost linearly with temperature in three different regimes. It is also found that the effective fractal dimension of the interfaces undergoes a sharp crossover around T c , between the values 1 and 1.75 at low and high temperatures, respectively. We also check the finite-size scaling hypothesis for the percolation probability, and the average mass of the largest spin cluster is in a good agreement with the theoretical predictions.
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.
A measure of cluster size heterogeneity (H), introduced by Lee et al. [Phys. Rev. E 84, 020101 (2011)] in the context of explosive percolation, was recently applied to random percolation and to domains of parallel spins in the Ising and Potts models. It is defined as the average number of different domain sizes in a given configuration and a new exponent was introduced to explain its scaling with the size of the system. In thermal spin models, however, physical clusters take into account the temperature-dependent correlation between neighboring spins and encode the critical properties of the phase transition. We here extend the measure of H to these clusters and, moreover, present new results for the geometric domains for both d=2 and 3. We show that the heterogeneity associated with geometric domains has a previously unnoticed double peak, thus being able to detect both the thermal and percolative transitions. An alternative interpretation for the scaling of H that does not introduce a new exponent is also proposed.
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.