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

Zero-temperature ordering dynamics in a two-dimensional biaxial next-nearest-neighbor Ising model

Soham Biswas, Mauricio Martin Saavedra Contreras

Published Oct 22, 2019• 2 citations
arXivHTMLPDF

Abstract

We investigate the dynamics of a two-dimensional biaxial next-nearest-neighbor Ising model following a quench to zero temperature. The Hamiltonian is given by H=-J_{0}∑_{i,j=1}^{L}[(S_{i,j}S_{i+1,j}+S_{i,j}S_{i,j+1})-κ(S_{i,j}S_{i+2,j}+S_{i,j}S_{i,j+2})]. For κ<1, the system does not reach the equilibrium ground state and keeps evolving in active states forever. For κ≥1, though, the system reaches a final state, but it does not reach the ground state always and freezes to a striped state with a finite probability like a two-dimensional ferromagnetic Ising model and an axial next-nearest-neighbor Ising (ANNNI) model. The overall dynamical behavior for κ>1 and κ=1 is quite different. The residual energy decays in a power law for both κ>1 and κ=1 from which the dynamical exponent z has been estimated. The persistence probability shows algebraic decay for κ>1 with an exponent θ=0.22±0.002 while the dynamical exponent for ordering z=2.33±0.01. For κ=1, the system belongs to a completely different dynamical class with θ=0.332±0.002 and z=2.47±0.04. We have computed the freezing probability for different values of κ. We have also studied the decay of autocorrelation function with time for a different regime of κ values. The results have been compared with those of the two-dimensional ANNNI model.

DOI
10.1103/physreve.100.042129

Related papers

Zero-temperature dynamics in the two-dimensional axial next-nearest-neighbor Ising model

Soham Biswas, Anjan Kumar Chandra, Parongama Sen

We investigate the dynamics of a two-dimensional axial next-nearest-neighbor Ising model following a quench to zero temperature. The Hamiltonian is given by H= -J_(0) summation operator(L)_(i,j=1)S_(i,j)S_(i+1,j)-J_(1)summation operator_(i,j=1)(S_{i,j}S_{i,j+1}-kappaS_{i,j}S_{i,j+2}) . For kappa<1 , the system does not reach the equilibrium ground state but slowly evolves to a metastable state. For kappa>1 , the system shows a behavior similar to that of the two-dimensional ferromagnetic Ising model in the sense that it freezes to a striped state with a finite probability. The persistence probability shows algebraic decay here with an exponent theta=0.235+/-0.001 while the dynamical exponent of growth z=2.08+/-0.01 . For kappa=1 , the system belongs to a completely different dynamical class; it always evolves to the true ground state with the persistence and dynamical exponent having unique values. Much of the dynamical phenomena can be understood by studying the dynamics and distribution of the number of domain walls. We also compare the dynamical behavior to that of a Ising model in which both the nearest and next-nearest-neighbor interactions are ferromagnetic.

arXiv:0810.2672•Physical Review E•Oct 21, 2008•10 citations

Kosterlitz-Thouless phase transition of the axial next-nearest-neighbor Ising model in two dimensions

Takayuki Shirakura, Fumitaka Matsubara, Naoki Suzuki

The spin structure of an axial next-nearest-neighbor Ising (ANNNI) model in two dimensions (2D) is a renewed problem because different Monte Carlo (MC) simulation methods predicted different spin orderings. The usual equilibrium simulation predicts the occurrence of a floating incommensurate (IC) Kosterlitz-Thouless (KT) type phase, which never emerges in non-equilibrium relaxation (NER) simulations. In this paper, we first examine previously published results of both methods, and then investigate a higher transition temperature ${T}_{c1}$ between the IC and paramagnetic phases. In the usual equilibrium simulation, we calculate the chain magnetization on larger lattices (up to $512\ifmmode\times\else\texttimes\fi{}512$ sites) and estimate ${T}_{c1}\ensuremath{\approx}1.16J$ with frustration ratio $\ensuremath{\kappa}(\ensuremath{\equiv}\ensuremath{-}{J}_{2}/{J}_{1})=0.6$. We examine the nature of the phase transition in terms of the Binder ratio ${g}_{L}$ of spin overlap functions and the correlation-length ratio $\ensuremath{\xi}/L$. In the NER simulation, we observe the spin dynamics in equilibrium states by means of an autocorrelation function and also observe the chain magnetization relaxations from the ground and disordered states. These quantities exhibit an algebraic decay at $T\ensuremath{\lesssim}1.17J$. We conclude that the two-dimensional ANNNI model actually admits an IC phase transition of the KT type.

arXiv:1406.6464•Physical Review B•Oct 8, 2014•8 citations

Nonequilibrium critical dynamics of the two-dimensional Ising model quenched from a correlated initial state

László Környei, Michel Pleimling, Ferenc Iglói

