Relativistic motion in a Coulomb field

Live Script exploring motion of a relativistic particle in an inverse-square law force field.
2 Descargas
Actualizado 2 oct 2025

Ver licencia

It is well known that, in non-relativistic mechanics, only the harmonic oscillator and inverse square force law potentials give rise to closed (elliptical) orbits while perturbations to these potentials result in precession of the orbital axes. This educational Live Script solves the relativistic form of Newton's 2nd Law for the motion of a particle in Coulomb force field and illustrates precession and other unusual behavior.
The script reproduces and enhances trajectory plots from the Master's Thesis by Angelo Di Fiore [2], which are based on the analytic solutions presented by T.H. Boyer [1]. The script visualizes six distinct cases of relativistic orbits, including scattering, bound, and spiral trajectories. For each case, the analytic time-dependence is also computed and validated using a numerical integration of the relativistic equations of motion.
The script also compares the relativistic and non-relativistic trajectories for a beam of particles for both attractive and repulsive Coulomb potentials. In the non-relativistic case, no trajectory enters a halo region behind the scattering center, a region bounded by a quadratic curve. In the relativistic case, the halo region shrinks considerably.
This script may interest students and instructors in physics and is suitable for a course in intermediate level mechanics. It contains 'Try this' suggestions, coding 'Challenges", hyperlinks, and references for further exploration. Additional educational Live Scripts by the author are here.

Citar como

Duncan Carlsmith (2025). Relativistic motion in a Coulomb field (https://es.mathworks.com/matlabcentral/fileexchange/182191-relativistic-motion-in-a-coulomb-field), MATLAB Central File Exchange. Recuperado .

Compatibilidad con la versión de MATLAB
Se creó con R2025b
Compatible con cualquier versión
Compatibilidad con las plataformas
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

RelativisticMotionInACoulombField folder

Versión Publicado Notas de la versión
1.0.0