Analysis of Three-Particle Elastic Collisions Using Newtonian Mechanics and Vector Geometry

Shuhei Kobayashi

REVIEW 2 major objections 3 minor 1 references

Analysis of Three-Particle Elastic Collisions Using Newtonian Mechanics and Vector Geometry

T0 review · 2 major / 3 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A rotation angle fixed by the masses determines the total collision count of three elastic point masses on a line, the paper argues.

desk verdict The headline closed-form collision count and velocity formulas are wrong on a simple example; the energy-transfer identity and eigenvalue checks are the only parts that hold up. read the letter →

arxiv 2509.02628 v2 pith:62C4YEUZ submitted 2025-09-01 physics.class-ph

Shuhei Kobayashi This is my paper · ORCID

classification physics.class-ph MSC 70F35

keywords elasticcollisionsthree-particlesystemcollisioncountingHouseholderreflectionmatrixdiagonalizationclosed-formvelocitiesenergytransfer

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper aims to give a complete closed-form description of a one-dimensional system of three point masses that collide elastically, one pair at a time. By writing the two-step collision map (A-B then B-C) as a matrix D and switching to a mass-weighted velocity vector, it finds that D is a rotation by an angle θ fixed only by the masses. That rotation yields explicit formulas for the velocities after any number of collisions, and it turns the total number of collisions into a geometric count: n = 1 + floor((Ω − φ_BC)/θ) + ceil((Ω − φ_BC)/θ), where Ω is the angle between the two collision-reflection lines and φ_BC is the phase after the first impact. If the identity is right, the whole future of the system is determined by masses and initial data without simulating each impact.

What carries the argument

The carrying object is the mass-weighted velocity vector u in s⊥, where s = (√m_A, √m_B, √m_C) is the conserved total-momentum direction. Collisions are Householder reflections—reflections of u⊥ across two lines L_AB and L_BC; the composite D is a rotation by θ fixed by the masses, so counting collisions is counting kθ rotations within a wedge of angle Ω. The phase φ_BC after the first impact is what anchors the count.

What would settle it

An event-driven simulation of a three-mass line—say masses (9,1,4) with initial velocities (0,5,0), A initially right of B right of C—would settle whether (3.22) is correct: compute Ω, φ_BC, and θ from the paper's definitions, count collisions by simulation, and compare. A mismatch where the post-first-collision vector has u·n_AB > 0 would show the trajectory left the wedge W and falsify the counting premise.

Watch

Extended reading notes

Core claim

The central claim is that the composite collision map D = D_BC D_AB, acting on the mass-weighted vector u = (√m_A v_A, √m_B v_B, √m_C v_C), has eigenvalue 1 with conserved direction s = (√m_A, √m_B, √m_C), and a pair of complex-conjugate eigenvalues on the unit circle. Thus, on the plane s⊥, D is exactly a rotation by θ = arg λ_+, with θ a function of the mass ratios alone. Since an A-B collision reflects u⊥ across one line L_AB and a B-C collision reflects across a second line L_BC, the sequence is a billiard in a wedge of angle Ω on s⊥. Counting complete AB→BC cycles is counting how many θ-rotations fit in the angular interval from the post-first-impact direction to the opposite side, givi

Load-bearing premise

The counting formula rests on the assumption that after the first collision the relative velocity stays inside the wedge formed by the two reflection lines and that each completed AB→BC cycle advances the phase by exactly θ; if the post-collision vector ever leaves that wedge, the floor/ceil expression no longer counts the collisions.

Editorial extensions

If this is right

