Astronomy & Celestial Mechanics
the sky, computed
The collection

Every essay

One idea per essay, ordered so that the earlier ones set up the later ones — but nothing here depends on being read in sequence. Essays 1–20 of 514.

Orbits

Kepler's three laws, and the family of curves a single force allows.

An orbit at eccentricity 0.6. An orbit of eccentricity 0.6. The primary sits at a focus, offset from the centre by 0.6 of the semi-major axis, and the closest and furthest points differ by a factor of 4.00.

The orbit is an ellipse, and the Sun is not in the middle of it

Kepler's first law is usually drawn wrong. The interesting content is not the ellipse — it is the focus, and the fact that one of the two is empty.

9 figures · The ellipse Equal areas in equal times, at eccentricity 0.65. Positions of an orbiting body at equal intervals of time, obtained by solving Kepler's equation. The two shaded sectors span the same interval and enclose the same area — a long thin one at the far end, a short fat one at the close approach.

Equal areas in equal times, which is angular momentum in disguise

Kepler's second law is a statement about the area a radius line sweeps. It looks like an odd thing to have noticed, and it turns out to be a conservation law arriving eighty years early.

10 figures · Angular momentum What an eccentricity does to the shape, and what it does to the offset. The fractional flattening 1 − b/a and the fractional focal offset c/a, against eccentricity. The offset is first order in e and the flattening is second, so at Earth's e = 0.0167 the outline is 0.014% from a circle while the Sun sits 1.67% of the semi-major axis off the centre — a factor of 120.

An orbit can look exactly like a circle and still not be one

Earth's orbit departs from a circle by fourteen parts in a hundred thousand. The Sun's offset from its centre is a hundred and twenty times larger, and everything interesting is in the offset.

8 figures · The ellipse A vector that does not change, drawn five times. The Laplace–Runge–Lenz vector A = v × L − GM r̂, constructed at five points of one orbit at e = 0.45 under a force ∝ r^−2, each drawn from the body rather than from the focus so that its length and direction can be compared point by point. Under the inverse square every one of them is the same vector: identical to machine precision in length and in direction. Its length is 0.4500, which is the eccentricity, and it points at pericentre — so the orbit's orientation is a conserved quantity and not a constant of integration, which is why the ellipse does not turn.

The third thing that is conserved

Energy fixes an orbit's size and angular momentum fixes its shape. Neither fixes which way it points — so a curve that closes needs a third conserved quantity, and only two force laws in the universe supply one.

8 figures · Conic sections The three anomalies at E = 1.15 rad. The auxiliary circle construction. The eccentric anomaly E is measured at the centre, the true anomaly ν at the focus, and the mean anomaly M is time expressed as an angle. Here E = 1.150, ν = 1.827 and M = 0.602 radians, related by M = E − e sin E.

The position that has no formula, and is computed anyway

Kepler's equation relates where a body is to when it is there. It cannot be solved in elementary functions, Kepler said so, and nobody has managed it since — which has stopped nothing.

8 figures · The ellipse Period against size for the planets, around the Sun. Orbital period against semi-major axis on logarithmic axes. The line has slope exactly three-halves — the harmonic law — and the measured bodies sit on it.

The law that links period to size, and weighs everything

Kepler found that the square of the period goes as the cube of the orbit. Newton found the constant of proportionality, and that constant is a mass — which is how every mass in astronomy has been obtained since.

9 figures · Harmonic law Every orbit one force allows. Circle, ellipse, parabola and hyperbola, all sharing a focus and a closest approach. The eccentricity alone decides which one a body is on, and whether it returns.

Every orbit one force allows, and the number that picks between them

Circle, ellipse, parabola, hyperbola. A single inverse-square force permits exactly these four, and one number decides which — including whether the body ever comes back.

10 figures · Conic sections An orbit at i = 42°, Ω = 35°, ω = 55°. The three orientation elements. The orbit is generated in its own plane and rotated by the standard sequence, so the inclination, the node line and the argument of periapsis are the angles that produced the drawn curve rather than labels applied to it.

Six numbers that fix an orbit for all time, and the sixth is the awkward one

