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:2406.18782

Numerical simulations for the Ising model on three dimensional lattices with coordination number equal 5: static and dynamic critical phenomena

Lourdes Bibiana Merino-Solís, Francisco Sastre

Published Jun 26, 2024• 0 citations
arXivHTMLPDF

Abstract

In this work we performed numerical simulations for the Ising model on three dimensional lattices with coordination number equal 5. With Monte Carlo simulations in the static case we evaluated the critical temperature and the static critical exponents $ν$, $γ$ and $β$. Once that we have the critical temperature value we investigated the dynamical critical behavior with Glauber dynamics, starting from disordered states. From our simulations we obtained the values for the dynamic exponent $z=2.037(8)$, the exponent for the autocorrelation $λ/z=1.364(5)$, the exponent of the critical initial increase $θ'=0.109(8)$ and the asymptotic value of the fluctuation-dissipation ratio $X^\infty=0.433(1)$. All of these results are in good agreement to the previous values reported for the 3D Ising model universality class.

DOI
10.48550/arxiv.2406.18782

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 static and dynamic simulations of a three-dimensional Ising critical model

F. Livet

The critical dynamics of `model A'' of Hohenberg and Halperin has been studied by the Monte Carlo method. Simulations have been carried out in the three-dimensional (3d) simple cubic Ising model for lattices of sizes $L=16$ to $L=512$. Using Wolff's cluster algorithm, the critical temperature is precisely found as ${\ensuremath{\beta}}_{c}=0.221\phantom{\rule{0.16em}{0ex}}654\phantom{\rule{0.16em}{0ex}}68(5)$. By Fourier transform of the lattice configurations, the critical scattering intensities $I(\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}})$ can be obtained. After circular averaging, the static simulations with $L=256$ and $L=512$ provide an estimate of the critical exponent $\ensuremath{\gamma}/\ensuremath{\nu}=2\ensuremath{-}\ensuremath{\eta}=1.9640(5)$. The $|\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}}|$-dependent distribution of $I(\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}})$ showed an exponential distribution, corresponding to a Gaussian distribution of the scattering amplitudes for a large $q$ domain. The time-dependent intensities were then used for the study of the critical dynamics of 3d lattices at the critical point. To simulate results of an x-ray photon correlation spectroscopy experiment, the time-dependent correlation function of the intensities was studied for each $|\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}}|$-value. In the $q$ region where $I(\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}})$ had an exponential distribution, the time correlations can be fit to a stretched exponential, where the exponent $\ensuremath{\mu}=\ensuremath{\gamma}/\ensuremath{\nu}z\ensuremath{\simeq}0.975$ provides an estimate of the dynamic exponent $z$. This corresponds to $z=2.0145$, in agreement with the observed variations of the characteristic fluctuation time of the intensity: $\ensuremath{\tau}(q)\ensuremath{\propto}{q}^{\ensuremath{-}z}$, which gives $z=2.015$. These results agree with the $\ensuremath{\epsilon}$ expansion of field-theoretical methods (2.017). In this paper, the need to take account of the anomalous time behavior ($\ensuremath{\mu}<1$) in the dynamics is exemplified. This dynamics reflects a nonlinear time behavior of model A, and its large time extension is discussed in detail.

Physical review. E•Feb 27, 2020•1 citation

Short-time critical dynamics of the three-dimensional Ising model

A Jaster, J. Mainville, L Schülke +1 more

Comprehensive Monte Carlo simulations of the short-time dynamic behaviour are reported for the three-dimensional Ising model at criticality. Besides the exponent of the critical initial increase and the dynamic exponent z , the static critical exponents and as well as the critical temperature are determined from the power-law scaling behaviour of observables at the beginning of the time evolution. States of very high temperature as well as of zero temperature are used as initial states for the simulations.

arXiv:cond-mat/9808131•Journal of Physics A Mathematical and General•Jan 1, 1999•76 citations

Dynamic critical exponent z of the three-dimensional Ising universality class: Monte Carlo simulations of the improved Blume-Capel model

Martin Hasenbusch B

We study purely dissipative relaxational dynamics in the three-dimensional Ising universality class. To this end, we simulate the improved Blume-Capel model on the simple cubic lattice by using local algorithms. We perform a finite size scaling analysis of the integrated autocorrelation time of the magnetic susceptibility in equilibrium at the critical point. We obtain z=2.0245(15) for the dynamic critical exponent. As a complement, fully magnetized configurations are suddenly quenched to the critical temperature, giving consistent results for the dynamic critical exponent. Furthermore, our estimate of z is fully consistent with recent field theoretic results.

arXiv:1908.01702•Physical review. E•Feb 24, 2020•39 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

Universal dynamic scaling in three-dimensional Ising spin glasses

Cheng-Wei Liu, Anatoli Polkovnikov, Anders W. Sandvik +1 more

We use a nonequilibrium Monte Carlo simulation method and dynamical scaling to study the phase transition in three-dimensional Ising spin glasses. The transition point is repeatedly approached at finite velocity v (temperature change versus time) in Monte Carlo simulations starting at a high temperature. This approach has the advantage that the equilibrium limit does not have to be strictly reached for a scaling analysis to yield critical exponents. For the dynamic exponent we obtain z=5.85(9) for bimodal couplings distribution and z=6.00(10) for the Gaussian case. Assuming universal dynamic scaling, we combine the two results and obtain z=5.93±0.07 for generic 3D Ising spin glasses.

arXiv:1411.6745•Physical Review E•Aug 19, 2015•25 citations
In this work we performed numerical simulations for the Ising model on three dimensional lattices with coordination number equal 5. With Monte Carlo simulations in the static case we evaluated the critical temperature and the static critical exponents ννν, γγγ and βββ. Once that we have the critical temperature value we investigated the dynamical critical behavior with Glauber dynamics, starting from disordered states. From our simulations we obtained the values for the dynamic exponent z=2.037(8)z=2.037(8)z=2.037(8), the exponent for the autocorrelation λ/z=1.364(5)λ/z=1.364(5)λ/z=1.364(5), the exponent of the critical initial increase θ′=0.109(8)θ'=0.109(8)θ′=0.109(8) and the asymptotic value of the fluctuation-dissipation ratio X∞=0.433(1)X^\infty=0.433(1)X∞=0.433(1). All of these results are in good agreement to the previous values reported for the 3D Ising model universality class.
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.
The critical dynamics of `model A'' of Hohenberg and Halperin has been studied by the Monte Carlo method. Simulations have been carried out in the three-dimensional (3d) simple cubic Ising model for lattices of sizes L=16L=16L=16 to L=512L=512L=512. Using Wolff's cluster algorithm, the critical temperature is precisely found as ${\ensuremath{\beta}}_{c}=0.221\phantom{\rule{0.16em}{0ex}}654\phantom{\rule{0.16em}{0ex}}68(5)$. By Fourier transform of the lattice configurations, the critical scattering intensities $I(\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}})$ can be obtained. After circular averaging, the static simulations with L=256L=256L=256 and L=512L=512L=512 provide an estimate of the critical exponent $\ensuremath{\gamma}/\ensuremath{\nu}=2\ensuremath{-}\ensuremath{\eta}=1.9640(5)$. The $|\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}}|$-dependent distribution of $I(\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}})$ showed an exponential distribution, corresponding to a Gaussian distribution of the scattering amplitudes for a large qqq domain. The time-dependent intensities were then used for the study of the critical dynamics of 3d lattices at the critical point. To simulate results of an x-ray photon correlation spectroscopy experiment, the time-dependent correlation function of the intensities was studied for each $|\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}}|$-value. In the qqq region where $I(\stackrel{P\vec}{q}\phantom{\rule{0.16em}{0ex}})$ had an exponential distribution, the time correlations can be fit to a stretched exponential, where the exponent $\ensuremath{\mu}=\ensuremath{\gamma}/\ensuremath{\nu}z\ensuremath{\simeq}0.975$ provides an estimate of the dynamic exponent zzz. This corresponds to z=2.0145z=2.0145z=2.0145, in agreement with the observed variations of the characteristic fluctuation time of the intensity: $\ensuremath{\tau}(q)\ensuremath{\propto}{q}^{\ensuremath{-}z}$, which gives z=2.015z=2.015z=2.015. These results agree with the $\ensuremath{\epsilon}$ expansion of field-theoretical methods (2.017). In this paper, the need to take account of the anomalous time behavior ($\ensuremath{\mu}<1$) in the dynamics is exemplified. This dynamics reflects a nonlinear time behavior of model A, and its large time extension is discussed in detail.
Comprehensive Monte Carlo simulations of the short-time dynamic behaviour are reported for the three-dimensional Ising model at criticality. Besides the exponent of the critical initial increase and the dynamic exponent z , the static critical exponents and as well as the critical temperature are determined from the power-law scaling behaviour of observables at the beginning of the time evolution. States of very high temperature as well as of zero temperature are used as initial states for the simulations.
We study purely dissipative relaxational dynamics in the three-dimensional Ising universality class. To this end, we simulate the improved Blume-Capel model on the simple cubic lattice by using local algorithms. We perform a finite size scaling analysis of the integrated autocorrelation time of the magnetic susceptibility in equilibrium at the critical point. We obtain z=2.0245(15) for the dynamic critical exponent. As a complement, fully magnetized configurations are suddenly quenched to the critical temperature, giving consistent results for the dynamic critical exponent. Furthermore, our estimate of z is fully consistent with recent field theoretic results.
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 use a nonequilibrium Monte Carlo simulation method and dynamical scaling to study the phase transition in three-dimensional Ising spin glasses. The transition point is repeatedly approached at finite velocity v (temperature change versus time) in Monte Carlo simulations starting at a high temperature. This approach has the advantage that the equilibrium limit does not have to be strictly reached for a scaling analysis to yield critical exponents. For the dynamic exponent we obtain z=5.85(9) for bimodal couplings distribution and z=6.00(10) for the Gaussian case. Assuming universal dynamic scaling, we combine the two results and obtain z=5.93±0.07 for generic 3D Ising spin glasses.