The energy, the specific heat, and the magnetic susceptibility of the bond-diluted S=(1/2) Heisenberg ferromagnet on a two-dimensional square lattice are studied using a Monte Carlo sampling of the terms contributing to the partition function. The full lattice is studied in detail for sizes up to 20\ifmmode\times\else\texttimes\fi{}20, and evidence is presented to suggest that the low-temperature correlation length behaves as \ensuremath{\xi}\ensuremath{\sim}${e}^{\ensuremath{\nu}\ensuremath{\beta}J}$ with \ensuremath{\nu}\ensuremath{\sim}1.91\ifmmode\pm\else\textpm\fi{}0.06. The temperature dependence of the susceptibility, however, is not quite a simple exponential in \ensuremath{\beta}J. Knowledge of the correlation length and finite-size scaling allows the determination of the prefactor. Our conclusion is that the low-temperature susceptibility behaves as ${T}^{\mathrm{\ensuremath{-}}1}$${e}^{2\ensuremath{\nu}\ensuremath{\beta}J}$. Square lattices of size ranging from 6\ifmmode\times\else\texttimes\fi{}6 to 12\ifmmode\times\else\texttimes\fi{}12 are studied for
kBT/J ranging from 0.1 to 3 and for several values of the bond probability p. Not surprisingly, the most drastic thermodynamic behavior near the critical bond-percolation probability is the sharp increase in the low-temperature susceptibility for p>
pc. This corresponds to much larger correlation effects in the percolating lattice than below the percolation threshold. We also compare the thermodynamics at p=
pc to the scaling theories of Lubensky and Stanley et al. which suggest that the correlations at the percolation threshold are predominantly one dimensional.
The two-dimensional random bond Ising model with two distinct ferromagnetic nearest-neighbour couplings, J and J ', distributed randomly over the lattice with concentration 1/2 is studied for various ratios J '/ J at criticality. The critical temperature is exactly known from self-duality. The size-dependent specific heat shows a crossover from the perfect Ising to the randomness-dominated critical behaviour of doubly logarithmic form. In the latter region, our estimates for the critical exponents of the size-dependent magnetization and susceptibility are in very good agreement with analytic results suggesting β/ν = 1/8 and γ/ν = 7/4.
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.
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. This provides a possible explanation for the apparent specific-heat rounding in certain experiments.