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.
Randomness is an important subject in the study of phase transition as defect and impurity may present in any real material. The pre-existing ordered phase of a pure system can be affected or even ruined by the presence of randomness. Here we study ferromagnetic Ising model on a square lattice with a presence of randomness in the form of bond dilution. The pure system of this model is known to experience second order phase transition, separating between the high temperature paramagnetic and low-temperature ferromagnetic phase. We used Wang-Landau algorithm of Monte Carlo method to obtain the density of states from which we extract the ensemble average of energy and the specific heat. We observed the signature of phase transition indicated by the diverging peak of the specific heat as system sizes increase. These peaks shift to the lower temperature side as the dilution increases. The lower temperature ordered phase preserves up to certain concentration of dilution and is totally ruined when the bonds no longer percolates.