On the asymptotic behavior of the Repulsive Pressureless Euler-Poisson System

The main objective of this paper is a study of the asymptotic behavior of distributional solutions to the one-dimensional repulsive pressureless Euler-Poisson system. The system is a model for the dynamics of a mass distribution evolving on \mathbb{R} whose masses exert outward forces on one another. A discrete (describing the evolution of finitely many particles) solution is called sticky if, upon collision, particles stick together and move as one for all subsequent time, according to the conservation of mass and momentum principles. We prove results on the total energy (Hamiltonian) of the system and demonstrate the existence and uniqueness of so-called “perfect” states, where the Hamiltonian is constant over all time and the solution converges to equilibrium, a single stationary particle. We provide a necessary and a sufficient condition for finite-time collapse, and present a quadratic envelope within which a solution must remain in order to collapse. We demonstrate various (counter)examples that illustrate the unique behavior of the repulsive scheme with the sticky condition, analytically and with a computer simulation.

💡 Research Summary

The paper investigates the long‑time behavior of distributional solutions to the one‑dimensional repulsive pressureless Euler‑Poisson (PEP) system, a model in which a mass distribution on the real line exerts outward forces on itself. The authors focus on “sticky” particle solutions: when particles collide they coalesce into a single particle whose mass is the sum of the colliding masses and whose velocity is determined by conservation of momentum.

First, the authors recall that for the attractive PEP system an explicit projection formula provides the exact sticky‑particle solution. They show that this formula fails in the repulsive case because particle trajectories can cross and then uncross, violating the sticky condition (Example 2.1). Consequently, a new analytical framework is required.

The paper’s main contributions are:

1. Energy analysis and perfect states – The total Hamiltonian
   (E(t)=\int\bigl(\tfrac12 v^{2}\rho+\tfrac14(\operatorname{sgn}*\rho)\rho\bigr)dy)
   is shown to be non‑decreasing in time. A “perfect” state is defined as a solution for which (E(t)) remains constant; such a state necessarily converges to the equilibrium ((\delta_{0},0)). For a finite number of particles the authors prove existence and uniqueness of perfect states, establishing that they represent the most energetically efficient way to reach equilibrium.

2. Finite‑time collapse criteria – Two complementary theorems give necessary and sufficient conditions for a finite‑time collision (collapse). The necessary condition involves a quadratic form of the initial positions and velocities; the sufficient condition requires a pointwise inequality for every pair of particles, essentially demanding that the relative velocity be sufficiently negative compared to the squared distance weighted by the combined mass.

3. Quadratic envelope – Theorem 4.1 introduces a sharp quadratic envelope (Q(t)=At^{2}+Bt+C) that any collapsing solution must stay inside. This envelope is derived from conserved quantities (energy and center‑of‑mass motion) and provides a geometric picture of the admissible region for collapse, a tool that is absent in the attractive case.

4. Generalized sticky solutions and ill‑posedness – The set (\mathcal S) of sticky‑particle solutions is not closed under limits of initial data. Examples 2.3 and 2.4 demonstrate sequences of sticky discrete solutions converging to a non‑sticky solution of the repulsive PEP equations. Even the closure (\overline{\mathcal S}) fails to guarantee uniqueness: Example 2.5 shows infinitely many generalized sticky solutions emanating from the same Dirac initial data. This non‑stability stems from the repulsive interaction, which can separate particles after they have merged, allowing repeated merging and splitting.

5. Numerical experiments – The authors implement a Python simulation for three‑particle and many‑particle configurations. The simulations confirm the theoretical collapse conditions and the quadratic envelope: initial data satisfying the sufficient condition collapse rapidly into a single Dirac mass, whereas data close to a perfect state exhibit extremely slow convergence.

Overall, the paper establishes that the repulsive pressureless Euler‑Poisson system behaves fundamentally differently from its attractive counterpart. The sticky‑particle rule, natural for attractive forces, is antagonistic to repulsion, leading to non‑uniqueness, lack of well‑posedness in the classical sense, and the necessity of new analytical tools such as the quadratic envelope. The results provide a solid foundation for future work on non‑local interaction models, kinetic formulations, and robust numerical schemes for systems where long‑range repulsion competes with mass aggregation.