Ising Model with Power Law Resetting

This paper investigates the nonequilibrium dynamics of the Ising model under power-law stochastic resetting, revealing that heavy-tailed reset times generate distinct quasi-ferromagnetic phases and rich phase diagrams in the $(T, \alpha)$ plane that differ significantly from those observed with exponential resetting.

The Gist

Imagine you are trying to find your way out of a giant, confusing maze. This is a bit like how particles in a physical system (like a magnet) try to find their "comfort zone" or equilibrium.

Usually, if you just let a system run on its own, it eventually settles down. But what happens if you keep interrupting it? What if, every now and then, you suddenly grab the system, throw it back to the starting line, and say, "Start over!"?

This is the concept of Stochastic Resetting. It's like a GPS that keeps resetting your route to the original destination every few minutes. Scientists have known for a while that doing this with a regular rhythm (like a clock ticking every 10 seconds) changes how the system behaves.

But this new paper asks a different question: What if the interruptions happen at random, unpredictable times, following a "Power Law"?

The "Power Law" Analogy: The Unpredictable Boss

Think of the difference between a strict boss and a chaotic one:

• Exponential Resetting (The Strict Boss): Your boss calls you back to the office every 10 minutes, no matter what. It's predictable. You know exactly when the reset will happen.
• Power Law Resetting (The Chaotic Boss): Your boss calls you back at random times. Sometimes they call after 1 minute. Sometimes they call after 1 hour. But here's the twist: sometimes they might not call for 10 years.

In a "Power Law" distribution, short breaks are common, but long, rare breaks are possible. This creates a "heavy tail" of time. The system gets to relax for a very long time, then suddenly gets yanked back to the start.

The Experiment: The Magnet (Ising Model)

The authors studied a classic model called the Ising Model. Imagine a grid of tiny magnets (spins) that can point Up or Down.

• Goal: They want to see how these magnets align (get a "magnetization") when they are constantly being reset to a specific starting pattern.
• The Variable: They changed the "temperature" (how chaotic the magnets are) and the "reset style" (how long the random breaks are).

The Surprising Results

The paper found that this "Chaotic Boss" style of resetting creates weird, new states of matter that you never see in normal physics or with the "Strict Boss" (exponential) resetting.

