🌀 nonlinear sciences

Breakdown of hydrodynamics in a one-dimensional cold gas

This paper analytically and numerically demonstrates that a one-dimensional cold gas of point particles with mass ratio $m/\mu$ exhibits a breakdown of hydrodynamics into ballistic splatter for most mass ratios, while identifying a specific discrete set of mass ratios where the system remains hydrodynamic-free with at most three moving particles, a distinction that vanishes when initial positions are randomized.

Original authors: Taras Holovatch, Yuri Kozitsky, Krzysztof Pilorz, Yurij Holovatch

Published 2026-07-28

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Taras Holovatch, Yuri Kozitsky, Krzysztof Pilorz, Yurij Holovatch

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 where science tries to understand how energy travels through a crowd. It's a bit like watching a wave of excitement ripple through a stadium, but instead of people cheering, we are tracking tiny, invisible particles bumping into each other. This field, known as statistical mechanics, asks a big question: How do the chaotic, random movements of individual atoms create the smooth, predictable flows we see in the real world, like wind or water? Usually, scientists assume that if you have enough particles, they will eventually settle into a predictable pattern called "hydrodynamics." Think of it as a crowd eventually moving together in a single, flowing direction. But what happens if the crowd is made of different types of people—some heavy, some light—and they are lined up perfectly? Does the smooth flow still happen, or does the perfect order break the rules?

This paper dives into that exact puzzle using a one-dimensional model, which is like a single-file line of marbles. The researchers set up a scenario where a heavy marble at the start of the line gets a sudden kick, sending it zooming into a line of resting marbles. As they collide, they bounce off each other, transferring energy down the line. The big question is: Does this energy spread out smoothly like a wave (hydrodynamics), or does it behave in a weird, unpredictable way? The authors found that the answer depends entirely on the ratio of the heavy marble's weight to the light ones. For some specific weight ratios, the system behaves exactly as the smooth-flow theory predicts. But for a very special, rare set of weight ratios, the system refuses to behave like a fluid at all. Instead, it acts like a perfect, mechanical domino chain where the energy never spreads out the way we expect.

The Great Marble Line-Up

Let's set the scene. Imagine an infinite line of marbles sitting on a track, stretching out to the right. They are perfectly spaced, like soldiers standing at attention. Most of these marbles are light, but every other one is a heavy, dense boulder. At the very start, the first marble (the heavy one) gets a sudden, powerful shove to the right. It zooms forward and smacks into the next marble. Because they are "elastic," they don't stick together; they bounce off like super-bouncy balls, trading speed and direction.

In a typical scenario, you might expect this chain reaction to create a "blast front"—a leading edge of moving particles that pushes outward. Usually, physics tells us this front should move in a specific, slowing-down way, spreading energy out like a ripple in a pond. This is the "hydrodynamic" way, and it's what scientists predicted would happen for many years, especially when the particles are arranged randomly.

The Surprise: When the Rules Break

The authors of this paper ran thousands of computer simulations to see what happens when they change the weight of the heavy marbles compared to the light ones. They discovered something fascinating: if the heavy marbles are exactly the right weight, the whole system changes its personality.

