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

This paper presents an alternative method for solving several important classical mechanics problems by factorizing total mechanical energy using complex numbers, offering new exact and approximate analytical solutions suitable for undergraduate teaching.

The Gist

Imagine you are trying to solve a puzzle. Usually, when physics students try to figure out how a swinging pendulum or a bouncing ball moves, they start with Newton's famous laws. They write down a complicated equation that describes how forces push and pull, and then they have to solve a difficult math problem (a second-order differential equation) to find the answer. For first-year students, this is like trying to climb a steep mountain without a map.

This paper proposes a different, much smoother path up the mountain. The authors, Karlo Lelas and Dario Jukić, suggest a method they call "Energy Factorization." Instead of wrestling with forces and acceleration, they start with the total energy of the system and use a little bit of complex numbers (imaginary numbers) to break the problem apart.

Here is how their approach works, using simple analogies:

The Core Idea: The Energy Split

Think of the total energy of a moving object like a fixed amount of money in a bank account. This money is split between two types of accounts:

1. Kinetic Energy: Money spent on speed (moving fast).
2. Potential Energy: Money saved up in position (like being high up on a hill).

In standard physics, you have to track how the money moves back and forth between these accounts by calculating the speed at every single moment.

The authors say: "Let's look at the total amount of money first." They take the equation for total energy and, using a trick with imaginary numbers (the square root of -1), they split it into two parts that look like a pair of complex conjugates.

The "Phasor" Analogy: A Rotating Clock Hand

