Pith. sign in

REVIEW 2 major objections 4 minor 20 references

An alternative approach to several important systems in classical mechanics: energy factorization

T0 review · 2 major / 4 minor · reviewed 2026-08-03 · deepseek-v4-flash

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

Pith's one-line read The paper shows that factoring total energy as a product of complex conjugates, with a time-dependent phase, yields exact solutions for the harmonic oscillator, vertical projectile, inverse-cube force, and linearly damped oscillator, plus a

desk verdict Core derivations hold up, but the new approximate solution (56) has a first-order residual and should be reframed; otherwise this is a solid teaching paper. read the letter →

arxiv 2601.18957 v2 pith:4AWQXR2Z submitted 2026-01-26 physics.class-ph

Karlo Lelas , Dario Jukić This is my paper · ORCID
classification physics.class-ph
keywords energyfactorizationcomplexnumberssimpleharmonicoscillatordampedweakdampingapproximationphaseintegralundergraduateteachingclassicalmechanics
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 attempts to establish that the conservation of mechanical energy, written as a complex factorization, is enough to derive the full dynamics of several standard one-dimensional systems. By introducing a phase φ(t) through v = √(2E/m) cosφ and √U = √E sinφ, the authors turn energy conservation into first-order phase equations that avoid solving Newton's second-order equation. The approach recovers the known solutions for the simple harmonic oscillator, vertical projectile motion, and repulsive inverse-cube motion, and it exactly reproduces the standard linearly damped oscillator solution. In the weak-damping limit it produces an approximate energy decay and a new approximate position formula that matches the exact solution better at turning points than the usual e^{-γt} cos(ω0t). If correct, the paper offers a pedagogically simpler route to these results and a new approximation for weakly damped motion.

What carries the argument

The central object is the complex factorization of total energy and the associated phase φ. Starting from E = mv²/2 + U(x) = (√(m/2)v + i√U)(√(m/2)v − i√U), the authors define φ by √(m/2)v(t) + i√U(x(t)) = √E e^{iφ(t)}, which yields v = √(2E/m) cosφ and √U = √E sinφ. These two equations convert energy conservation into a first-order differential equation for φ; whenever the resulting integral is elementary, x(t) follows. For the damped oscillator, the same phase carries the energy decay through dE/dt = −bv², and the phase integral (37) can be solved by the substitution u = tanφ, leading to the exact solution.

What would settle it

Take the weak-damping approximation (56) and the exact solution (48) for γ/ω0 = 0.5 and compare them over several periods near the turning points; if the fractional error of (56) exceeds a few percent while the usual e^{-γt}cos(ω0t) remains within 1%, the claim that (56) is an excellent weak-damping approximation would be refuted.

Watch

Extended reading notes

Core claim

The central claim is that for potentials U(x) ≥ 0, the identity E = (√(m/2)v + i√U)(√(m/2)v − i√U) lets one introduce a real phase φ(t) such that √(m/2)v(t) + i√U(x(t)) = √E e^{iφ(t)}. Splitting real and imaginary parts gives v = √(2E/m) cosφ and √U = √E sinφ. For any system where inserting the specific U and differentiating x(φ) leads to an elementary integral for dφ/dt, this yields exact closed-form x(t). The paper shows this for U ∝ x², U ∝ x, and U ∝ x^{-2}, and for the damped oscillator it keeps the same factorization with time-dependent E(t), obtaining the exact phase integral and solution, plus a weak-damping approximation x_approx(t) = x0 e^{-γt}(1 + γ/(2ω0) sin(2ω0t)) cos(ω0t) that

Load-bearing premise

The method requires the potential energy to remain non-negative along the entire trajectory, because the factorization uses the real square root √U; if a trajectory enters a region where U < 0, the phase equations no longer apply.

Editorial extensions

If this is right

  • The simple harmonic oscillator solution follows in a few steps from energy conservation alone, without invoking uniform circular motion, giving teachers a new way to present the topic.
  • Vertical projectile motion reduces to substituting the known constant-acceleration velocity into the energy equation, bypassing integral calculus.
  • For inverse-cube repulsive potentials, the phase integral is elementary and yields x(t) = sqrt((K/(m x0²)) t² + x0²) for a start from rest.
  • The damped oscillator's exact solution is recovered from the phase integral, and the same factorization yields the weak-damping energy approximation E0 e^{-2γt}(1 + (γ/ω0) sin(2ω0t)) and a new position approximation that improves turning-point matching.
  • The approach clarifies why generic power-law potentials U ∝ x^n are not solvable in closed form: the phase integral reduces to an incomplete elliptic integral or hypergeometric function for generic n.