For any prescribed pair of masses and initial velocities, the velocities after the k-th pair of collisions are explicit functions of sin kθ and cos kθ, so no step-by-step event simulation is needed.
The total number of collisions is given in closed form by masses and initial data, making a previously process-dependent quantity a formula.
Because det D = 1, the map is invertible; initial data can be recovered from any post-collision state, consistent with time-reversal symmetry.
Energy changes at each collision are proportional to the total momentum of the colliding pair and to the offset of their average velocity from the initial AB average, and these energy transfers shrink as velocities homogenise.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Beyond the paper, the same wedge-and-rotation picture should apply to any alternating two-boundary reflection process, suggesting a general lemma: the collision count of a two-reflector billiard is a floor+ceil rotation count.
The eigenvalue computation may serve as a prototype for longer one-dimensional elastic chains, where the collision operator could decompose into 1 plus several rotations and the total collision count could become a sum of floor/ceil terms.
Since θ is determined by masses alone, randomly sampled mass triples provide an immediate numerical test of (3.22); any deviation would pinpoint trajectories that leave the wedge W after the first impact.
The energy-transfer proportionality might give a macroscopic shortcut for granular or traffic models, where total collision energy loss could be estimated from pair momentum and mean velocity rather than from individual impacts.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague,

I read this paper expecting a complete closed-form solution to the three-particle 1D elastic collision problem. The abstract makes large claims, but the claims don't survive a direct check. I ran the numbers for masses (9,1,4) with initial velocities (0,5,0), using the paper's own collision equations. The true dynamics gives 3 collisions, with the state after two collisions being (1, 2.4, -1.6). The paper's velocity formula (2.23) at k=1 gives about (14.9, -9.9, -14.8). So the central computational result is simply false.

The same example kills the counting formula. Formula (3.22) is built on a wedge W where both u·n_AB and u·n_BC are negative. That region is never entered: the signs alternate between (AB<0, BC>0) and (AB>0, BC<0). The phase φ_BC cannot count cycles because the trajectory doesn't live where the paper thinks it does. Direct simulation gives 3 collisions; (3.22) does not.

The problems start earlier. The determinant of P is -4di, not -1/(4di), and the printed eigenvectors are not eigenvectors. I checked v+ from (2.6) against D from (1.8) numerically; D v+ ≠ λ+ v+. So the diagonalization that feeds (2.23)-(2.25) is wrong. The energy-transfer identity (4.6) does reproduce the exact algebra of an AB collision, and the eigenvalues γ ± iδ are correct. Those are the only parts I'd trust.

There are no references, so novelty is impossible to assess; the mass-weighted reflection and rotation picture is classical. The paper reads like a motivated but unchecked student derivation.

The stress-test note's concern is valid. Reject. If the author corrects the spectral decomposition and discards the wedge argument, the energy-transfer section and the eigenvalue facts could be the basis of a short exercise, but not a paper as is.

Recommend: desk reject or, if you want to give the author a detailed path, send to a friendly referee who can pin down the exact errors. Either way, not salvageable without rewriting the core.

Referee Report

2 major / 3 minor

Summary. The paper studies one-dimensional elastic collisions of three point masses, writing the collision laws as matrices D_AB and D_BC. It claims to compute D^k via an eigendecomposition, giving closed-form velocities after 2k and 2k+1 collisions. It then introduces a vector u in the plane s^⊥ normal to the momentum-like vector s, models AB and BC collisions as Householder reflections across two lines, and proposes a closed-form total collision count n = 1 + floor((Ω−φ_BC)/θ) + ceil((Ω−φ_BC)/θ). Finally, it derives energy-transfer identities for AB and BC collisions.

Significance. If correct, the paper would provide a complete closed-form solution for the three-particle one-dimensional elastic collision problem, including both post-collision velocities and the total number of collisions, with no fitted parameters. The eigenvalue computation and the energy-transfer identity for a single AB collision are correct and are useful building blocks. However, the two central claims — the explicit D^k formula and the collision-count formula — are invalid, and the latter is contradicted by a simple numerical simulation. Because these are the paper's main results, the contribution as it stands cannot be accepted.

major comments (2)