Once they split the energy, they introduce a concept called a phase (let's call it $\phi$). You can imagine this as a clock hand rotating on a dial.

• The length of the hand represents the total energy (which stays the same for a perfect, non-damping system).
• The angle of the hand tells you how the energy is currently split.
  • If the hand points straight up, all the energy is "saved" (Potential Energy).
  • If the hand points straight right, all the energy is "spent" on speed (Kinetic Energy).
  • If it's in between, the energy is shared.

By figuring out how fast this clock hand needs to rotate, the authors can instantly write down the position and speed of the object. It's like knowing the time on a clock tells you exactly where the sun is in the sky, without needing to calculate the sun's trajectory from scratch.

What They Solved

Using this "rotating clock hand" method, they derived exact solutions for several classic physics problems that are usually taught with much harder math:

1. The Simple Pendulum (Harmonic Oscillator): They showed how a spring or pendulum swings back and forth. Their method reveals that the "clock hand" rotates at a perfectly steady speed, which is a very intuitive way to understand why the motion is smooth and rhythmic.
2. Throwing a Ball Up (Vertical Projectile): They solved the motion of a ball thrown straight up against gravity. Here, the "clock hand" doesn't rotate at a steady speed; it speeds up and slows down, which perfectly matches how a ball slows down as it rises and speeds up as it falls.
3. Repulsive Forces: They solved a tricky case where a force pushes things away (like two magnets repelling each other), showing how the "clock hand" rotates in the opposite direction.
4. Damped Oscillators (The "Real World" Spring): This is the most impressive part. Real springs lose energy due to friction (air resistance). Usually, this makes the math very messy. The authors showed that even with friction, you can still use this clock-hand idea. The hand gets shorter over time (energy is lost) while it rotates. They found an exact formula for this and even created a simpler, highly accurate approximation for weak friction that is easier to understand than standard textbook methods.

The Limits of the Method

The authors are honest about where this trick doesn't work. It works beautifully for specific types of "energy landscapes" (like springs, gravity, and inverse-square forces). However, if the energy landscape is shaped in a very weird or complex way (like a jagged mountain range), the "clock hand" rotation becomes too complicated to solve with simple math. They note that this isn't a failure of their method; standard physics methods hit the exact same wall with these complex shapes.

They also mention that while they solved the case of "linear friction" (where drag increases steadily with speed), other types of friction (like sliding friction or drag that increases with the square of speed) are harder to solve exactly with this method, though they might still be able to find good approximations.

Why This Matters for Students

The main goal of this paper is educational. The authors argue that this method is perfect for undergraduate students because:

• It avoids the scary, complex calculus usually required to solve Newton's laws.
• It uses basic algebra and the concept of imaginary numbers, which students are already learning.
• It provides a visual, intuitive way to understand energy conservation: the "clock hand" rotating and changing length.

In short, the paper offers a new, elegant way to look at the motion of objects by treating energy not just as a number, but as a rotating vector in a complex plane, making difficult physics problems feel like simple geometry.

Technical Summary

Technical Summary: Energy Factorization in Classical Mechanics

Problem Statement
The paper addresses the pedagogical and analytical challenge of solving fundamental problems in classical mechanics, specifically for systems with positive-definite potential energy. Traditional undergraduate treatments often rely on solving second-order Newtonian differential equations, which can be mathematically demanding for first-year students. Alternatively, solutions for simple harmonic oscillators are frequently derived via geometric analogies to uniform circular motion, while damped systems are often presented without derivation. The authors propose an alternative framework that utilizes the factorization of total mechanical energy using complex numbers to derive exact analytical solutions using only elementary calculus and basic complex number theory.

Methodology
The core methodology involves factorizing the conservation of total mechanical energy equation, $E = \frac{1}{2}mv^2 + U(x)$, where $U(x) \geq 0$.

1. Auxiliary Function: An auxiliary real-valued function $u(x)$ is defined such that $U(x) = u^2(x)$. Note that $u(x)$ is not strictly the positive square root of $U(x)$ but can take negative values to preserve continuity.
2. Complex Factorization: The energy equation is rewritten as a product of complex conjugates:
   $$\left(\sqrt{\frac{m}{2}}v + i u(x)\right) \left(\sqrt{\frac{m}{2}}v - i u(x)\right) = E$$
3. Phase Introduction: The term in the first parenthesis is expressed in polar form as $\sqrt{E}e^{i\phi(t)}$, introducing a time-dependent phase $\phi(t)$. This leads to two coupled first-order relations:
   $$v(t) = \sqrt{\frac{2E}{m}} \cos \phi(t)$$
   $$u(x(t)) = \sqrt{E} \sin \phi(t)$$
4. Solution Strategy: By differentiating the position relation and equating it to the velocity relation, a first-order differential equation for the phase $\phi(t)$ is obtained. Solving this phase equation allows for the determination of $x(t)$ and $v(t)$ without directly integrating Newton's second law.

Key Results and Applications
The authors demonstrate the efficacy of this approach by deriving exact analytical solutions for several canonical systems:

• Simple Harmonic Oscillator (SHO): For $U(x) = \frac{1}{2}kx^2$, the method yields the standard solution $x(t) = \sqrt{2E/k}\sin(\omega_0 t + \phi_0)$ with $\omega_0 = \sqrt{k/m}$. The phase rotates at a constant angular velocity.
• Vertical Projectile Motion: For a homogeneous gravitational field $U(x) = mgx$, the method recovers the standard kinematic equations $x(t) = h + v_0 t - \frac{1}{2}gt^2$. The analysis reveals that the associated "phasor" rotates counter-clockwise with a time-varying angular velocity that diverges at the launch and landing points.
• Repulsive Inverse-Cube Force: For $U(x) = K/(2x^2)$, the method derives the trajectory $x(t) = \sqrt{(K/mx_0^2)t^2 + x_0^2}$. The phasor rotates clockwise with a monotonically decreasing angular velocity.
• Central Force Fields: The approach is extended to effective 1D radial problems. It successfully handles inverse-square effective potentials (including the centrifugal barrier). However, for the Kepler problem ($U(r) \propto 1/r$), while the phase integral is solvable, the resulting expression for $r(t)$ cannot be obtained in closed form, mirroring limitations in standard approaches.
• Damped Harmonic Oscillator: The method is extended to dissipative systems where energy $E(t)$ is time-dependent.
  • Exact Solution: By incorporating the energy dissipation rate $dE/dt = -bv^2$, the authors derive an exact analytical solution for the underdamped regime: $x(t) = x_0 e^{-\gamma t} (\cos \omega t + \frac{\gamma}{\omega} \sin \omega t)$, where $\omega = \sqrt{\omega_0^2 - \gamma^2}$.
  • Weak Damping Approximation: For weak damping ($\gamma \ll \omega_0$), the authors assume $\phi(t) \approx \omega_0 t$ to derive a new approximate analytical solution for position: $x_{approx}(t) = x_0 e^{-\gamma t} (1 + \frac{\gamma}{2\omega_0}\sin(2\omega_0 t))\cos(\omega_0 t)$. This approximation is shown to match the exact solution more closely at turning points than the standard approximation $x_0 e^{-\gamma t} \cos(\omega_0 t)$.

Limitations and Scope
The authors explicitly note the boundaries of their approach:

• Potential Constraints: The method requires positive-definite potential energy ($U(x) \geq 0$). While adaptable to negative-definite potentials using hyperbolic functions, the standard factorization relies on $U(x) \geq 0$.
• Solvability: The approach yields exact closed-form solutions only when the resulting phase integral is elementary. This restricts exact solvability to specific power-law potentials (linear, harmonic, and inverse-square). For generic power-law potentials, the phase integral reduces to elliptic or hypergeometric functions, preventing closed-form solutions for $x(t)$, a limitation shared with standard integration methods.
• Non-Linear Damping: The method provides exact solutions for linear damping ($F_d \propto -v$) but fails to provide exact solutions for Coulomb damping ($F_d \propto -\text{sgn}(v)$) or quadratic damping ($F_d \propto -v^2$) because the energy and phase remain coupled in a way that prevents separation. However, the authors suggest the approximation technique used for weak linear damping could be adapted for these cases.

Significance and Claims
The paper claims that this "energy factorization" approach offers a unified, mathematically simpler alternative to solving Newton's equations of motion for undergraduate physics education.

• Pedagogical Value: It derives exact solutions for important systems using only first-order differential equations and elementary complex number properties, avoiding the need for second-order calculus or geometric analogies to circular motion.
• Novelty in Damping: The derivation of the exact solution for the damped harmonic oscillator and the specific approximate analytical solution for weak damping are presented as novel contributions not typically found in standard literature.
• Physical Insight: The approach provides a "phasor-like" description of energy distribution between kinetic and potential forms, offering a visualizable interpretation of energy flow and phase evolution in both conservative and dissipative systems.

The authors conclude that while the method does not solve all classical mechanics problems (particularly those with non-elementary integrals), it provides a powerful, accessible tool for understanding the dynamics of a wide class of solvable systems and offers new approximate solutions for dissipative systems.