The universality class, even the order of the transition, of the two-dimensional Ising model depends on the range and the symmetry of the interactions (Onsager model, Baxter-Wu model, Turban model, etc.), but the critical temperature is generally the same due to self-duality. Here we consider a sudden change in the form of the interaction and study the nonequilibrium critical dynamical properties of the nearest-neighbor model. The relaxation of the magnetization and the decay of the autocorrelation function are found to display a power law behavior with characteristic exponents that depend on the universality class of the initial state.

arXiv:0710.2829•Physical Review E•Jan 30, 2008•2 citations

Finite-size-scaling estimation of critical exponents in the two-dimensional axial next-nearest-neighbor Ising model

Marcelo D. Grynberg, Horacio Ceva

The magnetic critical exponent of the two-dimensional axial next-nearest-neighbor Ising model is estimated by means of the finite-size scaling of the order parameter. For a finite system without ordering field, this magnitude is computed from the two leading eigenvectors of the transfer matrix. The ferromagnetic-paramagnetic transition temperatures estimated with this method exhibit a small lattice-size dependence. Moreover, they are in good agreement with those obtained using phenomenological renormalization, and also with analytical approximations. The correlation-length critical exponent is calculated by means of phenomenological renormalization. The results found with these numerical techniques indicate that the ferromagnetic-paramagnetic transition of the model seems to be continuous, in contradiction with a previous result obtained using the persistence-length criterion, which indicated that the transition is always first order. We found that the phase transition is first order only when the ratio of the next-nearest-neighbor and nearest-neighbor interactions equals one-half.

Physical review. B, Condensed matter•Nov 1, 1988•6 citations

Persistence and dynamics in the ANNNI chain

Parongama Sen, Subinay Dasgupta

We investigate both the local and global persistence behaviour in the ANNNI (axial next-nearest-neighbour Ising) model. We find that when the ratio κ of the second neighbour interaction to the first neighbour interaction is less than 1, P ( t ), the probability of a spin to remain in its original state up to time t shows a stretched exponential decay. For κ > 1, P ( t ) has an algebraic decay but the exponent is different from that of the nearest-neighbour Ising model. The global persistence behaviour shows similar features. We also conduct some deeper investigations in the dynamics of the ANNNI model and conclude that it has a different dynamical behaviour compared to the nearest-neighbour Ising model.

arXiv:cond-mat/0411504•Journal of Physics A Mathematical and General•Dec 2, 2004•24 citations

Persistence in the two dimensional ferromagnetic Ising model

Thibault Blanchard, Leticia F. Cugliandolo, Marco Picco

We present very accurate numerical estimates of the time and size dependence of the zero-temperature local persistence in the 2d ferromagnetic Ising model. We show that the effective exponent decays algebraically to an asymptotic value θ that depends upon the initial condition. More precisely, we find that θ takes one universal value 0.199(2) for initial conditions with short-range spatial correlations as in a paramagnetic state, and the value 0.033(1) for initial conditions with the long-range spatial correlations of the critical Ising state. We checked universality by working with a square and a triangular lattice, and by imposing free and periodic boundary conditions. We found that the effective exponent suffers from stronger finite size effects in the former case.

