Dynamics and Stability of the Two-Body Problem with Yukawa Correction to Newton’s Gravity, Revisited and Applied Numerically to the Solar System

Abstract

In this manuscript, we review the motion of a two-body celestial system (planet–sun) for a Yukawa-type correction on Newton’s gravitational potential using Hamilton’s formulation. We reexamine the stability using the corresponding linearization Jacobian matrix, and verify that the conditions of Bertrand’s theorem are met for radii ≪1015 m, meaning that bound closed orbits are expected. Applied to the solar system, we present the equation of motion of the planet, then solve it both analytically and numerically. Making use of the analytical expression of the orbit, we estimate the Yukawa strength α and find it to be larger than the nominal value (10−8) adopted in previous studies, in that it is of order (α=10−4−10−5) for the terrestrial planets (Mercury, Venus, earth, Mars, and Pluto) and even larger (α=10−3) for the giant planets (Jupiter, Saturn, Uranus, and Neptune). Taking the inputs (rmin,vmas,e) observed by NASA, we analyse the orbits analytically and numerically for both the estimated and nominal values of α and determine the corresponding trajectories. For each obtained orbit, we recalculate the characterizing parameters (rmin,rmax,a,b,e) and compare their values according to the potential (Newton with/without Yukawa correction) and method (analytical and/or numerical) used. When compared to the observational data, we conclude that the path correction due to Yukawa correction is on the order of up to 80 million km (20 million km) as the maximum deviation occurring for Neptune (Pluto) for a nominal (estimated) value of α.