🔬 condensed matter

Glauber dynamics phase transitions in athermal random field Blume-Capel and Blume-Emery-Grifitths models

This paper analytically solves the Glauber dynamics and equilibrium properties of athermal random field Blume-Capel and Blume-Emery-Griffiths models on a complete graph, revealing that the variance of the random field acts as a temperature-like control parameter determining the coincidence of steady states, the nature of phase transitions, and the complex shapes of hysteresis loops which depend non-trivially on initial conditions and crystal field values.

Original authors: Sumedha, Aldrin B E

Published 2026-07-21

📖 8 min read🧠 Deep dive

CC BY 4.0

Original authors: Sumedha, Aldrin B E

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a world made of tiny, invisible magnets, each capable of pointing up, down, or staying perfectly still. In the quiet, frozen world of physics, these magnets usually settle into a calm, predictable pattern, like soldiers standing at attention. This is called "equilibrium," a state where everything has settled down and the rules are simple. But what happens if you shake things up? What if you push these magnets with a magnetic field, or if you throw a chaotic storm of random noise at them? This is the realm of "non-equilibrium," a chaotic dance where the final pose depends entirely on how you started the dance. Scientists study this to understand everything from how magnets work in your hard drive to how opinions spread in a crowd. The key question is: if you push a system hard enough, does it snap into a new state suddenly, or does it slide there slowly? And does it matter if you started with everyone pointing up or everyone pointing down?

This paper dives into that chaotic dance using two specific models of these three-state magnets: the Blume-Capel model and the Blume-Emery-Griffiths model. The researchers, working on a theoretical "complete graph" where every magnet talks to every other magnet, asked a simple but tricky question: If we let these magnets evolve according to strict, energy-saving rules (called Glauber dynamics) at absolute zero temperature, do they end up in the same place as they would if they were allowed to relax slowly and perfectly (equilibrium)? They found that the answer is a resounding "it depends." For low levels of random noise, the magnets get stuck in a memory of their starting position. If you start them all pointing up, they stay up longer than if you start them all pointing down. However, once the random noise gets strong enough—past a specific critical threshold—the system forgets its past, and the chaotic dance finally matches the calm equilibrium. The paper also explores what happens when you add a magnetic field, revealing that the magnets can get stuck in "hysteresis loops," where they refuse to flip back until pushed much harder than it took to flip them forward. These loops can take on wild shapes, from rectangles to wasp-waisted hourglasses, depending on the internal rules of the magnets.

The Story of the Three-Choice Magnets

Imagine a giant ballroom filled with thousands of dancers. In this story, the dancers are our "spins." Unlike a normal dance where you can only face forward or backward, these dancers have a third option: they can stand perfectly still. This is the world of the Blume-Capel and Blume-Emery-Griffiths models. The dancers want to minimize their energy, which means they prefer to align with their neighbors (if one is up, the others want to be up too) or stay still if the music is too weird.

But there's a catch. The ballroom is filled with a chaotic, invisible wind—the random field. This wind blows differently on every single dancer, pushing some up, some down, and some into the middle. The strength of this wind is measured by a number called $R$ (the variance). Think of $R$ as the "loudness" of the chaos. In this paper, the scientists act as the choreographer, watching how the dancers move when the music stops (zero temperature) and they are only allowed to move if it makes them more comfortable (lower energy). This is called Glauber dynamics.

The Great "It Depends" Discovery

The researchers discovered something fascinating about how the dancers behave. If the chaotic wind ( $R$ ) is very weak, the dancers remember where they started.

Scenario A: If you start with everyone facing Up, they might stay facing Up even when the music suggests they should switch to Still.
Scenario B: If you start with everyone facing Down, they might stay Down.

In this low-noise world, the final state of the dance floor depends entirely on the initial state. It's like a game of "Red Light, Green Light" where the players get stuck in a specific pose because they started in a specific spot. The scientists call this a path-dependent system. The system is stuck in a local valley of comfort and can't climb out to find the true best spot because the rules of the game (Glauber dynamics) only allow moves that lower energy, never moves that require a temporary step up.

However, as the scientists turned up the volume on the chaotic wind (increasing $R$ ), something magical happened. Once the wind passed a critical threshold called $R_c$ , the dancers stopped caring about where they started. The wind was so strong that it shook them loose from their local valleys. Suddenly, whether you started them Up or Down, they all ended up in the exact same spot. The chaotic dance finally matched the calm, perfect equilibrium. The paper calculates this critical point exactly. For the simplest model (Blume-Capel), this happens when $R_c = \sqrt{2/\pi} \approx 0.8$ .

The Twist: Frustration and the "Zero" Threshold