arXiv:1410.4007•Journal of Statistical Mechanics Theory and Experiment•Dec 22, 2014•21 citations
We investigate the dynamics of a two-dimensional biaxial next-nearest-neighbor Ising model following a quench to zero temperature. The Hamiltonian is given by H=-J_{0}∑{i,j=1}^{L}[(S{i,j}S_{i+1,j}+S_{i,j}S_{i,j+1})-κ(S_{i,j}S_{i+2,j}+S_{i,j}S_{i,j+2})]. For κ<1, the system does not reach the equilibrium ground state and keeps evolving in active states forever. For κ≥1, though, the system reaches a final state, but it does not reach the ground state always and freezes to a striped state with a finite probability like a two-dimensional ferromagnetic Ising model and an axial next-nearest-neighbor Ising (ANNNI) model. The overall dynamical behavior for κ>1 and κ=1 is quite different. The residual energy decays in a power law for both κ>1 and κ=1 from which the dynamical exponent z has been estimated. The persistence probability shows algebraic decay for κ>1 with an exponent θ=0.22±0.002 while the dynamical exponent for ordering z=2.33±0.01. For κ=1, the system belongs to a completely different dynamical class with θ=0.332±0.002 and z=2.47±0.04. We have computed the freezing probability for different values of κ. We have also studied the decay of autocorrelation function with time for a different regime of κ values. The results have been compared with those of the two-dimensional ANNNI model.
We investigate the dynamics of a two-dimensional axial next-nearest-neighbor Ising model following a quench to zero temperature. The Hamiltonian is given by H= -J_(0) summation operator(L)_(i,j=1)S_(i,j)S_(i+1,j)-J_(1)summation operator_(i,j=1)(S_{i,j}S_{i,j+1}-kappaS_{i,j}S_{i,j+2}) . For kappa<1 , the system does not reach the equilibrium ground state but slowly evolves to a metastable state. For kappa>1 , the system shows a behavior similar to that of the two-dimensional ferromagnetic Ising model in the sense that it freezes to a striped state with a finite probability. The persistence probability shows algebraic decay here with an exponent theta=0.235+/-0.001 while the dynamical exponent of growth z=2.08+/-0.01 . For kappa=1 , the system belongs to a completely different dynamical class; it always evolves to the true ground state with the persistence and dynamical exponent having unique values. Much of the dynamical phenomena can be understood by studying the dynamics and distribution of the number of domain walls. We also compare the dynamical behavior to that of a Ising model in which both the nearest and next-nearest-neighbor interactions are ferromagnetic.
The spin structure of an axial next-nearest-neighbor Ising (ANNNI) model in two dimensions (2D) is a renewed problem because different Monte Carlo (MC) simulation methods predicted different spin orderings. The usual equilibrium simulation predicts the occurrence of a floating incommensurate (IC) Kosterlitz-Thouless (KT) type phase, which never emerges in non-equilibrium relaxation (NER) simulations. In this paper, we first examine previously published results of both methods, and then investigate a higher transition temperature Tc1{T}_{c1}Tc1​ between the IC and paramagnetic phases. In the usual equilibrium simulation, we calculate the chain magnetization on larger lattices (up to $512\ifmmode\times\else\texttimes\fi{}512$ sites) and estimate ${T}_{c1}\ensuremath{\approx}1.16J$ with frustration ratio $\ensuremath{\kappa}(\ensuremath{\equiv}\ensuremath{-}{J}_{2}/{J}_{1})=0.6$. We examine the nature of the phase transition in terms of the Binder ratio gL{g}_{L}gL​ of spin overlap functions and the correlation-length ratio $\ensuremath{\xi}/L$. In the NER simulation, we observe the spin dynamics in equilibrium states by means of an autocorrelation function and also observe the chain magnetization relaxations from the ground and disordered states. These quantities exhibit an algebraic decay at $T\ensuremath{\lesssim}1.17J$. We conclude that the two-dimensional ANNNI model actually admits an IC phase transition of the KT type.
The universality class, even the order of the transition, of the two-dimensional Ising model depends on the range and the symmetry of the interactions (Onsager model, Baxter-Wu model, Turban model, etc.), but the critical temperature is generally the same due to self-duality. Here we consider a sudden change in the form of the interaction and study the nonequilibrium critical dynamical properties of the nearest-neighbor model. The relaxation of the magnetization and the decay of the autocorrelation function are found to display a power law behavior with characteristic exponents that depend on the universality class of the initial state.
The magnetic critical exponent of the two-dimensional axial next-nearest-neighbor Ising model is estimated by means of the finite-size scaling of the order parameter. For a finite system without ordering field, this magnitude is computed from the two leading eigenvectors of the transfer matrix. The ferromagnetic-paramagnetic transition temperatures estimated with this method exhibit a small lattice-size dependence. Moreover, they are in good agreement with those obtained using phenomenological renormalization, and also with analytical approximations. The correlation-length critical exponent is calculated by means of phenomenological renormalization. The results found with these numerical techniques indicate that the ferromagnetic-paramagnetic transition of the model seems to be continuous, in contradiction with a previous result obtained using the persistence-length criterion, which indicated that the transition is always first order. We found that the phase transition is first order only when the ratio of the next-nearest-neighbor and nearest-neighbor interactions equals one-half.
We investigate both the local and global persistence behaviour in the ANNNI (axial next-nearest-neighbour Ising) model. We find that when the ratio κ of the second neighbour interaction to the first neighbour interaction is less than 1, P ( t ), the probability of a spin to remain in its original state up to time t shows a stretched exponential decay. For κ > 1, P ( t ) has an algebraic decay but the exponent is different from that of the nearest-neighbour Ising model. The global persistence behaviour shows similar features. We also conduct some deeper investigations in the dynamics of the ANNNI model and conclude that it has a different dynamical behaviour compared to the nearest-neighbour Ising model.
We present very accurate numerical estimates of the time and size dependence of the zero-temperature local persistence in the 2d ferromagnetic Ising model. We show that the effective exponent decays algebraically to an asymptotic value θ that depends upon the initial condition. More precisely, we find that θ takes one universal value 0.199(2) for initial conditions with short-range spatial correlations as in a paramagnetic state, and the value 0.033(1) for initial conditions with the long-range spatial correlations of the critical Ising state. We checked universality by working with a square and a triangular lattice, and by imposing free and periodic boundary conditions. We found that the effective exponent suffers from stronger finite size effects in the former case.