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

Dynamical phase transitions in the nonreciprocal Ising model

Yael Avni, Michel Fruchart, David Martin, Daniel S. Seara, Vincenzo Vitelli

Published Mar 21, 2025• 24 citations
arXivHTMLPDF

Abstract

Nonreciprocal interactions in many-body systems lead to time-dependent states, commonly observed in biological, chemical, and ecological systems. The stability of these states in the thermodynamic limit and the critical behavior of the phase transition from static to time-dependent states are not yet fully understood. To address these questions, we study a minimalistic system endowed with nonreciprocal interactions: an Ising model with two spin species having opposing goals. The mean-field equation predicts three stable phases: disorder, static order, and a time-dependent swap phase. Large-scale numerical simulations support the following: (i) in two dimensions, the swap phase is destabilized by defects; (ii) in three dimensions, the swap phase is stable and has the properties of a time crystal; (iii) the transition from disorder to swap in three dimensions is characterized by the critical exponents of the 3D XY model and corresponds to the breaking of a continuous symmetry, time translation invariance; (iv) when the two species have fully antisymmetric couplings, the static-order phase is unstable in any finite dimension due to droplet growth; and (v) in the general case of asymmetric couplings, static order can be restored by a droplet-capture mechanism preventing the droplets from growing indefinitely. We provide details on the full phase diagram, which includes first- and second-order-like phase transitions, and study how the system coarsens into swap and static-order states.

DOI
10.1103/physreve.111.034124

Related papers

Spin Systems far from Equilibrium: Aging and Dynamic Phase Transition

Hyunhang Park

Among the many non-equilibrium processes encountered in nature we deal with two different but related aspects. One is the non-equilibrium relaxation process that is at the origin of \'aging phenomena••, and the other one is a non-equilibrium phase transition, called ••dynamic phase transition••. One of the main purposes of our research is to explore more realistic situations than studied previously. Indeed, in the study of aging phenomena certain kinds of disorder effects are considered, and we introduce the ••surface•• as a spatial boundary to the system undergoing the dynamic phase transition. In order to observe these processes as clearly as possible, we study in both cases simple spin systems. Using Monte Carlo simulations we first investigate aging in three-dimensional Ising spin glasses as well as in two-dimensional Ising models with disorder quenched to low temperatures. The time-dependent dynamical correlation length L(t) is determined numerically and the scaling behavior of various two-time quantities as a function of L(t)/L(s) is discussed where t and s are two different times. For disordered Ising models deviations of L(t) from algebraic growth law show up. The generalized scaling forms as a function of L(t)/L(s) reveal a generic simple aging scenario for Ising spin glasses as well as for disordered Ising ferromagnets. We also study the local critical phenomena at a dynamic phase transition by means of numerical simulations of kinetic Ising models with surfaces subjected to a periodic oscillating field. We examine layer-dependent quantities, such as the period-averaged magnetization per layer Q(z) and the layer susceptibility ¥ö(z), and determine local critical exponents through finite size scaling. Both for two and three dimensions, we find that the values of the surface exponents differ from those of the equilibrium critical surface. It is revealed that the surface phase diagram of the non-equilibrium system is not identical to that of the equilibrium system in three dimensions.

•Mar 22, 2013•0 citations

Phase Transitions and Static Spin Correlations in Ising Models with Free Surfaces

Kurt Binder, Pierre C. Hohenberg

Phase transitions in Ising models with tree surfaces are studied from various points of view, including a phenomenological Landau theory, high-temperature series expansions, and a scaling theory for thermodynamic quantities and correlation functions. In the presence of a surface a number of new critical exponents must be defined. These arise because of the existence of "surface" terms in the thermodynamic functions, and because of the anisotropy of space and lack of translational symmetry introduced by the surface. The need for these new critical exponents already appears in the phenomenological theory, which is discussed in detail and related to the microscopic mean-field approximation. The essential new parameter appearing in this theory is an extrapolation length $\ensuremath{\lambda}$ which enters the boundary condition on the magnetization at the surface. For magnetic systems this length is of the order of the interaction range, in contrast to superconductors, where it is usually much larger. In order to go beyond the mean-field theory, high-temperature series expansions are carried out for the Ising half-space, to tenth order in two dimensions and to eighth order in three dimensions. A scaling theory is developed both for thermodynamic functions and for spin correlations near the surface, and relations are found among the exponents of the half-space. Both the scaling theory and the numerical calculations are compared with the exact solution of the Ising half-plane (two dimensions) by McCoy and Wu, and agreement is found wherever the theory is applicable. In analogy to the bulk situation, the scaling theory is found to agree with mean-field theory in four dimensions. The prediction of the present work which is most easily accessible to experiment is the temperature dependence of the magnetization at the surface, with critical exponent estimated to be ${\ensuremath{\beta}}_{1}=\frac{2}{3}$. The mean-field result, ${\ensuremath{\beta}}_{1}=1$, seems to agree more closely with presently available experiment, and more work is needed to clarify the situation.

