We study the purely relaxational dynamics (model A) at criticality in three-dimensional disordered Ising systems whose static critical behaviour belongs to the randomly diluted Ising universality class. We consider the site-diluted and bond-diluted Ising models, and the ± J Ising model along the paramagnetic–ferromagnetic transition line. We perform Monte Carlo simulations at the critical point using the Metropolis algorithm and study the dynamic behaviour in equilibrium at various values of the disorder parameter. The results provide a robust evidence of the existence of a unique model-A dynamic universality class which describes the relaxational critical dynamics in all considered models. In particular, the analysis of the size dependence of suitably defined autocorrelation times at the critical point provides the estimate z = 2.35(2) for the universal dynamic critical exponent. We also study the off-equilibrium relaxational dynamics following a quench from to T c . In agreement with the field-theory scenario, the analysis of the off-equilibrium dynamic critical behaviour gives an estimate of z that is perfectly consistent with the equilibrium estimate z = 2.35(2).
We apply numerical simulations to study the criticality of the 3D Ising model with random site quenched dilution. The emphasis is given to the issues not being discussed in detail before. In particular, we attempt a comparison of different Monte Carlo techniques, discussing regions of their applicability and advantages/disadvantages depending on the aim of a particular simulation set. Moreover, besides evaluation of the critical indices we estimate the universal ratio
Γ+/Γ− for the magnetic susceptibility critical amplitudes. Our estimate
Γ+/Γ−=1.67±0.15 is in a good agreement with the recent MC analysis of the random-bond Ising model giving further support that both random-site and random-bond dilutions lead to the same universality class.
We study a short-range site-random Ising model in three dimensions (3D) using a Monte Carlo (MC) method and find a phase transition to a spin glass (SG) phase at a very high temperature ( T c ≃2.83 J). Values of the critical exponents estimated from a finite-size scaling analysis are different from those in both the 3D
∓J Ising model and the 3D Ising ferromagnet. We suggest that the SG phase in this model is a new SG phase which is characterized by a nonzero SG order parameter q ≠0.