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:cond-mat/9503073 · Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics

Method to study relaxation of metastable phases: Macroscopic mean-field dynamics

Jooyoung Lee, M. A. Novotny, Per Arne Rikvold

Published Jul 1, 1995• 43 citations
arXivHTMLPDF

Abstract

We propose two different macroscopic dynamics to describe the decay of metastable phases in many-particle systems with local interactions. These dynamics depend on the macroscopic order parameter m through the restricted free energy F(m) and are designed to give the correct equilibrium distribution for m. The connection between macroscopic dynamics and the underlying microscopic dynamic is considered in the context of a projection-operator formalism. Application to the square-lattice nearest-neighbor Ising ferromagnet gives good agreement with droplet theory and Monte Carlo simulations of the underlying microscopic dynamic. This includes quantitative agreement for the exponential dependence of the lifetime 〈\ensuremath{\tau}〉 on the inverse of the applied field H, and the observation of distinct field regions in which \ensuremath{\Lambda}\ensuremath{\equiv}d ln〈\ensuremath{\tau}〉/d\ensuremath{\Vert}H${\mathrm{\ensuremath{\Vert}}}^{1\mathrm{\ensuremath{-}}\mathit{d}}$ depends differently on \ensuremath{\Vert}H\ensuremath{\Vert}. In addition, at very low temperatures we observe oscillatory behavior of \ensuremath{\Lambda} with respect to \ensuremath{\Vert}H\ensuremath{\Vert}, which is due to the discreteness of the lattice and in agreement with rigorous results. Similarities and differences between this work and earlier works on finite Ising models in the fixed-magnetization ensemble are discussed.

DOI
10.1103/physreve.52.356

Related papers

Critical phenomena: General considerations. Mean-field theory (MFT)

Jean Zinn‐Justin

