Staggering domino-like blast front motion in a one-dimensional cold gas
Pith reviewed 2026-05-19 18:42 UTC · model grok-4.3
The pith
For an infinite family of special mass ratios, a one-dimensional alternating particle chain with equidistant spacing admits an exact solution where the blast front advances ballistically because only one triplet of particles moves at any时刻.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For an infinite family of mass ratios M_k the dynamics admit an exact solution in which at each moment only a single triplet m, μ, m is in motion, all other particles are at rest, and the shock front moves ballistically with average velocity equal to the initial one. The solution yields explicit formulas for the sequence M_k together with the time-dependent velocities and positions of every particle.
What carries the argument
The staggering triplet mechanism in which collisions propagate sequentially through isolated groups of three particles while all others remain motionless.
Load-bearing premise
The initial positions must be exactly equidistant and the mass ratio must be tuned precisely to one of the discrete values M_k; any deviation from these exact conditions destroys the staggered-triplet regime.
What would settle it
A direct numerical integration with mass ratio exactly equal to one computed M_k and perfectly equidistant initial positions should show only one triplet moving at any time and a front advancing at constant average speed; introducing even tiny random displacements in the initial positions should immediately restore the usual decelerating hydrodynamic front.
Figures
read the original abstract
One-dimensional alternating particle systems are widely used to study interconnections between the hydrodynamics of blast waves in a gas-like medium and the Newtonian dynamics of its corpuscular constituents. We study the model in which point particles with masses $m,\mu, m,\mu,\dots, (m\geq\mu)$ are distributed on the positive half-line $\mathbb{R}_{+}$. Their dynamics are initiated by giving a positive velocity to the leftmost particle; in its course, the particles undergo elastic collisions. For this model with $m/\mu=2$, it has previously been established that the dynamics that start from random initial positions are consistent with predictions based on Euler's hydrodynamic equation. In particular, they have the following properties: (i) the position of the rightmost particle (shock front) evolves as $t^\delta$ with $\delta<1$; (ii) recoiled particles behind the front enter the negative half-axis; (iii) particles with locations $x\leq0$ move ballistically and eventually take over the total energy of the system. In this paper, we present numerical and analytical results for the dynamics of this model with nonrandom (typically equidistant) initial positions and various values of $m/\mu$. For $m/\mu=2$ and equidistant initial positions, our results qualitatively agree with those just mentioned. At the same time, we found an infinite family of numbers $\{\mathcal{M}_k\}$ such that, for $m/\mu=\mathcal{M}_k$, the hydrodynamic behavior mentioned changes drastically to the following. At each moment, only a single triplet $m,\mu, m$ is in motion, whereas all other particles are at rest. As a result, the shock front moves ballistically with an average velocity equal to the initial one. Such a `staggering domino-like' picture is obtained as an exact solution, which yields, in particular, explicit formulas for $\mathcal{M}_k$ and the particle velocities and positions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper examines one-dimensional elastic collisions in an alternating-mass system (m, μ, m, μ, ...) on the positive half-line, initiated by imparting positive velocity to the leftmost particle. For generic mass ratios and random initial positions the dynamics are consistent with Euler hydrodynamics (sub-ballistic shock front, recoils into x<0, energy takeover by leftward particles). The central claim is that for a countable family of special ratios M_k the dynamics admit an exact closed-form solution in which motion is confined at all times to a single isolated triplet m-μ-m while every other particle remains at rest; the resulting shock front therefore advances ballistically with the same average speed as the initial impulse. Explicit algebraic expressions for the sequence M_k together with the successive velocities and positions are derived from the two-body collision map.
Significance. If the construction is correct, the result supplies a rare, fully explicit solvable case inside a model whose generic behavior is hydrodynamic and non-integrable. It demonstrates that the staggered-triplet regime is an exact, non-perturbative solution for discrete mass ratios under perfectly equidistant initial data, thereby furnishing a concrete benchmark against which numerical schemes and hydrodynamic approximations can be tested. The explicit formulas also make the dependence of the front velocity on the mass ratio immediately falsifiable.
major comments (2)
- [derivation of M_k and propagation argument] The manuscript asserts that the staggered pattern propagates indefinitely without out-of-sequence collisions. This claim is load-bearing for the exact-solution statement; an explicit inductive step or closed-form verification that the arrival time of the next triplet always precedes any possible crossing collision should be supplied (e.g., in the section deriving the recursion for M_k).
- [explicit formulas for M_k] The explicit formula for M_k is obtained by solving the finite collision map at each step. It is not immediately clear whether the resulting algebraic expression for general k remains positive and greater than unity for all k, which is required for physical masses; a short table or asymptotic analysis of M_k for large k would confirm this.
minor comments (2)
- [numerical results for m/μ = 2] The abstract states that for m/μ = 2 and equidistant positions the results 'qualitatively agree' with the hydrodynamic picture; a brief quantitative comparison (e.g., measured exponent δ versus the expected hydrodynamic value) would strengthen the contrast with the M_k cases.
- [velocity and position formulas] Notation for the successive velocities after each triplet collision is introduced without a compact summary table; adding such a table would improve readability of the exact solution.
Simulated Author's Rebuttal
We thank the referee for the positive assessment and the recommendation of minor revision. The two major comments identify points where additional explicit verification will improve clarity; we address each below and will incorporate the requested material in the revised manuscript.
read point-by-point responses
-
Referee: [derivation of M_k and propagation argument] The manuscript asserts that the staggered pattern propagates indefinitely without out-of-sequence collisions. This claim is load-bearing for the exact-solution statement; an explicit inductive step or closed-form verification that the arrival time of the next triplet always precedes any possible crossing collision should be supplied (e.g., in the section deriving the recursion for M_k).
Authors: We agree that an explicit inductive argument strengthens the presentation. In the revised version we will insert a new subsection immediately after the recursion for M_k. The base case (k=1) is verified by direct computation of the three-particle collision sequence. For the inductive step we show that if the pattern holds up to triplet k with mass ratio M_k, then the arrival time of the (k+1) triplet at the next stationary particle is strictly earlier than any possible crossing collision; the ordering follows from the velocity inequalities preserved by the two-body map and the recursive definition of M_{k+1}. revision: yes
-
Referee: [explicit formulas for M_k] The explicit formula for M_k is obtained by solving the finite collision map at each step. It is not immediately clear whether the resulting algebraic expression for general k remains positive and greater than unity for all k, which is required for physical masses; a short table or asymptotic analysis of M_k for large k would confirm this.
Authors: We thank the referee for highlighting this point. Although the recursive construction ensures M_k > 1 by design, we will add both a short table listing the first ten explicit values of M_k (all >1) and a brief asymptotic analysis showing that M_k = 1 + O(1/k) as k → ∞, thereby confirming positivity and physical admissibility for every k. revision: yes
Circularity Check
No significant circularity; exact solution constructed directly from collision rules
full rationale
The paper derives an exact staggered-triplet solution for a countable family of mass ratios M_k by iteratively solving the finite elastic collision map under equidistant initial positions. Explicit algebraic formulas for M_k, velocities, and positions are obtained step-by-step from the Newtonian rules without parameter fitting, self-referential predictions, or load-bearing self-citations. The construction is shown to propagate consistently without generating out-of-sequence collisions, making the derivation self-contained against the stated initial conditions and hypotheses. No reduction of the central claim to its own inputs occurs.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Particles interact only via instantaneous elastic collisions; no other forces act between collisions.
- domain assumption Initial positions are exactly equidistant (or otherwise deterministic and non-random).
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel echoes?
echoesECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.
For m/μ=M_k ... only a single triplet m,μ,m is in motion ... explicit formulas for M_k ... v0,s = -q_{2s-1}(ω_k) q_{2s+1}(ω_k) ... α_k = kπ/(2k+1)
-
IndisputableMonolith/Foundation/AlphaDerivationExplicit.leanphi_golden_ratio echoes?
echoesECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.
M2=2+√5 ... staggering domino-like motion
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
-
[1]
using scaling arguments applied to the corresponding hydrodynam ic equation and then confirmed by molecular (Newtonian) dynamics calculations with random initial positions of the particles. The exponents corresponding to other observables, obtained in a similar way, are also listed in Table I. Their physical significance is discussed below. In this work, we...
-
[2]
Note that H(t) = 0 for all t > 0 for the degenerate domino dynamics
to secure equality ∑ l pl(t) = 1, which is standard in this case. Note that H(t) = 0 for all t > 0 for the degenerate domino dynamics. Hence, the higher the value of H(t), the more uniform the distribution of the initial energy between the particles. IV. MOLECULAR DYNAMICS SIMULATIONS Here, we present and discuss the numerical calculation of the afor emen...
-
[3]
By combining the latter with (
and ( 6) H(t) = − 2 m ∑ l≥ 0 mlu2 l (t) log |ul(t)|+ 2 m ∑ l:odd mlu2 l (t) ∼ tα . By combining the latter with (
-
[4]
we obtain ∑ l≥ 0 mlu2 l (t) = m, ∑ l≥ 0 ml|ul(t)| ≥ P x≥ 0(t) ∼ tγ , (11) ∑ l≥ 0 mlu2 l (t) log |ul(t)| ∼ tα . Note that for m = 2, α is much smaller than γ, which means that the increase of the latter sum is much slower than that of ∑ l≥ 0 ml|ul(t)|. 10 3 10 4 10 5 /u1D635 5 × 10 0 6 × 10 0 7 × 10 0 FIG. 3. Time dependence of the Shannon entropy H(t); ...
-
[5]
and ( 12), the velocities after the next round are v0,k +1 v1,k +1 v2,k +1 = 1 0 0 0 − θ 1 + θ 0 1 − θ θ θ 1 − θ 0 1 + θ − θ 0 0 0 1 v0,k v1,k v2,k , which we abbreviate to the following formula Vk+1 = A1A0Vk, V k := v0,k v1,k v2,k , (15) where A0 and A1 describe collisions 0 ↔ 1 and 1 ↔ 2, respe...
-
[6]
Then, in the Appendix B, we solve the recursion relations in ( 32) – ( 34)
– ( 23). Then, in the Appendix B, we solve the recursion relations in ( 32) – ( 34). The first line of ( 13) allows one to exclude v2,k without losing the linearity of the recursion relation in ( 15). To this end, we introduce w0,k = v0,k , w 1,k = v1,k θ + 1 , (A1) and then get v2,k = 1 − w0,k − (1 − θ)w1,k . (A2) Then by ( 16) and ( A1), ( A2), we obtain...
-
[7]
(A24) Now, for a fixed k ≥ 2, let us find α ∈ [π/ 3, π/ 2) (hence m, see ( A24)) such that the following holds sin(2k + 1)α = 0, sin(2s − 1)α sin(2s + 1)α < 0, for all 2 ≤ s ≤ k − 1. (A25) By ( A23), for such α , v0,k = v1,k = 0, v 2,k = 1, and v0,s > 0, for all 2 ≤ s ≤ k − 1, (A26) cf. item ( A) of the staggering domino scenario. By the first condition in (...
-
[8]
Spohn, Kinetic equations from Hamiltonian dynamics: Markovian limits, Rev
H. Spohn, Kinetic equations from Hamiltonian dynamics: Markovian limits, Rev. Mod. Phys. 52, 569 (1980)
work page 1980
- [9]
-
[10]
A. N. Gorban, Hilbert’s sixth problem: the endless road t o rigour, Philosophical Transactions of the Royal Society A: Mathema tical, Physical
- [11]
-
[12]
S. Chakraborti, A. Dhar, and P. L. Krapivsky, A splash in a one-dimensional cold gas, SciPost Phys. 13, 074 (2022)
work page 2022
-
[13]
T. Holovatch, Y. Kozitsky, K. Pilorz, and Y. Holovatch, B reakdown of hydrodynamics in a one-dimensional cold gas, Phys. Rev. E 112, L052101 (2025)
work page 2025
-
[14]
G. I. Taylor, The formation of a blast wave by a very intens e explosion I. Theoretical discussion, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences 201, 159 (1950)
work page 1950
-
[15]
Sakurai, Blast Wave Theory, in Basic Developments in Fluid Dynamics , Vol
A. Sakurai, Blast Wave Theory, in Basic Developments in Fluid Dynamics , Vol. I, edited by M. Holt (Academic Press, New York and London, 1965) Chap. 4, p. 101
work page 1965
-
[16]
Sedov, Similarity and Dimensional Methods in Mechanics (Academic Press, New York and London, 1959) p
L. Sedov, Similarity and Dimensional Methods in Mechanics (Academic Press, New York and London, 1959) p. 363
work page 1959
-
[17]
J. v. Neumann, Collected Works, Volume VI: Theory of Games, Astrophysics, Hydrodynamics and Meteorology (Perga- mon Press, Oxford, 1963) p. 548
work page 1963
-
[18]
Y. B. Zel’dovich and Y. P. Raizer, Physics of Shock Waves and High-Temperature Hydrodynamic P henomena, edited by W. D. Hayes and R. F. Probstein (Dover Publications, Inc., Mi neola, New York, 2002) p. 916
work page 2002
-
[19]
L. Landau and E. Lifshitz, Course of Theoretical Physics, Volume VI: Fluid Mechanics (Academic Press, New York and London, 1987) p. 539
work page 1987
-
[20]
G. I. Taylor, The formation of a blast wave by a very inten se explosion. - II. The atomic explosion of 1945, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences 201, 175 (1950)
work page 1945
-
[21]
M. J. Edwards, A. J. MacKinnon, J. Zweiback, K. Shigemor i, D. Ryutov, A. M. Rubenchik, K. A. Keilty, E. Liang, B. A. Remington, and T. Ditmire, Investigation of u ltrafast laser-driven radiative blast waves, Phys. Rev. Lett. 87, 085004 (2001)
work page 2001
-
[22]
A. D. Edens, T. Ditmire, J. F. Hansen, M. J. Edwards, R. G. Adams, P. Rambo, L. Ruggles, I. C. Smith, and J. L. Porter, Study of high Mach number laser driven blast waves, Physics of Plasmas 11, 4968 (2004)
work page 2004
-
[23]
A. S. Moore, D. R. Symes, and R. A. Smith, Tailored blast w ave formation: Developing experiments pertinent to labora tory astrophysics, Physics of Plasmas 12, 052707 (2005)
work page 2005
-
[24]
J. P. Ostriker and C. F. McKee, Astrophysical blastwave s, Rev. Mod. Phys. 60, 1 (1988)
work page 1988
-
[25]
X. Tang and R. A. Chevalier, Shock evolution in non-radi ative supernova remnants, Monthly Notices of the Royal Astronomical Society 465, 3793 (2017)
work page 2017
-
[26]
Y. Du, H. Li, and L. P. Kadanoff, Breakdown of hydrodynami cs in a one-dimensional system of inelastic particles, Phys. Rev. Lett. 74, 1268 (1995)
work page 1995
- [27]
-
[28]
M. Barbier, D. Villamaina, and E. Trizac, Blast dynamic s in a dissipative gas, Phys. Rev. Lett. 115, 214301 (2015)
work page 2015
-
[29]
J. Joy, S. Pathak, and R. Rajesh, Shock propagation foll owing an intense explosion: Comparison between hydrodynam ics and simulations, Jounal of Statistical Physics 182, 34 (2021)
work page 2021
- [30]
-
[31]
A. Kumar and R. Rajesh, Shock propagation following an i ntense explosion in an inhomogeneous gas: Core scaling and hydrodynamics, Journal of Statistical Physics 192, 17, 18 (2025)
work page 2025
- [32]
-
[33]
Y. G. Sinai, Dynamical systems with elastic reflections , Russian Mathematical Surveys 25, 137 (1970)
work page 1970
-
[34]
R. MacKay and J. Meiss, eds., Hamiltonian Dynamical Systems (Adam Hilger, Bristol and Philadelphia, 1987)
work page 1987
-
[35]
Tabachnikov, Geometry and Billiards , Student Mathematical Library, Vol
S. Tabachnikov, Geometry and Billiards , Student Mathematical Library, Vol. 30 (Americal Mathemat ical Society, 2005) p. 178
work page 2005
-
[36]
L. Bunimovich, D. Burago, N. Chernov, E. Cohen, C. Dettm ann, S. Dorfman, J.R. Ferleger, R. Hirschl, A. Kononenko, J. Lebowitz, C. Liverani, T. Murphy, J. Piasecki, H. Posch, N . Sim´ anyi, Y. Sinai, D. Sz´ asz, T. T´ el, H. Beijeren, R. Zon, J. Vollmer, and L. Young, Hard Ball Systems and the Lorentz Gas , edited by D. Sz´ asz, Encyclopaedia of Mathemati...
work page 2000
-
[37]
P. Grassberger, W. Nadler, and L. Yang, Heat conduction and entropy production in a one-dimensional hard-particle gas, Phys. Rev. Lett. 89, 180601 (2002)
work page 2002
-
[38]
P. I. Hurtado and P. L. Garrido, A violation of universal ity in anomalous Fourier’s law, Scientific Reports 6, 38823, 10 (2016)
work page 2016
- [39]
-
[40]
P. Grassberger, D. Dhar, and P. K. Mohanty, Many univers ality classes in an interface model restricted to non-negat ive heights, Phys. Rev. E 107, 044112 (2023)
work page 2023
-
[41]
P. I. Hurtado, Breakdown of hydrodynamics in a simple on e-dimensional fluid, Phys. Rev. Lett. 96, 010601 (2006)
work page 2006
- [42]
-
[43]
S. Chakraborti, S. Ganapa, P. L. Krapivsky, and A. Dhar, Blast in a one-dimensional cold gas: From Newtonian dynamic s to hydrodynamics, Phys. Rev. Lett. 126, 244503 (2021)
work page 2021
-
[44]
A. Dhar, Heat conduction in a one-dimensional gas of ela stically colliding particles of unequal masses, Phys. Rev. Lett. 86, 3554 (2001)
work page 2001
-
[45]
P. L. Garrido, P. I. Hurtado, and B. Nadrowski, Simple on e-dimensional model of heat conduction which obeys Fourier ’s law, Phys. Rev. Lett. 86, 5486 (2001)
work page 2001
-
[46]
G. Casati and T. Prosen, Anomalous heat conduction in a o ne-dimensional ideal gas, Phys. Rev. E 67, 015203 (2003)
work page 2003
-
[47]
P. Cipriani, S. Denisov, and A. Politi, From anomalous e nergy diffusion to L´ evy walks and heat conductivity in one- dimensional systems, Phys. Rev. Lett. 94, 244301 (2005)
work page 2005
-
[48]
Williams, Algorithm 232: Heapsort, Communications of the ACM 7, 347 (1964)
J. Williams, Algorithm 232: Heapsort, Communications of the ACM 7, 347 (1964)
work page 1964
-
[49]
T. Holovatch, Y. Kozitsky, K. Pilorz, and Y. Holovatch, Effective and asymptotic scaling in a one-dimensional billi ard problem, Condens. Matter Phys. 28, 23401 (2025)
work page 2025
-
[50]
R. Baxter, Exactly Solved Models in Statistical Mechanics , Dover books on physics (Dover Publications, 2007)
work page 2007
-
[51]
W. Sheng and T. Zhang, The Riemann problem for the transp ortation equations in gas dynamics, Memoirs of the American Mathematical Society 137, 1 (1999)
work page 1999
-
[52]
M. Swamy, Generalized Fibonacci and Lucas polynomials , and their associated diagonal polynomials, The Fibonacci Quarterly 37, 213 (1999) . 21
work page 1999
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.