🔢 mathematics

Equilibria of the Pressureless Euler System in Dimension 1

This paper establishes a necessary and sufficient condition based solely on initial data for sticky particles solutions of the one-dimensional pressureless Euler equations to reach equilibrium, while also characterizing the equilibrium state and providing a formula for the collapse time.

Original authors: Nicholas Biglin, Joseph Crachiola, Thomas Kunz, Omkar Maralappanavar

Published 2026-06-16

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Nicholas Biglin, Joseph Crachiola, Thomas Kunz, Omkar Maralappanavar

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 long, straight highway stretching infinitely in both directions. On this highway, there are millions of tiny cars, each representing a tiny piece of "mass." These cars are driving according to two simple rules:

They never disappear: If two cars merge, they don't vanish; they become one bigger car.
They stick together: If a faster car catches up to a slower car in front of it, they don't crash and bounce off. Instead, they lock bumpers and stick together, moving forward as a single, heavier unit at a new, shared speed.

This is the world of the Pressureless Euler System described in this paper. The authors are trying to answer a very specific question about these sticky cars: Will they eventually stop moving relative to each other and form a stable, unchanging pattern?

Here is the breakdown of their findings using simple analogies.

The "Traffic Report" (The Main Problem)

In physics, we usually study how things move. But here, the authors are looking at the "end game." If you start with a specific arrangement of cars and speeds, will the traffic eventually settle down into a calm state (equilibrium), or will it keep changing forever?

The paper provides a crystal-clear "Yes or No" test based entirely on the starting conditions (where the cars are and how fast they are going). You don't need to simulate the traffic for hours; you just need to look at the initial "traffic report."

The "Momentum Balance" Test (The Condition)

