Violation of local equilibrium thermodynamics in one-dimensional Hamiltonian-Potts model

By numerically studying a one-dimensional Hamiltonian-Potts model with fractional spatial derivatives under steady heat conduction, this paper demonstrates that a stationary interface between coexisting phases exhibits a temperature deviation from equilibrium values, thereby confirming the violation of local equilibrium and the stabilization of metastable states in nonequilibrium first-order phase transitions.

The Gist

Imagine a long, thin rope made of tiny, interconnected beads. In the world of physics, this rope represents a material that can exist in two different "moods": an ordered mood (where the beads are neatly lined up, like soldiers) and a disordered mood (where the beads are jumbled and chaotic, like a crowd at a concert).

Usually, there is a specific "switching point" temperature where the rope changes from soldiers to a crowd. If you heat it up, it snaps to chaotic; if you cool it down, it snaps back to orderly.

The Problem: The "Local" Rulebook

For a long time, scientists believed in a rule called Local Equilibrium. Think of this rule as a strict traffic cop who says: "No matter what is happening elsewhere, every single inch of this rope must obey the standard temperature rules right where it is standing."

According to this old rule, if you have a hot end and a cold end, and a "battle line" (an interface) forms between the orderly soldiers and the chaotic crowd, that battle line should sit exactly at the standard switching temperature.

The Experiment: A One-Dimensional "Magic" Rope

The authors of this paper wanted to test if this traffic cop is actually right. They built a computer simulation of a one-dimensional rope (a line of beads).

Here's the catch: In real life, a simple one-dimensional line usually can't hold a battle between order and chaos; it's too weak. To fix this, the scientists gave the rope a special "magic" property. They used a mathematical trick called a fractional derivative.

The Analogy: Imagine the beads on the rope can't just talk to their immediate neighbors; they can also "feel" the vibe of beads far away, but in a very specific, long-range way. This trick makes the one-dimensional rope behave exactly like a much more complex, two-dimensional sheet of material, allowing the battle between order and chaos to happen.

They attached a cold bath to the left end and a hot bath to the right end, creating a steady flow of heat through the rope.

The Discovery: The "Rebel" Interface

When they watched the simulation, something surprising happened. The battle line (the interface) did not sit at the standard switching temperature.

• The Old Rule said: The interface should be at temperature $T_c$ (the standard switch point).
• The Reality showed: The interface was actually hotter than $T_c$.

It's as if the soldiers on the left side were being pushed by the heat current to stay in their orderly formation even though the local temperature was high enough that they should have turned into a chaotic crowd. The steady flow of heat acted like a "glue," stabilizing a state that should have been unstable.

The New Theory: "Global" Thermodynamics

The paper confirms that a newer theory, called Global Thermodynamics, predicted this exactly.

The Analogy:

• Local Thermodynamics is like judging a whole city's weather by looking at just one street corner. It assumes the corner knows nothing about the rest of the city.
• Global Thermodynamics is like looking at the city as one giant, connected organism. It realizes that the heat flowing from the hot side to the cold side changes the rules for everyone in the system, including the battle line.

The authors found that the temperature of the interface matched the "Global" prediction perfectly. This proves that when a system is under a steady heat current, the old "Local" rulebook breaks down. The system doesn't just follow local rules; it follows the rules of the whole system working together.

The Bottom Line

This study didn't just find a glitch; it found a fundamental truth about how nature works when things are out of balance.

1. Local rules fail: You cannot always assume that a small part of a system behaves like it's in a calm, isolated environment.
2. Heat currents stabilize the unstable: A steady flow of heat can lock a system into a "metastable" state (like superheated ice or supercooled water) that would normally disappear instantly.
3. It's universal: This happens even in a simple one-dimensional line, proving it's a fundamental feature of nature, not just a quirk of complex 3D shapes.

The paper concludes that this one-dimensional model is a perfect, simplified "laboratory" to study these complex thermodynamic limits, showing that the universe is more interconnected than the old "local" rules suggested.

Technical Summary

Technical Summary: Violation of Local Equilibrium Thermodynamics in a One-Dimensional Hamiltonian-Potts Model

