Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2603.21303

The phase boundary of the random site Ising model

Riccardo Ben Alì Zinati, Giacomo Gori, Alessandro Codello

Published Mar 22, 2026• 0 citations
arXivHTMLPDF

Abstract

We introduce a new approach to disordered two-dimensional Ising models based on the extension of the combinatorial solution to randomized supercells. Applying it to the site-diluted Ising model on the square lattice, we resolve the full phase boundary $T_c(p)$ from the pure-Ising point to the percolation limit $T_c(p_c)=0$ with, in principle, arbitrary precision. The critical eigenvalue governing the transition is found to follow a remarkably accurate linear interpolation between the Ising and percolation endpoints, whose small but systematic deviations reveal the nontrivial fine structure of the phase boundary. Near the percolation threshold, we confirm the crossover exponent $ϕ_{\rm RSIM}=1$ and extract the nonuniversal amplitude ${α_{\rm RSIM}\simeq 1.616}$.

DOI
10.48550/arxiv.2603.21303

Related papers

Transfer-matrix scaling for diluted Ising systems

S. L. A. de Queiroz, Robin B. Stinchcombe

A transfer-matrix scaling technique is developed for randomly diluted systems and applied to the site-diluted Ising model on a square lattice. For each connected configuration between adjacent columns, the contribution of the respective transfer matrix to the decay of correlations is considered only as far as the ratio of the two largest eigenvalues, allowing an economical incorporation of configurational averages. Predictions for the phase boundary at and near the percolation and pure Ising limits, and for the correlation exponent \ensuremath{\eta} at those limits, agree with exactly known results to within 1% error with largest strip widths of only L=5. The exponent \ensuremath{\eta} remains near the pure value (1/4) for all intermediate concentrations until it turns over to the percolation value at the threshold.

Physical review. B, Condensed matter•Sep 1, 1992•16 citations

Criticality and crossover in the bond-diluted random Ising model

Eytan Domany

By using the replica trick and a duality transformation, the bond-diluted random Ising model is mapped onto a new Hamiltonian. This demonstrates the higher order critical nature of the percolation point and identifies the appropriate crossover scaling variables. Taking the n to 0 limit of the replica method near the percolation point is shown to be equivalent to the q to 1 limit of the Potts model. The critical line near p c is calculated, yielding for a square lattice exp(-2K c )=2ln2(p- 1 / 2 )+0(p- 1 / 2 ) 2 .

Journal of Physics C Solid State Physics•Apr 28, 1978•52 citations

Crossover and self-averaging in the two-dimensional site-diluted Ising model: Application of probability-changing cluster algorithm

Yusuke Tomita, Yutaka Okabe

Using the newly proposed probability-changing cluster (PCC) Monte Carlo algorithm, we simulate the two-dimensional (2D) site-diluted Ising model. Since we can tune the critical point of each random sample automatically with the PCC algorithm, we succeed in studying the sample-dependent T(c)(L) and the sample average of physical quantities at each T(c)(L) systematically. Using the finite-size scaling (FSS) analysis for T(c)(L), we discuss the importance of corrections to FSS both in the strong-dilution and weak-dilution regions. The critical phenomena of the 2D site-diluted Ising model are shown to be controlled by the pure fixed point. The crossover from the percolation fixed point to the pure Ising fixed point with the system size is explicitly demonstrated by the study of the Binder parameter. We also study the distribution of critical temperature T(c)(L). Its variance shows the power-law L dependence, L(-n), and the estimate of the exponent n is consistent with the prediction of Aharony and Harris [Phys. Rev. Lett. 77, 3700 (1996)]. Calculating the relative variance of critical magnetization at the sample-dependent T(c)(L), we show that the 2D site-diluted Ising model exhibits weak self-averaging.

arXiv:cond-mat/0106154•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Aug 27, 2001•32 citations

Percolation in the two-dimensional Ising model

Tao Chen, Jinhong Zhu, Wei Zhong +2 more

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.

arXiv:2504.18861•Physical review. E•Jul 22, 2025•2 citations

Transfer-matrix scaling from disorder-averaged correlation lengths for diluted Ising systems

S. L. A. de Queiroz, Robin B. Stinchcombe

A transfer-matrix-scaling technique is developed for randomly diluted systems, and applied to the site-diluted Ising model on a square lattice in two dimensions. For each allowed disorder configuration between two adjacent columns, the contribution of the respective transfer matrix to the decay of correlations is considered only as far as the ratio of its two largest eigenvalues, allowing an economical calculation of a configuration-averaged correlation length. Standard phenomenological-renormalization procedures are then used to analyze aspects of the phase boundary which are difficult to assess accurately by alternative methods. For magnetic site concentration p close to ${\mathit{p}}_{\mathit{c}}$, the extent of exponential behavior of the ${\mathit{T}}_{\mathit{c}}$\ifmmode\times\else\texttimes\fi{}p curve is clearly seen for over two decades of variation of p-${\mathit{p}}_{\mathit{c}}$. Close to the pure-system limit, the exactly known reduced slope is reproduced to a very good approximation, though with nonmonotonic convergence. The averaged correlation lengths are inserted into the exponent-amplitude relationship predicted by conformal invariance to hold at criticality. The resulting exponent \ensuremath{\eta} remains near the pure value (1/4) for all intermediate concentrations until it crosses over to the percolation value at the threshold.

arXiv:cond-mat/9407075•Physical review. B, Condensed matter•Oct 1, 1994•25 citations

Critical and thermodynamic properties of the randomly dilute Ising model

