Virtual walks in the Ising model: finite time scaling

This paper analyzes the non-equilibrium dynamics of the Ising model in one and two dimensions using a virtual walk approach associated with spin states and local energy, demonstrating that finite-time scaling of these walks yields critical exponents consistent with known theoretical values and reveals a distinct non-equilibrium region and time-dependent critical point.

The Gist

Imagine a giant, crowded dance floor where thousands of dancers (the "spins") are moving to music. This is the Ising Model, a famous mathematical way physicists study how things like magnets or fluids change state (like ice melting into water).

Usually, to understand how this dance floor behaves near a critical moment (the "phase transition" where order turns to chaos), physicists have to build many different-sized dance floors and watch them for a very long time. It's like trying to predict the weather by building a hundred different model cities.

This paper introduces a clever shortcut: "Virtual Walks."

Instead of just watching the dancers, the authors imagine that every single dancer is carrying a tiny, invisible walker on their shoulder. Here is how the story unfolds:

1. The Setup: The Great Quench

The researchers start with the dance floor in a state of pure chaos. The music is so loud and fast (high temperature) that everyone is spinning wildly in random directions.

• The "Quench": Suddenly, they turn the music down to a slow, steady beat (lower temperature).
• The Goal: They want to see how the dancers organize themselves. Do they form neat lines (magnetism)? Or do they stay chaotic?

2. The Virtual Walk: A New Perspective

Here is the magic trick. For every dancer, they imagine a "walker" moving on a straight line.

• If a dancer is facing Right (Spin = +1), their walker takes a step Forward.
• If a dancer is facing Left (Spin = -1), their walker takes a step Backward.

As time goes on, the walker's path is just the sum of all the steps the dancer took.

• In the Chaos (High Temp): The dancer flips back and forth randomly. The walker zig-zags like a drunkard. The path looks like a messy, wiggly line that stays near the center.
• In the Order (Low Temp): The dancer gets stuck facing one way for a long time. The walker marches steadily in one direction, creating a long, straight line away from the center.

3. The "Shape" of the Crowd

The authors looked at the distribution of where all these walkers ended up.

• Below the Critical Point (Order): The walkers are mostly far away from the center (either far left or far right). The graph looks like a camel with two humps.
• Above the Critical Point (Chaos): The walkers are clustered right in the middle. The graph looks like a single hill (a bell curve).

The Discovery: By watching how the shape of this graph changes from a "two-humped camel" to a "single hill," they can pinpoint the exact moment the system changes from order to chaos. This is much easier than the old methods!

4. The "Energy Walk" (The Second Story)

The authors didn't stop there. They created a second type of walker.

• Instead of tracking the dancer's direction, this walker tracks the energy of the dancer's relationship with their neighbors.
• If neighbors are getting along (same direction), the energy is low. If they are fighting (opposite directions), the energy is high.
• This "Energy Walker" also changes its behavior at the critical point, giving the researchers a second, independent way to confirm their findings.

5. Why This Matters: The "Time Machine"

The biggest breakthrough is Finite Time Scaling.

• The Old Way: To find the "Critical Temperature" (the exact moment of change), you usually need to simulate huge systems for a long time and compare different sizes. It's slow and computationally expensive.
• The New Way: The authors show that by using these virtual walks, you can figure out the critical temperature and the "rules of the game" (critical exponents) using just one system size and looking at how things evolve over time.

It's like being able to predict the entire history of a storm just by watching a single raindrop fall for a few seconds, rather than needing to map the whole sky.

Summary

This paper is about a new, clever way to watch a crowd of magnetic spins. By turning their chaotic movements into "virtual walks," the researchers found a simple visual signal (the shape of the walk distribution) that tells us exactly when the system is about to change its state. It's a faster, smarter way to understand the physics of change, using the power of "time" instead of just "size."

Technical Summary

Here is a detailed technical summary of the paper "Virtual walks in the Ising model: finite time scaling" by Amit Pradhan, Parongama Sen, and Sagnik Seth.

1. Problem Statement

The paper addresses the challenge of detecting critical points and estimating critical exponents in the Ising model, particularly in non-equilibrium scenarios following a temperature quench.

• Context: While equilibrium behavior of the Ising model is well-understood, non-equilibrium dynamics (specifically domain coarsening after a quench from high temperature) are complex.
• Limitation of Standard Methods: The standard approach, Finite Size Scaling (FSS), requires numerical simulations of multiple system sizes ($L$) evolved over long times to extract critical temperatures ($T_c$) and exponents. This is computationally expensive.
• Goal: The authors propose an alternative method using Virtual Walks combined with Finite Time Scaling (FTS). The objective is to extract the complete set of static and dynamic critical exponents and $T_c$ using data from a single system size evolved over time, thereby avoiding the computational cost of varying system sizes.

2. Methodology

The study employs the Glauber dynamics scheme to simulate the Ising model in one dimension ($d=1$) and two dimensions ($d=2$). The system is quenched from a random high-temperature state to a lower temperature $T$.

A. Virtual Walk Definitions

Two types of virtual walks are defined for each spin (or site) $i$:

1. Spin Walk ($x_i$):
   • The displacement is the cumulative sum of the spin state over time.
   • Equation: $x_i(t+1) = x_i(t) + S_i(t+1)$, where $S_i = \pm 1$.
   • This maps the spin history to a 1D random walk trajectory.
