pith. sign in

arxiv: 2601.19733 · v2 · submitted 2026-01-27 · 🧮 math.AP

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

Nicholas Biglin , Joseph Crachiola , Jack Curtis , Thomas Kunz , Omkar Maralappanavar , Adrian Tudorascu This is my paper

Pith reviewed 2026-05-16 10:30 UTC · model grok-4.3

classification 🧮 math.AP
keywords Euler-Poisson systempressureless fluidssticky collisionsasymptotic behaviorfinite-time collapseHamiltonian energyrepulsive forcesone-dimensional dynamics
0
0 comments X

The pith

The repulsive pressureless Euler-Poisson system in one dimension admits unique perfect states in which total energy stays constant and the mass distribution converges to a single stationary particle.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies the long-term evolution of mass distributions whose particles repel one another according to the one-dimensional pressureless Euler-Poisson equations. It restricts attention to distributional solutions that obey a sticky collision rule: when particles meet they merge and thereafter travel as one unit while conserving total mass and momentum. Under this rule the authors prove that the system's total energy, or Hamiltonian, is nonincreasing, and they establish the existence and uniqueness of perfect states in which the Hamiltonian is exactly constant and the configuration collapses to equilibrium. They supply a necessary and sufficient condition for the entire mass to reach a single point in finite time and show that any such collapsing solution must remain inside a quadratic envelope in space-time. These statements are illustrated by explicit examples and by direct numerical simulation of the discrete particle system.

Core claim

For the one-dimensional repulsive pressureless Euler-Poisson system with the sticky collision rule, there exists a unique perfect state associated with any given initial data. In this state the Hamiltonian remains constant for all future times and the solution converges to the equilibrium consisting of one stationary particle whose mass equals the total mass. A necessary and sufficient condition for finite-time collapse is identified, and every collapsing solution is shown to lie inside an explicit quadratic envelope.

What carries the argument

The sticky collision rule, under which colliding particles merge into a single entity that conserves mass and momentum and continues to evolve under the repulsive force field.

If this is right

  • The Hamiltonian is conserved if and only if the solution is perfect and strictly decreases otherwise until equilibrium is reached.
  • Finite-time collapse occurs precisely when a computable integral condition on the initial positions, velocities, and masses is satisfied.
  • Any solution that collapses must remain inside the quadratic envelope for all time up to the collapse instant.
  • The uniqueness of the perfect state determines the long-time asymptotics completely from the initial data.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The quadratic envelope supplies a direct numerical test that can be applied to any discrete initial data to decide whether collapse will occur.
  • Generalizing the sticky rule to two or three dimensions would likely produce analogous perfect states but with different collapse criteria.
  • The energy-constancy property of perfect states may serve as a stability benchmark for numerical schemes applied to related repulsive fluid models.

Load-bearing premise

That distributional solutions to the system remain well-defined for every initial datum once the sticky collision rule is imposed in one space dimension.

What would settle it

An explicit initial configuration of two or more particles whose computed trajectory collapses in finite time yet exits the predicted quadratic envelope at some moment before collapse.

Figures

Figures reproduced from arXiv: 2601.19733 by Adrian Tudorascu, Jack Curtis, Joseph Crachiola, Nicholas Biglin, Omkar Maralappanavar, Thomas Kunz.

Figure 2.1
Figure 2.1. Figure 2.1: The trajectories in Example 2.1. Left is the attractive case, right is the repulsive case. [PITH_FULL_IMAGE:figures/full_fig_p006_2_1.png] view at source ↗
Figure 2.2
Figure 2.2. Figure 2.2: Left: The discrete approximation 𝜌 3 𝑡 from Example 2.4. Right: 𝜌𝑡 as in Example 2.4. Example 2.5. Consider the initial condition (𝛿0 , 0). We claim that, while there is only one solution in  satisfying this condition, there are infinitely many in ̄. Indeed, let 𝑐 > 0 and consider the sticky particles solution to the initial conditions 𝜌 𝑛 0 = 1 2 𝛿− 1 𝑛 + 1 2 𝛿 1 𝑛 , 𝑣0 = 𝑐 ⋅ sgn . Taking 𝑛 → ∞, we se… view at source ↗
Figure 2.3
Figure 2.3. Figure 2.3: The distribution described in Example 2.7. [PITH_FULL_IMAGE:figures/full_fig_p013_2_3.png] view at source ↗
Figure 2.4
Figure 2.4. Figure 2.4: Glancing and non-glancing sticky collisions. [PITH_FULL_IMAGE:figures/full_fig_p019_2_4.png] view at source ↗
Figure 4.1
Figure 4.1. Figure 4.1: The quadratic bound described in Theorem 4.1 is shown in dashed lines. Left: the solution [PITH_FULL_IMAGE:figures/full_fig_p036_4_1.png] view at source ↗
Figure 5.1
Figure 5.1. Figure 5.1: An example DSPS and the state of the solution as the position and velocity of a particle varies. [PITH_FULL_IMAGE:figures/full_fig_p044_5_1.png] view at source ↗
read the original abstract

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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper studies the asymptotic behavior of distributional solutions to the one-dimensional repulsive pressureless Euler-Poisson system. It defines sticky discrete solutions in which particles merge upon collision while conserving mass and momentum. The authors establish results on the total Hamiltonian energy, prove existence and uniqueness of 'perfect' states in which the Hamiltonian remains constant and the solution converges to a single stationary particle, give necessary and sufficient conditions for finite-time collapse together with a quadratic envelope bound, and illustrate the dynamics through analytical examples and numerical simulations.