Five of the orbital elements describe a curve that never changes. The sixth says where on it the body is, and it is the only one that has to keep being measured.

10 figures · Orbital elements The energy budget of an orbit at e = 0.7. Kinetic, potential and total energy per unit mass against distance from the primary, in units where GM = 1. The total is a horizontal line — it depends only on the semi-major axis — and the vis-viva relation is that statement solved for the speed.

One equation for the speed anywhere, and the eccentricity is not in it

The vis-viva relation gives the speed at any point of any orbit from two numbers. What it leaves out is the surprise — the shape of the orbit does not appear at all.

10 figures · Vis-viva One series, two eccentricities, and a limit between them. The Lagrange series for E − M, summed to 1, 2, 3, 6, 22 terms, against the exact solution of Kepler's equation, over half a revolution. Left, at e = 0.5: the partial sums close on the exact curve and the last two are indistinguishable from it. Right, at e = 0.8: they do not, and the 22-term sum is worse than the 1-term one, missing the exact value by 0.24 radians at M = π/2. Nothing about the orbit changes between the two panels — an eccentricity of 0.8 is an ordinary comet — and nothing about the equation changes either. What changes is that a singularity of E as a function of complex e has come inside the circle of radius e, at the Laplace limit 0.6627434, which is the root of e·exp√(1+e²) = 1 + √(1+e²) and has no astronomical meaning whatever. The coefficients are computed from the Bessel expansion in logarithms; the first two are sin M and ½sin 2M exactly, which is what the generator asserts before drawing.

The formula that exists, and is not used

Kepler's equation does have a closed-form solution — an infinite series in the eccentricity, written down by Lagrange. It converges up to e = 0.6627434 and not one part beyond, and that number has nothing to do with astronomy.

8 figures · The ellipse Three rotations, applied in order. An orbit of eccentricity 0.45 carried into space by the three orientation angles, one panel per rotation: the argument of periapsis ω = 40°, then the inclination i = 42°, then the longitude of the ascending node Ω = 55°. The last panel applies the same three angles in the reverse order and arrives somewhere else, because rotations about different axes do not commute.

Three rotations that put an orbit in space, and they do not commute

An orbit's orientation takes three angles. Give the same three angles in a different order and the orbit ends up somewhere else — which is why the convention is part of the data.

8 figures · Orbital elements The effective potential, for three angular momenta. The radial motion of an orbiting body is one-dimensional motion in an effective potential: the attraction −GM/r plus the centrifugal term L²/2r² that the angular momentum contributes. The barrier at small radius is what stops a body with any angular momentum at all from reaching the centre, and the bottom of each well is the circular orbit.

The wall that angular momentum builds

A body falling toward a star almost never arrives. Sideways motion, which looks like a detail of the initial conditions, turns the attraction into a well with a wall around the middle of it — and the wall is why it costs more fuel to hit the Sun than to leave the solar system.

8 figures · Effective potential The osculating semi-major axis of a perturbed orbit. The semi-major axis a test particle would have if the perturber vanished, computed from its position and velocity at every step of an integration over 26 orbits of the perturber. It is not constant: a short-period ripple rides on a slow trend, and only the trend accumulates.

Elements that do not stay constant

Six numbers fix an orbit for all time, and the phrase is only true in a universe containing two bodies. Add a third and the six start moving — some of them wandering and returning, one or two of them drifting in one direction forever, and the difference between those two behaviours is the whole of celestial mechanics after Newton.

8 figures · Perturbations The mass correction against the mass. Each planet's departure from the massless harmonic law, against its own mass in solar units, on logarithmic axes. The exact law puts every point on the diagonal. Jupiter and Saturn are the only planets whose mass correction is larger than the perturbations from everything else, and Saturn's measured departure has the opposite sign.

The third law is wrong by the mass of the planet

Kepler's harmonic law says the square of the period goes as the cube of the size. Newton's version has one more term in it, and the term is the orbiting body's own mass — negligible for a planet, decisive for a binary star, and the reason a period can be converted into a mass at all.