Reading between the lines

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

  • A natural extension, which the paper only gestures at, is to apply the same phase-ansatz approximation to sliding friction and quadratic drag; the paper notes the energy-phase coupling prevents exact solutions, but the weak-damping approximation strategy could yield new formulas.
  • The method is essentially a factorization of the Hamiltonian into action-angle-like variables; in systems with more degrees of freedom, a similar complex factorization might connect to normal-mode decomposition, though the paper does not pursue this.
  • The improved turning-point match of the approximate position suggests that the amplitude envelope should be corrected by the factor (1 + γ/(2ω0) sin(2ω0t)); this correction is a prediction that could be verified in a laboratory damped spring-mass experiment with γ/ω0 ≲ 0.1.
  • The pedagogical claim could be tested directly: if the phase-based derivation is taught to an undergraduate cohort, time-to-solution and error rates could be compared against the standard circular-motion derivation.
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.

Referee Report

2 major / 4 minor

Summary. The paper proposes an alternative method for solving one-dimensional classical mechanics problems by factorizing the total mechanical energy as (sqrt(m/2)v + i sqrt(U))(sqrt(m/2)v - i sqrt(U)) and introducing a phase phi(t) via equations (5) and (6). The method is applied to the simple harmonic oscillator, vertical projectile motion, the repulsive inverse-cube force, and the linearly damped harmonic oscillator. For the weakly damped oscillator, the paper derives an approximate energy decay (54) and an approximate position (56), which is claimed to be a new approximate analytical solution. The paper also discusses limitations, including power-law potentials and nonlinear damping.

Significance. The exact derivations in Sections II, III, IV, and VI A are sound and reproduce standard textbook results through an elementary and potentially useful pedagogical route. The derivation of the approximate energy (54) is simpler than in previous work and is correct. If the weak-damping position (56) were a genuine first-order approximate solution, the paper would offer a modest but useful contribution to undergraduate teaching. However, the main novelty claim rests on (56), and that claim is not supported by the analysis presented. The paper's exact results are correct and its limitations section is honest, but the central new approximate solution requires either correction or a clear reframing.

major comments (2)
  1. [Section VI B, Eqs. (51)-(56)] The function (56) is not an approximate solution of the damped oscillator at first order in gamma/omega0. Writing (56) as x_ap = x0 e^{-gamma t}[cos(omega0 t) + (gamma/4 omega0)(sin 3 omega0 t + sin omega0 t)], one obtains L[x_ap] = -gamma^2 x0 e^{-gamma t} cos(omega0 t) - 2 gamma omega0 x0 e^{-gamma t} sin 3 omega0 t + O(gamma^2). The sin 3 term is O(gamma omega0), not O(gamma^2); for gamma/omega0 = 0.1 this residual is about 20% of omega0^2 x at t = pi/(4 omega0). The standard first-order truncation of the exact solution (48) is x0 e^{-gamma t}[cos(omega0 t) + (gamma/omega0) sin(omega0 t)], which has residual O(gamma^2). The difference between (56) and that first-order solution is -(gamma/omega0) x0 e^{-gamma t} sin^3(omega0 t), i.e. O(gamma/omega0). Therefore calling (56) an 'approximate analytical solution' is unsupported; at best it is a curve matched to the exact turning-point enve
  2. [Section VI B, Eqs. (51)-(56)] The paper explicitly omits a validity analysis of the energy approximation (54) and refers to Ref. [12]; that is acceptable for (54). However, the position approximation (56) is new, and its only evidence is Fig. 1, which compares curves rather than testing the equation of motion. The phase approximation phi(t) ≈ omega0 t + phi0 is introduced without any quantitative error bound. A curve comparison does not establish that (56) is an approximate solution. Please provide an analytic error estimate (e.g. a bound on the equation residual or on |x_ap - x_exact|) and state the range of gamma/omega0 for which (56) is intended to be accurate. Without this, the abstract's 'new approximate analytical solution' and the conclusion's 'excellent approximation' are not established.
