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:1907.10296 · Physical review. E

Approximate dynamical eigenmodes of the Ising model with local spin-exchange moves

Wei Zhong, Debabrata Panja, G. T. Barkema

Published Jul 23, 2019• 4 citations
arXivHTMLPDF

Abstract

We establish that the Fourier modes of the magnetization serve as the dynamical eigenmodes for the two-dimensional Ising model at the critical temperature with local spin-exchange moves, i.e., Kawasaki dynamics. We obtain the dynamical scaling properties for these modes and use them to calculate the time evolution of two dynamical quantities for the system, namely, the autocorrelation function and the mean-square deviation of the line magnetizations. At intermediate times 1≲t≲L^{z_{c}}, where z_{c}=4-η=15/4 is the dynamical critical exponent of the model, we find that the line magnetization undergoes anomalous diffusion. Following our recent work on anomalous diffusion in spin models, we demonstrate that the generalized Langevin equation with a memory kernel consistently describes the anomalous diffusion, verifying the corresponding fluctuation-dissipation theorem with the calculation of the force autocorrelation function.

DOI
10.1103/physreve.100.012132

Related papers

Critical dynamical behavior of the Ising model

Zihua Liu, Erol Vatansever, G. T. Barkema +1 more

We investigate the dynamical critical behavior of the two- and three-dimensional Ising models with Glauber dynamics in equilibrium. In contrast to the usual standing, we focus on the mean-squared deviation of the magnetization $M, {\mathrm{MSD}}_{M}$, as a function of time, as well as on the autocorrelation function of $M$. These two functions are distinct but closely related. We find that ${\mathrm{MSD}}_{M}$ features a first crossover at time ${\ensuremath{\tau}}_{1}\ensuremath{\sim}{L}^{{z}_{1}}$, from ordinary diffusion with ${\mathrm{MSD}}_{M}\ensuremath{\sim}t$, to anomalous diffusion with ${\mathrm{MSD}}_{M}\ensuremath{\sim}{t}^{\ensuremath{\alpha}}$. Purely on numerical grounds, we obtain the values ${z}_{1}=0.45(5)$ and $\ensuremath{\alpha}=0.752(5)$ for the two-dimensional Ising ferromagnet. Related to this, the magnetization autocorrelation function crosses over from an exponential decay to a stretched-exponential decay. At later times, we find a second crossover at time ${\ensuremath{\tau}}_{2}\ensuremath{\sim}{L}^{{z}_{2}}$. Here, ${\mathrm{MSD}}_{M}$ saturates to its late-time value $\ensuremath{\sim}{L}^{2+\ensuremath{\gamma}/\ensuremath{\nu}}$, while the autocorrelation function crosses over from stretched-exponential decay to simple exponential one. We also confirm numerically the value ${z}_{2}=2.1665(12)$, earlier reported as the single dynamic exponent. Continuity of ${\mathrm{MSD}}_{M}$ requires that $\ensuremath{\alpha}({z}_{2}\ensuremath{-}{z}_{1})=\ensuremath{\gamma}/\ensuremath{\nu}\ensuremath{-}{z}_{1}$. We speculate that ${z}_{1}=1/2$ and $\ensuremath{\alpha}=3/4$, values that indeed lead to the expected ${z}_{2}=13/6$ result. A complementary analysis for the three-dimensional Ising model provides the estimates ${z}_{1}=1.35(2), \ensuremath{\alpha}=0.90(2)$, and ${z}_{2}=2.032(3)$. While ${z}_{2}$ has attracted significant attention in the literature, we argue that for all practical purposes ${z}_{1}$ is more important, as it determines the number of statistically independent measurements during a long simulation.

arXiv:2307.01837•Physical review. E•Sep 21, 2023•14 citations

Dynamical critical behavior of the two-dimensional three-state Potts model

Erol Vatansever, G. T. Barkema, Nikolaos G. Fytas

