Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:0709.1450 · Journal of Experimental and Theoretical Physics

Computer simulation of the critical behavior of 3D disordered ising model

V. V. Prudnikov, Pavel Vladimirovich Prudnikov, A. N. Vakilov, A. S. Krinitsyn

Published Aug 1, 2007• 35 citations
arXivHTMLPDF

Abstract

The critical behavior of the disordered ferromagnetic Ising model is studied numerically by the Monte Carlo method in a wide range of variation of concentration of nonmagnetic impurity atoms. The temperature dependences of correlation length and magnetic susceptibility are determined for samples with various spin concentrations and various linear sizes. The finite-size scaling technique is used for obtaining scaling functions for these quantities, which exhibit a universal behavior in the critical region; the critical temperatures and static critical exponents are also determined using scaling corrections. On the basis of variation of the scaling functions and values of critical exponents upon a change in the concentration, the conclusion is drawn concerning the existence of two universal classes of the critical behavior of the diluted Ising model with different characteristics for weakly and strongly disordered systems.

DOI
10.1134/s1063776107080092

Related papers

Critical behavior of a cubic-lattice 3D Ising model for systems with quenched disorder

Akai Kurbanovich Murtazaev, Ibragimkhan K. Kamilov, Albert Babaevich Babaev

A Monte Carlo method is applied to simulate the static critical behavior of a cubic-lattice 3D Ising model for systems with quenched disorder. Numerical results are presented for the spin concentrations of p = 1.0, 0.95, 0.9, 0.8, 0.6 on L × L × L lattices with L = 20–60 under periodic boundary conditions. The critical temperature is determined by the Binder cumulant method. A finite-size scaling technique is used to calculate the static critical exponents α, β, γ, and ν (for specific heat, susceptibility, magnetization, and correlation length, respectively) in the range of p under study. Universality classes of critical behavior are discussed for three-dimensional diluted systems.

Journal of Experimental and Theoretical Physics•Dec 1, 2004•35 citations

Monte Carlo simulation of strongly disordered Ising ferromagnets

Hans-Otto Heuer

We report extensive Monte Carlo simulations of disordered Ising systems in the ferromagnetic region with concentrations of magnetic sites between p=1.0 and 0.5. The magnetization, the susceptibility, and the caloric properties have been studied in the critical region. The critical exponents \ensuremath{\beta} and \ensuremath{\gamma} as well as the universal amplitude A of the magnetization and the ratio ${\mathit{C}}_{+}$/${\mathit{C}}_{\mathrm{\ensuremath{-}}}$ of the susceptibility amplitudes have been determined with high precision. In addition to the cusplike specific heat, we have also measured the magnetization-energy correlation function \ensuremath{\Gamma}, which is a divergent thermal quantity. The corresponding critical exponent \ensuremath{\zeta}, which is related to the other exponents by \ensuremath{\zeta}=1-\ensuremath{\beta}=(\ensuremath{\alpha}+\ensuremath{\gamma})/2, has been determined. We have found that all quantities show power-law behavior within the temperature range of our simulation. All critical exponents change continuously with dilution. Even in the range of weak dilution (p\ensuremath{\ge}0.8), the effective critical exponents are concentration dependent and are clearly different from their pure system values. In the strongly diluted regime the critical exponents gradually reach new asymptotic behavior at p=0.5--0.6 with values of \ensuremath{\beta}=0.335\ifmmode\pm\else\textpm\fi{}0.01 and \ensuremath{\gamma}=1.49\ifmmode\pm\else\textpm\fi{}0.02. The exponent \ensuremath{\alpha} of the specific heat becomes -0.17\ifmmode\pm\else\textpm\fi{}0.04, which corresponds to a cusplike singularity. We conclude that disorder profoundly changes the critical behavior for weakly as well as strongly disordered spin systems.

Physical review. B, Condensed matter•Oct 1, 1990•55 citations

Short-time dynamics and critical behavior of the three-dimensional site-diluted Ising model

Vladimir V. Prudnikov, Pavel Vladimirovich Prudnikov, Aleksandr Krinitsyn +3 more

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.