minor comments (4)
  1. [Sections III and VIII] The statement that the approach 'completely bypasses solving Newton's equations of motion' is somewhat overstated: Eq. (13) is the constant-acceleration kinematic solution, and the energy dissipation rate dE/dt = -b v^2 used in Section VI A is a consequence of Newton's second law. The wording could be softened.
  2. [Section III] The factorization (2) requires U(x) >= 0. For the vertical projectile with U(x) = mgx, the derivation as written assumes x >= 0 (or that the solution is used only above the reference level). This restriction is not stated and could confuse students applying the result to negative heights.
  3. [Section V, Eq. (29)] The typeset form of Eq. (29) is nearly unreadable. The identity is U_eff = [L/(sqrt(2m) r) + (k/L) sqrt(m/2)]^2 - m k^2/(2 L^2); please typeset it cleanly so that dimensions and terms are clear.
  4. [Section VI B] In the text after Eq. (53), the expansion exp[(gamma/omega0) sin(2 omega0 t)] ≈ 1 + (gamma/omega0) sin(2 omega0 t) is used. It would be helpful to state explicitly that this is valid for gamma/omega0 << 1 uniformly in t, and to mention that this is a separate approximation from the phase approximation.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the exact and approximate derivations are self-contained; the only minor issue is a self-citation used to justify omitting a validity check for the energy approximation.

full rationale

The paper's derivation chain starts from the energy identity (1), factorizes it as (2), and introduces the phase through (3)-(4). For each conservative example, x(t) is expressed via U(x)=E sin^2 phi, and comparing the time derivative of that expression with v = sqrt(2E/m) cos phi yields a first-order ODE for phi (e.g. (9), (21)) that is then integrated without importing the final solution. The damped-oscillator section is likewise self-contained: Eqs. (31)-(33) are just the energy-conservation identities combined with the independent dissipation law dE/dt = -b v^2; solving the phase integral (37) gives (38) and leads to the standard solution (48). No parameter is fitted to data, and no target result appears in an assumption. The weak-damping result (56) follows from an explicit analytic ansatz phi ≈ omega0 t + phi0, not from the solution one is trying to predict. The only self-citation, Ref. [12] by co-author Lelas, is used to support the claim that (54) overlaps the exact energy and to justify omitting a validity analysis; however, (54) is re-derived here in Eqs. (51)-(54), and the new solution (56) is checked against the exact solution in Fig. 1. Thus the self-citation is not load-bearing for the central derivation. The skeptic's residual/O(gamma omega0) concern about Eq. (56) failing the damped equation of motion is a correctness/accuracy issue, not a circularity, and does not change this assessment.

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

The paper introduces no new free parameters fitted to data: all constants (m, k, g, K, b, E) are physical inputs. The phase function φ(t) is a mathematical reparameterization, not a new physical entity. The main burden on the reader is accepting the standard energy and dissipation laws plus the ad hoc weak-damping phase assumption for the approximate section.

assumptions (5)
  • domain assumption Total mechanical energy is conserved for conservative forces and is given by Eq. (1).
    Starting point for all conservative examples; the paper cites the work-energy theorem as the derivation.
  • domain assumption The potential energy U(x) is positive definite so that sqrt(U(x)) is real.
    Required for the complex factorization (2); explicitly restricted in the introduction and limitations.
  • domain assumption For the damped oscillator, the mechanical energy decays as dE/dt = -b v^2.
    Standard consequence of the linear damping force F_d = -bv; used to close the phase ODE in Section VI A.
  • ad hoc to paper For weak damping, the phase can be approximated as φ(t) ≈ ω0t + φ0.
    Basis for the approximate energy (53)-(54) and the new approximate solution (56); validated only numerically, not by an error bound.
  • standard math Elementary calculus operations (separation of variables, substitution u = tanφ, inverse trigonometric integration) are valid.
    Used to solve the phase equations in Sections II, IV, VI A and Appendix A.

how reviews work

0 comments
Cite this review

Pith. "Pith review of An alternative approach to several important systems in classical mechanics: energy factorization." pith.science (2026). https://pith.science/paper/4AWQXR2Z

@misc{pith2026260118957,
  author       = {Pith},
  title        = {Pith review of: An alternative approach to several important systems in classical mechanics: energy factorization},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4AWQXR2Z}},
  note         = {Machine review of arXiv:2601.18957}
}
read the original abstract

We show how several important classical problems, with positive definite potential energy, can be solved by starting from the factorization of the total mechanical energy using complex numbers. In particular, we derive in a new way exact analytical solutions for: simple harmonic oscillator, vertical projectile motion, motion under a repulsive inverse cube force, and damped harmonic oscillator (with linear damping). We also show how this approach easily yields an excellent approximation of the energy decay and a new approximate analytical solution in the case of a weakly damped harmonic oscillator. Our derivations are suitable for undergraduate physics teaching as an alternative to solving Newton's equations of motion. In addition, we comment on the limitations of our approach, but also on the insights it provides and opportunities for further research.

Figures

Figures reproduced from arXiv: 2601.18957 by the authors.