1. The "Quasi-Ferro" State (When it's Hot)

Normally, if a magnet is hot (above a critical temperature), the tiny magnets are chaotic and point in random directions. The average magnetism is zero.

• With Power Law Resetting: Even though it's hot, the system gets stuck in a weird "Quasi-Ferro" state.
• The Analogy: Imagine a crowd of people trying to decide on a dance move. Usually, they just flail around randomly. But because the "Chaotic Boss" sometimes leaves them alone for years, they manage to get into a rhythm. Then, the boss suddenly yells "Start over!" and they jump back to the original move.
• The Result: The crowd ends up spending a lot of time in two places: either flailing randomly (zero magnetism) or doing the original move (initial magnetism). They never settle in the middle. The paper calls this a Double-Peaked state because the probability of finding the system in these two states is huge.

2. The "Single vs. Double Peak" Switch (When it's Cold)

When the magnet is cold, it naturally wants to align (all Up or all Down).

• The Crossover: The authors found a magical switch point (called $\alpha^*$).
  • If the breaks are "short" enough (Low $\alpha$): The system behaves like a normal magnet. It settles into its natural cold state (Single Peak).
  • If the breaks are "long" enough (High $\alpha$): The system gets confused. It spends so much time relaxing that it almost reaches the natural state, but then the "Chaotic Boss" yells "Start over!" just as it's about to finish.
  • The Result: The system gets stuck in a tug-of-war. It spends time in its natural cold state AND the starting state. This creates a Double-Peaked state where the system is split between two identities.

Why Does This Matter?

This isn't just about magnets. The "Power Law" pattern of interruptions is everywhere in real life:

• Earthquakes: They happen in bursts with long quiet periods.
• Stock Markets: Prices jump wildly after long periods of calm.
• Human Behavior: We work for a while, then take a break, then work for a very long time, then take a vacation.

The paper shows that when you have a system that interacts with itself (like magnets, or people in a crowd), and you interrupt it with these "bursty" patterns, the system doesn't just behave "normally." It creates entirely new phases of matter that are a mix of order and chaos.

The Takeaway

If you treat a complex system like a clockwork machine (regular resets), it behaves predictably. But if you treat it like a real-life system with unpredictable, "bursty" interruptions (Power Law resets), it develops a rich, complex personality. It can be stuck in two places at once, or refuse to settle down, creating a "Quasi-Ferro" state that defies the usual rules of physics.

In short: Unpredictable interruptions don't just delay progress; they fundamentally change the destination.

Technical Summary

1. Problem Statement

The paper investigates the nonequilibrium dynamics of the nearest-neighbour Ising model subjected to stochastic resetting. While previous studies (e.g., Magoni et al.) have analyzed this system with exponentially distributed reset times (Poissonian process), this work explores the effects of power-law distributed inter-reset times.

• The Core Question: How does the heavy-tailed nature of power-law resetting ($\rho(\tau) \propto \tau^{-(\alpha+1)}$) alter the phase behavior, magnetization distributions, and steady-state properties of the Ising model compared to exponential resetting or standard equilibrium dynamics?
• Context: Power-law waiting times are ubiquitous in natural and human-made systems (earthquakes, neuron firing, financial markets) and introduce "heavy tails," allowing for rare but arbitrarily long intervals between resets. This creates a complex interplay between the system's intrinsic relaxation and the resetting mechanism.

2. Methodology

The authors employ a combination of analytical renewal theory and numerical simulations for both 1D and 2D Ising models.

• Model Setup:

  • System: $d$-dimensional Ising model with $N$ sites and periodic boundary conditions.
  • Dynamics: Glauber dynamics governs the evolution between resets.
  • Reset Protocol: The system is reset to an initial configuration with magnetization $m_0$ at random times $\tau$ drawn from a power-law distribution:
    $$\rho(\tau) = \frac{\alpha \tau_0^\alpha}{\tau^{\alpha+1}}, \quad \tau \in [\tau_0, \infty)$$
  • Observables: The probability distribution of the magnetization $P_r(m, t)$ at time $t$.

• Analytical Framework:

  • The authors use a renewal approach. The distribution at time $t$ is determined by the evolution following the last reset prior to $t$.
  • Renewal Equation:
    $$P_r(m, t) = \int_0^t d\tau f(t, t-\tau) P_0(m, \tau)$$
    where $P_0(m, \tau)$ is the magnetization distribution without resetting (deterministic evolution $\delta(m - m(\tau))$), and $f(t, t-\tau)$ is the probability density of the last reset occurring at $t-\tau$.
  • Regimes: The analysis distinguishes between $\alpha < 1$ (no steady state, time-dependent) and $\alpha > 1$ (steady state exists).
  • Approximations:
    • 1D: Exact analytical solution using $m(t) = m_0 e^{-(1-\gamma)t}$.
    • 2D ($T > T_c$): Exponential decay approximation $m(t) \approx m_0 e^{-\lambda t}$.
    • 2D ($T < T_c$): Stretched exponential relaxation $m(t) \approx m_{eq} \pm a e^{-b t^c}$.
    • 2D ($T = T_c$): Critical power-law decay $m(t) \approx b_c t^{-\phi}$.

• Numerical Validation:

  • Monte Carlo simulations were performed on lattices (e.g., $256 \times 256$ for 2D) to validate the analytical predictions across different $\alpha$ and temperature regimes.

3. Key Contributions

1. Discovery of New Nonequilibrium Phases: The study reveals that power-law resetting generates magnetization distributions fundamentally different from those produced by exponential resetting. Specifically, it introduces "quasi-ferro" and "double-peak ferro" states that do not exist in the exponential case.
2. Crossover Exponent ($\alpha^*$): For the 2D ferromagnetic phase ($T < T_c$), the authors identify a critical crossover exponent $\alpha^* = 1 - c$ (where $c$ is the Glauber dynamical exponent). This separates a single-peak phase from a double-peak phase.
3. Analytical Expressions for Heavy-Tailed Resetting: The paper derives explicit closed-form expressions for the magnetization probability density functions (PDFs) in 1D and 2D for both $\alpha < 1$ and $\alpha > 1$, incorporating the specific relaxation dynamics of the Ising model.
4. Phase Diagram Construction: A comprehensive phase diagram in the $(T, \alpha)$ plane is constructed, mapping out the transitions between paramagnetic, quasi-ferro, single-peak ferro, and double-peak ferro states.

4. Key Results

A. One-Dimensional Ising Model

• Equilibrium: Paramagnetic ($m_{eq} = 0$) at all finite temperatures.
• With Power-Law Resetting:
  • $\alpha < 1$: No steady state. The distribution $P_r(m, t)$ is time-dependent and exhibits a double-peaked structure with divergences at $m=0$ and $m=m_0$.
  • $\alpha > 1$: A stationary state emerges. The distribution remains double-peaked with divergences at $m=0$ and $m=m_0$.
  • Classification: This is termed a "Quasi-Ferro" (QF) state. Despite the underlying paramagnetic nature, the competition between short bursts of resetting (favoring $m_0$) and long relaxation tails (favoring $m=0$) creates a finite average magnetization.

B. Two-Dimensional Ising Model ($T > T_c$)

• Equilibrium: Paramagnetic ($m_{eq} = 0$).
• With Power-Law Resetting:
  • Similar to 1D, the system enters a Quasi-Ferro (QF) state for all $\alpha$.
  • $\alpha < 1$: Time-dependent double-peak distribution diverging at $m=0$ and $m=m_0$.
  • $\alpha > 1$: Stationary double-peak distribution with the same divergences.
  • Significance: This differs from the "pseudo-ferro" phase seen in exponential resetting (which has a gap at $m=0$). Here, the distribution diverges at $m=0$.

C. Two-Dimensional Ising Model ($T < T_c$)

• Equilibrium: Ferromagnetic ($m_{eq} > 0$).
• With Power-Law Resetting: Two distinct regimes separated by $\alpha^* = 1 - c$:
  1. Single-Peak Ferro (SPF) for $\alpha < \alpha^*$:
     • The distribution peaks only at the equilibrium magnetization $m_{eq}$.
     • The resetting is not frequent enough (or the tail is too heavy) to sustain a peak at $m_0$.
  2. Double-Peak Ferro (DPF) for $\alpha > \alpha^*$:
     • The distribution develops two peaks (divergences) at both $m_{eq}$ and $m_0$.
     • $\alpha < 1$: Time-dependent DPF.
     • $\alpha > 1$: Stationary DPF.
• Significance: This crossover is unique to power-law resetting. In exponential resetting, the system typically remains in a single ferro phase or transitions to a pseudo-ferro phase without this specific $\alpha$-dependent structural change in the ferromagnetic regime.

D. Critical Temperature ($T = T_c$)

• The distribution exhibits distinct power-law regimes determined by the interplay of the Ising critical exponent $\phi$ and the reset exponent $\alpha$.
• For $\alpha > 1$, a crossover occurs at $\alpha = \phi + 1$, shifting the distribution from decaying to growing power laws.

5. Significance and Implications

• Universality of Resetting Statistics: The work demonstrates that the statistics of the resetting process (exponential vs. power-law) are as critical as the system dynamics themselves in determining the nonequilibrium phase. Heavy tails introduce memory effects that prevent the system from settling into simple steady states or create new phases with divergent probabilities.
• Control of Nonequilibrium States: The findings suggest that power-law resetting can be used as a tool to "engineer" specific nonequilibrium phases (e.g., inducing a quasi-ferro state in a paramagnet or creating a double-peak ferromagnet).
• Beyond Ising: The authors argue that these features (heavy-tailed induced phases, crossover exponents) are likely generic to interacting many-body systems, offering a new framework for understanding complex systems with bursty dynamics.
• Theoretical Advancement: It bridges the gap between renewal theory and statistical mechanics of interacting spin systems, providing a rigorous analytical framework for non-Markovian resetting protocols.

In conclusion, the paper establishes that power-law resetting fundamentally rewrites the phase diagram of the Ising model, creating a rich landscape of non-equilibrium phases characterized by divergent magnetization distributions and novel crossover behaviors that are absent in both equilibrium and exponential-reset scenarios.