Abstract This chapter is devoted to a brief review of general properties of phase transitions in macroscopic physics and, in particular in lattice models. Some of these lattice models actually appear as lattice regularizations of Euclidean (imaginary time) quantum physics theory (QFT). Most of the transitions considered in this work have the following character: spins on the lattice, or macroscopic particles in the continuum, interact through short-range forces, assumed, for simplicity, to decay exponentially. For simple systems, it is possible to find a local observable, called order parameter, whose expectation values depend on the phase in the several phase region, for example, the spin in ferromagnetic systems. In the disordered phase, the connected two-point function decreases exponentially at large distance, at a rate characterized by the correlation length (the inverse of the smallest physical mass in particle physics). In continuous transitions, the correlation length diverges at the critical temperature. Within the mean-field approximation (consistent with Landau's theory of critical phenomena), it can be shown that the singular behaviour of thermodynamic quantities at the critical temperature is universal. These properties can also be reproduced by calculating correlation functions with a perturbed Gaussian measure. It is then shown that the leading corrections to the mean-field approximation, in Ising-like systems, diverge at the critical temperature for dimensions smaller than or equal to $4$.

•Apr 15, 2021•0 citations

Dynamics of three-dimensional Ising spin glasses in thermal equilibrium

Andrew T. Ogielski

I present an analysis of the dynamic behavior of short-range Ising spin glasses observed in stochastic simulations. The time dependence of the order parameter q(t)=〈${S}_{x}$(0)${S}_{x}$(t)〉\ifmmode\bar\else\textasciimacron\fi{}---which is the same as that of the structure factor---and the time dependence of the related dynamic correlation functions have been recorded with good statistics and very long observation times. The spin-glass model with a symmetric distribution of discrete nearest-neighbor \ifmmode\pm\else\textpm\fi{}J interactions on a simple-cubic lattice was used. Simulations were performed with a special fast computer, allowing for the first-time investigation of the equilibrium dynamics for a wide range of temperatures (0.7\ensuremath{\le}kT/J\ensuremath{\le}5.0) and lattice sizes (${8}^{3}$, ${16}^{3}$, ${32}^{3}$, and ${64}^{3}$). I have found that the empirical formula q(t)=${\mathrm{ct}}^{\mathrm{\ensuremath{-}}x}$exp(-\ensuremath{\omega}${t}^{\ensuremath{\beta}}$) with temperature-dependent exponents x(T) and \ensuremath{\beta}(T) describes the decay very well at all temperatures above the spin-glass transition. In the spin-glass phase, only the algebraic decay q(t)=${\mathrm{ct}}^{\mathrm{\ensuremath{-}}x}$ could be observed, with different temperature dependences of the exponent x(T). The dynamic scaling hypothesis and finite-size scaling explain well the observed temperature and size dependence of the data, and the functional form of the correlation functions is com- patible with the scaling form if corrections to scaling are taken into account. The scaling behavior and the dynamic and static critical exponents found in my simulations are in reasonable agreement with recent experiments performed on insulating spin glasses, showing that despite its simplicity the discrete model of spin glasses analyzed in this work displays behavior similar to that seen in nature.

Physical review. B, Condensed matter•Dec 1, 1985•773 citations

Global Speed Limit for Finite-Time Dynamical Phase Transition and Nonequilibrium Relaxation

Kristian Blom, Aljaž Godec

Recent works unraveled an intriguing finite-time dynamical phase transition in the thermal relaxation of the mean field Curie-Weiss model. The phase transition reflects a sudden switch in the dynamics. Its existence in systems with a finite range of interaction, however, remained unclear. Here we demonstrate the dynamical phase transition for nearest-neighbor Ising systems on the square and Bethe lattices through extensive computer simulations and by analytical results. Combining large-deviation techniques and Bethe-Guggenheim theory we prove the existence of the dynamical phase transition for arbitrary quenches, including those within the two-phase region. Strikingly, for any given initial condition we prove and explain the existence of non-trivial speed limits for the dynamical phase transition and the relaxation of magnetization, which are fully corroborated by simulations of the microscopic Ising model but are absent in the mean field setting. Pair correlations, which are neglected in mean field theory and trivial in the Curie-Weiss model, account for kinetic constraints due to frustrated local configurations that give rise to a global speed limit.

arXiv:2209.14287•Sep 28, 2022•1 citation

Critical fluctuations at a finite-time dynamical phase transition

Nalina Vadakkayil, Massimiliano Esposito, Jan Meibohm

We explore the critical properties of the recently discovered finite-time dynamical phase transition in the nonequilibrium relaxation of Ising magnets after a temperature quench. The transition is characterized by a sudden switch in the relaxation dynamics and it occurs at a sharp critical time. While previous works have focused either on mean-field interactions or on investigating the critical time, we analyze the critical fluctuations at the phase transition in the nearest-neighbor Ising model on a square lattice using Monte Carlo simulations. By means of a finite-size scaling analysis, we extract the critical exponents for the transition. In two spatial dimensions, we find that the exponents are consistent with those of the two-dimensional Ising universality class when the system is initially in the vicinity of the critical point. For initial temperatures below the critical one, however, the critical exponents differ from the Ising-exponents, indicating a distinct dynamical critical phenomenon.

arXiv:2403.10505•Physical review. E•Dec 30, 2024•6 citations

Macroscopic Fokker-Planck equation from microscopic Glauber dynamics for the Nagle-Kardar model

Jean-François de Kemmeter, Stefano Ruffo, Stefano Gherardini

In the companion Letter we have highlighted the dynamical universality class of the Nagle-Kardar model, where a mean-field interaction is added to the one-dimensional nearest-neighbor Ising model. Starting from the microscopic Glauber dynamics, this paper provides a complete derivation of the Fokker-Planck equation that describes the time evolution of the macroscopic variables (magnetization and defect density) appearing in the model's Hamiltonian. The study of the Langevin equation, associated with the Fokker-Planck equation, allowed us to prove that the model belongs to the universal class of systems with diffusive dynamics and non-conserved order parameter (model A). To this goal, we used several features of the model at equilibrium, including the phase diagram and the fluctuations of macroscopic variables, which are here discussed for completeness. The derivation of the Fokker-Planck equation requires the solution of some combinatorial problems that appear in the counting of configurations at an intermediate level between the microscopic Glauber dynamics and the macroscopic Fokker-Planck one. Finally, we apply both the Glauber and the Langevin dynamics to the study of the average first passage time between local equilibrium states. We confirm that this time obeys an exponential Arrhenius law in terms of system's size, offering a direct link between microscopic energy landscapes and macroscopic relaxation mechanisms.

arXiv:2606.25919•Jun 24, 2026•0 citations

Theory for the dynamics of clusters near the critical point. I. Relaxation of the Glauber kinetic Ising model

Kurt Binder, Dietrich Stauffer, Heiner Muller-Krumbhaar

A semiphenomenological cluster theory is developed for dynamic critical properties, which is not limited to small deviations from equilibrium. Explicit numerical expressions are derived for linear and nonlinear response functions of the kinetic Ising model, which are compatible with dynamic scaling and recent calculations of dynamic critical exponents. An uniaxial magnet (or Ising system) and its critical fluctuations are approximated by a "droplet" or "cluster" model, the critical part of the free energy being the sum of the contributions of clusters each containing $l$ spins. These clusters are assumed to grow and shrink with a phenomenological rate $\ensuremath{\propto}{l}^{r}$, $r\ensuremath{\approx}1$. These rates are compatible with direct Monte Carlo simulations in two and three dimensions. The resulting cluster reaction and diffusion equation is an approximation to the master equation of the Glauber kinetic Ising model (in which the magnetization is not conserved), and has the same structure as used in nucleation theories. In linear response to space- and time-dependent fields the relaxation times of energy and magnetization are expressed as triple integrals. All relaxation times diverge as ${|1\ensuremath{-}\frac{T}{{T}_{c}}|}^{\ensuremath{-}(2\ensuremath{-}r)\ensuremath{\beta}\ensuremath{\delta}}$. This treatment is consistent with dynamic scaling and universality, and is more general than the conventional (Van Hove) theory of critical slowing down. Also the (smaller) exponents found for the response to localized variations ("autocorrelation functions") agree with dynamic scaling. The wave-vector dependence of the relaxation times arises from the static correlation function only, and not from the cluster diffusion term. If the equilibrium cluster distribution is assumed to be that of the Fisher droplet model, and for special values of $r$, the (discrete) eigenvalue spectrum and eigenfunctions (generalized Laguerre polynomials) of the cluster reaction equation are found explicitly. Relaxation functions and frequency-dependent susceptibilities are expressed by hypergeometric series, even in the case of nonlinear response. At the critical point, the spectrum becomes continuous and the exponential decay for large times is replaced by power-law behavior. Similarly, it is found that dynamic scaling applies also to the nucleation rate for this model, and the scaling exponent of the nucleation rate is $2\ensuremath{-}\ensuremath{\alpha}+(2\ensuremath{-}r)\ensuremath{\beta}\ensuremath{\delta}$.

Physical review. B, Solid state•Dec 1, 1975•98 citations
We propose two different macroscopic dynamics to describe the decay of metastable phases in many-particle systems with local interactions. These dynamics depend on the macroscopic order parameter m through the restricted free energy F(m) and are designed to give the correct equilibrium distribution for m. The connection between macroscopic dynamics and the underlying microscopic dynamic is considered in the context of a projection-operator formalism. Application to the square-lattice nearest-neighbor Ising ferromagnet gives good agreement with droplet theory and Monte Carlo simulations of the underlying microscopic dynamic. This includes quantitative agreement for the exponential dependence of the lifetime 〈\ensuremath{\tau}〉 on the inverse of the applied field H, and the observation of distinct field regions in which \ensuremath{\Lambda}\ensuremath{\equiv}d ln〈\ensuremath{\tau}〉/d\ensuremath{\Vert}H\ensuremath∥1\ensuremath−d{\mathrm{\ensuremath{\Vert}}}^{1\mathrm{\ensuremath{-}}\mathit{d}}\ensuremath∥1\ensuremath−d depends differently on \ensuremath{\Vert}H\ensuremath{\Vert}. In addition, at very low temperatures we observe oscillatory behavior of \ensuremath{\Lambda} with respect to \ensuremath{\Vert}H\ensuremath{\Vert}, which is due to the discreteness of the lattice and in agreement with rigorous results. Similarities and differences between this work and earlier works on finite Ising models in the fixed-magnetization ensemble are discussed.
Abstract This chapter is devoted to a brief review of general properties of phase transitions in macroscopic physics and, in particular in lattice models. Some of these lattice models actually appear as lattice regularizations of Euclidean (imaginary time) quantum physics theory (QFT). Most of the transitions considered in this work have the following character: spins on the lattice, or macroscopic particles in the continuum, interact through short-range forces, assumed, for simplicity, to decay exponentially. For simple systems, it is possible to find a local observable, called order parameter, whose expectation values depend on the phase in the several phase region, for example, the spin in ferromagnetic systems. In the disordered phase, the connected two-point function decreases exponentially at large distance, at a rate characterized by the correlation length (the inverse of the smallest physical mass in particle physics). In continuous transitions, the correlation length diverges at the critical temperature. Within the mean-field approximation (consistent with Landau's theory of critical phenomena), it can be shown that the singular behaviour of thermodynamic quantities at the critical temperature is universal. These properties can also be reproduced by calculating correlation functions with a perturbed Gaussian measure. It is then shown that the leading corrections to the mean-field approximation, in Ising-like systems, diverge at the critical temperature for dimensions smaller than or equal to 444.
I present an analysis of the dynamic behavior of short-range Ising spin glasses observed in stochastic simulations. The time dependence of the order parameter q(t)=〈Sx{S}_{x}Sx​(0)Sx{S}_{x}Sx​(t)〉\ifmmode\bar\else\textasciimacron\fi{}---which is the same as that of the structure factor---and the time dependence of the related dynamic correlation functions have been recorded with good statistics and very long observation times. The spin-glass model with a symmetric distribution of discrete nearest-neighbor \ifmmode\pm\else\textpm\fi{}J interactions on a simple-cubic lattice was used. Simulations were performed with a special fast computer, allowing for the first-time investigation of the equilibrium dynamics for a wide range of temperatures (0.7\ensuremath{\le}kT/J\ensuremath{\le}5.0) and lattice sizes (83{8}^{3}83, 163{16}^{3}163, 323{32}^{3}323, and 643{64}^{3}643). I have found that the empirical formula q(t)=${\mathrm{ct}}^{\mathrm{\ensuremath{-}}x}$exp(-\ensuremath{\omega}${t}^{\ensuremath{\beta}}$) with temperature-dependent exponents x(T) and \ensuremath{\beta}(T) describes the decay very well at all temperatures above the spin-glass transition. In the spin-glass phase, only the algebraic decay q(t)=${\mathrm{ct}}^{\mathrm{\ensuremath{-}}x}$ could be observed, with different temperature dependences of the exponent x(T). The dynamic scaling hypothesis and finite-size scaling explain well the observed temperature and size dependence of the data, and the functional form of the correlation functions is com- patible with the scaling form if corrections to scaling are taken into account. The scaling behavior and the dynamic and static critical exponents found in my simulations are in reasonable agreement with recent experiments performed on insulating spin glasses, showing that despite its simplicity the discrete model of spin glasses analyzed in this work displays behavior similar to that seen in nature.
Recent works unraveled an intriguing finite-time dynamical phase transition in the thermal relaxation of the mean field Curie-Weiss model. The phase transition reflects a sudden switch in the dynamics. Its existence in systems with a finite range of interaction, however, remained unclear. Here we demonstrate the dynamical phase transition for nearest-neighbor Ising systems on the square and Bethe lattices through extensive computer simulations and by analytical results. Combining large-deviation techniques and Bethe-Guggenheim theory we prove the existence of the dynamical phase transition for arbitrary quenches, including those within the two-phase region. Strikingly, for any given initial condition we prove and explain the existence of non-trivial speed limits for the dynamical phase transition and the relaxation of magnetization, which are fully corroborated by simulations of the microscopic Ising model but are absent in the mean field setting. Pair correlations, which are neglected in mean field theory and trivial in the Curie-Weiss model, account for kinetic constraints due to frustrated local configurations that give rise to a global speed limit.
We explore the critical properties of the recently discovered finite-time dynamical phase transition in the nonequilibrium relaxation of Ising magnets after a temperature quench. The transition is characterized by a sudden switch in the relaxation dynamics and it occurs at a sharp critical time. While previous works have focused either on mean-field interactions or on investigating the critical time, we analyze the critical fluctuations at the phase transition in the nearest-neighbor Ising model on a square lattice using Monte Carlo simulations. By means of a finite-size scaling analysis, we extract the critical exponents for the transition. In two spatial dimensions, we find that the exponents are consistent with those of the two-dimensional Ising universality class when the system is initially in the vicinity of the critical point. For initial temperatures below the critical one, however, the critical exponents differ from the Ising-exponents, indicating a distinct dynamical critical phenomenon.
In the companion Letter we have highlighted the dynamical universality class of the Nagle-Kardar model, where a mean-field interaction is added to the one-dimensional nearest-neighbor Ising model. Starting from the microscopic Glauber dynamics, this paper provides a complete derivation of the Fokker-Planck equation that describes the time evolution of the macroscopic variables (magnetization and defect density) appearing in the model's Hamiltonian. The study of the Langevin equation, associated with the Fokker-Planck equation, allowed us to prove that the model belongs to the universal class of systems with diffusive dynamics and non-conserved order parameter (model A). To this goal, we used several features of the model at equilibrium, including the phase diagram and the fluctuations of macroscopic variables, which are here discussed for completeness. The derivation of the Fokker-Planck equation requires the solution of some combinatorial problems that appear in the counting of configurations at an intermediate level between the microscopic Glauber dynamics and the macroscopic Fokker-Planck one. Finally, we apply both the Glauber and the Langevin dynamics to the study of the average first passage time between local equilibrium states. We confirm that this time obeys an exponential Arrhenius law in terms of system's size, offering a direct link between microscopic energy landscapes and macroscopic relaxation mechanisms.
A semiphenomenological cluster theory is developed for dynamic critical properties, which is not limited to small deviations from equilibrium. Explicit numerical expressions are derived for linear and nonlinear response functions of the kinetic Ising model, which are compatible with dynamic scaling and recent calculations of dynamic critical exponents. An uniaxial magnet (or Ising system) and its critical fluctuations are approximated by a "droplet" or "cluster" model, the critical part of the free energy being the sum of the contributions of clusters each containing lll spins. These clusters are assumed to grow and shrink with a phenomenological rate $\ensuremath{\propto}{l}^{r}$, $r\ensuremath{\approx}1$. These rates are compatible with direct Monte Carlo simulations in two and three dimensions. The resulting cluster reaction and diffusion equation is an approximation to the master equation of the Glauber kinetic Ising model (in which the magnetization is not conserved), and has the same structure as used in nucleation theories. In linear response to space- and time-dependent fields the relaxation times of energy and magnetization are expressed as triple integrals. All relaxation times diverge as ${|1\ensuremath{-}\frac{T}{{T}_{c}}|}^{\ensuremath{-}(2\ensuremath{-}r)\ensuremath{\beta}\ensuremath{\delta}}$. This treatment is consistent with dynamic scaling and universality, and is more general than the conventional (Van Hove) theory of critical slowing down. Also the (smaller) exponents found for the response to localized variations ("autocorrelation functions") agree with dynamic scaling. The wave-vector dependence of the relaxation times arises from the static correlation function only, and not from the cluster diffusion term. If the equilibrium cluster distribution is assumed to be that of the Fisher droplet model, and for special values of rrr, the (discrete) eigenvalue spectrum and eigenfunctions (generalized Laguerre polynomials) of the cluster reaction equation are found explicitly. Relaxation functions and frequency-dependent susceptibilities are expressed by hypergeometric series, even in the case of nonlinear response. At the critical point, the spectrum becomes continuous and the exponential decay for large times is replaced by power-law behavior. Similarly, it is found that dynamic scaling applies also to the nucleation rate for this model, and the scaling exponent of the nucleation rate is $2\ensuremath{-}\ensuremath{\alpha}+(2\ensuremath{-}r)\ensuremath{\beta}\ensuremath{\delta}$.