2. Energy Walk ($y_i$):
   • Introduced specifically for the 2D model.
   • The displacement is the cumulative sum of the local interaction energy.
   • Equation: $y_i(t+1) = y_i(t) + E_i(t+1)$, where $E_i(t) = -S_i(t) \sum_{j \in NN} S_j(t)$.
   • Step sizes vary between $\pm 4$.

B. Analysis Techniques

• Probability Distribution Analysis: The authors analyze the probability distribution $P(x, t)$ of the walker displacements. They observe a crossover from a bimodal (two-peaked) distribution (indicating ordered/ballistic behavior) to a unimodal Gaussian (disordered/diffusive behavior) as temperature increases or time evolves.
• Ratio Method: They define a ratio $r = P(x=0, t) / P(x=x_m, t)$, where $x_m$ is the peak position. This ratio helps identify the crossover time $t_c$ and the critical temperature $T_c$.
• Binder Cumulant: A modified Binder cumulant $U(T, t) = 1 - \langle x^4 \rangle / (3\langle x^2 \rangle^2)$ is constructed using walker displacements. The crossing point of curves for different times $t$ yields $T_c$.
• Finite Time Scaling (FTS): Instead of varying system size $L$, the authors scale time $t$. They assume quantities follow scaling forms like $Q(T, t) = t^{-\kappa} f((T-T_c)t^{1/\nu z_c})$, where $\nu$ is the correlation length exponent and $z_c$ is the dynamic exponent.

3. Key Contributions

1. Virtual Walk Framework: The paper rigorously extends the concept of virtual walks (previously used in econophysics and opinion dynamics) to the Ising model, establishing a direct link between spin dynamics and random walk statistics.
2. Single-System Criticality Detection: The authors demonstrate that $T_c$ and critical exponents can be accurately determined using Finite Time Scaling on a single system size, bypassing the need for the computationally intensive Finite Size Scaling method.
3. Introduction of Energy Walk: The novel introduction of the "Energy Walk" provides an independent channel to estimate critical exponents, specifically the dynamic exponent $z_c$ and the anomalous dimension $\eta$.
4. Analytical and Numerical Consistency: The paper provides detailed analytical derivations (Appendices A-C) for the scaling forms of fluctuations in both $d=1$ and $d=2$, which are then validated against numerical simulations.

4. Key Results

One Dimension ($d=1$)

• Critical Point: $T_c = 0$.
• Behavior: At $T=0$, the walk is ballistic ($x \sim t$). The probability distribution is U-shaped.
• Persistence: The persistence exponent $\theta$ is extracted from the tips of the distribution, yielding $\theta_{1d} \approx 0.375$, consistent with known exact results.
• Fluctuations: The variance of the spin walk $\sigma_x^2(t)$ shows an exponential divergence as $T \to 0$, consistent with the analytical form $\sigma_x^2/t \sim \exp(4J/k_B T)$.

Two Dimensions ($d=2$)

• Critical Temperature:
  • Using the Ratio Method and Binder Cumulant, $T_c$ is estimated to be $\approx 2.28$ (slightly overestimated due to finite size effects, close to the exact $2.269$).
  • The crossing of Binder cumulant curves for different times $t$ confirms the critical point.
• Order Parameter: The scaled most probable displacement $x_m/t$ behaves exactly like the bulk magnetization $m$, vanishing as $(T_c - T)^\beta$ with $\beta = 1/8$.
• Finite Time Scaling Collapse:
  • Binder Cumulant: Collapses onto a universal curve using the scaling variable $(T-T_c)t^{1/\nu z_c}$ with $1/\nu z_c \approx 0.46$.
  • Spin Walk Fluctuation: The variance $\sigma_x^2$ scales as $t^{2-(d-2+\eta)/z_c}$. The data collapses perfectly, confirming the exponent relation.
  • Energy Walk Fluctuation: The variance $\sigma_y^2$ of the energy walk scales as $t^{2-2(d-1/\nu)/z_c}$. This provides an independent check for the exponents.
• Exponent Extraction: By combining the scaling of the order parameter, spin fluctuations, and energy fluctuations, the authors successfully extract:
  • Static exponents: $\beta, \nu, \eta$.
  • Dynamic exponent: $z_c$.
  • All values are consistent with the known 2D Ising universality class.

5. Significance

• Computational Efficiency: The proposed method offers a significant advantage over traditional Finite Size Scaling. It allows researchers to determine critical properties from a single simulation run of a fixed system size, reducing computational overhead.
• Unified Approach: It unifies the study of order parameters, fluctuations, and persistence under a single "virtual walk" framework.
• Completeness: The inclusion of the Energy Walk is a crucial innovation. While spin walks are related to magnetization, the energy walk (related to local energy density) allows for the independent estimation of the dynamic exponent $z_c$ and the anomalous dimension $\eta$, making the method capable of extracting the complete set of critical exponents.
• Theoretical Validation: The work provides strong numerical and analytical evidence that non-equilibrium time evolution in a single system contains sufficient information to characterize equilibrium critical phenomena, reinforcing the validity of Finite Time Scaling as a robust tool in statistical physics.

In conclusion, the paper establishes "Virtual Walks" as a powerful, efficient, and theoretically sound alternative to Finite Size Scaling for analyzing phase transitions in the Ising model.