Kepler's Orbits and Special Relativity in Introductory Classical Mechanics

Kepler's orbits with corrections due to Special Relativity are explored using the Lagrangian formalism. A very simple model includes only relativistic kinetic energy by defining a Lagrangian that is consistent with both the relativistic momentum of Special Relativity and Newtonian gravity. The corresponding equations of motion are solved in a Keplerian limit, resulting in an approximate relativistic orbit equation that has the same form as that derived from General Relativity in the same limit and clearly describes three characteristics of relativistic Keplerian orbits: precession of perihelion; reduced radius of circular orbit; and increased eccentricity. The prediction for the rate of precession of perihelion is in agreement with established calculations using only Special Relativity. All three characteristics are qualitatively correct, though suppressed when compared to more accurate general-relativistic calculations. This model is improved upon by including relativistic gravitational potential energy. The resulting approximate relativistic orbit equation has the same form and symmetry as that derived using the very simple model, and more accurately describes characteristics of relativistic orbits. For example, the prediction for the rate of precession of perihelion of Mercury is one-third that derived from General Relativity. These Lagrangian formulations of the special-relativistic Kepler problem are equivalent to the familiar vector calculus formulations. In this Keplerian limit, these models are supposed to be physical based on the likeness of the equations of motion to those derived using General Relativity. The derivation of this orbit equation is approachable by undergraduate physics majors and nonspecialists whom have not had a course dedicated to relativity.