arXiv:1005.5218•Physical Review E•Jan 26, 2010•44 citations

Zeros of Partition Function and Critical Exponents of the 3D Diluted Ising Model

Andrey N. Vakilov

We used a Monte Carlo simulation of the structurally disordered three dimensional Ising model. For the systems with spin concentrations p = 0.95 ,0.8, 0.6 and 0.5 we calculated the correlation-length critical exponent ν by finite-size scaling. Extrapolations to the thermodynamic limit yield ν(0.95) = 0.705(5) ,ν(0.8) = 0.711(6),ν(0.6) = 0.736(6) and ν(0.5) = 0.744(6). These results are compatible with some previous estimates from a variety of sources. The analysis of the results demonstrates the nonuniversality of the critical behavior in the disordered Ising model.

Materials science forum•Mar 1, 2016•0 citations

Weak Universality in the Disordered Two-Dimensional Antiferromagnetic Potts Model on a Triangular Lattice

Albert Babaevich Babaev, Akai Kurbanovich Murtazaev

The critical behavior of the disordered two-dimensional antiferromagnetic Potts model with the number of spin states q= 3 on a triangular lattice with disorder in the form of nonmagnetic impurities is studied by the Monte Carlo method. The critical exponents for the susceptibility γ, magnetization β, specific heat α, and correlation radius ν are calculated in the framework of the finite-size scaling theory at spin concentrations p = 0.90, 0.80, 0.70, and 0.65. It is found that the critical exponents increase with the degree of disorder, whereas the ratios and do not change, thus holding the scaling equality $$\frac{{2\beta }}{\nu } + \frac{\gamma }{\nu } = d$$ . Such behavior of the critical exponents is related to the weak universality of the critical behavior characteristic of disordered systems. All results are obtained using independent Monte Carlo algorithms, such as the Metropolis and Wolff algorithms.

Journal of Experimental and Theoretical Physics Letters•May 1, 2018•5 citations

Monte Carlo Simulation of Strongly Disordered Ising Ferromagnets

Hans-O Heuer

We report about extensive Monte Carlo simulations of disordered Ising systems with concentrations of magnetic sites between p = 1.0 and p = 0.5. We have measured the magnetization, the susceptibility and the caloric properties in the critical region. In addition the magnetization-energy correlation function has been investigated and the critical exponent ζ = 1 - β = (α + γ)/2 has been determined. We find that power law behaviour and scaling is fulfilled in the whole concentration range. The effective critical exponents γ, β and ζ are shown to be concentration dependent. In the range of strong disorder γ and β increase until they reach a kind of asymptotic behaviour at β = 0.335 ± 0.01 and γ = 1.49 ± 0.02 for p = 0.5. In accord with the exponent relation ζ = 1 - β the critical exponent ζ decreases with dilution to ζ = 0.675 ± 0.01.