The story gets even wilder with the Blume-Emery-Griffiths model. Here, the dancers have an extra rule: a "bi-quadratic" interaction. Imagine that if two neighbors are both "still," they get a bonus, but if they are both "moving," they get a penalty. If this rule is negative (repulsive), it creates frustration. It's like a game where the dancers are told, "If you stand still, you're happy, but if your neighbor stands still, you must move!" This creates a conflict that makes it hard for the group to agree.

The paper found a surprising result here. In some cases of this frustration, the critical threshold $R_c$ drops to zero. This means that even with the tiniest bit of chaotic wind, the dancers forget their starting position immediately. There is no "memory" phase; the system is chaotic from the very first second.

But wait, there's a crossover! The authors found that for certain settings of the "crystal field" (a rule that prefers the "still" state), the system starts with $R_c = 0$ (no memory) but as the wind gets stronger, it suddenly switches to behaving like a normal system with a non-zero $R_c$ . It's as if the dancers, initially too confused to remember anything, suddenly gain a collective memory as the chaos increases. This is a unique "crossover" behavior that the paper predicts and confirms with computer simulations.

The Hysteresis Loops: The Dance of the Wasp-Waist

Finally, the scientists added a magnetic field, which is like a loudspeaker shouting "UP!" or "DOWN!" to the whole room. When you slowly turn up the volume to shout "UP!", the dancers eventually flip. But when you turn the volume down to shout "DOWN!", they don't flip back immediately. They get stuck. This lag is called hysteresis.

In the simplest models, this loop is a perfect rectangle. But in these complex models, the shape of the loop changes based on the internal rules ( $K$ and $\Delta$ ). The paper maps out a zoo of loop shapes:

Rectangular: The standard, sharp flip.
Wasp-waisted: The loop pinches in the middle, looking like an hourglass.
Double Hysteresis: Two separate loops, meaning the dancers flip in two distinct stages.
Hexagonal and Parallelogram: Weird, slanted shapes.

The paper shows that while the size of the loop (the area) shrinks as the chaotic wind ( $R$ ) gets stronger, the shape of the loop is determined by the rules of the game when the wind is zero. Even more interestingly, in the frustrated models, the transition isn't always a sharp snap. Sometimes, the dancers flip one by one in a smooth, continuous flow, creating a loop with no sharp corners at all. This happens when the frustration is high enough to break the "snap" behavior.

The Bottom Line

This paper is a masterclass in how the history of a system matters. It proves that for these three-state magnets, the "equilibrium" state (the perfect, calm solution) is not always the state the system actually reaches if you push it with Glauber dynamics. If the noise is low, the system gets stuck in a memory of its past. If the noise is high, the memory is wiped clean.

The authors didn't just guess this; they derived exact mathematical equations for where these transitions happen and confirmed them with massive computer simulations of 1,000 dancers. They showed that the "critical noise" needed to reset the system's memory can be zero, or it can jump from zero to a specific value depending on the internal rules. They also provided a complete catalog of the strange, beautiful shapes that hysteresis loops can take, from the simple rectangle to the complex wasp-waist.

In short, the paper tells us that in the chaotic world of magnets and noise, how you start matters, unless the noise is loud enough to make you forget everything. And when you do remember, the path you take can look like a rectangle, a wasp, or a hexagon, depending on the hidden rules of the dance.

Technical Summary: Glauber Dynamics Phase Transitions in Athermal Random Field Blume-Capel and Blume-Emery-Griffiths Models

Problem Statement
The paper addresses the lack of a universal framework for understanding non-equilibrium steady states (NESS) compared to equilibrium phase transitions. Specifically, it investigates the zero-temperature ( $T=0$ ) behavior of three-state spin-1 models—the Random Field Blume-Capel Model (RFBCM) and the Random Field Blume-Emery-Griffiths Model (RFBEGM)—under Glauber dynamics. Unlike finite-temperature Glauber dynamics which satisfies detailed balance and reaches equilibrium, the $T=0$ limit breaks detailed balance, leading to non-equilibrium steady states that may depend on initial conditions. The study aims to compare these non-equilibrium steady states with equilibrium states in the presence of Gaussian random fields (variance $R$ ) and external magnetic fields ( $H$ ), focusing on how disorder ( $R$ ) acts as a tuning parameter analogous to temperature.

Methodology
The authors employ a mean-field approach on a fully connected (complete) graph, allowing for exact analytical solutions.