8 figures · Harmonic law 1I/ʻOumuamua: an orbit at e = 1.201, and the angle it turned through. The open branch of a conic at eccentricity 1.201 and periapsis 0.2559 AU, with the Sun at the occupied focus. Two numbers fix everything else: |a| = q/(e−1) = 1.273 AU, an impact parameter b = |a|√(e²−1) = 0.847 AU, an asymptote at ν∞ = arccos(−1/e) = 146.37° from periapsis, and a deflection of δ = 2ν∞ − 180° = 112.74° — which is the same statement as sin(δ/2) = 1/e. The asymptotes cross at |a|e = 1.529 AU from the Sun on the periapsis side, where a bound orbit's centre would be on the other. The speed left over at infinity is √(μ/|a|) = 26.40 km/s. Schematic in one respect: the drawing runs at about 119 px to the AU, so the Sun's own disc would be far smaller than its marker.

The orbit that has no period

Above an eccentricity of 1 the conic is open, the energy is positive and the semi-major axis is negative — and the vis-viva relation survives the sign change without a single alteration. What replaces the period is a speed, and that speed is what says where a visitor came from.

10 figures · Hyperbolic orbits The aim point and the miss distance are the same number far out and nothing like it close in. Periapsis distance against aim point for a hyperbolic approach to Jupiter at v∞ = 5.6 km/s, both in planet radii. The diagonal is where the two would be equal — where gravity did nothing — and the curve falls below it everywhere, by more the closer in the aim is. A trajectory aimed at 10.68 radii grazes the surface, because gravitational focusing means the planet's effective size is √(1 + v_esc²/v∞²) times its radius. The slope of this curve is what a navigation team cares about: it is 0.208 at an aim of 12 radii and 0.936 at 150, so the same correction manoeuvre changes the periapsis distance by 4.5 times as much at one end of the range as at the other, while changing B by exactly the same amount at both; far outside this plot, where focusing has run out, it reaches 1.000 and the two numbers become the same one. That is why the aim point is the coordinate a manoeuvre is quoted in, why an error ellipse is published in the B-plane, and why the turn angle — 2 arctan(μ/Bv∞²), 156.0° at 12 radii and 77.8° at 70 — is thought of as a function of B rather than of anything the spacecraft does.

Aiming at a plane instead of at a planet

A spacecraft arriving at a planet is not aimed at a periapsis distance. It is aimed at a point in a plane perpendicular to its own incoming asymptote, because that is the one coordinate in which the miss distance responds linearly to a correction — and every navigation product ever published for a flyby is written in it.

8 figures · Hyperbolic orbits Ceres from five directions and no distance. Ceres seen five times over 41 days, from an Earth on a circular orbit, reduced in the plane. Each sighting gives a direction and no range, so the object is somewhere on its sight line; the five lines here span 1.37° of geocentric arc altogether, and Earth's own motion supplies the only baseline there is — 0.403 AU of its 0.691 AU of travel lies across the sight lines. A planar orbit is four numbers, so five angles over-determine it and one orbit comes out: a = 2.7658 AU, e = 0.0785. That is the answer and not an input — the sightings were generated at a = 2.7658 AU and e = 0.0785, and the solve, which sees only the directions and the dates, returns them to 3e-12. The two shaded sectors are what closes the determination: between the first and middle sightings the radius vector sweeps 0.2993 AU² in 21.0 days and between the middle and last 0.2853 AU² in 20.0 days, a ratio of 1.04918 against a time ratio of 1.04918. Slide all three crossings out along their sight lines together and that equality fails at once, so it fixes the distance by itself, with no propagation anywhere in the argument — and it gives a = 2.7658 AU over again. The two dashed curves are candidates that thread the same three sight lines at 85% and 108% of the recovered distance: a = 1.86 AU at e = 0.35, sweeping its areas in a ratio 4.5% wrong; and a = 5.61 AU at e = 0.49, sweeping its areas in a ratio 2.9% wrong. A few per cent in the distance is an orbit of another kind, which is the same fact the conditioning panel measures: one arcsecond of angle error moves a by 0.60% on this arc.

Five directions and no distance among them