We investigate the dynamical critical behavior of the two-dimensional three-state Potts model with single spin-flip dynamics in equilibrium. We focus on the mean-squared deviation of the magnetization M (MSD_{M}) as a function of time, as well as on the autocorrelation function of M. Our simulations reveal the existence of two crossover behaviors at times τ_{1}∼L^{z_{1}} and τ_{2}∼L^{z_{2}}, separating three dynamical regimes. MSD_{M} appears to shift from ordinary diffusion in the first regime, to anomalous diffusion in the second, and finally to be constant in the third regime. The magnetization autocorrelation function, on the other hand, is found to fluctuate between exponential decay, stretched-exponential decay, and then again exponential decay along these three regimes. This behavior is in agreement with the one reported recently for the two-dimensional Ising ferromagnet [Phys. Rev. E 108, 034118 (2023)2470-004510.1103/PhysRevE.108.034118], indicating that the existence of two dynamic critical exponents is not a peculiarity of the Ising model itself. A comparison of both MSD_{M} and the magnetization's autocorrelation function suggests that within our numerical accuracy the exponents z_{1} and z_{2} are shared between the Ising and three-state Potts models at least for the particular case of single spin-flip dynamics studied here, even though their equilibrium universality classes are clearly distinct. Continuity of MSD_{M} requires that α(z_{2}-z_{1})=γ/ν-z_{1}, in which α is the anomalous exponent in the intermediate regime. Since the ratio γ/ν is not shared between the two models, it follows that α is not shared either, an aspect well verified in our simulations. Finally, we also discuss the relevance of our main findings using another useful observable, namely the line magnetization M_{l}.

arXiv:2407.08198•Physical review. E•Jul 23, 2024•2 citations

Non-equilibrium critical dynamics of the ferromagnetic Ising model with Kawasaki dynamics

Claude Godrèche, Florent Krzaka a, Federico Ricci‐Tersenghi

We investigate the temporal evolution of a ferromagnetic system of Ising spins evolving under Kawasaki dynamics from a random initial condition, in spatial dimensions one and two. We examine in detail the asymptotic behaviour of the two-time correlation and response functions. The linear response is measured without applying a field, using a recently proposed algorithm. For the chain at vanishingly small temperature, we introduce an accelerated dynamics which has the virtue of projecting the system into the asymptotic scaling regime. This allows us to revisit critically previous works on the behaviour at large time of the two-time autocorrelation and response functions. We also analyse the case of the two-dimensional system at criticality. A comparison with Glauber dynamics is performed in both dimensionalities, in order to underline the similarities and differences in the phenomenology of the two dynamics.

arXiv:cond-mat/0401334•Journal of Statistical Mechanics Theory and Experiment•May 1, 2004•41 citations

Critical dynamical exponent of the two-dimensional scalar ϕ4 model with local moves

Wei Zhong, G. T. Barkema, Debabrata Panja +1 more

We study the scalar one-component two-dimensional (2D) ${\ensuremath{\phi}}^{4}$ model by computer simulations, with local Metropolis moves. The equilibrium exponents of this model are well established, e.g., for the 2D ${\ensuremath{\phi}}^{4}$ model $\ensuremath{\gamma}=1.75$ and $\ensuremath{\nu}=1$. The model has also been conjectured to belong to the Ising universality class. However, the value of the critical dynamical exponent ${z}_{c}$ is not settled. In this paper, we obtain ${z}_{c}$ for the 2D ${\ensuremath{\phi}}^{4}$ model using two independent methods: (a) by calculating the relative terminal exponential decay time $\ensuremath{\tau}$ for the correlation function $\ensuremath{\langle}\mathrm{\ensuremath{\Phi}}(t)\mathrm{\ensuremath{\Phi}}(0)\ensuremath{\rangle}$, and thereafter fitting the data as $\ensuremath{\tau}\ensuremath{\sim}{L}^{{z}_{c}}$, where $L$ is the system size, and (b) by measuring the anomalous diffusion exponent for the order parameter, viz., the mean-square displacement $\ensuremath{\langle}\mathrm{\ensuremath{\Delta}}{\mathrm{\ensuremath{\Phi}}}^{2}(t)\ensuremath{\rangle}\ensuremath{\sim}{t}^{c}$ as $c=\ensuremath{\gamma}/(\ensuremath{\nu}{z}_{c})$, and from the numerically obtained value $c\ensuremath{\approx}0.80$, we calculate ${z}_{c}$. For different values of the coupling constant $\ensuremath{\lambda}$, we report that ${z}_{c}=2.17\ifmmode\pm\else\textpm\fi{}0.03$ and ${z}_{c}=2.19\ifmmode\pm\else\textpm\fi{}0.03$ for the two methods, respectively. Our results indicate that ${z}_{c}$ is independent of $\ensuremath{\lambda}$, and is likely identical to that for the 2D Ising model. Additionally, we demonstrate that the generalized Langevin equation formulation with a memory kernel, identical to those applicable for the Ising model and polymeric systems, consistently captures the observed anomalous diffusion behavior.