Problem Statement
The dynamics of Hamiltonian systems are typically described by hydrodynamic equations assuming local equilibrium thermodynamics. However, this approximation is known to fail in the presence of strong heat currents and steep temperature gradients, particularly during first-order phase transitions. A key question is whether the temperature at the interface between coexisting ordered and disordered phases in a nonequilibrium steady state (NESS) equals the equilibrium transition temperature ($T_c$). While global thermodynamics predicts a deviation from $T_c$ due to the stabilization of metastable states by heat currents, previous numerical verifications using two-dimensional Hamiltonian-Potts models suffered from significant finite-size and finite-roughness effects. These effects, arising from coarse spatial discretization and limited system sizes, often mask the violation of local equilibrium, necessitating computationally prohibitive simulations to reach the thermodynamic limit.

Methodology
To overcome these limitations, the authors propose a minimal, controlled numerical model: a one-dimensional Hamiltonian-Potts model incorporating fractional-order spatial derivatives.

• Model Construction: The standard one-dimensional model with short-range interactions does not support phase transitions. To induce a first-order phase transition (specifically for $n=5$ states), the authors introduce a nonlocal spatial derivative operator $D$. This operator is defined via a Fourier representation where the spectral density at low wavenumbers scales as $\rho(k) = \text{const}$, effectively reproducing the low-wave-number density of states of a standard two-dimensional model. This dimensional reduction allows for the study of phase coexistence in a 1D setting free from transverse geometrical constraints.
• Hamiltonian: The system is governed by a Hamiltonian $H = \int dx \, h$, where the kinetic and gradient terms involve the fractional derivative $D$. The potential term is constructed to have $n$ equivalent minima arranged as vertices of a regular simplex.
• Simulation Setup: Nonequilibrium steady states are realized by attaching heat baths at temperatures $T_1$ and $T_2$ (satisfying $T_1 < T_c < T_2$) to the boundaries of a central Hamiltonian region. The boundary regions are modeled using Langevin dynamics, while the central region evolves via energy-conserving Hamiltonian equations.
• Numerical Implementation: The equations of motion are integrated using a second-order symplectic scheme for Hamiltonian dynamics and a second-order stochastic Runge-Kutta scheme for Langevin dynamics. The system size ($L_x = 256$) and lattice spacing ($\Delta x = 1/4$) are chosen to minimize discretization effects while maintaining computational feasibility.

Key Contributions and Results

1. Realization of Phase Coexistence in 1D: The study successfully demonstrates that the introduction of fractional derivatives enables a one-dimensional system to exhibit a first-order phase transition and stable phase coexistence under steady heat conduction, a phenomenon typically restricted to higher dimensions or long-range interactions.
2. Violation of Local Equilibrium: Numerical simulations reveal that the temperature at the stationary interface ($\theta$) systematically deviates from the equilibrium transition temperature ($T_c \approx 0.634$). Specifically, the interface temperature is found to be higher than $T_c$, indicating the stabilization of a superheated ordered state. This deviation persists well beyond statistical uncertainty, providing clear evidence that the local equilibrium description breaks down at the interface.
3. Quantitative Agreement with Global Thermodynamics: The observed interface temperatures are in quantitative agreement with the predictions of global thermodynamics. The global temperature $\tilde{T}$, defined as the density-weighted spatial average of the local temperature, and the theoretical formula for the interface temperature deviation (derived from a variational principle) accurately describe the simulation data.
4. Thermal Conductivity Asymmetry: The results indicate that the thermal conductivity in the ordered phase ($\kappa_o$) is smaller than that in the disordered phase ($\kappa_d$), consistent with the theoretical prediction that the interface temperature shifts toward the phase with lower thermal conductivity.

Significance
The paper claims that the stabilization of metastable states and the violation of local equilibrium are intrinsic, universal features of nonequilibrium first-order phase transitions, independent of spatial dimensionality. By isolating the essential mechanisms in a minimal one-dimensional framework, the study confirms that these phenomena are governed by macroscopic constraints (as described by global thermodynamics) rather than microscopic details or specific dimensionalities. The work establishes a robust, computationally efficient model for exploring the fundamental limits of thermodynamic descriptions in nonequilibrium steady states, offering a controlled alternative to the computationally expensive and effect-prone simulations of higher-dimensional systems.