Glauber Dynamics: The system evolves via single spin-flip dynamics where spins flip only if the energy decreases. The dynamics are governed by local fields $L_1$ (linear in spin) and $L_2$ (quadratic in spin). The authors derive self-consistent equations for the order parameters: magnetization ( $m = \langle s \rangle$ ) and quadrupole moment ( $q = \langle s^2 \rangle$ ). The steady state is determined by the stability of fixed points in these equations.
Equilibrium Analysis: To obtain the equilibrium phase diagram, the authors use the large deviation principle to derive a generalized free energy functional (rate function) $I(m, q)$ . The zero-temperature limit of this functional is analyzed to find global minima, which define the equilibrium state.
Comparison: The study contrasts the equilibrium global minima with the non-equilibrium steady states, which are determined by the local stability of fixed points (basins of attraction) and the initial state of the system.
Hysteresis: In the presence of an external field $H$ , the authors analyze quasi-static evolution to map out hysteresis loops, deriving analytical expressions for coercive fields and loop shapes.

Key Contributions and Results

Dependence on Initial State and Disorder Variance ( $R$ ):
- For low values of disorder variance $R$ , the steady state of the Glauber dynamics depends on the initial state (e.g., starting with $m=0$ vs. $m=1$ ). This leads to path-dependent phase transitions.
- Beyond a critical variance $R_c$ , the steady state becomes independent of the initial state, and the non-equilibrium steady state coincides with the equilibrium state.
- Unlike equilibrium where first-order transition points often require numerical determination, the authors derive closed-form equations for the transition coordinates in the Glauber dynamics steady state.
RFBCM Phase Diagram:
- The model exhibits a tricritical point at $R_{tcp} = \sqrt{2/e\pi} \approx 0.48394$ .
- For $R > R_{tcp}$ , the transition is continuous and matches the equilibrium phase diagram.
- For $R < R_{tcp}$ , the system exhibits first-order transitions. The transition point for Glauber dynamics is determined by the condition where a local minimum of the free energy functional becomes a local maximum (loss of local stability), rather than the global minimum condition used in equilibrium. This results in different transition boundaries for different initial states.
RFBEGM and Frustration:
- The RFBEGM includes a bi-quadratic coupling $K$ . Positive $K$ behaves similarly to RFBCM, while negative $K$ introduces frustration.
- Crossover Phenomenon: A significant finding is that for repulsive bi-quadratic interactions ( $K < 0$ ) and specific crystal fields ( $\Delta$ ), the system can undergo a crossover as $R$ increases. Initially, for small $R$ , the system behaves as if $R_c = 0$ (no transition). However, as $R$ increases and the quadrupole moment $q$ drops sufficiently, the system crosses over to a regime with a finite critical variance $R_c = \sqrt{2/\pi} \approx 0.8$ , exhibiting an order-disorder transition.
- In regimes with repulsive bi-quadratic interaction and positive crystal fields, the model can exhibit $R_c = 0$ , meaning no phase transition occurs in the steady state for certain parameter ranges.
Hysteresis Loops:
- In the presence of an external field $H$ , the models exhibit a rich variety of hysteresis loop shapes, including rectangular, hexagonal, parallelogram, wasp-waisted, and double hysteresis loops.
- The authors derive that while the area of the hysteresis loop depends on the disorder variance $R$ , the shape of the loop is determined by the deterministic behavior of the model at $R=0$ .
- Analytical equations for the coercive field are derived as a function of $R$ for both models.
- For RFBEGM with repulsive bi-quadratic coupling ( $K < -1$ ), the dynamics can become non-abelian, leading to hysteresis loops that contain both continuous and first-order transition regions.

Significance and Claims
The paper claims to provide an exact analytical framework for understanding non-equilibrium steady states in driven disordered spin systems at zero temperature. By solving these models on a complete graph, the authors demonstrate that:

The variance of the random field $R$ acts effectively as a temperature, restoring equilibrium behavior above a critical threshold.
Non-equilibrium phase transitions are fundamentally distinct from equilibrium transitions in the low-disorder regime, characterized by path dependence and local stability criteria rather than global free energy minimization.
The RFBEGM serves as a minimal, analytically tractable model to study the interplay between frustration (introduced via repulsive bi-quadratic coupling) and disorder, revealing novel crossover behaviors and complex hysteresis shapes (such as double and wasp-waisted loops) that are relevant to understanding material properties like magnetic rocks and ferroelectrics.
The work bridges the gap between equilibrium statistical mechanics and kinetic studies, offering exact solutions for transition points and coercive fields that are typically obtained only numerically in equilibrium studies.

The Story of the Three-Choice Magnets

The Great "It Depends" Discovery

The Twist: Frustration and the "Zero" Threshold

The Hysteresis Loops: The Dance of the Wasp-Waist

The Bottom Line

More like this