arXiv:1807.05123•Physical review. E•Dec 19, 2018•8 citations

Dynamical scaling and decay of correlations for spinodal decomposition at Tc

Francis J. Alexander, David A. Huse, Steven A. Janowsky

We have investigated the nonequilibrium critical dynamics of a two-dimensional, spin-exchange (i.e., conserved magnetization), kinetic Ising model following quenches to ${\mathit{T}}_{\mathit{c}}$ from random initial states. As in the case of subcritical quenches, we observe dynamical scaling of the equal time order parameter correlation function G(r,t,t), but here the nonequilibrium correlation length grows as L(t)\ensuremath{\sim}${\mathit{t}}^{1/\mathit{z}}$, where z=4/15 is the usual dynamic critical exponent for this system. Moreover, we find dynamical scaling for correlations with the initial condition G(r,0,t). The temporal decay of autocorrelations at ${\mathit{T}}_{\mathit{c}}$ yields a critical exponent ${\ensuremath{\lambda}}_{\mathit{c}}$=1.96\ifmmode\pm\else\textpm\fi{}0.05 for the nonequilibrium dynamics: G(r,0,t)\ensuremath{\approxeq}${\mathit{L}}_{\mathit{c}}^{\mathrm{\ensuremath{-}}\ensuremath{\lambda}}$(t)G\ifmmode \tilde{}\else \~{}\fi{}(r/L(t)).

Physical review. B, Condensed matter•Jul 1, 1994•16 citations

Crossover of the dynamical behavior in two-dimensional random transverse Ising model

Shu-Xia Chen, Yin-Yang Shen, Xiang-Mu Kong

The dynamics of the two-dimensional spin-$\frac{1}{2}$ random transverse Ising model is studied by means of the recurrence relations method. Both the autocorrelation functions and the corresponding spectral densities are calculated in the cases that the exchange couplings or fields satisfy three typical distributions, respectively. For the cases of Gaussian and double-Gaussian distribution, when the standard deviation of the random variable is small, the dynamics of the system undergoes a crossover from a central-peak behavior to a collective-mode one, the result is similar to that of the bimodal distribution. While the standard deviation is large enough, the crossover disappears and the dynamics of the system shows a central-peak behavior or a disordered behavior.

