We evaluate by Monte Carlo simulations various singular thermodynamic quantities
X for ensembles of quenched random Ising and Ashkin-Teller models. The measurements are taken at
Tc and we study how the distributions
P(X) (and, in particular, their relative squared width,
RX) over the ensemble depend on the system size
l. The Ashkin-Teller model was studied in the regime where bond randomness is irrelevant and we found weak self-averaging; ${R}_{X}\ensuremath{\sim}{l}^{\ensuremath{\alpha}/\ensuremath{\nu}}\ensuremath{\rightarrow}0,$ where $\ensuremath{\alpha}<0$ and $\ensuremath{\nu}$ are the exponents (of the pure model fixed point) governing the transition. For the site-dilute Ising model on a cubic lattice, known to be governed by a random fixed point, we find that
RX tends to a constant, as predicted by Aharony and Harris. We tested whether this constant is universal. However, this constant is different for canonical and grand canonical disorder. We identify the pseudocritical temperature of each sample
i, Tc(i,l), as the temperature at which the susceptibility reaches its maximal value ${\ensuremath{\chi}}^{\mathrm{max}}.$ The distribution of these sample dependent
Tc(i,l) was investigated; we found that its variance scales as $[\ensuremath{\delta}{T}_{c}(l){]}^{2}\ensuremath{\sim}{l}^{\ensuremath{-}2/\ensuremath{\nu}}.$ Our previously proposed finite size scaling ansatz for disordered systems was tested and found to hold. We did observe deviations from a single function, which imply that sample dependent scaling functions are needed. These deviations are, however, relatively small and hence to obtain a fixed statistical error it may be more computationally efficient to measure ${\ensuremath{\chi}}^{\mathrm{max}}$ than the commonly used $\ensuremath{\chi}{(T}_{c}^{\ensuremath{\infty}}).$
We simulated site dilute Ising models in
d=3 dimensions for several lattice sizes
L. For each
L singular thermodynamic quantities
X were measured at criticality and their distributions
P(X) were determined for ensembles of several thousand random samples. For $L\ensuremath{\rightarrow}\ensuremath{\infty}$ the relative width of
P(X) tends to a universal constant: there is no self-averaging. The width of the distribution of the sample
(i) dependent pseudocritical temperatures
Tc(i,L) scales as $\ensuremath{\delta}{T}_{c}(L)\ensuremath{\sim}{L}^{\ensuremath{-}1/\ensuremath{\nu}}$ and not as $\mathrm{\ensuremath{\sim}}{L}^{\ensuremath{-}d/2}$. The sample dependence of
Xi(T,L) enters dominantly, but not exclusively, via
Tc(i,L).
The universality class of thermally diluted Ising systems, in which the realization of the disposition of magnetic atoms and vacancies is taken from the local distribution of spins in the pure original Ising model at criticality, is investigated by finite size scaling techniques using the Monte Carlo method. We find that the critical temperature, the critical exponents, and therefore the universality class of these thermally diluted Ising systems depart markedly from the ones of short-range correlated disordered systems. Our results agree fairly well with theoretical predictions previously made by Weinrib and Halperin for systems with long-range correlated disorder.
We simulated site dilute Ising models in
d=3 dimensions for several lattice sizes
L. For each
L singular thermodynamic quantities
X were measured at criticality and their distributions
P(X) were determined, for ensembles of several thousand random samples. For
L→∞ the width of
P(X) tends to a universal constant, i.e. there is no self averaging. The width of the distribution of the sample dependent pseudocritical temperatures
Tc(i,L) scales as
δTc(L)∼L−1/ν and NOT as
∼L−d/2. Finite size scaling holds; the sample dependence of
Xi(Tc) enters predominantly through
Tc(i,L).
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.
We study by Monte Carlo simulations the influence of bond dilution on the three-dimensional Ising model. This paradigmatic model in its pure version displays a second-order phase transition with a positive specific heat critical exponent
α . According to the Harris criterion disorder should hence lead to a new fixed point characterized by new critical exponents. We have determined the phase diagram of the diluted model, starting from the pure model limit down to the neighbourhood of the percolation threshold. For the estimation of critical exponents, we have first performed a finite-size scaling study, where we concentrated on three different dilutions to check the stability of the disorder fixed point. We emphasize in this work the great influence of the cross-over phenomena between the pure, disorder and percolation fixed points which lead to effective critical exponents dependent on the concentration. In a second set of simulations, the temperature behaviour of physical quantities has been studied in order to characterize the disorder fixed point more accurately. In particular this allowed us to estimate ratios of some critical amplitudes. In accord with previous observations for other models this provides stronger evidence for the existence of the disorder fixed point since the amplitude ratios are more sensitive to the universality class than the critical exponents. Moreover, the question of non-self-averaging at the disorder fixed point is investigated and compared with recent results for the bond-diluted q = 4 Potts model. Overall our numerical results provide evidence that, as expected on theoretical grounds, the critical behaviour of the bond-diluted model is indeed governed by the same universality class as the site-diluted model.