Physical review. B, Solid state•Nov 1, 1972•452 citations

Dynamical phase transitions in the X Y model: A Monte Carlo and mean-field-theory study

Mainak Pal, William Baez, Pushan Majumdar +2 more

We investigate the dynamical phases and phase transitions arising in a classical two-dimensional anisotropic XY model under the influence of a periodically driven temporal external magnetic field in the form of a symmetric square wave. We use a combination of finite temperature classical Monte Carlo simulation, implemented within a CPU+GPU paradigm, utilizing local dynamics provided by the Glauber algorithm and a phenomenological equation-of-motion approach based on relaxational dynamics governed by the time-dependent free energy within a mean-field approximation to study the model. We investigate several parameter regimes of the variables (magnetic field, anisotropy, and the external drive frequency) that influence the anisotropic XY system. We identify four possible dynamical phases: Ising-SBO, Ising-SRO, XY-SBO, and XY-SRO. Both techniques indicate that only three of them (Ising-SRO, Ising-SBO, and XY-SRO) are stable dynamical phases in the thermodynamic sense. Within the Monte Carlo framework, a finite-size scaling analysis, shows that XY-SBO does not survive in the thermodynamic limit giving way to either an Ising-SBO or a XY-SRO regime. The finite-size scaling analysis further shows that the transitions between the three remaining dynamical phases either belong to the two-dimensional Ising universality class or are first-order in nature. Within the mean-field calculations yield three stable dynamical phases, i.e., Ising-SRO, Ising-SBO and XY-SRO, where the final steady state is independent of the initial condition chosen to evolve the equations of motion, as well as a region of bistability where the system flows to either Ising-SBO or XY-SRO (Ising-SRO) depending on the initial condition. Unlike the stable dynamical phases, the XY-SBO represents a transient feature that is eventually lost to either Ising-SBO or XY-SRO. Our mean-field analysis highlights the importance of the competition between switching of the stationary point(s) of the free energy after each half cycle of the external field and the two-dimensional nature of the phase space for the equations of motion.

arXiv:2402.07505•Physical review. E•Nov 12, 2024•4 citations

Nonequilibrium critical dynamics of two dimensional interacting monomer-dimer model: non-Ising criticality

Keekwon Nam, Bongsoo Kim, Sung Jong Lee

We investigate the nonequilibrium relaxation dynamics of an interacting monomer-dimer model with nearest neighbor repulsion on a square lattice, which possesses two symmetric absorbing states. The model is known to exhibit two nearby continuous transitions: the Z 2 symmetry-breaking order-disorder transition and the absorbing transition with directed percolation criticality. We performed a more detailed analysis of our extensive simulations on bigger lattice systems which reaffirms that the symmetry-breaking transition exhibits a non-Ising critical behavior with β ≃ 0.149(2) and η ≃ 0.30(1) that are distinct from those values of a pure two dimensional Ising model. Finite size scaling of dimer density near the symmetry breaking transition gives logarithmic scaling ( α = 0.0) which is consistent with the hyperscaling relation but the corresponding exponent of ν B ≃ 1.37(2) exhibits a conspicuous deviation from the pure Ising value of 1. The value of dynamic critical exponent z , however, is found to be close to that of the kinetic Ising model as 1/ z ≃ 0.466(5) from the relaxation of staggered magnetization (and also similar but slightly smaller values from coarsening).

Journal of Statistical Mechanics Theory and Experiment•Aug 14, 2014•2 citations

Thermal melting of discrete time crystals: A dynamical phase transition induced by thermal fluctuations

Mingxi Yue, Xiaoqin Yang, Zi Cai

