Kepler's laws of planetary motion

Kepler's laws describe three connected features of an elliptic orbit: its shape, how fast it is traversed, and how its size determines its orbital period. Kepler's first law gives an ellipse with the attracting centre at one focus of a conic; Kepler's second law gives equal swept areas in equal times; Kepler's third law relates the squared orbital period to the cubed semi-major axis. Together they turn a drawing of a Kepler orbit into a prediction of where the bodies will be and when they will return. ^[ref]

Here the mathematical model is an isolated two-body problem in Newtonian gravity, with point masses $M,m>0$ and a bound, noncolliding elliptic orbit. Let $\mathbf{r}={\mathbf{r}}_m-{\mathbf{r}}_M$ be the relative displacement, $r=|\mathbf{r}|$ the separation, and $\mu =G(M+m)$ the gravitational parameter, where $G$ is the gravitational constant. Newton's law of universal gravitation then givesThe relative Kepler orbit has semi-major axis $a$ and orbital eccentricity $0\le e<1$. Both bodies move around their shared astronomical barycenter, the system's centre of mass. In its rest frame their displacements are ${\mathbf{r}}_m=M\mathbf{r}/(M+m)$ and ${\mathbf{r}}_M=-m\mathbf{r}/(M+m)$. Their individual semi-major axes therefore add to $a$; the relative ellipse's attracting focus of a conic represents zero separation, rather than a body fixed in space.

One example will connect all three laws. After nondimensionalization, choose $a=1$, $e=0.6$ and $\mu =1$. The semi-minor axis is $b=a\sqrt{1-e^2}=0.8$, while the periapsis distance and apoapsis distance areThus $a=(0.4+1.6)/2=1$ measures the ellipse's size even though the separation changes fourfold. The specific angular momentum will be $0.8$, the two endpoint speeds $2$ and $0.5$, and the orbital period $2\pi$. These are values in the chosen scaled units, not measurements of a particular planet.

Johannes Kepler published Kepler's first law and Kepler's second law in Astronomia Nova in 1609 and Kepler's third law in Harmonices Mundi in 1619. These publication dates are documented by ^[ref]. Their appeal is that three short statements organize both geometry and timing; Newtonian gravity supplies a common dynamical explanation.

The bound relative Kepler orbit is an ellipse with the attracting centre at one focus of a conic. In polar coordinates, its separation is the displayed function of true anomaly $f$, the angle measured from periapsis. For $a=1$ and $e=0.6$, it becomes $r=0.64/(1+0.6\cos f)$, giving $0.4$ at $f=0$ and $1.6$ at $f=\pi$. Setting orbital eccentricity to zero gives a circular orbit. This specifies the path; Kepler's second law specifies how it is traversed.

For an isolated two-body problem in Newtonian gravity, the relative elliptic orbit's orbital period $T$ obeysHere $a$ is the relative semi-major axis, and $M+m$ is the sum of the two masses. The separation $r$ varies during an elliptic orbit, so replacing $a$ by $r$ works only for a circular orbit, where $r=a$. Likewise, either body's individual semi-major axis about the centre of mass must first be converted to the relative one. When $m$ is negligible compared with $M$, the denominator can be approximated by $GM$.

For fixed gravitational parameter, the orbital period depends on the semi-major axis and is independent of orbital eccentricity. Our $a=1$, $\mu =1$ example gives $T=2\pi$ both for $e=0.6$ and for a circular orbit with $e=0$. The eccentric body's rapidly changing speed compensates for its different path. A compact way to see the cancellation uses Kepler's second law: divide the ellipse's area $\pi ab$ by its areal velocity $h/2$. With $b=a\sqrt{1-e^2}$ and specific angular momentum $h=\sqrt{\mu a(1-e^2)}$, the result is $2\pi \sqrt{a^3/\mu }$.

Compare two separate Kepler orbits with the same gravitational parameter, $a_1=1$ and $a_2=4$:With $\mu =1$, their orbital periods are $2\pi$ and $16\pi$; measured in units of $T_1$, these are $1$ and $8$. Quadrupling the size therefore takes eight times as long to complete a circuit. This power law makes comparisons possible without first choosing physical units.

The same relation can be inverted to weigh a system: $M+m=4{\pi }^2{a}^3/(G{T}^2)$. An observed relative semi-major axis and orbital period determine the total mass, while additional information is needed to separate the two masses. This is useful for binary stars, as discussed in NASA, Orbits and Kepler's Laws: measuring masses. Geometry and a clock together become a weighing instrument.

As a heuristic approximation, inserting $r=6GM/c^2$ into the circular orbit form of Kepler's third law, with $M$ here denoting the total mass, givesThe dominant quadrupole gravitational-wave frequency is twice the orbital frequency, $f_{\mathrm{GW}}=2/T$.

The relative displacement of a Kepler orbit sweeps equal areas about the attracting focus of a conic in equal times. A short interval near periapsis covers a wide angle at small separation; the same interval near apoapsis covers a narrow angle at large separation. The swept areas are equal even though the travelled distances are different. In particular, the sectors have curved ellipse boundaries: triangles formed by joining their endpoints only approximate the swept areas for short intervals.

Let $h=|\mathbf{r}\times \dot{\mathbf{r}}|$ be the relative specific angular momentum. The areal velocity isA central force exerts zero torque about its centre, so conservation of angular momentum makes this rate constant. This explanation also applies to other central forces; the ellipse shape specifically comes from Newtonian gravity.

In the $a=1$, $e=0.6$, $\mu =1$ example, specific angular momentum is $h=\sqrt{\mu a(1-e^2)}=\sqrt{0.64}=0.8$, giving areal velocity $0.4$. At periapsis and apoapsis, the velocity is perpendicular to the relative displacement, so $h=rv$. The endpoint speeds are consequentlyThe fourfold speed difference is consistent with the fourfold separation difference. Away from these endpoints, $h/r$ gives only the transverse component of velocity; the full speed follows from the vis-viva equation.

To actually place equal-time markers, use Kepler's equation. Write mean anomaly as $\mathcal{M}=nt$, with $t=0$ at periapsis and mean motion $n=\sqrt{\mu /a^3}=1$. The different symbol $\mathcal{M}$ avoids confusing this angle with the mass $M$. Solve for eccentric anomaly $E$ in radians, then compute the relative position:This timing procedure is explained in NASA educational archive, How Orbital Motion is Calculated.

For the first interval $0\le t\le \pi /6$, the final eccentric anomaly is $E\approx 1.041495$. For the equal interval $\pi \le t\le 7\pi /6$ starting at apoapsis, it ends at $E\approx 3.471065$. Both swept areas are $\Delta A=(h/2)\Delta t=0.4\pi /6=\pi /15\approx 0.209440$, one twelfth of the entire ellipse's area $0.8\pi$. The diagram computes its markers from Kepler's equation, making visible how uniform time steps become uneven spatial steps.