Significance. If the proofs are complete and the sticky solutions are verified to remain distributional, the results would clarify long-time behavior and collapse criteria for repulsive pressureless systems with inelastic collisions. The characterization of perfect states and the envelope bound are potentially useful for related models in aggregation or plasma dynamics. The combination of analysis and simulation is a positive feature, though the absence of explicit function spaces, stability estimates, and post-collision weak-form verification limits immediate applicability.

major comments (3)
  1. [Definition of sticky solutions and weak formulation] The central claim that sticky collisions preserve the distributional formulation of the momentum equation under the nonlocal repulsive force is load-bearing for all subsequent results on energy and collapse, yet no explicit verification or orthogonality condition is supplied after a merger; the weak form may fail to hold across collision times because the internal repulsion term is discarded by the sticky rule.
  2. [Results on perfect states and Hamiltonian] Theorem on existence and uniqueness of perfect states (where Hamiltonian is constant and convergence to a single stationary particle occurs) assumes global-in-time distributional solutions exist for arbitrary initial data, but provides no stability estimate or function-space setting (e.g., measures with bounded variation) to justify that the nonlocal force remains well-defined post-collision.
  3. [Finite-time collapse criteria] The necessary and sufficient condition for finite-time collapse and the quadratic envelope bound rest on the same global existence assumption; without an a-priori estimate controlling the nonlocal integral term immediately after sticking, the collapse criterion may only hold under additional regularity that is not stated.
minor comments (2)
  1. [Abstract] The abstract refers to a 'quadratic envelope' without an equation number or explicit statement of the bound; this should be labeled for cross-reference.
  2. [Examples and simulations] Numerical simulation details (time-step size, number of particles, discretization of the nonlocal force) are not provided, making reproducibility difficult.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We agree that additional clarifications are needed regarding the preservation of the weak formulation under sticky collisions, the precise functional setting, and a-priori estimates. We will incorporate these into a revised version. Below we respond point by point to the major comments.

read point-by-point responses

  1. Referee: The central claim that sticky collisions preserve the distributional formulation of the momentum equation under the nonlocal repulsive force is load-bearing for all subsequent results on energy and collapse, yet no explicit verification or orthogonality condition is supplied after a merger; the weak form may fail to hold across collision times because the internal repulsion term is discarded by the sticky rule.

    Authors: We acknowledge this point and the need for explicit verification. The sticky rule is defined via conservation of mass and momentum at collision instants, which ensures the post-merger velocity is consistent with the pre-collision data. Because the repulsive force is nonlocal and acts on the entire measure, the net force on the merged particle remains well-defined in the distributional sense; the discarded internal component is precisely the one that would violate inelasticity. Nevertheless, to address the concern directly, we will add a dedicated lemma in the revised manuscript verifying that the weak form of the momentum equation holds across collision times, including an orthogonality check for compactly supported test functions at merger instants. revision: yes

  2. Referee: Theorem on existence and uniqueness of perfect states (where Hamiltonian is constant and convergence to a single stationary particle occurs) assumes global-in-time distributional solutions exist for arbitrary initial data, but provides no stability estimate or function-space setting (e.g., measures with bounded variation) to justify that the nonlocal force remains well-defined post-collision.

    Authors: The referee correctly identifies the missing functional framework. In the current manuscript the solutions are constructed as Radon measures with finite second moment, but stability estimates are indeed not stated explicitly. We will revise the statement of the theorem on perfect states to work in the space of measures with bounded total variation, and we will insert a stability estimate showing that the nonlocal force operator remains continuous with respect to weak convergence after each collision. This estimate follows from the repulsive character of the kernel and the conservation of the Hamiltonian. revision: yes

  3. Referee: The necessary and sufficient condition for finite-time collapse and the quadratic envelope bound rest on the same global existence assumption; without an a-priori estimate controlling the nonlocal integral term immediately after sticking, the collapse criterion may only hold under additional regularity that is not stated.

    Authors: We agree that an a-priori control on the nonlocal term immediately after sticking is required for the collapse criteria to be rigorous. We will add such an estimate in the revised proofs, derived directly from the non-increasing property of the Hamiltonian and the positivity of the repulsive kernel. The resulting bound will be uniform in time between collisions and will be used to justify both the necessary and sufficient conditions as well as the quadratic envelope without invoking extra regularity. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results follow from model equations and sticky collision definition

full rationale

