Johannes Kepler

Johannes Kepler initiated the transformation of observational cosmology into quantitative celestial mechanics. Before Kepler, astronomy was predominantly kinematic and instrumentalist: Ptolemaic and Copernican models sought to "save the appearances" through geometric constructions (deferents, epicycles, equants) without asserting physical causality. Kepler established that planetary trajectories are the necessary consequences of an underlying physical force (vis motrix) centered in the Sun.

His architecture rests upon three structural interventions: 1. The Rupture with Uniform Circular Motion: Demonstrating that planetary orbits are non-circular, breaking an unbroken two-millennium Aristotelian paradigm by isolating an 8-arcminute residual in Tycho Brahe's telemetry. 2. Conservation of Areal Velocity: Discovering that the vector connecting the Sun to a planet sweeps equal areas in equal times, which Newton later formalized as the conservation of orbital angular momentum under a central force field. 3. The Harmonic Invariant: Discovering the scale-invariant power-law link ($T^2 \propto a^3$) coupling the orbital period and semi-major axis across all planetary bodies.

2.1 First Law: The Elliptic Orbit

In polar coordinates $(r, \theta)$ with the Sun situated at the origin (one focus), an elliptical trajectory is governed by: $$r(\theta) = \frac{p}{1 + e \cos \theta} = \frac{a(1 - e^2)}{1 + e \cos \theta}$$ where: - $a$ is the semi-major axis. - $b = a\sqrt{1 - e^2}$ is the semi-minor axis. - $e = \frac{\sqrt{a^2 - b^2}}{a} \in [0, 1)$ is the orbital eccentricity. - $p = a(1 - e^2)$ is the semi-latus rectum. - $\theta$ is the true anomaly (polar angle measured from perihelion).

Distance at perihelion: $r_{\text{min}} = a(1 - e)$.
Distance at aphelion: $r_{\text{max}} = a(1 + e)$.

2.2 Second Law: Area Law & Angular Momentum Conservation

The differential area $dA$ swept out by the position vector $\mathbf{r}(t)$ in differential time $dt$ is given by: $$dA = \frac{1}{2} \|\mathbf{r} \times d\mathbf{r}\| = \frac{1}{2} r^2 \frac{d\theta}{dt} \, dt$$

The rate of areal sweep is strictly constant: $$\frac{dA}{dt} = \frac{1}{2} r^2 \dot{\theta} = C = \text{constant}$$

In classical mechanics, the specific angular momentum $\mathbf{h} = \mathbf{r} \times \mathbf{v}$ has magnitude: $$h = \|\mathbf{r} \times \mathbf{v}\| = r^2 \dot{\theta} = 2 \frac{dA}{dt}$$ Because the gravitational force $\mathbf{F} = -\frac{GMm}{r^2} \hat{\mathbf{r}}$ is strictly central: $$\boldsymbol{\tau} = \mathbf{r} \times \mathbf{F} = \mathbf{0} \implies \frac{d\mathbf{h}}{dt} = \mathbf{0}$$ Kepler's Second Law is therefore mathematically identical to the conservation of orbital angular momentum.

2.3 Kepler’s Equation & Orbital Time-Stepping

To find a planet's position as an explicit function of time $t$, Kepler introduced the Eccentric Anomaly $E$. The relationship between mean anomaly $M$ (proportional to elapsed time) and $E$ is given by the transcendental Kepler Equation: $$M = \frac{2\pi}{T}(t - t_0) = E - e \sin E$$

The true anomaly $\theta$ is recovered via the half-angle relation: $$\tan\left(\frac{\theta}{2}\right) = \sqrt{\frac{1 + e}{1 - e}} \tan\left(\frac{E}{2}\right)$$

Because $M = E - e \sin E$ cannot be inverted analytically for $E$, Kepler devised the first numerical root-finding iterative scheme (precursor to Newton-Raphson): $$E_{k+1} = E_k - \frac{E_k - e \sin E_k - M}{1 - e \cos E_k}$$

2.4 Third Law: The Harmonic Scale Invariant

