Damage spreading method in studies of the critical behavior of disordered systems

The computer simulations of the critical dynamics of the structurally disordered three-dimensional Ising model are performed by the damage spreading method. For the systems with spin densities $p = 0.6$ and $0.8$, we calculate the critical temperature and dynamic critical exponent $z$ characterizing the behavior of relaxation properties near the critical point. The analysis of the results demonstrates the nonuniversality of the critical dynamics in the disordered Ising model. To interpret such dynamics, we introduce the concept of two universal types of critical behavior corresponding to weak and strong disorder, respectively.