Using the Monte Carlo technique the authors have studied a two-dimensional Ising model with quenched non-magnetic impurities. In particular they calculated the temperature dependence of the order parameter, of the zero-field susceptibility and of the specific heat. These results indicate that a second-order phase transition is preserved for p>or=0.9 where 1-p denotes the concentration of randomly distributed non-magnetic sites. Moreover, they yield beta approximately=1/8 and gamma approximately=7/4, indicating that these exponents equal the corresponding ones for the pure Ising model.
Using the Monte Carlo technique the authors have studied a two-dimensional Ising model with quenched non-magnetic impurities. In particular they calculated the temperature dependence of the order parameter, of the zero-field susceptibility and of the specific heat. These results indicate that a second-order phase transition is preserved for p>or=0.9 where 1-p denotes the concentration of randomly distributed non-magnetic sites. Moreover, they yield beta approximately=1/8 and gamma approximately=7/4, indicating that these exponents equal the corresponding ones for the pure Ising model.
A review is given of Monte Carlo simulations for the phase transition in simple cubic Ising ferro-, ferri-, and antiferromagnets with quenched non-magnetic impurities of concentration x. Special emphasis is given on the description of a fully vectorized multispin coding program and on the problems related to the analysis of the Monte Carlo data. The following three subjects are discussed: 1) For ferro- and antiferromagnets, the critical behaviour of the order parameter is different from the one of the corresponding pure system, in agreement with the prediction of the Harris criterion. The possibly effective exponent β depends on x, whereas for ferrimagnets no such concentration dependence is revealed by the present study for x=0, 0.1, and 0.2.2) For strongly diluted systems, there is a maximum for the Kouvel-Fisher exponent γ(T) of the order-parameter susceptibility outside the critical regime, in accordance with the predictions of the correlated molecular field theory. 3) For the diluted antiferromagnet (x=0.5), a cusp emerges at T c for the uniform susceptibility, in agreement with the theory of Fishman and Aharony.
Monte Carlo simulations are performed to analyze phase transitions in three-dimensional systems described by the 3-state Potts model with nonmagnetic impurities. Numerical results are presented for systems with spin concentrations p = 1.00, 0.95, 0.90, 0.80, 0.70, and 0.65 on lattices of size L varying between 20 and 44. Binder’s cumulant analysis shows that the introduction of quenched disorder in the form of non-magnetic impurities induces a crossover from first-order to second-order phase transition. The finite-size scaling method is used to calculate the static critical exponents α, γ, β, and ν for specific heat, susceptibility, magnetization, and correlation length, respectively.
The thermal and magnetic properties of the simple cubic Ising model with random, quenched, site impurities have been determined by Monte Carlo calculations. The transition appears to remain sharp, and the critical temperature decreases monotonically towards zero as the impurity concentration is increased towards the percolation limit. No clear evidence is found for any change in the critical exponents away from the pure lattice values. These results are discussed in terms of Harris's predictions for the critical behavior of impure systems.
The influence of quenched nonmagnetic impurities on phase transitions and critical phenomena in the 3D Potts model with the spin state number q = 4 is studied using the Monte Carlo method. Systems with the linear size L = 20–32 and spin concentrations p = 1.00, 0.90, 0.65 are considered. The fourth order Binder cumulant method is used to demonstrate that in the strongly diluted regime, a phase transition of the second kind is observed in this model for the spin concentration p = 0.65, and a phase transition of the first kind is observed for the pure ( p = 1.00) and weakly diluted ( p = 0.90) models. The theory of finite-dimensional scaling is used to calculate the static critical parameters of heat capacity α, susceptibility γ, magnetization β, and correlation radius ν.
The effect of quenched nonmagnetic impurities on phase transitions in the three-dimensional Potts model with the number of spin states q=3 is studied using the Wolff single-cluster algorithm of the Monte Carlo method. By the method of fourth-order Binder cumulants, it is demonstrated that the second-order phase transition occurs in the model under study at spin concentrations p=0.9, 0.8, 0.7, and 0.65, while the first-order phase transition is observed in the pure model (p=1.0). The static critical exponents (CEs) α (heat capacity), γ (susceptibility), β (magnetization), and ν (correlation length) are calculated based on the finite-size scaling theory.
The influence of quenched nonmagnetic impurities on phase transitions in the three-dimensional Potts model with the number of spin states q = 3 is investigated using the Wolff single-cluster algorithm of the Monte Carlo method. The systems with linear sizes L = 20–44 at the spin concentrations p = 1.0, 0.9, 0.8, and 0.7 are analyzed. It is demonstrated with the use of the method of fourth-order Binder cumulants that the second-order phase transition occurs in the model under consideration at the spin concentrations p = 0.9, 0.8, and 0.7 and that the first-order phase transition is observed in the pure model ( p = 1.0). The static critical exponents α (heat capacity), γ (susceptibility), β (magnetization), and ν (correlation length) are calculated in the framework of the finite-size scaling theory. The problem regarding the universality classes of the critical behavior of weakly diluted systems is discussed.