They found a specific list of "magic numbers" for the weight ratio (let's call them $M_i$ ). If the heavy marbles weigh exactly these amounts compared to the light ones, and the marbles are lined up perfectly evenly, the smooth, fluid-like behavior completely disappears. Instead of a wave spreading out, the system acts like a perfectly tuned machine.

Here is what happens in these "magic" cases:

No Splatter: In normal cases, some particles get knocked backward, flying into the negative space behind the starting line. This is called a "splatter." But in the magic cases, nothing flies backward. The energy stays strictly on the forward track.
The Three-Marble Dance: Usually, many marbles are moving at once. In the magic cases, the authors found that at any given moment, there are never more than three marbles moving at the same time. It's like a relay race where the baton is passed so efficiently that only three runners are ever on the track.
The Bullet-Proof Front: The leading edge of the moving particles (the blast front) doesn't slow down or spread out. It moves in a straight, "ballistic" line, like a bullet fired from a gun, maintaining a constant speed. This is the opposite of the hydrodynamic "ripple" effect.

The paper explicitly shows that for these specific weights (up to the 700th "magic number" they calculated), the system is not behaving hydrodynamically. It is behaving in a rigid, predictable, mechanical way that defies the usual fluid rules.

The Twist: Order vs. Chaos

The most playful part of this discovery is how sensitive the system is. The authors showed that if you take those same "magic" weights but shuffle the starting positions of the marbles so they aren't perfectly evenly spaced (making them random), the system suddenly snaps back to the normal, fluid-like behavior. The "magic" only works when the setup is perfectly ordered.

This suggests that the system is "non-ergodic," which is a fancy way of saying that the history of how you set it up matters forever. If you start with a perfect line, it stays perfect and mechanical. If you start with a messy line, it becomes a fluid. The paper proves that for these specific weights, the "perfect line" setup creates a unique state of matter that refuses to follow the standard rules of energy spreading.

What They Found (and What They Didn't)

The researchers didn't just guess this; they calculated the exact formula for these magic weights and ran massive computer simulations to confirm it. They found that for weights like 2, 7, 12, and 13 (which are not the magic numbers), the system behaves normally, with the blast front slowing down and spreading out. But for the magic numbers (like $M_3 \approx 9.1$ ), the system stays ballistic.

They also noted that if you get close to a magic number but not exactly on it, the system might look ballistic for a while, but eventually, it will cross over to the normal fluid behavior. The closer you get to the magic number, the longer it takes for the system to "wake up" and start behaving like a fluid. But if you hit the magic number exactly, it never wakes up; it stays in that rigid, mechanical state forever.

In short, this paper reveals that in a world of colliding particles, perfect order can create a "frozen" flow that ignores the usual laws of fluid dynamics. It's a reminder that sometimes, the most predictable setups lead to the most surprising exceptions.

Technical Summary: Breakdown of Hydrodynamics in a One-Dimensional Cold Gas

Problem Statement
The paper investigates the interconnection between microscopic dynamics and macroscopic hydrodynamic descriptions in a one-dimensional cold gas system following an instantaneous energy injection (an "explosion"). Specifically, the authors study an asymmetric Alternating Hard Particle (AHP) gas model consisting of point particles with alternating masses $m$ and $\mu$ (where $m \ge \mu$ ) distributed on the positive half-axis. The dynamics are initiated by imparting a positive velocity to the particle at the origin.

Previous work (e.g., Chakraborti et al., 2022) established that for specific mass ratios (such as $m/\mu = 2$ ) and random initial particle positions, the system evolves into a hydrodynamic regime characterized by:

A blast front (rightmost particle) moving sub-ballistically as $R(t) \sim t^\delta$ with $\delta < 1$ .
The formation of a "splatter" (particles recoiling into the negative half-axis).
A specific decay of energy in the positive half-axis and a power-law growth in the number of collisions.

The central question addressed here is whether this hydrodynamic behavior is universal or if specific mass ratios and initial configurations can lead to a breakdown of this regime, preserving ballistic dynamics.

Methodology
The authors employ a dual approach combining analytical derivation and numerical simulation:

Analytical Derivation: The study utilizes the exact elastic collision rules for point particles to trace the velocity evolution of particle triplets. By solving for conditions where the first two particles in a triplet come to a complete stop after a finite number of collision rounds, the authors derive a specific set of mass ratios $\{M_i\}$ .
Numerical Simulation: Molecular Dynamics (MD) simulations are performed using a heap data structure to manage collision events efficiently. The simulations track observables including the blast front position $R(t)$ , the total number of collisions $C(t)$ , and the energy of particles in the positive half-axis $E_{x \ge 0}(t)$ .
Scaling Analysis: Time dependencies of observables are fitted to power laws ( $O(t) \sim t^\alpha$ ) to determine exponents $\delta$ , $\beta$ , and $\eta$ . The relationship between these exponents is tested against hydrodynamic predictions ( $\eta = 2\delta$ and $\delta = (2-\beta)/3$ ).

Key Contributions and Results

Identification of Special Mass Ratios ( $M_i$ ):
The authors explicitly derive a collection of positive numbers $\{M_i : i \in \mathbb{N}\}$ defined by the formula:
$M_i = \cot\left(\frac{\pi}{2(2i + 1)}\right) \cot\left(\frac{\pi}{2i + 1}\right)$
For mass ratios $m/\mu = M_i$ (verified for $i \le 700$ ), the system exhibits fundamentally different behavior compared to the hydrodynamic regime.
Breakdown of Hydrodynamics for Equidistant Initial Conditions:
When particles are initially placed at equidistant positions ( $x_k = k$ ) and $m/\mu = M_i$ :
- Absence of Splatter: No particles enter the negative half-axis; the normalized energy $\epsilon(m)$ remains 1, indicating no energy loss to the "splatter."
- Limited Active Particles: The number of simultaneously moving particles never exceeds three.
- Ballistic Propagation: The blast front moves ballistically ( $R(t) \sim t^1$ ), deviating from the hydrodynamic prediction ( $\delta \approx 0.628$ ).
- Mechanism: The dynamics involve a recursive process where a triplet of particles undergoes exactly $i$ rounds of mutual collisions, after which the first two particles stop completely, transferring the motion to the next triplet.
Dependence on Initial Conditions (Non-ergodicity):
The paper demonstrates that the system's behavior is highly sensitive to initial conditions.
- For $m/\mu = M_i$ with equidistant initial positions, the system remains ballistic.
- For $m/\mu = M_i$ with randomly distributed initial positions (sampled from a uniform ensemble), the system recovers the hydrodynamic behavior described in previous literature (sub-ballistic front, splatter formation).
- This dependence on the initial state suggests the model is non-ergodic for these specific mass ratios.
Transition Regimes:
For mass ratios close to, but not exactly equal to, $M_i$ (e.g., $m=9$ near $M_3 \approx 9.09$ ), the system exhibits a crossover. Initially, the dynamics appear ballistic-like, but over longer timescales (which increase as $m$ approaches $M_i$ ), the system transitions to the hydrodynamic regime.

Significance
The paper claims to identify a specific set of parameters where the standard hydrodynamic description of energy transport in a 1D cold gas fails, even in the presence of elastic collisions. The significance lies in:

Refuting Universality: Demonstrating that the hydrodynamic scaling laws observed in previous studies are not universal for all mass ratios and initial configurations.
Exact Solvability: Providing an exact analytical condition ( $M_i$ ) under which the complex many-body collision problem simplifies to a sequence of independent triplet interactions, resulting in ballistic motion.
Non-ergodicity: Highlighting that for specific mass ratios, the long-time behavior of the system is determined by the precise initial configuration, challenging assumptions of ergodicity in such 1D systems.

The authors note that while these results apply to the asymmetric AHP gas, the findings also extend to the symmetric AHP gas model studied in prior literature. The work concludes that the hydrodynamic regime is a specific outcome of the interplay between mass ratios and initial disorder, rather than an inevitable feature of 1D elastic particle systems.

The Great Marble Line-Up

The Surprise: When the Rules Break

The Twist: Order vs. Chaos

What They Found (and What They Didn't)

Technical Summary: Breakdown of Hydrodynamics in a One-Dimensional Cold Gas

More like this