Relativistic Orbits

Published:

ME275 is the first grad class I have taken at Cal and the third class I have taken with Professor O’Reilly. ME275’s course title is Advanced Dynamics and delves into Lagrangian and Hamiltonian equations of motion. The course culminates into an individual project that further explores a concept previously learned in class.

After being exposed to the ideas of relativity in the class, Professor O’Reilly suggested I compare the Mercury’s orbit in the relativistic and non-relativistic case. (I didn’t really prepare another idea :P) For the uninitiated, non-relativistic is grounded in Newtonian gravity and that it is a force of attraction between two massive objects. The relativistic case is rooted in that planetary orbits are a result of the geometry of space-time, warped by gravity.

After setting forth the Lagrangian equations of motion for the two body system and applying a coordinate transformation, I plotted the system’s relativistic and non-relativistic phase portrait using pplane to analyze their orbits.

Standard Newtonian orbit phase portrait
Relativistic orbits near the origin
Relativistic orbits far away from the origin

From left to right: Newtonian phase portrait, relativistic orbits near and far away from the origin.

The non-relativistic orbits form closed loops, which is expected in order to produce the physical orbits (from circular to hyperbolic) we observed in the universe. However, the relativistic phase portraits tell us a slightly different message. In the far view, the phase portraits look similar in the sense that there are closed loops encircling a red dot, the equilibrium distance. However, for some initial conditions in the relativistic case, such as being close to the origin (0,0), the phase portrait goes haywire and demonstrates a behavior to converge on (1,0), as seen in the left top portrait. This strange behavior is actually a manifestation of a black hole, as predicted by Einstein’s theories. For initial conditions that pass beyond the vertical line at x = 1, it would implicate that the object has slipped beyond the event horizon of the black hole and would be lost forever. How’s that for cool?

Integrating the results of the aforementioned phase portraits, I reconstructed the Mercury’s physical orbits for the non-relativistic and relativistic case. In contrast, the non-relativistic orbits, depending on your initial conditions with origin at (0,0), circular orbits (blue) can transform into more eccentric orbits (cyan) until finally it becomes parabolic and hyperbolic orbits (green and magenta, respectively). We can observe this in Earth, which has an orbit with low eccentricity and for comets, which has orbits that are highly hyperbolic.

Integrated orbits in the relativistic case
Integrated orbits in the standard Newtonian case

From left to right: Representative integrated orbits for the relativistic case and standard Newtonian case.

For the relativistic case, I observed that changing the initial conditions generated a response characterized above. Irresistibly pretty, Mercury’s orbit precesses around the origin (0,0) as time goes on. In green, Mercury’s orbit for one Mercury year is shown. Six Mercury years later, Mercury’s orbit (magenta) will have precessed and be physically in a different place in space. This behavior was one of the supporting evidence for Einstein’s theories because Newtonian gravity couldn’t have explained this phenomenon.

For a complete analysis of this work, which I have summarized above, you may download my paper . All figures were created in MATLAB.