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:1005.1678 · Physical Review E

Short-time critical dynamics of damage spreading in the two-dimensional Ising model

M. Leticia Rubio Puzzo, Ezequiel Vicente Albano

Published May 12, 2010• 5 citations
arXivHTMLPDF

Abstract

The short-time critical dynamics of propagation of damage in the Ising ferromagnet in two dimensions is studied by means of Monte Carlo simulations. Starting with equilibrium configurations at $T=\ensuremath{\infty}$ and magnetization $M=0$, an initial damage is created by flipping a small amount of spins in one of the two replicas studied. In this way, the initial damage is proportional to the initial magnetization ${M}_{0}$ in one of the configurations upon quenching the system at ${T}_{C}$, the Onsager critical temperature of the ferromagnetic-paramagnetic transition. It is found that, at short times, the damage increases with an exponent ${\ensuremath{\theta}}_{D}=1.915(3)$, which is much larger than the exponent $\ensuremath{\theta}=0.197$ characteristic of the initial increase of the magnetization $M(t)$. Also, an epidemic study was performed. It is found that the average distance from the origin of the epidemic $(⟨{R}^{2}(t)⟩)$ grows with an exponent ${z}^{\ensuremath{\ast}}\ensuremath{\approx}\ensuremath{\eta}\ensuremath{\approx}1.9$, which is the same, within error bars, as the exponent ${\ensuremath{\theta}}_{D}$. However, the survival probability of the epidemics reaches a plateau so that $\ensuremath{\delta}=0$. On the other hand, by quenching the system to lower temperatures one observes the critical spreading of the damage at ${T}_{D}\ensuremath{\simeq}0.51{T}_{C}$, where all the measured observables exhibit power laws with exponents ${\ensuremath{\theta}}_{D}=1.026(3)$, $\ensuremath{\delta}=0.133(1)$, and ${z}^{\ensuremath{\ast}}=1.74(3)$.

DOI
10.1103/physreve.81.051116

Related papers

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

The Damage Spreading Method in Monte Carlo Simulations: A Brief Overview and Applications to Confined Magnetic Materials

M. Leticia Rubio Puzzo, Ezequiel Vicente Albano

The Damage Spreading (DS) method allows the investigation of the effect caused by tiny perturbations, in the initial conditions of physical systems, on their final stationary or equilibrium states. The damage (D(t)) is determined during the dynamic evolution of a physical system and measures the time dependence of the difference between a reference (unperturbed) configuration and an initially perturbed one. In this paper we first give a brief overview of Monte Carlo simulation results obtained by applying the DS method. Different model systems under study often exhibit a transition between a state where the damage becomes healed (the frozen phase) and a regime where the damage spreads arriving at a finite (stationary) value (the damaged phase), when a control parameter is finely tuned. These kinds of transitions are actually true irreversible phase transitions themselves, and the issue of their universality class is also discussed. Subsequently, the attention is focused on the propagation of damage in magnetic systems placed in confined geometries. The influence of interfaces between magnetic domains of different orientation on the spreading of the perturbation is also discussed, showing that the presence of interfaces enhances the propagation of the damage. Furthermore, the critical transition between propagation and nonpropagation of the damage is discussed. In all cases, the determined critical exponents suggest that the DS transition does not belong to the universality class of Directed Percolation, unlike many other systems exhibiting irreversible phase transitions. This result reflects the dramatic influence of interfaces on the propagation of perturbations in magnetic systems.

arXiv:0803.1602•Communications in Computational Physics•Jan 1, 2008•6 citations

On the dynamic exponent of the two-dimensional Ising model

K MacIsaac, Naeem Jan

The authors use the method of 'damage spreading' to measure the time taken for the equilibrium damage (the magnetization) to disappear in finite systems (10<or=L<or=103) at the critical temperature of the infinite system. This time is found to scale as L z with an exponent z of 2.16+or-0.02. This value is compatible with the expected value of the dynamic exponent for the 2D Ising model, but they argue that this method is unable to measure fluctuations properly which they show are responsible for the higher values recently reported in the literature. They suggest that the slowest mode of relaxation has an exponent of z>or=2.27.

Journal of Physics A Mathematical and General•Apr 21, 1992•24 citations

Damage spreading in the Ising model with a microcanonical constraint

F. G. Brady Moreira, Adauto J. F. de Souza, Ananias M. Mariz

We study damage spreading in the ferromagnetic Ising model, on the square lattice, by performing numerical simulations in the microcanonical ensemble. Although the damage spreads for all temperature values, a dynamical critical phenomenon is observed at ${\mathit{T}}_{\mathit{d}}$=1.020(2)${\mathit{T}}_{\mathit{C}}$ (${\mathit{T}}_{\mathit{C}}$ is the Curie temperature): For T\ensuremath{\ge}${\mathit{T}}_{\mathit{d}}$ the damage D(t) is insensitive to its initial value D(0), while for T${\mathit{T}}_{\mathit{d}}$ a weak dependence on D(0) is present in such a way that for each T there are two values of D(t) corresponding to the initial conditions D(0)1/2 and D(0)>1/2. We have also identified a convenient order parameter to characterize the transition occurring at ${\mathit{T}}_{\mathit{d}}$. \textcopyright{} 1996 The American Physical Society.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jan 1, 1996•11 citations

