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

Physical review. B, Condensed matter

Nonlinear-response theory of Ising-like systems

James H. Luscombe

Published May 1, 1984• 13 citations
Publisher

Abstract

The nonequilibrium response of a two-dimensional, paramagnetic system to a rapid and arbitrarily large change (quench) in a symmetry-breaking field is investigated. The complete time dependence of the conjugate order parameter is calculated with a real-space renormalization-group (RG) method. The nonlinear relaxation "rate" is found to be time dependent. Its approach to the corresponding prediction of linear-response theory is characterized. The nonlinear response is found to go over to linear response via an early, exponential regime followed by a crossover to an asymptotic, algebraic approach to linear response at long times. The initial transient decay is quench-strength dependent in contrast with the final behavior, which is largely "universal" and quench insensitive. We study a nearest-neighbor Ising model in a magnetic field in the paramagnetic, "disordered" phase. The dynamics are controlled by a kinetic Ising model with simultaneous spin-flip (adsorption-desorption) and spin-exchange (diffusion) kinetics. The model is applicable to non-equilibrium "dosing" experiments on adsorbed systems. Diffusion is shown not to affect the growth rate of the total order parameter (magnetization or coverage). General recursion relations for the temperature and magnetic field are derived for the entire thermodynamic plane. The RG analysis reproduces linear-response theory for arbitrary magnetic field and is consistent with known scaling laws concerning nonlinear relaxation. The time-dependent magnetization in the critical region has a power-law regime before crossing over into its final exponential decay.

DOI
10.1103/physrevb.29.5128

Related papers

Nonequilibrium dynamics of random field Ising spin chains: Exact results via real space renormalization group

Daniel S. Fisher, Pierre Le Doussal, Cécile Monthus

