orbital-mechanics

# Orbital Mechanics

Orbital mechanics is the branch of celestial mechanics that answers: where will this object be at time t, given its position and velocity now? The foundational discoveries were Kepler's three laws (1609, 1619) and Newton's derivation of them from universal gravitation (1687). Modern practice adds perturbation theory for multi-body problems and relativistic corrections for high-precision ephemerides. This skill covers the classical core: the six orbital elements, Kepler's laws, the vis-viva equation, orbit types, transfer orbits, gravity assists, and a handful of strategies for computing positions in practice.

**Agent affinity:** hubble (catalog cross-reference), payne-gaposchkin (binary star orbits and dynamical masses)

**Concept IDs:** astro-keplers-laws, astro-earth-moon-sun-geometry, astro-planetary-motion

## The Two-Body Problem

Two point masses interacting only via gravity produce a relative motion that is exactly solvable. The relative orbit is a conic section — circle, ellipse, parabola, or hyperbola — with the total mass concentrated at one focus.

**Assumptions:**

- Only two bodies (all others negligible)
- Point masses (or spherically symmetric)
- No radiation pressure, atmosphere, or non-gravitational forces
- Newtonian gravity (relativistic correction needed for Mercury's perihelion, binary pulsars)

When these assumptions hold, the motion is fully determined by six constants — the **orbital elements** — and an epoch.

## The Six Orbital Elements

An orbit in three-dimensional space has six degrees of freedom. There are many choices of six numbers; the classical set is:

| Element | Symbol | Meaning |
|---|---|---|
| Semi-major axis | a | Size of the orbit |
| Eccentricity | e | Shape — how elongated (0 = circle, <1 = ellipse, 1 = parabola, >1 = hyperbola) |
| Inclination | i | Tilt of the orbit plane relative to a reference plane |
| Longitude of ascending node | Omega | Where the orbit crosses the reference plane going up |
| Argument of periapsis | omega | Angle from ascending node to closest approach |
| True anomaly | nu | Current angular position from periapsis |

For solar-system work the reference plane is usually the ecliptic. For Earth satellites it is the equator. For exoplanets around a host star, the plane is defined by observation geometry.

**An alternative sixth element.** The true anomaly is time-dependent. For cataloging a fixed orbit you can substitute the **time of perihelion passage** (T_0) or the **mean anomaly at epoch** (M_0), and derive nu from them at any later time using Kepler's equation.

## Kepler's Three Laws

### First law (1609)

Planets move in elliptical orbits with the Sun at one focus.

**Consequence:** Circular orbits are a special case (e = 0). Most orbits are slightly elliptical — Earth's eccentricity is 0.017. Mercury is 0.206. Pluto is 0.249. Halley's comet is 0.967.

### Second law (1609) — Equal areas in equal times

The line from the Sun to a planet sweeps out equal areas in equal time intervals.

**Consequence:** A planet moves fastest at perihelion (closest to the Sun) and slowest at aphelion. This is equivalent to conservation of angular momentum.

**Formula:** dA/dt = L / (2m) where L is the orbital angular momentum. For Earth the ratio of perihelion to aphelion speeds is about 1.034, matching the eccentricity.

### Third law (1619) — Period and size

The square of the orbital period is proportional to the cube of the semi-major axis:

    T^2 = (4 pi^2 / G * M_total) * a^3

For planets around the Sun, using years and AU:

    T^2 = a^3

where T is in years and a is in AU. Jupiter has a = 5.2 AU, so T = 5.2^1.5 ~ 11.86 years. Checks.

**Binary stars.** The same law applies with M_total = M_1 + M_2. Measuring a and T gives the total mass. Measuring the two stars' individual orbits around the barycenter gives the mass ratio, separating the two masses.

## The Vis-Viva Equation

The workhorse of orbital mechanics:

    v^2 = G * M * (2/r - 1/a)

where v is orbital speed, r is current distance from the focus, and a is the semi-major axis. **Vis-viva** is Latin for "living force" — an old name for kinetic energy.

**Corollaries:**

- **Circular orbit** (r = a): v_circ = sqrt(G M / r). For Earth at 1 AU, this gives 29.78 km/s.
- **Escape velocity** (a = infinity): v_esc = sqrt(2 G M / r) = v_circ * sqrt(2). For Earth surface, 11.2 km/s.
- **Parabolic trajectory** (e = 1): v_esc marginally exceeded, approaches zero at infinity.
- **Hyperbolic trajectory** (a negative in the vis-viva convention): v at infinity is nonzero — you leave and keep going.

**Use.** Given a spacecraft's position and desired speed, vis-viva tells you what semi-major axis you are on, which fixes period and future positions.

## Orbit Types

### Circular (e = 0)

Constant speed, constant distance. Low Earth orbit, geostationary orbit, most spacecraft parking orbits. Not strictly realized in nature but a useful idealization.

### Elliptical (0 < e < 1)

Closed orbit. All planets, most asteroids, most comets, most moons, most binary stars. The two special points are **perihelion** (closest to focus) and **aphelion** (farthest).

For orbits around bodies other than the Sun, the terminology changes: perigee/apogee (Earth), perijove/apojove (Jupiter), periastron/apastron (star), perihelion/aphelion (Sun). The root is the body name.

### Parabolic (e = 1)

Marginally unbound. Never repeats. Speed at infinity equals zero. A theoretical boundary case — real "parabolic" comets are typically slightly elliptical or slightly hyperbolic.

### Hyperbolic (e > 1)

Unbound. The object visits once and leaves forever. Interstellar objects like 1I/Oumuamua (2017) and 2I/Borisov (2019) follow hyperbolic orbits through the Solar System.

## Kepler's Equation

To get position as a function of time in an elliptical orbit, you need to solve Kepler's transcendental equation:

    M = E - e * sin(E)

where M is the **mean anomaly** (linear in time, M = 2 pi * t /