The paper defines sticky discrete solutions directly from conservation of mass and momentum upon collision. Hamiltonian constancy for perfect states and convergence to a single stationary particle are then proved from the resulting dynamics and the repulsive force law. Necessary/sufficient conditions for finite-time collapse and the quadratic envelope are stated as consequences of the same equations, not obtained by fitting or by redefining the target quantity in terms of itself. No step reduces a claimed prediction to a self-citation chain or to an ansatz imported from prior work by the same authors. The derivation remains self-contained against the stated 1D distributional framework.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work relies on standard mathematical assumptions for distributional solutions of hyperbolic systems and the definition of the sticky collision rule; no free parameters or invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption Distributional solutions to the Euler-Poisson system satisfy conservation of mass and momentum under the sticky collision rule.
    Invoked throughout the abstract as the basis for defining sticky solutions and deriving energy properties.

pith-pipeline@v0.9.0 · 5494 in / 1234 out tokens · 35687 ms · 2026-05-16T10:30:58.765662+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages

  1. [1]

    BECERRA, J

    B. BECERRA, J. LINDEROTH, H. PESIN, A. TUDORASCU, R. WASSINK,Sticky Particles solutions 45 for the attractive pressureless Euler-Poisson system; a projection formula and asymptotic behavior, J. Math. Anal. Appl.,550, No. 2 (2025)

    work page 2025

  2. [2]

    BRENIER, E

    Y. BRENIER, E. GRENIER,Stickyparticlesandscalarconservationlaws,SIAMJ.Numer.Anal.,35 (1998)

    work page 1998

  3. [3]

    BRENIER, W

    Y. BRENIER, W. GANGBO, G. SAVARÉ, M. WESTDICKENBERG,Sticky particle dynamics with interactions, J. Math. Pures Appl.99(2013), 577–617

    work page 2013

  4. [4]

    W. E, YU.G. RYKOV, YA.G.SINAI,Generalized variational principles, global weak solutions and behavior with random initial data for systems of conservation laws arising in adhesion particle dynamics, Comm. Math. Phys,177, No. 2 (1996), 349–380

    work page 1996

  5. [5]

    F.HUANG,Z.WANG,WellPosednessforPressurelessFlow,Comm.Math.Phys.,222,No.1(2001), 117–146

    work page 2001

  6. [6]

    HYND,Lagrangian Coordinates for the Sticky Particle System, SIAM J

    R. HYND,Lagrangian Coordinates for the Sticky Particle System, SIAM J. Math. Anal.51, No. 5 (2019), 3769–3795

    work page 2019

  7. [7]

    HYND,AtrajectorymapforthepressurelessEulerequations,Trans.AMSNo.10(2020),6777– 6815

    R. HYND,AtrajectorymapforthepressurelessEulerequations,Trans.AMSNo.10(2020),6777– 6815

    work page 2020

  8. [8]

    R. HYND, A. TUDORASCU,Asymptotics of the Sticky Particles evolution, Nonlinearity37, No. 7 (2024), 1–22

    work page 2024

  9. [9]

    NATILE, G

    L. NATILE, G. SAVARÉ,A Wasserstein Approach to the One-Dimensional Sticky Particle System, SIAM J. Math. Anal.41, No. 4 (2009), 1340–1365

    work page 2009

  10. [10]

    NGUYEN, A

    T. NGUYEN, A. TUDORASCU,Pressureless Euler/Euler-Poisson Systems via Adhesion Dynamics and Scalar Conservation Laws, SIAM J. Math. Anal.40(2008), 754–775

    work page 2008

  11. [11]

    NGUYEN, A

    T. NGUYEN, A. TUDORASCU,One-Dimensional Pressureless Gas Systems with/without Viscosity, Comm. Part. Diff. Eq.,40, No. 9 (2015), 1619–1665. 46

    work page 2015

  12. [12]

    SUDER, A

    M. SUDER, A. TUDORASCU,On the Lagrangian description and uniqueness for the one- dimensional Pressureless Euler system, Indiana Univ. Math. J., https://www.iumj.indiana.edu/ IUMJ/Preprints/9354.pdf (2021)

    work page 2021

  13. [13]

    TUDORASCU,Sticky Particles with Sticky Boundary, Meth

    A. TUDORASCU,Sticky Particles with Sticky Boundary, Meth. Appl. Anal.30, No. 4 (2023), 129– 148

    work page 2023

  14. [14]

    A.TUDORASCU,R.WASSINK,CompactlysupportedStickyParticlessolutionsdecaytoequilibrium arbitrarily slowly, submitted, 2024

    work page 2024

  15. [15]

    VILLANI, Topics in Optimal Transportation, American Mathematical Society, 2003

    C. VILLANI, Topics in Optimal Transportation, American Mathematical Society, 2003

    work page 2003

  16. [16]

    YA. B. ZELDOVICH,Gravitational instability: an approximate theory for large density perturba- tions. Astro. & Astrophys.5(1970), 84–89. 47

    work page 1970