[Section 3, Eqs. (3.19)–(3.22)] The wedge W = {u·n_AB < 0, u·n_BC < 0} is not the region visited by the trajectory after the first collision. The physical collision conditions alternate: before an AB collision one has u·n_AB < 0 and u·n_BC > 0, while after the AB reflection one has u·n_AB > 0 and u·n_BC < 0; after the subsequent BC reflection the signs are again reversed. No post-first-collision state has both dot products negative. Consequently the phase φ_BC defined in Eq. (3.21) is measured outside the wedge used for counting, and (Ω−φ_BC)/θ does not count completed cycles. Concrete counterexample: for masses (m_A,m_B,m_C)=(9,1,4) and initial velocities (0,5,0), direct iteration of Eqs. (1.1)–(1.4) gives exactly 3 collisions (the sequence is (1,-4,0), then (1,12/5,-8/5), then (32/25,-3/25,-8/5)). However, the post-first-collision vector is outside W, and substituting the resulting φ_BC into Eq. (3.22) does not produ
[Section 2, Eqs. (2.6)–(2.15) and (2.23)] The claimed eigenvectors are not eigenvectors of D, so the spectral decomposition and the closed-form velocity formulas are invalid. For equal masses, D is the 3-cycle permutation matrix with eigenvalues 1, ω, ω^2, and the eigenvector for λ=ω is (1, ω, ω^2) up to scale. Equation (2.6) instead gives (−1+i√3/2, 1−i√3/2, 1), which is not proportional to any eigenvector. As a direct check, take equal masses and initial velocities (0,2,0): D v_0 should be (2,0,0), but Eq. (2.23) with k=1 yields v_A2 = 17/4. This invalidates the claimed ability to compute velocities after any number of collisions, which is used in the counting argument and elsewhere.

minor comments (3)

[Abstract and Section 1] There are numerous grammatical and typographical issues: 'under vacuum' is unnecessary in one-dimensional mechanics, 'Through this concept' is unclear, and the sentence 'all terms letters initial velocities rely on only masses' is ungrammatical. The paper also has no references; it should compare with known treatments of one-dimensional multi-particle collisions.
[Section 3, Figure 1] Figure 1 is referenced but not included in the text I received. If it is present in the actual submission, it should be marked with labels for L_AB, L_BC, Ω, and φ_BC clearly; the phase convention is central and must be unambiguous.
[Section 2, notation] The notation for the index k is inconsistent: n=2k and n=2k+1 are used, but k is said to be a non-negative integer, and later the same letter k appears in v_k and v'_k in Section 4 with a different meaning. This should be cleaned up.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the paper contains no fitted inputs, no self-citations, and no derivation step that reduces to its own conclusion.

full rationale

The paper is entirely self-contained and analytic. Collision velocities are derived directly from the elastic collision laws (1.1)-(1.4); the matrix powers D^k are obtained by explicit diagonalization with stated eigenvalues (2.1)-(2.3) and eigenvectors (2.8), and the resulting closed-form formulas (2.23)/(2.25) are expressed in terms of masses and initial data only. No parameter is fitted to any subset of output data, and no prediction is a renamed input. The collision-count formula (3.22) is derived from a geometric wedge construction and phase definitions (3.19)-(3.21). Even if that construction is mathematically flawed—as the skeptical reading argues—that is an incorrectness or unsupported premise, not circularity. The paper contains no references, so there is no self-citation chain and no imported uniqueness theorem. The central derivations do not assume the results they claim to prove. Therefore the circularity score is 0.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The paper is self-contained analytically: no parameters are fitted to data, and all constants (theta, Omega, phi_BC, nu^2) are computed from the masses and initial velocities. The central claim rests on the physical collision model of Section 1, the classical mass-weighted reflection picture of Section 3 (correct in outline), and two ad-hoc premises specific to this paper: the wedge-region and counting premise of (3.19)-(3.22), which fails on a generic example, and the unproven homogenization assertion of Section 4. No invented entities.

assumptions (5)

