Monte Carlo simulations of the short-time dynamic behavior are reported for three-dimensional weakly site-diluted Ising model with spin concentrations p=0.95 and 0.8 at criticality. In contrast to studies of the critical behavior of the pure systems by the short-time dynamics method, our investigations of site-diluted Ising model have revealed three stages of the dynamic evolution characterizing a crossover phenomenon from the critical behavior typical for the pure systems to behavior determined by the influence of disorder. The static and dynamic critical exponents are determined with the use of the corrections to scaling for systems starting separately from ordered and disordered initial states. The obtained values of the exponents demonstrate a universal behavior of weakly site-diluted Ising model in the critical region. The values of the exponents are compared to results of numerical simulations which have been obtained in various works and, also, with results of the renormalization-group description of this model.
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).
The time evolution of the three-dimensional critical Ising model relaxing from a non-equilibrium initial state is studied by means of Monte Carlo simulation. We observe the characteristic initial increase of the (spatially) averaged magnetization predicted by Janssen et al. The exponent theta ' that governs the initial behaviour is determined, and the dependence of the long-time linear decay on the initial magnetization, m 0 , is analysed. Our simulation corroborates earlier results derived from continuum models.
Exact analytic results are presented for the critical properties of the dilute Ising model on a Bethe lattice, where a fraction p of the sites are occupied at random by spins. It is found that the critical behaviour is the same as that of the pure system (p=1) except for transitions which occur by varying p at zero temperature, where the critical behaviour is different and corresponds to the percolation problem. For T to 0 and p to p c , the percolation concentration, the zero-field magnetization shows crossover behaviour of the form expected from scaling theories. Outside the critical region, a rapidly convergent numerical method is developed and used to calculate the spontaneous magnetization.
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 We have performed a numerical investigation of the influence of disorder on the dynamical non-equilibrium evolution of a 3D site-diluted Ising model from a low-temperature initial state with magnetization m 0 = 1. It is shown that two-time dependences of the autocorrelation and integrated response functions for systems with spin concentrations p = 1.0, 0.95, 0.8, 0.6 and 0.5 demonstrate ageing properties with anomalous slowing-down relaxation and violation of the fluctuation-dissipation ratio. It was revealed that during non-equilibrium critical dynamics in the long-time regime t − t w ≫ t w ≫ 1 the autocorrelation functions for diluted systems are extremely slow due to the pinning of domain walls on impurity sites. We have found that the fluctuation-dissipation ratio X ∞ = 0 for diluted systems with spin concentration p < 1 while the pure system is characterized by X ∞ = 0.784 ( 7 ) . The autocorrelation function power-law delay becomes the same as for the time dependence of the magnetization in the critical point and is characterized by exponent − β / z ν . Also, for diluted systems we reveal memory effects for critical evolution in the ageing regime with realization of cyclic temperature change and quenching at T < T c .