Europhysics Letters (EPL)•Jul 15, 1990•27 citations
The critical behavior of the disordered ferromagnetic Ising model is studied numerically by the Monte Carlo method in a wide range of variation of concentration of nonmagnetic impurity atoms. The temperature dependences of correlation length and magnetic susceptibility are determined for samples with various spin concentrations and various linear sizes. The finite-size scaling technique is used for obtaining scaling functions for these quantities, which exhibit a universal behavior in the critical region; the critical temperatures and static critical exponents are also determined using scaling corrections. On the basis of variation of the scaling functions and values of critical exponents upon a change in the concentration, the conclusion is drawn concerning the existence of two universal classes of the critical behavior of the diluted Ising model with different characteristics for weakly and strongly disordered systems.
A Monte Carlo method is applied to simulate the static critical behavior of a cubic-lattice 3D Ising model for systems with quenched disorder. Numerical results are presented for the spin concentrations of p = 1.0, 0.95, 0.9, 0.8, 0.6 on L × L × L lattices with L = 20–60 under periodic boundary conditions. The critical temperature is determined by the Binder cumulant method. A finite-size scaling technique is used to calculate the static critical exponents α, β, γ, and ν (for specific heat, susceptibility, magnetization, and correlation length, respectively) in the range of p under study. Universality classes of critical behavior are discussed for three-dimensional diluted systems.
We report extensive Monte Carlo simulations of disordered Ising systems in the ferromagnetic region with concentrations of magnetic sites between p=1.0 and 0.5. The magnetization, the susceptibility, and the caloric properties have been studied in the critical region. The critical exponents \ensuremath{\beta} and \ensuremath{\gamma} as well as the universal amplitude A of the magnetization and the ratio C+{\mathit{C}}_{+}C+​/${\mathit{C}}_{\mathrm{\ensuremath{-}}}$ of the susceptibility amplitudes have been determined with high precision. In addition to the cusplike specific heat, we have also measured the magnetization-energy correlation function \ensuremath{\Gamma}, which is a divergent thermal quantity. The corresponding critical exponent \ensuremath{\zeta}, which is related to the other exponents by \ensuremath{\zeta}=1-\ensuremath{\beta}=(\ensuremath{\alpha}+\ensuremath{\gamma})/2, has been determined. We have found that all quantities show power-law behavior within the temperature range of our simulation. All critical exponents change continuously with dilution. Even in the range of weak dilution (p\ensuremath{\ge}0.8), the effective critical exponents are concentration dependent and are clearly different from their pure system values. In the strongly diluted regime the critical exponents gradually reach new asymptotic behavior at p=0.5--0.6 with values of \ensuremath{\beta}=0.335\ifmmode\pm\else\textpm\fi{}0.01 and \ensuremath{\gamma}=1.49\ifmmode\pm\else\textpm\fi{}0.02. The exponent \ensuremath{\alpha} of the specific heat becomes -0.17\ifmmode\pm\else\textpm\fi{}0.04, which corresponds to a cusplike singularity. We conclude that disorder profoundly changes the critical behavior for weakly as well as strongly disordered spin systems.
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 used a Monte Carlo simulation of the structurally disordered three dimensional Ising model. For the systems with spin concentrations p = 0.95 ,0.8, 0.6 and 0.5 we calculated the correlation-length critical exponent ν by finite-size scaling. Extrapolations to the thermodynamic limit yield ν(0.95) = 0.705(5) ,ν(0.8) = 0.711(6),ν(0.6) = 0.736(6) and ν(0.5) = 0.744(6). These results are compatible with some previous estimates from a variety of sources. The analysis of the results demonstrates the nonuniversality of the critical behavior in the disordered Ising model.
The critical behavior of the disordered two-dimensional antiferromagnetic Potts model with the number of spin states q= 3 on a triangular lattice with disorder in the form of nonmagnetic impurities is studied by the Monte Carlo method. The critical exponents for the susceptibility γ, magnetization β, specific heat α, and correlation radius ν are calculated in the framework of the finite-size scaling theory at spin concentrations p = 0.90, 0.80, 0.70, and 0.65. It is found that the critical exponents increase with the degree of disorder, whereas the ratios and do not change, thus holding the scaling equality 2βν+γν=d\frac{{2\beta }}{\nu } + \frac{\gamma }{\nu } = dν2β​+νγ​=d . Such behavior of the critical exponents is related to the weak universality of the critical behavior characteristic of disordered systems. All results are obtained using independent Monte Carlo algorithms, such as the Metropolis and Wolff algorithms.
We report about extensive Monte Carlo simulations of disordered Ising systems with concentrations of magnetic sites between p = 1.0 and p = 0.5. We have measured the magnetization, the susceptibility and the caloric properties in the critical region. In addition the magnetization-energy correlation function has been investigated and the critical exponent ζ = 1 - β = (α + γ)/2 has been determined. We find that power law behaviour and scaling is fulfilled in the whole concentration range. The effective critical exponents γ, β and ζ are shown to be concentration dependent. In the range of strong disorder γ and β increase until they reach a kind of asymptotic behaviour at β = 0.335 ± 0.01 and γ = 1.49 ± 0.02 for p = 0.5. In accord with the exponent relation ζ = 1 - β the critical exponent ζ decreases with dilution to ζ = 0.675 ± 0.01.