[Physics] Mathematically calculate if a Planet is in Retrograde

astronomyorbital-motionplanets

How do you mathematically calculate if a Planet is in Retrograde when you are given both the Earth and the Planet Keplerian Elements value and rate?

Keplerian Elements are:
a – Semi-Major Axis;
e – Eccentricity;
I – Inclination;
N – Longitude of the Ascending Node;
LP – Longitude of the Perihelion;
ML – Mean Longitude

From these you can find the sun-centered positions and orbits with velocities.

Just looking for a simple equation for checking Retrograde.

Best Answer

What you are asking for is not simple at all. Retrograde motion occurs when the line joining two planets rotates with respect to a fixed coordinate system (or the fixed stars) in the opposite direction as the planets. With both planets in motion in orbits that are not nicely aligned with each other, the times between retrograde motion will only be described by a complicated expression.

If you are numerically simulating the planets' orbits, here's what you can do:

  1. Choose a coordinate system so that the Earth lies in the x-y plane. It's motion will be described by some function of time: $$\vec{x}_{earth} = \begin{bmatrix}x_{earth}(t)\\y_{earth}(t)\\0\end{bmatrix}$$
  2. Define the other planet's motion in the same coordinate system. $$\vec{x}_{planet} = \begin{bmatrix}x_{planet}(t)\\y_{planet}(t)\\z_{planet}(t)\end{bmatrix}$$
  3. The angle between the line joining the planets and the fixed coordinate system is found by subtracting the two vectors and using the two-argument inverse tangent function on the x- and y-coordinates. $$\vec{x}_{planet} - \vec{x}_{earth} = \begin{bmatrix}\Delta x(t)\\\Delta y(t)\\\Delta z(t)\end{bmatrix}$$ $$\theta(t) = \textrm{atan2}(\Delta y(t),\,\Delta x(t))$$
  4. The planet is in retrograde motion with respect to Earth when $$\frac{\textrm{d}\theta}{\textrm{d}t} < 0$$ assuming the planets orbit anticlockwise in the coordinate system.
Related Question