The critical behavior of the disordered ferromagnetic Ising model is studied numerically by the Monte Carlo method in a wide range of variation of concentration of nonmagnetic impurity atoms. The temperature dependences of correlation length and magnetic susceptibility are determined for samples with various spin concentrations and various linear sizes. The finite-size scaling technique is used for obtaining scaling functions for these quantities, which exhibit a universal behavior in the critical region; the critical temperatures and static critical exponents are also determined using scaling corrections. On the basis of variation of the scaling functions and values of critical exponents upon a change in the concentration, the conclusion is drawn concerning the existence of two universal classes of the critical behavior of the diluted Ising model with different characteristics for weakly and strongly disordered systems.
The critical behavior of the disordered two-dimensional antiferromagnetic Potts model with the number of spin states q= 3 on a triangular lattice with disorder in the form of nonmagnetic impurities is studied by the Monte Carlo method. The critical exponents for the susceptibility γ, magnetization β, specific heat α, and correlation radius ν are calculated in the framework of the finite-size scaling theory at spin concentrations p = 0.90, 0.80, 0.70, and 0.65. It is found that the critical exponents increase with the degree of disorder, whereas the ratios and do not change, thus holding the scaling equality
ν2β+νγ=d . Such behavior of the critical exponents is related to the weak universality of the critical behavior characteristic of disordered systems. All results are obtained using independent Monte Carlo algorithms, such as the Metropolis and Wolff algorithms.
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.
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.
Abstract. The critical properties of the randomly site-diluted two-dimensional Ising model were studied using the Wang-Landau algorithm. The concentration of nonmagnetic sites was q = 0.1; the remaining sites were occupied by magnetic particles. The study was carried out in the appro-priate restricted but dominant energy subspaces. The main effort was focused on the specific heat and magnetic susceptibility by using the density of states for the model for lattices with linear size L = 20-120. The finite size scaling behavior of the specific heat and susceptibility was studied. The ratio of the critical exponents (α/ν) is negative (consistent with the Harris criterion), while the ratio (γ/ν) appears to assume its pure Ising model value (weak universality). 1.
The critical relaxation of the magnetization in a two-dimensional Ising model with quenched random non-magnetic impurities has been studied by numerical simulation. A squared lattice of size 400 2 with spin concentrations p=1.0; 0.95, 0,9, 0.85, 0.8, 0.75 and 0.7 was considered. The dynamical critical exponent z was determined by the Monte Carlo method combined with the dynamical renormalization group method. The following values were received for z(p): z(1)=2.24+or-0.07, z(0.95)=2.24+or-0.06, z(0.9)=2.24+or-0.06, z(0.85)=2.38+or-0.05, z(0.8)=2.51+or-0.06, z(0.75)=2.66+or-0.07 and z(0.7)=2.88+or-0.06. A singular dynamic scaling behaviour with z=A' mod ln(p-p c ) mod +B' was found for systems with spin concentrations p<or=0.85 near the percolation threshold p c , where A'=0.56+or-0.07 and B'=1.62+or-0.07 are temperature-independent constants.