Ciriyam Jayaprakash, Eberhard K. Riedel, Michael Wortis

The randomly bond-dilute two-dimensional nearest-neighbor Ising model on the square lattice is studied by renormalization-group methods based on the Migdal-Kadanoff approximate recursion relations. Calculations give both thermal and magnetic exponents associated with the percolative fixed point. Differential recursion relations yield a phase diagram which is in quantitative agreement with all known results. Curves for the specific heat, percolation probability, and magnetization are displayed. The critical region of the specific heat becomes unobservably narrow well above the percolation threshold ${p}_{c}$. This provides a possible explanation for the apparent specific-heat rounding in certain experiments.

Physical review. B, Condensed matter•Sep 1, 1978•114 citations
We introduce a new approach to disordered two-dimensional Ising models based on the extension of the combinatorial solution to randomized supercells. Applying it to the site-diluted Ising model on the square lattice, we resolve the full phase boundary Tc(p)T_c(p)Tc​(p) from the pure-Ising point to the percolation limit Tc(pc)=0T_c(p_c)=0Tc​(pc​)=0 with, in principle, arbitrary precision. The critical eigenvalue governing the transition is found to follow a remarkably accurate linear interpolation between the Ising and percolation endpoints, whose small but systematic deviations reveal the nontrivial fine structure of the phase boundary. Near the percolation threshold, we confirm the crossover exponent ϕRSIM=1ϕ_{\rm RSIM}=1ϕRSIM​=1 and extract the nonuniversal amplitude αRSIM≃1.616{α_{\rm RSIM}\simeq 1.616}αRSIM​≃1.616.
A transfer-matrix scaling technique is developed for randomly diluted systems and applied to the site-diluted Ising model on a square lattice. For each connected configuration between adjacent columns, the contribution of the respective transfer matrix to the decay of correlations is considered only as far as the ratio of the two largest eigenvalues, allowing an economical incorporation of configurational averages. Predictions for the phase boundary at and near the percolation and pure Ising limits, and for the correlation exponent \ensuremath{\eta} at those limits, agree with exactly known results to within 1% error with largest strip widths of only L=5. The exponent \ensuremath{\eta} remains near the pure value (1/4) for all intermediate concentrations until it turns over to the percolation value at the threshold.
By using the replica trick and a duality transformation, the bond-diluted random Ising model is mapped onto a new Hamiltonian. This demonstrates the higher order critical nature of the percolation point and identifies the appropriate crossover scaling variables. Taking the n to 0 limit of the replica method near the percolation point is shown to be equivalent to the q to 1 limit of the Potts model. The critical line near p c is calculated, yielding for a square lattice exp(-2K c )=2ln2(p- 1 / 2 )+0(p- 1 / 2 ) 2 .
Using the newly proposed probability-changing cluster (PCC) Monte Carlo algorithm, we simulate the two-dimensional (2D) site-diluted Ising model. Since we can tune the critical point of each random sample automatically with the PCC algorithm, we succeed in studying the sample-dependent T(c)(L) and the sample average of physical quantities at each T(c)(L) systematically. Using the finite-size scaling (FSS) analysis for T(c)(L), we discuss the importance of corrections to FSS both in the strong-dilution and weak-dilution regions. The critical phenomena of the 2D site-diluted Ising model are shown to be controlled by the pure fixed point. The crossover from the percolation fixed point to the pure Ising fixed point with the system size is explicitly demonstrated by the study of the Binder parameter. We also study the distribution of critical temperature T(c)(L). Its variance shows the power-law L dependence, L(-n), and the estimate of the exponent n is consistent with the prediction of Aharony and Harris [Phys. Rev. Lett. 77, 3700 (1996)]. Calculating the relative variance of critical magnetization at the sample-dependent T(c)(L), we show that the 2D site-diluted Ising model exhibits weak self-averaging.
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.
A transfer-matrix-scaling technique is developed for randomly diluted systems, and applied to the site-diluted Ising model on a square lattice in two dimensions. For each allowed disorder configuration between two adjacent columns, the contribution of the respective transfer matrix to the decay of correlations is considered only as far as the ratio of its two largest eigenvalues, allowing an economical calculation of a configuration-averaged correlation length. Standard phenomenological-renormalization procedures are then used to analyze aspects of the phase boundary which are difficult to assess accurately by alternative methods. For magnetic site concentration p close to pc{\mathit{p}}_{\mathit{c}}pc​, the extent of exponential behavior of the Tc{\mathit{T}}_{\mathit{c}}Tc​\ifmmode\times\else\texttimes\fi{}p curve is clearly seen for over two decades of variation of p-pc{\mathit{p}}_{\mathit{c}}pc​. Close to the pure-system limit, the exactly known reduced slope is reproduced to a very good approximation, though with nonmonotonic convergence. The averaged correlation lengths are inserted into the exponent-amplitude relationship predicted by conformal invariance to hold at criticality. The resulting exponent \ensuremath{\eta} remains near the pure value (1/4) for all intermediate concentrations until it crosses over to the percolation value at the threshold.
The randomly bond-dilute two-dimensional nearest-neighbor Ising model on the square lattice is studied by renormalization-group methods based on the Migdal-Kadanoff approximate recursion relations. Calculations give both thermal and magnetic exponents associated with the percolative fixed point. Differential recursion relations yield a phase diagram which is in quantitative agreement with all known results. Curves for the specific heat, percolation probability, and magnetization are displayed. The critical region of the specific heat becomes unobservably narrow well above the percolation threshold pc{p}_{c}pc​. This provides a possible explanation for the apparent specific-heat rounding in certain experiments.