The nonequilibrium dynamics of classical random Ising spin chains with nonconserved magnetization are studied using an asymptotically exact real space renormalization group (RSRG). We focus on random field Ising model (RFIM) spin chains with and without a uniform applied field, as well as on Ising spin glass chains in an applied field. For the RFIM we consider a universal regime where the random field and the temperature are both much smaller than the exchange coupling. In this regime, the Imry-Ma length that sets the scale of the equilibrium correlations is large and the coarsening of domains from random initial conditions (e.g., a quench from high temperature) occurs over a wide range of length scales. The two types of domain walls that occur diffuse in opposite random potentials, of the form studied by Sinai, and domain walls annihilate when they meet. Using the RSRG we compute many universal asymptotic properties of both the nonequilibrium dynamics and the equilibrium limit. We find that the configurations of the domain walls converge rapidly toward a set of system-specific time-dependent positions that are independent of the initial conditions. Thus the behavior of this nonequilibrium system is pseudodeterministic at long times because of the broad distributions of barriers that occur on the long length scales involved. Specifically, we obtain the time dependence of the energy, the magnetization, and the distribution of domain sizes (found to be statistically independent). The equilibrium limits agree with known exact results. We obtain the exact scaling form of the two-point equal time correlation function <S0(t)Sx(t)> and the two-time autocorrelations <S0(t('))S0(t)>. We also compute the persistence properties of a single spin, of local magnetization, and of domains. The analogous quantities for the +/-J Ising spin glass in an applied field are obtained from the RFIM via a gauge transformation. In addition to these we compute the two-point two-time correlation function <S0(t)Sx(t)><S0(t('))Sx(t('))> which can in principle be measured by experiments on spin-glass-like systems. The thermal fluctuations are studied and found to be dominated by rare events; in particular all moments of truncated equal time correlations are computed. Physical properties which are typically measured in aging experiments are also studied, focusing on the response to a small magnetic field which is applied after waiting for the system to equilibrate for a time t(w). The nonequilibrium fluctuation-dissipation ratio X(t,t(w)) is computed. We find that for (t-t(w)) approximately t(alpha)(w) with alpha;<1, the ratio equal to its equilibrium value X=1, although time translational invariance does not hold in this regime. For t-t(w) approximately t(w) the ratio exhibits an aging regime with a nontrivial X=X(t/t(w)) not equal to 1, but the behavior is markedly different from mean field theory. Finally the distribution of the total magnetization and of the number of domains is computed for large finite size systems. General issues about convergence toward equilibrium and the possibilities of weakly history-dependent evolution in other random systems are discussed.

arXiv:cond-mat/0012290•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Nov 14, 2001•77 citations

Temperature quenches in the disordered phase of a spin system

Gene F. Mazenko

The problem of a sudden temperature quench of a spin system is studied for the case of a two-dimensional spin-flip kinetic Ising model in zero external field. The analysis is carried out for the case where both the initial and final temperatures are in the disordered phase. The problem is analyzed with the use of a generalization of the real-space dynamic renormalization-group method. The interesting features found in studying the time-dependent structure factor are the following: (i) a nonmonotonic time dependence at finite wave numbers corresponding to "structure pulses"; (ii) a slow shift and decay in the peak of the structure factor for fixed wave number as a function of temperature; (iii) a slow approach to the final exponential relaxation to equilibrium.

Physical review. B, Condensed matter•Nov 1, 1982•13 citations

Nonequilibrium Kondo model: Crossover from weak to strong coupling

Mikhail Pletyukhov, Herbert Schoeller

We analyze the nonequilibrium Kondo model at finite voltage and temperature by using a new formulation of the real-time renormalization group method with the Laplace variable as the flow parameter. We evaluate the energy-dependent spin relaxation rate and nonlinear conductance, and derive an approximate form for the universal line shape for the latter in the whole crossover regime from weak to strong coupling. The results are shown to agree well with exact methods in equilibrium, Fermi-liquid theory, weak-coupling expansions, and recent experiments. For the transient spin dynamics we find a universal exponential decay in the long-time limit along with a truncation-dependent pre-exponential power law. For multichannel models a pure power-law decay typical for non-Fermi-liquid behavior is predicted.

arXiv:1201.6295•Physical Review Letters•Jun 26, 2012•89 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

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

Critical behavior of the order parameter at the nonequilibrium phase transition of the Ising model

Bin Li, Chao Gao, Gao Xianlong +1 more

After a quench of transverse field, the asymptotic long-time state of the Ising model displays a transition from a ferromagnetic phase to a paramagnetic phase as the post-quench field strength increases, which is revealed by the vanishing of the order parameter defined as the averaged magnetization over time. We estimate the critical behavior of the magnetization at this nonequilibrium phase transition by using the mean-field approximation. In the vicinity of the critical field, the magnetization vanishes as the inverse of a logarithmic function, which is significantly distinguished from the critical behavior of order parameter at the corresponding equilibrium phase transition, i.e. a power-law function.

arXiv:1808.07395•Journal of Physics Condensed Matter•Dec 6, 2018•2 citations
The nonequilibrium response of a two-dimensional, paramagnetic system to a rapid and arbitrarily large change (quench) in a symmetry-breaking field is investigated. The complete time dependence of the conjugate order parameter is calculated with a real-space renormalization-group (RG) method. The nonlinear relaxation "rate" is found to be time dependent. Its approach to the corresponding prediction of linear-response theory is characterized. The nonlinear response is found to go over to linear response via an early, exponential regime followed by a crossover to an asymptotic, algebraic approach to linear response at long times. The initial transient decay is quench-strength dependent in contrast with the final behavior, which is largely "universal" and quench insensitive. We study a nearest-neighbor Ising model in a magnetic field in the paramagnetic, "disordered" phase. The dynamics are controlled by a kinetic Ising model with simultaneous spin-flip (adsorption-desorption) and spin-exchange (diffusion) kinetics. The model is applicable to non-equilibrium "dosing" experiments on adsorbed systems. Diffusion is shown not to affect the growth rate of the total order parameter (magnetization or coverage). General recursion relations for the temperature and magnetic field are derived for the entire thermodynamic plane. The RG analysis reproduces linear-response theory for arbitrary magnetic field and is consistent with known scaling laws concerning nonlinear relaxation. The time-dependent magnetization in the critical region has a power-law regime before crossing over into its final exponential decay.
The nonequilibrium dynamics of classical random Ising spin chains with nonconserved magnetization are studied using an asymptotically exact real space renormalization group (RSRG). We focus on random field Ising model (RFIM) spin chains with and without a uniform applied field, as well as on Ising spin glass chains in an applied field. For the RFIM we consider a universal regime where the random field and the temperature are both much smaller than the exchange coupling. In this regime, the Imry-Ma length that sets the scale of the equilibrium correlations is large and the coarsening of domains from random initial conditions (e.g., a quench from high temperature) occurs over a wide range of length scales. The two types of domain walls that occur diffuse in opposite random potentials, of the form studied by Sinai, and domain walls annihilate when they meet. Using the RSRG we compute many universal asymptotic properties of both the nonequilibrium dynamics and the equilibrium limit. We find that the configurations of the domain walls converge rapidly toward a set of system-specific time-dependent positions that are independent of the initial conditions. Thus the behavior of this nonequilibrium system is pseudodeterministic at long times because of the broad distributions of barriers that occur on the long length scales involved. Specifically, we obtain the time dependence of the energy, the magnetization, and the distribution of domain sizes (found to be statistically independent). The equilibrium limits agree with known exact results. We obtain the exact scaling form of the two-point equal time correlation function <S0(t)Sx(t)> and the two-time autocorrelations <S0(t('))S0(t)>. We also compute the persistence properties of a single spin, of local magnetization, and of domains. The analogous quantities for the +/-J Ising spin glass in an applied field are obtained from the RFIM via a gauge transformation. In addition to these we compute the two-point two-time correlation function <S0(t)Sx(t)><S0(t('))Sx(t('))> which can in principle be measured by experiments on spin-glass-like systems. The thermal fluctuations are studied and found to be dominated by rare events; in particular all moments of truncated equal time correlations are computed. Physical properties which are typically measured in aging experiments are also studied, focusing on the response to a small magnetic field which is applied after waiting for the system to equilibrate for a time t(w). The nonequilibrium fluctuation-dissipation ratio X(t,t(w)) is computed. We find that for (t-t(w)) approximately t(alpha)(w) with alpha;<1, the ratio equal to its equilibrium value X=1, although time translational invariance does not hold in this regime. For t-t(w) approximately t(w) the ratio exhibits an aging regime with a nontrivial X=X(t/t(w)) not equal to 1, but the behavior is markedly different from mean field theory. Finally the distribution of the total magnetization and of the number of domains is computed for large finite size systems. General issues about convergence toward equilibrium and the possibilities of weakly history-dependent evolution in other random systems are discussed.
The problem of a sudden temperature quench of a spin system is studied for the case of a two-dimensional spin-flip kinetic Ising model in zero external field. The analysis is carried out for the case where both the initial and final temperatures are in the disordered phase. The problem is analyzed with the use of a generalization of the real-space dynamic renormalization-group method. The interesting features found in studying the time-dependent structure factor are the following: (i) a nonmonotonic time dependence at finite wave numbers corresponding to "structure pulses"; (ii) a slow shift and decay in the peak of the structure factor for fixed wave number as a function of temperature; (iii) a slow approach to the final exponential relaxation to equilibrium.
We analyze the nonequilibrium Kondo model at finite voltage and temperature by using a new formulation of the real-time renormalization group method with the Laplace variable as the flow parameter. We evaluate the energy-dependent spin relaxation rate and nonlinear conductance, and derive an approximate form for the universal line shape for the latter in the whole crossover regime from weak to strong coupling. The results are shown to agree well with exact methods in equilibrium, Fermi-liquid theory, weak-coupling expansions, and recent experiments. For the transient spin dynamics we find a universal exponential decay in the long-time limit along with a truncation-dependent pre-exponential power law. For multichannel models a pure power-law decay typical for non-Fermi-liquid behavior is predicted.
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}$.
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.
After a quench of transverse field, the asymptotic long-time state of the Ising model displays a transition from a ferromagnetic phase to a paramagnetic phase as the post-quench field strength increases, which is revealed by the vanishing of the order parameter defined as the averaged magnetization over time. We estimate the critical behavior of the magnetization at this nonequilibrium phase transition by using the mean-field approximation. In the vicinity of the critical field, the magnetization vanishes as the inverse of a logarithmic function, which is significantly distinguished from the critical behavior of order parameter at the corresponding equilibrium phase transition, i.e. a power-law function.