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.
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.