Persistence in Ferromagnetic Ordering: Dependence upon initial configuration

Saikat Chakraborty, Subir K. Das

We study the dynamics of ordering in ferromagnets via Monte Carlo simulations of the Ising model, employing the Glauber spin-flip mechanism, in space dimensions $d=2$ and $3$. Results for the persistence probability and the domain growth are discussed for quenches to various temperatures ($T_f$) below the critical one ($T_{c}$), from different initial temperatures $T_{i} \geq T_{c}$. In long time limit, for $T_{i} > T_{c}$, the persistence probability exhibits power-law decay with exponents $θ\simeq 0.22$ and $\simeq 0.18$ in $d=2$ and $3$, respectively. For finite $T_i$, the early time behavior is a different power-law whose life-time diverges and exponent decreases as $T_{i} \rightarrow T_{c}$. The crossover length between the two steps diverges as the equilibrium correlation length. $T_i=T_c$ is expected to provide a {\it{new universality class}} for which we obtain $θ\simeq 0.035$ in $d=2$ and $\simeq 0.10$ in $d=3$. The time dependence of the average domain size $\ell$, however, is observed to be rather insensitive to the choice of $T_i$.

arXiv:1501.05178•Jan 21, 2015•0 citations

Spin-1 Ising model: Exact damage-spreading relations and numerical simulations

Aécia S. D. dos Anjos, Ananias M. Mariz, Fernando Dantas Nobre +1 more

The nearest-neighbor-interaction spin-1 Ising model is investigated within the damage-spreading approach. Exact relations involving quantities computable through damage-spreading simulations and thermodynamic properties are derived for such a model, defined in terms of a very general Hamiltonian that covers several spin-1 models of interest in the literature. Such relations presuppose translational invariance and hold for any ergodic dynamical procedure, leading to an efficient tool for obtaining thermodynamic properties. The implementation of the method is illustrated through damage-spreading simulations for the ferromagnetic spin-1 Ising model on a square lattice. The two-spin correlation function and the magnetization are obtained, with precise estimates of their associated critical exponents and of the critical temperature of the model, in spite of the small lattice sizes considered. These results are in good agreement with the universality hypothesis, with critical exponents in the same universality class of the spin- 12 Ising model. The advantage of the present method is shown through a significant reduction of finite-size effects by comparing its results with those obtained from standard Monte Carlo simulations.

