We study the effects of dilution to the critical properties of site-diluted Ising model in two dimensions using Monte Carlo simulations. Quenched disorder from the dilution is incorporated into the Ising model via random empty sites on the square lattice of Ising spins. Thermodynamic quantities such as the magnetization
M per spin, energy
E per spin, magnetic susceptibility
χ per spin, and specific heat
C per spin are then calculated after the system has equilibrated. At small dilution concentrations
d<0.1, we find that the value of the critical exponent
β does not deviate from its pure Ising value. At higher dilution concentrations
d>0.1, however, we find
β to strongly depend on the value of
d. We are able to locate a critical temperature
Tc and a critical dilution concentration
dc where the phase transition occurs. We find
Tc to depend linearly on
d. In the phase diagrams of
M,
E,
χ, and
C, we find that the phase transition line eventually disappears at high dilutions. Our results suggest that there is a transition from Strong Universality at low dilution to Weak Universality at high dilution. Lastly, we find a wide and nonsingular increase in the magnetic susceptibility
χ at the low temperature and high dilution region.
We consider the two-dimensional randomly site diluted Ising model and the random-bond +/-J Ising model (also called the Edwards-Anderson model), and study their critical behavior at the paramagnetic-ferromagnetic transition. The critical behavior of thermodynamic quantities can be derived from a set of renormalization-group equations, in which disorder is a marginally irrelevant perturbation at the two-dimensional Ising fixed point. We discuss their solutions, focusing in particular on the universality of the logarithmic corrections arising from the presence of disorder. Then, we present a finite-size scaling analysis of high-statistics Monte Carlo simulations. The numerical results confirm the renormalization-group predictions, and in particular the universality of the logarithmic corrections to the Ising behavior due to quenched dilution.
We present a finite-size scaling analysis of high-statistics Monte Carlo simulations of the three-dimensional randomly site-diluted and bond-diluted Ising model. The critical behaviour of these systems is affected by slowly decaying scaling corrections which make the accurate determination of their universal asymptotic behaviour quite hard, requiring an effective control of the scaling corrections. For this purpose we exploit improved Hamiltonians, for which the leading scaling corrections are suppressed for any thermodynamic quantity, and improved observables, for which the leading scaling corrections are suppressed for any model belonging to the same universality class. The results of the finite-size scaling analysis provide strong numerical evidence that phase transitions in three-dimensional randomly site-diluted and bond-diluted Ising models belong to the same randomly dilute Ising universality class. We obtain accurate estimates of the critical exponents, ν = 0.683(2), η = 0.036(1), α = −0.049(6), γ = 1.341(4), β = 0.354(1), δ = 4.792(6), and of the leading and next-to-leading correction-to-scaling exponents, ω = 0.33(3) and ω 2 = 0.82(8).
We present a different way of probing the universality class of the site-diluted two-dimensional Ising model. We analyze Monte Carlo data for the magnetic susceptibility, introducing a fitting procedure in the critical region applicable even for a single sample with quenched disorder. This gives us the possibility to fit simultaneously the critical exponent, the critical amplitude, and the sample-dependent pseudocritical temperature. The critical amplitude ratio of the magnetic susceptibility is seen to be independent of the concentration q of the empty sites for all investigated values of q < or =0.25. At the same time the average effective exponent gamma(eff) is found to vary with the concentration q, which may be argued to be due to logarithmic corrections to the power law of the pure system. These corrections are canceled in the susceptibility amplitude ratio as predicted by theory. The central charge of the corresponding field theory was computed and compared well with the theoretical predictions.