Physical Review B•Nov 3, 2010•25 citations
We establish that the Fourier modes of the magnetization serve as the dynamical eigenmodes for the two-dimensional Ising model at the critical temperature with local spin-exchange moves, i.e., Kawasaki dynamics. We obtain the dynamical scaling properties for these modes and use them to calculate the time evolution of two dynamical quantities for the system, namely, the autocorrelation function and the mean-square deviation of the line magnetizations. At intermediate times 1≲t≲L^{z_{c}}, where z_{c}=4-η=15/4 is the dynamical critical exponent of the model, we find that the line magnetization undergoes anomalous diffusion. Following our recent work on anomalous diffusion in spin models, we demonstrate that the generalized Langevin equation with a memory kernel consistently describes the anomalous diffusion, verifying the corresponding fluctuation-dissipation theorem with the calculation of the force autocorrelation function.
We investigate the dynamical critical behavior of the two- and three-dimensional Ising models with Glauber dynamics in equilibrium. In contrast to the usual standing, we focus on the mean-squared deviation of the magnetization M,MSDMM, {\mathrm{MSD}}_{M}M,MSDM​, as a function of time, as well as on the autocorrelation function of MMM. These two functions are distinct but closely related. We find that MSDM{\mathrm{MSD}}_{M}MSDM​ features a first crossover at time ${\ensuremath{\tau}}_{1}\ensuremath{\sim}{L}^{{z}_{1}}$, from ordinary diffusion with ${\mathrm{MSD}}_{M}\ensuremath{\sim}t$, to anomalous diffusion with ${\mathrm{MSD}}_{M}\ensuremath{\sim}{t}^{\ensuremath{\alpha}}$. Purely on numerical grounds, we obtain the values z1=0.45(5){z}_{1}=0.45(5)z1​=0.45(5) and $\ensuremath{\alpha}=0.752(5)$ for the two-dimensional Ising ferromagnet. Related to this, the magnetization autocorrelation function crosses over from an exponential decay to a stretched-exponential decay. At later times, we find a second crossover at time ${\ensuremath{\tau}}_{2}\ensuremath{\sim}{L}^{{z}_{2}}$. Here, MSDM{\mathrm{MSD}}_{M}MSDM​ saturates to its late-time value $\ensuremath{\sim}{L}^{2+\ensuremath{\gamma}/\ensuremath{\nu}}$, while the autocorrelation function crosses over from stretched-exponential decay to simple exponential one. We also confirm numerically the value z2=2.1665(12){z}_{2}=2.1665(12)z2​=2.1665(12), earlier reported as the single dynamic exponent. Continuity of MSDM{\mathrm{MSD}}_{M}MSDM​ requires that $\ensuremath{\alpha}({z}_{2}\ensuremath{-}{z}_{1})=\ensuremath{\gamma}/\ensuremath{\nu}\ensuremath{-}{z}_{1}$. We speculate that z1=1/2{z}_{1}=1/2z1​=1/2 and $\ensuremath{\alpha}=3/4$, values that indeed lead to the expected z2=13/6{z}_{2}=13/6z2​=13/6 result. A complementary analysis for the three-dimensional Ising model provides the estimates ${z}_{1}=1.35(2), \ensuremath{\alpha}=0.90(2)$, and z2=2.032(3){z}_{2}=2.032(3)z2​=2.032(3). While z2{z}_{2}z2​ has attracted significant attention in the literature, we argue that for all practical purposes z1{z}_{1}z1​ is more important, as it determines the number of statistically independent measurements during a long simulation.
We investigate the dynamical critical behavior of the two-dimensional three-state Potts model with single spin-flip dynamics in equilibrium. We focus on the mean-squared deviation of the magnetization M (MSD_{M}) as a function of time, as well as on the autocorrelation function of M. Our simulations reveal the existence of two crossover behaviors at times τ_{1}∼L^{z_{1}} and τ_{2}∼L^{z_{2}}, separating three dynamical regimes. MSD_{M} appears to shift from ordinary diffusion in the first regime, to anomalous diffusion in the second, and finally to be constant in the third regime. The magnetization autocorrelation function, on the other hand, is found to fluctuate between exponential decay, stretched-exponential decay, and then again exponential decay along these three regimes. This behavior is in agreement with the one reported recently for the two-dimensional Ising ferromagnet [Phys. Rev. E 108, 034118 (2023)2470-004510.1103/PhysRevE.108.034118], indicating that the existence of two dynamic critical exponents is not a peculiarity of the Ising model itself. A comparison of both MSD_{M} and the magnetization's autocorrelation function suggests that within our numerical accuracy the exponents z_{1} and z_{2} are shared between the Ising and three-state Potts models at least for the particular case of single spin-flip dynamics studied here, even though their equilibrium universality classes are clearly distinct. Continuity of MSD_{M} requires that α(z_{2}-z_{1})=γ/ν-z_{1}, in which α is the anomalous exponent in the intermediate regime. Since the ratio γ/ν is not shared between the two models, it follows that α is not shared either, an aspect well verified in our simulations. Finally, we also discuss the relevance of our main findings using another useful observable, namely the line magnetization M_{l}.
We investigate the temporal evolution of a ferromagnetic system of Ising spins evolving under Kawasaki dynamics from a random initial condition, in spatial dimensions one and two. We examine in detail the asymptotic behaviour of the two-time correlation and response functions. The linear response is measured without applying a field, using a recently proposed algorithm. For the chain at vanishingly small temperature, we introduce an accelerated dynamics which has the virtue of projecting the system into the asymptotic scaling regime. This allows us to revisit critically previous works on the behaviour at large time of the two-time autocorrelation and response functions. We also analyse the case of the two-dimensional system at criticality. A comparison with Glauber dynamics is performed in both dimensionalities, in order to underline the similarities and differences in the phenomenology of the two dynamics.
We study the scalar one-component two-dimensional (2D) ${\ensuremath{\phi}}^{4}$ model by computer simulations, with local Metropolis moves. The equilibrium exponents of this model are well established, e.g., for the 2D ${\ensuremath{\phi}}^{4}$ model $\ensuremath{\gamma}=1.75$ and $\ensuremath{\nu}=1$. The model has also been conjectured to belong to the Ising universality class. However, the value of the critical dynamical exponent zc{z}_{c}zc​ is not settled. In this paper, we obtain zc{z}_{c}zc​ for the 2D ${\ensuremath{\phi}}^{4}$ model using two independent methods: (a) by calculating the relative terminal exponential decay time $\ensuremath{\tau}$ for the correlation function $\ensuremath{\langle}\mathrm{\ensuremath{\Phi}}(t)\mathrm{\ensuremath{\Phi}}(0)\ensuremath{\rangle}$, and thereafter fitting the data as $\ensuremath{\tau}\ensuremath{\sim}{L}^{{z}_{c}}$, where LLL is the system size, and (b) by measuring the anomalous diffusion exponent for the order parameter, viz., the mean-square displacement $\ensuremath{\langle}\mathrm{\ensuremath{\Delta}}{\mathrm{\ensuremath{\Phi}}}^{2}(t)\ensuremath{\rangle}\ensuremath{\sim}{t}^{c}$ as $c=\ensuremath{\gamma}/(\ensuremath{\nu}{z}_{c})$, and from the numerically obtained value $c\ensuremath{\approx}0.80$, we calculate zc{z}_{c}zc​. For different values of the coupling constant $\ensuremath{\lambda}$, we report that ${z}_{c}=2.17\ifmmode\pm\else\textpm\fi{}0.03$ and ${z}_{c}=2.19\ifmmode\pm\else\textpm\fi{}0.03$ for the two methods, respectively. Our results indicate that zc{z}_{c}zc​ is independent of $\ensuremath{\lambda}$, and is likely identical to that for the 2D Ising model. Additionally, we demonstrate that the generalized Langevin equation formulation with a memory kernel, identical to those applicable for the Ising model and polymeric systems, consistently captures the observed anomalous diffusion behavior.
We have investigated the nonequilibrium critical dynamics of a two-dimensional, spin-exchange (i.e., conserved magnetization), kinetic Ising model following quenches to Tc{\mathit{T}}_{\mathit{c}}Tc​ from random initial states. As in the case of subcritical quenches, we observe dynamical scaling of the equal time order parameter correlation function G(r,t,t), but here the nonequilibrium correlation length grows as L(t)\ensuremath{\sim}t1/z{\mathit{t}}^{1/\mathit{z}}t1/z, where z=4/15 is the usual dynamic critical exponent for this system. Moreover, we find dynamical scaling for correlations with the initial condition G(r,0,t). The temporal decay of autocorrelations at Tc{\mathit{T}}_{\mathit{c}}Tc​ yields a critical exponent ${\ensuremath{\lambda}}_{\mathit{c}}$=1.96\ifmmode\pm\else\textpm\fi{}0.05 for the nonequilibrium dynamics: G(r,0,t)\ensuremath{\approxeq}${\mathit{L}}_{\mathit{c}}^{\mathrm{\ensuremath{-}}\ensuremath{\lambda}}$(t)G\ifmmode \tilde{}\else \~{}\fi{}(r/L(t)).
The dynamics of the two-dimensional spin-12\frac{1}{2}21​ random transverse Ising model is studied by means of the recurrence relations method. Both the autocorrelation functions and the corresponding spectral densities are calculated in the cases that the exchange couplings or fields satisfy three typical distributions, respectively. For the cases of Gaussian and double-Gaussian distribution, when the standard deviation of the random variable is small, the dynamics of the system undergoes a crossover from a central-peak behavior to a collective-mode one, the result is similar to that of the bimodal distribution. While the standard deviation is large enough, the crossover disappears and the dynamics of the system shows a central-peak behavior or a disordered behavior.