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 report results of large Monte Carlo simulations of an experimentally realizable random-field system: the three-dimensional dilute Ising antiferromagnetic in a field. We find that the correlation time diverges dramatically as $T\ensuremath{\rightarrow}{T}_{c}$; the results are consistent with a proposed new type of activated dynamic scaling. The transition appears to be continuous, with effective critical exponents $\ensuremath{\eta}=0.5\ifmmode\pm\else\textpm\fi{}0.1$ and $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\eta}}=\ensuremath{-}1.0\ifmmode\pm\else\textpm\fi{}0.3$ for strong fields, away from the weak-field regime where crossover effects distort exponent estimates. These strong-field exponents satisfy recently derived inequalities.
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.
The phase transitions and critical phenomena in three-dimensional (3D) site-diluted 3-and 4-state Potts models is investigated by Monte-Carlo method based on the highly efficient Wolff algorithm. The systems with linear sizesL=20-44 at spin concentrationsp=1.00, 0.95, 0.90, 0.80, 0.70, 0.65 are explored. The second-order phase transition is shown to occur in the three-dimensional 3-state Potts model with nonmagnetic impurities. In the 4-state Potts model there are observed first-order phase transitions in weakly diluted state, when the model is strongly diluted the first-order phase transitions change to the second-order one. On the basis of the finite-size scaling theory static critical exponents of specific heatα, susceptibilityγ, magnetizationβ, and exponent of correlation radiusνfor the systems under study are calculated.
Field-induced phase transitions in the bcc Ising antiferromagnet with nearest-neighbor interactions have been studied using a Monte Carlo method. The properties of $N\ifmmode\times\else\texttimes\fi{}N\ifmmode\times\else\texttimes\fi{}N$ lattices with periodic boundary conditions have been determined for $N\ensuremath{\le}16$ and the size effects analyzed using finite-size scaling theory. We find no field dependence of the critical exponents and only a slow variation of critical amplitudes in small fields in agreement with the "smoothness postulate." At low temperatures the critical-field curve shows a weak, but clear maximum which is about 2% greater than the zero-temperature value
Hc. Crossover to noncritical exponential behavior due to thermally induced spin flips is observed as we move away from the low-temperature phase boundary.