domain assumption The dynamics is exactly the alternating sequence of 1D elastic collision formulas (1.1)-(1.4), with no triple collisions and no simultaneous or tangent events.
Section 1. This is the physical model; the paper asserts the AB, BC, AB alternation without proof and never discusses degeneracies such as zero relative velocity at contact, which occur for example for equal masses with initial velocities (0,v,0).
domain assumption The mass-weighted vector u and the plane s-perp fully represent the collision dynamics, with each collision a Householder reflection across L_AB or L_BC.
Section 3, eqs. (3.1)-(3.12). This reduction is classical and in outline correct; an independent reflection computation reproduces the physical collision sequence.
ad hoc to paper The region of motion is the wedge W = {u dot n_AB < 0, u dot n_BC < 0}, and the phase phi_BC defined in (3.21) measures the trajectory position so that (Omega - phi_BC)/theta counts completed cycles.
Section 3, eqs. (3.19)-(3.22). Load-bearing and incorrect: for masses (9,1,4) and velocities (0,5,0) the trajectory never enters W, and the formula yields 0 or 4 collisions instead of the actual 3.
ad hoc to paper Velocities homogenise as collisions proceed, so the energy transferred per collision must decrease.
Section 4, last paragraph. With |beta| = 1, the quantity v_k is oscillatory in k; monotone decrease of the energy transfer does not follow from (4.6)-(4.9) as stated.
standard math D is diagonalizable with the displayed eigenvalues and eigenvectors, and formulas (2.15)-(2.25) are the powers of D.
Section 2. The claim |lambda_plus| = 1 is analytically correct, but the displayed det P = -1/(4di) conflicts with P times P^-1 = I (direct computation gives det P = -4di), so the asserted matrix products are not internally consistent as printed.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Analysis of Three-Particle Elastic Collisions Using Newtonian Mechanics and Vector Geometry." pith.science (2026). https://pith.science/paper/62C4YEUZ

@misc{pith2026250902628,
  author       = {Pith},
  title        = {Pith review of: Analysis of Three-Particle Elastic Collisions Using Newtonian Mechanics and Vector Geometry},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/62C4YEUZ}},
  note         = {Machine review of arXiv:2509.02628}
}

abstract

We study one-dimensional elastic collisions of three point masses on a line under vacuum, with no triple collisions. We express momentum conservation in matrix form and analyze the composite map $D=D_{BC}D_{AB}$ and its powers $D^k$, which yield the velocities after any prescribed number of collisions for arbitrary mass ratios and initial data. After that, using vector $u$ on a plane $s^\perp$, the total number of collisions is \[ n\;=\;1+\Big\lfloor\tfrac{\Omega-\phi_{BC}}{\theta}\Big\rfloor+\Big\lceil\tfrac{\Omega-\phi_{BC}}{\theta}\Big\rceil, \] Through this concept, $D$ is recognised as giving $u$ a rotation with angle $\theta$ which is determined by only mass ratios. And, we calculated energy transfer through collisions. With the work, we find that the change of energy is proportional to total momentum of two particles and average velocity of particles based on initial average velocity of A and B before collision.

Figures

Figures reproduced from arXiv: 2509.02628 by the authors.

**Figure 1.** Wedge in s ⊥ showing LBC , LAB, Ω, and ϕBC . Each AB (resp. BC) collision reflects u⊥ across LAB (resp. LBC). Defining one cycle as AB → BC, applying D = DBCDAB to the initial state advances the phase by θ = θAB +θBC, and the number of completed cycles is Kmax = (Ω − ϕBC)/θ . where ϕBC is a phase of u after first collision of A and B measured from LBC. To calculate this, we introduce a unit vector along a new axi… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

1 extracted references · 1 canonical work pages

[1]

Analysis of Three-Particle Elastic Collisions Using Newtonian Mechanics and Vector Geometry

Analysis of Three-Particle Elastic Collisions Using Newtonian Mechanics and Vector Geometry Shuhei Kobayashi∗ August 31, 2025 Department of Physics and Astronomy, Faculty of Science and Technology, Tokyo University of Science, Noda City, Chiba, Japan 6222041@ed.tus.ac.jp(long-term:xiaolinxiuping1@gmail.com) Abstract We study one-dimensional elastic collis...

work page Pith review arXiv 2025

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.

[1] [1]

Analysis of Three-Particle Elastic Collisions Using Newtonian Mechanics and Vector Geometry

Analysis of Three-Particle Elastic Collisions Using Newtonian Mechanics and Vector Geometry Shuhei Kobayashi∗ August 31, 2025 Department of Physics and Astronomy, Faculty of Science and Technology, Tokyo University of Science, Noda City, Chiba, Japan 6222041@ed.tus.ac.jp(long-term:xiaolinxiuping1@gmail.com) Abstract We study one-dimensional elastic collis...

work page Pith review arXiv 2025