In Harmonices Mundi (1619), Kepler derived the empirical relation: $$\frac{T^2}{a^3} = K_{\odot} \approx \text{constant for all planets}$$

Under Newtonian gravitation with solar mass $M_\odot$ and planetary mass $m$: $$T = \frac{2\pi a b}{dA/dt} = \frac{2\pi a (a\sqrt{1-e^2})}{\frac{1}{2}\sqrt{G(M_\odot + m) a(1-e^2)}} = \frac{2\pi a^{3/2}}{\sqrt{G(M_\odot + m)}}$$ $$\implies \frac{a^3}{T^2} = \frac{G(M_\odot + m)}{4\pi^2} \approx \frac{G M_\odot}{4\pi^2} \quad (\text{since } m \ll M_\odot)$$

3.1 Geometric Structure of Keplerian Orbit

                     Semi-minor axis b
                           ^
                           |
             . - - - - - - + - - - - - - .
         '                 |                 '
       /                   |                   \
      |    Sun (Focus 1)   |       Focus 2      |
  ----+---------*----------+----------*---------+----> Major Axis
     Perihelion |<- c = ae ->|                  Aphelion
      |  r_min             |                     |
       \                   |                    /
         .                 |                 .
             ' - - - - - - + - - - - - - '
                           |
                     <- -  a  - ->

3.2 Numerical Orbit Inversion Flow (Kepler Solver)

[Target Time t]
       |
       v
[Compute Mean Anomaly: M = n(t - t_0)]
       |
       v
[Iterative Newton-Raphson Solver: E_{k+1} = E_k - (E_k - e sin E_k - M)/(1 - e cos E_k)]
       |
       v
[Compute True Anomaly: tan(theta/2) = sqrt((1+e)/(1-e)) tan(E/2)]
       |
       v
[Calculate Radial Distance: r = a(1 - e cos E)]
       |
       v
[State Coordinates: x = r cos(theta), y = r sin(theta)]

5.1 Failure under Multi-Body Perturbations

Kepler's laws assume a strictly unperturbed two-body system ($N = 2$). In the actual solar system, gravitational interactions between massive outer planets (e.g., Jupiter-Saturn 5:2 mean-motion resonance) cause secular and periodic drift in orbital elements ($a, e, i, \Omega, \omega$).

5.2 Singularity at Parabolic / Hyperbolic Escape

For parabolic ($e = 1$) and hyperbolic ($e > 1$) orbits, Kepler's Equation breaks down ($M = E - e \sin E$ diverges or becomes unphysical). Barker's equation ($e=1$) or hyperbolic anomaly $H$ ($M_h = e \sinh H - H$) must be substituted.

5.3 Relativistic Perihelion Precession

In General Relativity, Schwarzschild curvature introduces an effective $-\frac{GM \ell^2}{c^2 r^3}$ potential term, generating a perihelion shift per orbit: $$\Delta \varpi \approx \frac{6\pi GM}{c^2 a (1 - e^2)}$$ For Mercury, this produces the famous $43''$ arcseconds/century advance unexplained by Keplerian dynamics.

1. Kepler, J. (1596). Mysterium Cosmographicum (The Cosmographic Mystery). Gruppenbach, Tübingen.
2. Kepler, J. (1604). Ad Vitellionem Paralipomena, quibus astronomiae pars optica traditur (Optical Part of Astronomy). Claudius Marnius & Heirs of Jean Aubry, Frankfurt.
3. Kepler, J. (1609). Astronomia Nova AITIOLOGETOS, seu physica coelestis, tradita commentariis de motibus stellae Martis (New Astronomy). Voegelin, Heidelberg.
4. Kepler, J. (1611). Dioptrice (Dioptrics). David Franck, Augsburg.
5. Kepler, J. (1619). Harmonices Mundi Libri V (The Harmony of the World in Five Books). Gottfried Tampach, Linz.
6. Kepler, J. (1618–1621). Epitome Astronomiae Copernicanae (Epitome of Copernican Astronomy, Vols. I–III). Linz & Frankfurt.
7. Kepler, J. (1627). Tabulae Rudolphinae (The Rudolphine Tables). Jonas Saur, Ulm.