Figure 1
Figure 1. FIG. 1: Solid red curve corresponds to the approximate solution [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

20 extracted references · 7 canonical work pages

  1. [12]

    F. S. Crawford,Waves : Berkeley physics course / volume 3(McGraw-Hill, New York, USA, 1968)

    1968

  2. [1]

    Implementation of cutting-edge research and its application as part of the Scientific Center of Excellence for Quantum and Complex Systems, and Representations of Lie Algebras

    andϕ 0 =π/2. Thus, cotϕ 0 = 0 in (23). Relation (18) becomes x(t) =x 0 sin arctan x2 0 r m K t−1 −1 .(24) We can use identity sin[arctan(α)] =α/ √ 1 +α 2 in (24), and obtain x(t) = s K mx2 0 t2 +x 2 0 .(25) V. CENTRAL FORCE FIELD: EFFECTIVE POTENTIAL ENERGY It is well known that for a particle of massmmoving in a three-dimensional central force field, i.e...

    2021

  3. [2]

    Cutnell and K

    J. Cutnell and K. Johnson,Physics(John Wiley & Sons, 2009), ISBN 9780470223550, URLhttps://books.google.hr/ books?id=en1sBgAAQBAJ

    2009

  4. [3]

    H. D. Young, R. A. Freedman, A. L. Ford, and F. W. Sears,Sears and Zemansky’s university physics: with modern physics (Pearson Addison Wesley, San Francisco, 2004)

    2004

  5. [4]

    Halliday, R

    D. Halliday, R. Resnick, and J. Walker,Fundamentals of Physics(John Wiley & Sons, 2013)

    2013

  6. [5]

    Gauthier, Int

    N. Gauthier, Int. J. Math. Educ. Sci. Technol.35, 446–452 (2004), URLhttps://doi.org/10.1080/ 00207390410001686580

    2004

  7. [6]

    C. C. Tisdell, Int. J. Math. Educ. Sci. Technol.50, 950–959 (2019), URLhttps://doi.org/10.1080/0020739X.2018. 1516826

    work page doi:10.1080/0020739x.2018 2019

  8. [7]

    J. Tran, L. Doughty, and J. K. Freericks, American Journal of Physics93, 437–440 (2025), URLhttps://doi.org/10. 1119/5.0220797

    2025

Show all 20 references
  1. [8]

    D. J. Griffiths,Introduction to Electrodynamics(Cambridge University Press, 2017), 4th ed

    2017

  2. [9]

    Goldstein, C

    H. Goldstein, C. P. Poole, and J. L. Safko,Classical Mechanics(Addison-Wesley, 2002), 3rd ed

    2002

  3. [10]

    D. J. Morin,Introduction to classical mechanics: with problems and solutions(Cambridge University Press, Cambridge, UK, 2008)

    2008

  4. [11]

    P. A. Dourmashkin,Classical Mechanics: MIT 8.01 Course Notes(Wiley Custom Learning Solutions, 2014)

    2014

  5. [13]

    Lelas and R

    K. Lelas and R. Pezer, European Journal of Physics46, 015004 (2024), URLhttps://doi.org/10.1088/1361-6404/ ad9559

    2024 doi

  6. [14]

    I. R. Lapidus, American Journal of Physics38, 1360–1361 (1970), URLhttps://doi.org/10.1119/1.1976111. 12

    1970 doi

  7. [15]

    Marchewka, D

    A. Marchewka, D. S. Abbott, and R. J. BeichnerKamela, American Journal of Physics72, 477–483 (2004), URLhttps: //doi.org/10.1119/1.1624113

    2004 doi

  8. [16]

    Anastasios Adamopoulosa and N

    A. Anastasios Adamopoulosa and N. Adamopoulos, International Journal of Mathematical Education in Science and Technology53, 3151 (2022), URLhttps://doi.org/10.1080/0020739X.2021.1954253

    2022

  9. [17]

    Kamela, The Physics Teacher45, 110 (2007), URLhttp://dx.doi.org/10.1119/1.2432089

    M. Kamela, The Physics Teacher45, 110 (2007), URLhttp://dx.doi.org/10.1119/1.2432089

    2007 doi

  10. [18]

    B. R. J. Smith, American Journal of Physics80(2012), URLhttps://doi.org/10.1119/1.4729440

    2012 doi

  11. [19]

    C. E. Mungan and T. C. Lipscombe, European Journal of Physics34, 1243–1253 (2013), URLhttps://doi.org/10.1088/ 0143-0807/34/5/1243

    2013

  12. [20]

    X. Wang, C. Schmitt, and M. Payne, European Journal of Physics23, 155 (2002), URLhttps://iopscience.iop.org/ article/10.1088/0143-0807/23/2/309

    2002 doi

Pith tools

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