We have performed Monte Carlo calculations of the thermodynamic properties of the bond-dilute Ising model on the square lattice. The results are obtained for the magnetization, internal energy, specific heat, and susceptibility for several concentrations of missing bonds. The dependence of the Curie temperature on dilution is in good agreement with the results obtained by other methods. The phase transition is found to remain sharp with the same value of the susceptibility exponent as for the pure Ising model.
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 report extensive Monte Carlo simulations of disordered Ising systems in the ferromagnetic region with concentrations of magnetic sites between p=1.0 and 0.5. The magnetization, the susceptibility, and the caloric properties have been studied in the critical region. The critical exponents \ensuremath{\beta} and \ensuremath{\gamma} as well as the universal amplitude A of the magnetization and the ratio
C+/${\mathit{C}}_{\mathrm{\ensuremath{-}}}$ of the susceptibility amplitudes have been determined with high precision. In addition to the cusplike specific heat, we have also measured the magnetization-energy correlation function \ensuremath{\Gamma}, which is a divergent thermal quantity. The corresponding critical exponent \ensuremath{\zeta}, which is related to the other exponents by \ensuremath{\zeta}=1-\ensuremath{\beta}=(\ensuremath{\alpha}+\ensuremath{\gamma})/2, has been determined. We have found that all quantities show power-law behavior within the temperature range of our simulation. All critical exponents change continuously with dilution. Even in the range of weak dilution (p\ensuremath{\ge}0.8), the effective critical exponents are concentration dependent and are clearly different from their pure system values. In the strongly diluted regime the critical exponents gradually reach new asymptotic behavior at p=0.5--0.6 with values of \ensuremath{\beta}=0.335\ifmmode\pm\else\textpm\fi{}0.01 and \ensuremath{\gamma}=1.49\ifmmode\pm\else\textpm\fi{}0.02. The exponent \ensuremath{\alpha} of the specific heat becomes -0.17\ifmmode\pm\else\textpm\fi{}0.04, which corresponds to a cusplike singularity. We conclude that disorder profoundly changes the critical behavior for weakly as well as strongly disordered spin systems.
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.