The authors discovered a specific mathematical condition (let's call it the "No-Left-Leak" rule) that determines if the system will settle down.

Imagine you are standing at the start of the highway (the left side). You look at all the cars to your left and calculate their total "push" (momentum).

The Rule: For the traffic to eventually settle, the total "push" of any group of cars starting from the far left must never be negative.
The Analogy: Think of the cars as a team pulling a rope. If the cars on the far left are pulling the rope to the left with more force than the cars on the right can pull back, that chunk of cars will just keep sliding left forever, never stopping.
The Verdict: If, at any point, a group of cars on the left has a net "leftward pull," the system will never reach equilibrium. It will keep drifting apart. However, if the "leftward pull" is never negative (meaning the system is balanced or pushing right), the traffic will eventually settle down.

What Does "Settling Down" Look Like?

If the system passes the test, the paper describes exactly what the final traffic jam looks like. It's a mix of two behaviors:

The "Clumping" Zones: In areas where the cars were moving in a way that created tension (like a fast car chasing a slow one), they will all crash into each other and merge into a single, giant "super-car" (a point mass). This is like a traffic jam where everyone stops and becomes one big block.
The "Stagnant" Zones: In areas where the cars were already perfectly balanced (no one was trying to overtake anyone), they will just stay exactly where they started, moving in perfect unison forever.

The authors provide a map that tells you exactly which parts of the highway will turn into a single "super-car" and which parts will remain as a steady stream of individual cars.

The "Time to Freeze" (Collapse Time)

The paper also answers: How long does it take for the traffic to stop changing?

They found a formula to calculate the exact moment when the "clumping" stops.

Imagine the cars are sliding toward each other. The formula calculates the precise second when the last fast car finally catches the last slow car and they lock bumpers.
After this specific time, the traffic pattern is frozen. Nothing changes anymore. The paper gives a recipe to calculate this "freeze time" just by looking at the starting positions and speeds.

Why This Matters (According to the Paper)

The authors note that while other versions of this problem (like cars that repel each other) are chaotic and hard to predict, this specific "sticky" version is surprisingly orderly.

Predictability: You don't need to know the future to know the outcome; the starting conditions tell the whole story.
Universality: This applies whether the cars are spread out over a small neighborhood or scattered across the entire universe.

Summary

In short, this paper is a guidebook for a sticky traffic jam. It says:

Check the balance: If the left side of the traffic ever has too much "leftward push," the jam will never settle.
If it balances: The traffic will eventually freeze into a pattern of giant clumps and steady streams.
The timing: We can calculate the exact second the chaos ends and the order begins, using only the starting data.

The paper does not discuss using this for real-world traffic control, weather prediction, or medical applications. It is purely a mathematical description of how these specific "sticky particles" behave over time.

Technical Summary: Equilibria of the Pressureless Euler System in Dimension 1

Problem Statement
This paper investigates the long-time behavior of the one-dimensional pressureless Euler system (PE), originally proposed by Zeldovich to model large-scale structure formation. The system is given by the continuity and momentum equations:
$\partial_t \rho + \partial_y(\rho v) = 0, \quad \partial_t(\rho v) + \partial_y(\rho v^2) = 0$
where $\rho$ is a mass distribution (probability measure with finite second moment) and $v$ is the velocity field. The authors focus specifically on Sticky Particle Solutions (SPS), a class of distributional solutions where mass particles stick together upon collision. While the existence of such solutions is established via a projection formula (Natile-Savaré), the conditions under which these solutions converge to an equilibrium state, the nature of that equilibrium, and the time required for collapse have remained open questions. Unlike the attractive Euler-Poisson system, where equilibrium is typically a single point mass, the pressureless Euler system allows for a rich variety of equilibrium distributions.

Methodology
The analysis relies heavily on the projection formula derived in [10], which translates the problem from physical "position-space" to "mass-space." By utilizing the optimal transport map $N_t: (0,1) \to \mathbb{R}$ (where $N_t \# \mathcal{L} = \rho_t$ ), the evolution of the system is characterized by projecting the linear trajectory $N_0 + tV_0$ onto the convex cone $K$ of non-decreasing functions in $L^2(0,1)$ :
$N_t = \text{proj}_K(N_0 + tV_0)$
Here, $V_0 = v_0 \circ N_0$ . The authors employ tools from convex analysis, specifically properties of the convex envelope ( $\text{conv}$ ) and the right derivative, to analyze the asymptotic behavior of $N_t$ . They define a primitive function $F(x) = \int_0^x V_0(w) dw$ and utilize the Galilean invariance of the system to normalize the total momentum to zero.

Key Contributions and Results

Necessary and Sufficient Condition for Equilibrium:
The paper establishes a precise condition, denoted as $(\star)$ , based solely on the initial data $(N_0, V_0)$ that determines whether the SPS converges to an equilibrium. The condition requires:
- Zero total momentum: $\int_0^1 V_0(w) dw = 0$ (equivalent to $F(1)=0$ ).
- Non-negative partial momentum: $\int_0^x V_0(w) dw \geq 0$ for all $x \in (0,1)$ (equivalent to $F(x) \geq 0$ ).
The authors prove that the SPS converges to an equilibrium in the 1-Wasserstein metric if and only if condition $(\star)$ holds. If $F(x) < 0$ for any $x$ , the solution diverges (mass spreads to infinity).
Characterization of the Equilibrium State:
When condition $(\star)$ is satisfied, the paper provides a complete description of the limiting distribution $\rho_\infty$ . The equilibrium is determined by the zeros of the function $F$ :
- Stationary Regions: On any interval where $F$ is identically zero, the mass distribution remains stationary for all time ( $N_t(x) = N_0(x)$ ).
- Collapse Regions: On any interval $(a, b)$ where $F(a)=F(b)=0$ and $F > 0$ strictly inside, the mass within that interval collapses to a single Dirac measure (point mass) located at the center of mass of that initial "chunk" of mass.
  This implies that the equilibrium is a superposition of stationary mass and collapsed point masses.
Time of Collapse:
The authors derive an explicit formula for the time of collapse ( $T_c$ ), defined as the minimal time $T$ such that $\rho_t = \rho_\infty$ for all $t \geq T$ . If the solution converges, $T_c$ is given by the supremum over all intervals $(a,b)$ where $F$ is positive:
$T_c = \sup_{[a,b] \in I} \sup_{x \in (a,b)} \frac{\int_x^b N_0(w) dw + \frac{x-b}{b-a}\int_a^b N_0(w) dw}{\int_0^x V_0(w) dw}$
If this value is infinite, the solution approaches equilibrium asymptotically but never reaches it in finite time.
Convergence Metrics:
The results are primarily framed in terms of the 1-Wasserstein distance ( $W_1$ ). The paper notes that if the support of the solution remains bounded (which occurs if the initial support is bounded and $(\star)$ holds), convergence in $W_p$ for any $p$ follows. The authors discuss the extension to $W_2$ convergence, noting that while it holds under specific monotonicity conditions near the boundaries of the mass interval, a general condition depending only on initial data for $W_2$ convergence remains an open difficulty due to the global nature of the convex envelope operation.

Significance and Claims
The paper claims to provide a complete characterization of the equilibrium behavior for sticky particle solutions to the 1D pressureless Euler system using only initial conditions. It resolves the three central questions regarding existence of equilibrium, the specific form of that equilibrium, and the time of collapse. The authors emphasize that their approach avoids the need for assumptions regarding compact support over time, which were previously required in related literature (e.g., [8]). The work distinguishes itself from the repulsive Euler-Poisson system, where such predictive methods fail due to pathological behavior, highlighting that the momentum-driven nature of the pressureless Euler system allows for this precise determination. The results rely on the projection formula and convex analysis, offering a rigorous link between the geometry of the initial velocity distribution and the asymptotic state of the system.