Physical Review E•Sep 3, 2008•1 citation
The short-time critical dynamics of propagation of damage in the Ising ferromagnet in two dimensions is studied by means of Monte Carlo simulations. Starting with equilibrium configurations at T=\ensuremath∞T=\ensuremath{\infty}T=\ensuremath∞ and magnetization M=0M=0M=0, an initial damage is created by flipping a small amount of spins in one of the two replicas studied. In this way, the initial damage is proportional to the initial magnetization M0{M}_{0}M0​ in one of the configurations upon quenching the system at TC{T}_{C}TC​, the Onsager critical temperature of the ferromagnetic-paramagnetic transition. It is found that, at short times, the damage increases with an exponent \ensuremathθD=1.915(3){\ensuremath{\theta}}_{D}=1.915(3)\ensuremathθD​=1.915(3), which is much larger than the exponent \ensuremathθ=0.197\ensuremath{\theta}=0.197\ensuremathθ=0.197 characteristic of the initial increase of the magnetization M(t)M(t)M(t). Also, an epidemic study was performed. It is found that the average distance from the origin of the epidemic (⟨R2(t)⟩)(⟨{R}^{2}(t)⟩)(⟨R2(t)⟩) grows with an exponent z\ensuremath∗\ensuremath≈\ensuremathη\ensuremath≈1.9{z}^{\ensuremath{\ast}}\ensuremath{\approx}\ensuremath{\eta}\ensuremath{\approx}1.9z\ensuremath∗\ensuremath≈\ensuremathη\ensuremath≈1.9, which is the same, within error bars, as the exponent \ensuremathθD{\ensuremath{\theta}}_{D}\ensuremathθD​. However, the survival probability of the epidemics reaches a plateau so that \ensuremathδ=0\ensuremath{\delta}=0\ensuremathδ=0. On the other hand, by quenching the system to lower temperatures one observes the critical spreading of the damage at TD\ensuremath≃0.51TC{T}_{D}\ensuremath{\simeq}0.51{T}_{C}TD​\ensuremath≃0.51TC​, where all the measured observables exhibit power laws with exponents \ensuremathθD=1.026(3){\ensuremath{\theta}}_{D}=1.026(3)\ensuremathθD​=1.026(3), \ensuremathδ=0.133(1)\ensuremath{\delta}=0.133(1)\ensuremathδ=0.133(1), and z\ensuremath∗=1.74(3){z}^{\ensuremath{\ast}}=1.74(3)z\ensuremath∗=1.74(3).
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.
The Damage Spreading (DS) method allows the investigation of the effect caused by tiny perturbations, in the initial conditions of physical systems, on their final stationary or equilibrium states. The damage (D(t)) is determined during the dynamic evolution of a physical system and measures the time dependence of the difference between a reference (unperturbed) configuration and an initially perturbed one. In this paper we first give a brief overview of Monte Carlo simulation results obtained by applying the DS method. Different model systems under study often exhibit a transition between a state where the damage becomes healed (the frozen phase) and a regime where the damage spreads arriving at a finite (stationary) value (the damaged phase), when a control parameter is finely tuned. These kinds of transitions are actually true irreversible phase transitions themselves, and the issue of their universality class is also discussed. Subsequently, the attention is focused on the propagation of damage in magnetic systems placed in confined geometries. The influence of interfaces between magnetic domains of different orientation on the spreading of the perturbation is also discussed, showing that the presence of interfaces enhances the propagation of the damage. Furthermore, the critical transition between propagation and nonpropagation of the damage is discussed. In all cases, the determined critical exponents suggest that the DS transition does not belong to the universality class of Directed Percolation, unlike many other systems exhibiting irreversible phase transitions. This result reflects the dramatic influence of interfaces on the propagation of perturbations in magnetic systems.
The authors use the method of 'damage spreading' to measure the time taken for the equilibrium damage (the magnetization) to disappear in finite systems (10<or=L<or=103) at the critical temperature of the infinite system. This time is found to scale as L z with an exponent z of 2.16+or-0.02. This value is compatible with the expected value of the dynamic exponent for the 2D Ising model, but they argue that this method is unable to measure fluctuations properly which they show are responsible for the higher values recently reported in the literature. They suggest that the slowest mode of relaxation has an exponent of z>or=2.27.
We study damage spreading in the ferromagnetic Ising model, on the square lattice, by performing numerical simulations in the microcanonical ensemble. Although the damage spreads for all temperature values, a dynamical critical phenomenon is observed at Td{\mathit{T}}_{\mathit{d}}Td​=1.020(2)TC{\mathit{T}}_{\mathit{C}}TC​ (TC{\mathit{T}}_{\mathit{C}}TC​ is the Curie temperature): For T\ensuremath{\ge}Td{\mathit{T}}_{\mathit{d}}Td​ the damage D(t) is insensitive to its initial value D(0), while for TTd{\mathit{T}}_{\mathit{d}}Td​ a weak dependence on D(0) is present in such a way that for each T there are two values of D(t) corresponding to the initial conditions D(0)1/2 and D(0)>1/2. We have also identified a convenient order parameter to characterize the transition occurring at Td{\mathit{T}}_{\mathit{d}}Td​. \textcopyright{} 1996 The American Physical Society.
We study the dynamics of ordering in ferromagnets via Monte Carlo simulations of the Ising model, employing the Glauber spin-flip mechanism, in space dimensions d=2d=2d=2 and 333. Results for the persistence probability and the domain growth are discussed for quenches to various temperatures (TfT_fTf​) below the critical one (TcT_{c}Tc​), from different initial temperatures Ti≥TcT_{i} \geq T_{c}Ti​≥Tc​. In long time limit, for Ti>TcT_{i} > T_{c}Ti​>Tc​, the persistence probability exhibits power-law decay with exponents θ≃0.22θ\simeq 0.22θ≃0.22 and ≃0.18\simeq 0.18≃0.18 in d=2d=2d=2 and 333, respectively. For finite TiT_iTi​, the early time behavior is a different power-law whose life-time diverges and exponent decreases as Ti→TcT_{i} \rightarrow T_{c}Ti​→Tc​. The crossover length between the two steps diverges as the equilibrium correlation length. Ti=TcT_i=T_cTi​=Tc​ is expected to provide a {\it{new universality class}} for which we obtain θ≃0.035θ\simeq 0.035θ≃0.035 in d=2d=2d=2 and ≃0.10\simeq 0.10≃0.10 in d=3d=3d=3. The time dependence of the average domain size ℓ\ellℓ, however, is observed to be rather insensitive to the choice of TiT_iTi​.
The nearest-neighbor-interaction spin-1 Ising model is investigated within the damage-spreading approach. Exact relations involving quantities computable through damage-spreading simulations and thermodynamic properties are derived for such a model, defined in terms of a very general Hamiltonian that covers several spin-1 models of interest in the literature. Such relations presuppose translational invariance and hold for any ergodic dynamical procedure, leading to an efficient tool for obtaining thermodynamic properties. The implementation of the method is illustrated through damage-spreading simulations for the ferromagnetic spin-1 Ising model on a square lattice. The two-spin correlation function and the magnetization are obtained, with precise estimates of their associated critical exponents and of the critical temperature of the model, in spite of the small lattice sizes considered. These results are in good agreement with the universality hypothesis, with critical exponents in the same universality class of the spin- 12 Ising model. The advantage of the present method is shown through a significant reduction of finite-size effects by comparing its results with those obtained from standard Monte Carlo simulations.