An image of a moving point records an angle and throws the range away, so an orbit has to be assembled out of angles alone. How many angles are needed is not a detail of the method — it is the whole of what a determination is.

10 figures · Orbit determination Four conics through one periapsis, drawn by one solve. Distance from the Sun against time for four orbits sharing a periapsis of 0.5 AU, at e = 0.6, e = 1, e = 1.0001, e = 1.4, over 900 days. Every point on every curve came from the same universal Kepler solve — no branch on the conic class anywhere in it — and each curve was then checked against the classical solution of its own kind at the midpoint: Kepler's equation at e = 0.6 agrees to machine precision; Barker's cubic at e = 1 agrees to machine precision; e sinh H − H at e = 1.0001 agrees to machine precision; e sinh H − H at e = 1.4 agrees to machine precision. The curve to read twice is e = 1.0001: over this arc it lies within 0.02% of the parabola and is indistinguishable from it, and it is the only one of the four whose fate the drawing cannot show. What the classical parameterisation costs there is arithmetic rather than impossibility: at the midpoint of this arc, e sinh H − H throws away 3.3 of its sixteen digits to cancellation, against 0.2 at e = 1.4 — enough to matter to an ephemeris and not enough to stop a plot.

The one solve that does not ask which conic it is

Kepler's equation is for ellipses, Barker's cubic for parabolas, and a hyperbolic sine for the rest — three parameterisations of one motion, each worst exactly where its neighbour takes over. The universal variable removes the question, and the removal is not a convenience.

10 figures · Universal variables Four averages of one distance, and the two of them that are the semi-major axis. The average distance of a body from its primary, against eccentricity and in units of the semi-major axis, computed four ways: averaged over time, over true anomaly, over eccentric anomaly, and as the harmonic mean in time. Every curve is a quadrature over the orbit — 2,048 panels uniform in eccentric anomaly, with Kepler's equation supplying the time weight — and not a closed form. Two of the four are exactly a at every eccentricity, which is why they are drawn as one line: the eccentric-anomaly average, because the mean of cos E over a turn is zero, and the harmonic mean in time, because the time weight cancels 1/r at every node before the sum begins. The other two are not: the time average is a(1 + e²/2), which rises to 1.4050 a at e = 0.9, and the true-anomaly average is a√(1−e²) — the semi-minor axis — which falls to 0.4359 a there. So a is the average distance in two senses out of four, and the ordering b ≤ a ≤ ⟨r⟩ₜ holds at every eccentricity with equality only on the circle. At Earth's e = 0.0167 the four agree to 0.014%, and at Mercury's e = 0.2056 the spread is 2.14%. The distinction is invisible for the planets and unavoidable for a comet, and it is the reason a quoted "mean distance" has to say which mean.

The average depends on what is being averaged

Four ways of averaging one orbit's distance from its primary give four different numbers, and only two of them are the semi-major axis. Which two is not a matter of convention, and the same arithmetic decides how much sunlight a planet receives in a year.

8 figures · Orbital averages One squaring, and the singularity is gone. Left, a harmonic oscillator: an ellipse centred on the origin, marked at 24 equal steps of its own phase. Right, the same points after squaring as complex numbers, u ↦ u². The image is an ellipse with the origin at a focus — checked here by the focal property, r₁ + r₂ = 2a to nine figures at four points — with semi-axes (A²+B²)/2 = 2.5000 and AB = 1.5000 and the focus at (A²−B²)/2 = 2.0000, which is ae exactly. Three things follow at once. The angle doubles, so one turn of the oscillator is one whole orbit traversed twice as fast in phase; the equal phase steps on the left arrive as the eccentric anomaly on the right, which is why that anomaly and not the true one is what the equations want; and the collision at r = 0, where the inverse square is infinite, is the point u = 0, where the oscillator has a perfectly ordinary velocity. The singularity was a property of the coordinates.

The singularity that is a change of variable

The Kepler problem blows up at zero separation, and a fixed-step integrator falls apart long before it gets there. Divide time by the radius and the equations become a harmonic oscillator — exactly, for every conic at once.

8 figures · Universal variables

Every sequence · Every named object · Search