The stability of a discrete time crystal against thermal fluctuations has been studied numerically by solving a stochastic Landau-Lifshitz-Gilbert equation of a periodically driven classical system composed of interacting spins, each of which couples to a thermal bath. It is shown that in the thermodynamic limit, even though the long-range temporal crystalline order is stable at low temperature, it is melting above a critical temperature, at which the system experiences a nonequilibrium phase transition. The critical behaviors of the continuous phase transition have been systematically investigated, and it is shown that despite the genuine nonequilibrium feature of such a periodically driven system, its critical properties fall into the three-dimensional Ising universality class with a dynamical exponent ($z=2$) identical to that in the critical dynamics of a kinetic Ising model without driving.

arXiv:2110.15506•Physical review. B./Physical review. B•Mar 22, 2022•11 citations

Critical Dynamics of Nonconserved Ising-Like Systems

Kevin E. Bassler, B. Schmittmann

We show that the dynamical fixed point of Ising-like models, characterized by a single scalar, nonconserved ordering field, is stable near four dimensions with respect to all dynamic perturbations, including those of a nonequilibrium nature.

Physical Review Letters•Dec 19, 1994•50 citations
Nonreciprocal interactions in many-body systems lead to time-dependent states, commonly observed in biological, chemical, and ecological systems. The stability of these states in the thermodynamic limit and the critical behavior of the phase transition from static to time-dependent states are not yet fully understood. To address these questions, we study a minimalistic system endowed with nonreciprocal interactions: an Ising model with two spin species having opposing goals. The mean-field equation predicts three stable phases: disorder, static order, and a time-dependent swap phase. Large-scale numerical simulations support the following: (i) in two dimensions, the swap phase is destabilized by defects; (ii) in three dimensions, the swap phase is stable and has the properties of a time crystal; (iii) the transition from disorder to swap in three dimensions is characterized by the critical exponents of the 3D XY model and corresponds to the breaking of a continuous symmetry, time translation invariance; (iv) when the two species have fully antisymmetric couplings, the static-order phase is unstable in any finite dimension due to droplet growth; and (v) in the general case of asymmetric couplings, static order can be restored by a droplet-capture mechanism preventing the droplets from growing indefinitely. We provide details on the full phase diagram, which includes first- and second-order-like phase transitions, and study how the system coarsens into swap and static-order states.
Among the many non-equilibrium processes encountered in nature we deal with two different but related aspects. One is the non-equilibrium relaxation process that is at the origin of \'aging phenomena••, and the other one is a non-equilibrium phase transition, called ••dynamic phase transition••. One of the main purposes of our research is to explore more realistic situations than studied previously. Indeed, in the study of aging phenomena certain kinds of disorder effects are considered, and we introduce the ••surface•• as a spatial boundary to the system undergoing the dynamic phase transition. In order to observe these processes as clearly as possible, we study in both cases simple spin systems. Using Monte Carlo simulations we first investigate aging in three-dimensional Ising spin glasses as well as in two-dimensional Ising models with disorder quenched to low temperatures. The time-dependent dynamical correlation length L(t) is determined numerically and the scaling behavior of various two-time quantities as a function of L(t)/L(s) is discussed where t and s are two different times. For disordered Ising models deviations of L(t) from algebraic growth law show up. The generalized scaling forms as a function of L(t)/L(s) reveal a generic simple aging scenario for Ising spin glasses as well as for disordered Ising ferromagnets. We also study the local critical phenomena at a dynamic phase transition by means of numerical simulations of kinetic Ising models with surfaces subjected to a periodic oscillating field. We examine layer-dependent quantities, such as the period-averaged magnetization per layer Q(z) and the layer susceptibility ¥ö(z), and determine local critical exponents through finite size scaling. Both for two and three dimensions, we find that the values of the surface exponents differ from those of the equilibrium critical surface. It is revealed that the surface phase diagram of the non-equilibrium system is not identical to that of the equilibrium system in three dimensions.
Phase transitions in Ising models with tree surfaces are studied from various points of view, including a phenomenological Landau theory, high-temperature series expansions, and a scaling theory for thermodynamic quantities and correlation functions. In the presence of a surface a number of new critical exponents must be defined. These arise because of the existence of "surface" terms in the thermodynamic functions, and because of the anisotropy of space and lack of translational symmetry introduced by the surface. The need for these new critical exponents already appears in the phenomenological theory, which is discussed in detail and related to the microscopic mean-field approximation. The essential new parameter appearing in this theory is an extrapolation length $\ensuremath{\lambda}$ which enters the boundary condition on the magnetization at the surface. For magnetic systems this length is of the order of the interaction range, in contrast to superconductors, where it is usually much larger. In order to go beyond the mean-field theory, high-temperature series expansions are carried out for the Ising half-space, to tenth order in two dimensions and to eighth order in three dimensions. A scaling theory is developed both for thermodynamic functions and for spin correlations near the surface, and relations are found among the exponents of the half-space. Both the scaling theory and the numerical calculations are compared with the exact solution of the Ising half-plane (two dimensions) by McCoy and Wu, and agreement is found wherever the theory is applicable. In analogy to the bulk situation, the scaling theory is found to agree with mean-field theory in four dimensions. The prediction of the present work which is most easily accessible to experiment is the temperature dependence of the magnetization at the surface, with critical exponent estimated to be ${\ensuremath{\beta}}_{1}=\frac{2}{3}$. The mean-field result, ${\ensuremath{\beta}}_{1}=1$, seems to agree more closely with presently available experiment, and more work is needed to clarify the situation.
We investigate the dynamical phases and phase transitions arising in a classical two-dimensional anisotropic XY model under the influence of a periodically driven temporal external magnetic field in the form of a symmetric square wave. We use a combination of finite temperature classical Monte Carlo simulation, implemented within a CPU+GPU paradigm, utilizing local dynamics provided by the Glauber algorithm and a phenomenological equation-of-motion approach based on relaxational dynamics governed by the time-dependent free energy within a mean-field approximation to study the model. We investigate several parameter regimes of the variables (magnetic field, anisotropy, and the external drive frequency) that influence the anisotropic XY system. We identify four possible dynamical phases: Ising-SBO, Ising-SRO, XY-SBO, and XY-SRO. Both techniques indicate that only three of them (Ising-SRO, Ising-SBO, and XY-SRO) are stable dynamical phases in the thermodynamic sense. Within the Monte Carlo framework, a finite-size scaling analysis, shows that XY-SBO does not survive in the thermodynamic limit giving way to either an Ising-SBO or a XY-SRO regime. The finite-size scaling analysis further shows that the transitions between the three remaining dynamical phases either belong to the two-dimensional Ising universality class or are first-order in nature. Within the mean-field calculations yield three stable dynamical phases, i.e., Ising-SRO, Ising-SBO and XY-SRO, where the final steady state is independent of the initial condition chosen to evolve the equations of motion, as well as a region of bistability where the system flows to either Ising-SBO or XY-SRO (Ising-SRO) depending on the initial condition. Unlike the stable dynamical phases, the XY-SBO represents a transient feature that is eventually lost to either Ising-SBO or XY-SRO. Our mean-field analysis highlights the importance of the competition between switching of the stationary point(s) of the free energy after each half cycle of the external field and the two-dimensional nature of the phase space for the equations of motion.
We investigate the nonequilibrium relaxation dynamics of an interacting monomer-dimer model with nearest neighbor repulsion on a square lattice, which possesses two symmetric absorbing states. The model is known to exhibit two nearby continuous transitions: the Z 2 symmetry-breaking order-disorder transition and the absorbing transition with directed percolation criticality. We performed a more detailed analysis of our extensive simulations on bigger lattice systems which reaffirms that the symmetry-breaking transition exhibits a non-Ising critical behavior with β ≃ 0.149(2) and η ≃ 0.30(1) that are distinct from those values of a pure two dimensional Ising model. Finite size scaling of dimer density near the symmetry breaking transition gives logarithmic scaling ( α = 0.0) which is consistent with the hyperscaling relation but the corresponding exponent of ν B ≃ 1.37(2) exhibits a conspicuous deviation from the pure Ising value of 1. The value of dynamic critical exponent z , however, is found to be close to that of the kinetic Ising model as 1/ z ≃ 0.466(5) from the relaxation of staggered magnetization (and also similar but slightly smaller values from coarsening).
The stability of a discrete time crystal against thermal fluctuations has been studied numerically by solving a stochastic Landau-Lifshitz-Gilbert equation of a periodically driven classical system composed of interacting spins, each of which couples to a thermal bath. It is shown that in the thermodynamic limit, even though the long-range temporal crystalline order is stable at low temperature, it is melting above a critical temperature, at which the system experiences a nonequilibrium phase transition. The critical behaviors of the continuous phase transition have been systematically investigated, and it is shown that despite the genuine nonequilibrium feature of such a periodically driven system, its critical properties fall into the three-dimensional Ising universality class with a dynamical exponent (z=2z=2z=2) identical to that in the critical dynamics of a kinetic Ising model without driving.