The collection

Every essay — page 3

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 41–60 of 514.

Gravitation

Two bodies pulling on each other, and everything that goes wrong at three.

Three bodies, integrated. Three equal masses integrated forward under mutual gravity: the figure-eight choreography — all three bodies on one closed curve. Every point is a step of the equations of motion, and the total energy is conserved to 2.5e-11 across the run.

Three bodies, and what "no solution" actually means

The three-body problem is routinely called unsolvable. Trajectories are computed for it every day, exact periodic solutions are known, and both statements are true — the word is doing more work than it looks.

8 figures · The three-body problem $The Kirkwood gaps. Asteroid numbers against semi-major axis, with the resonant radii marked. Each gap sits where the orbital period is a simple fraction of Jupiter's, and each of those radii is computed from the harmonic law rather than placed by eye.$

Resonance clears a gap in one place and locks a moon in another

When two orbital periods are in a simple ratio, small tugs stop averaging away and start accumulating. Sometimes that empties a region entirely. Sometimes it holds three moons together for the age of the solar system.

8 figures · Resonance $Zero-velocity curves, mass fraction 0.15. Level sets of the Jacobi constant in the frame that rotates with the two bodies. A particle with a given value of C is confined to the side of its own curve where the kinetic energy would be positive; as C falls the curves open, first at the inner point between the bodies, then behind the smaller one, and finally around the whole system.$

The curve that says where a body cannot go

The restricted three-body problem has no solution and one conserved quantity. That quantity is enough to draw a boundary the body can never cross — without integrating anything, without knowing where it started, and for all time.

8 figures · The three-body problem $The surface a star stops at. The equipotential through L₁ — the Roche lobe — at mass fractions 0.50, 0.20, 0.05, in the frame that rotates with the pair. Each is a level set of exactly the same function the zero-velocity curves are level sets of, at exactly the critical value, so nothing here is a new construction: the Roche lobe is the last closed equipotential, and it is closed only because the two lobes touch at a single point. Material that reaches that point is no longer bound to the star it came from, and it leaves through an opening of zero area. The lobes are drawn in the orbital plane; in three dimensions each is a teardrop, and its volume-equivalent radius is what "the size of a Roche lobe" means. As the mass ratio becomes extreme the smaller star's lobe shrinks towards it, which is why a white dwarf accreting from a companion has a lobe smaller than the Sun.$

The surface a star stops at

Around each star of a close pair there is a last closed equipotential, and the two touch at a single point. A star that swells to reach it hands its outer layers to its companion through an opening of zero area — and the transfer, once started, makes itself worse.

8 figures · Lagrange points

Mercury's perihelion, term by term. The observed advance of Mercury's perihelion is 5600 arcseconds per century against the equinox. Almost all of it is the equinox: the coordinate frame itself turns, and removing it leaves 574. Subtracting the perturbations of the other planets — computed by Le Verrier in 1859 and refined many times since — leaves 42.98 arcseconds a century that nothing in Newtonian gravitation accounts for. The bars are logarithmic in nothing; they are the real proportions, which is why the residual is barely visible beside the frame term.

Forty-three arcseconds, after everything else

Mercury's perihelion moves through 5,600 arcseconds a century. Nearly all of that is the coordinate system turning, and almost all of the rest is the other planets. What was left over was 43 — under one per cent of the raw number, and the most consequential residual in the history of the subject.

8 figures · Relativistic orbits

G: 14 determinations in two families. Published determinations of G, each with its quoted one-sigma interval, sorted into two families — torsion balance, in one form or another, against beam balance, pendulum, atom interferometry. The shaded band behind each family is that family's inverse-variance weighted mean: 6.67435 ± 0.00004 across 11 of them, against 6.67343 ± 0.00009 across 3. The difference is 0.00092 ± 0.00010 10⁻¹¹ m³ kg⁻¹ s⁻², which is 9.3 standard deviations, computed here from the quoted errors alone. The arithmetic is the same one the Hubble figure uses and here it should be distrusted, because the scatter inside each family already exceeds what the intervals allow: eleven torsion-balance determinations spread over 500 parts per million with quoted intervals of 12 to 130 cannot all be right, whatever the difference between the families comes to. That is why the recommended value's uncertainty is expanded far beyond any single experiment's rather than being the weighted combination drawn here — the disagreement is between laboratories using the same method, not between methods.

Nothing in the sky is weighed in kilograms

The Sun's gravitational parameter is known to eleven significant figures. The Sun's mass is known to five. The two statements are about the same object and the difference between them is a constant measured in basements, which is the worst-determined fundamental constant in physics.

8 figures · Gravitational constant

A day five hours long

The tidal bulge leads, so the Earth's spin is being paid into the Moon's orbit. Run the measured payment backwards and two curves come out of one integration — a timeline that is refuted by the Moon's own age, and a day length that is refuted by nothing.

8 figures · Tides

PSR B1913+16: 45 years of periastron arriving 88 seconds early. The cumulative shift of periastron passage for PSR B1913+16, Δt = ½(Ṗ_b/P_b)T², over 45 years from its discovery in 1975. The orbital period is 0.322997449 days and is measured to be shortening at 2.423e-12 seconds per second — a change in the twelfth decimal place, which over a career accumulates to 87.5 seconds. That is the whole reason the effect is measurable. Ṗ_b is taken here as a measured quantity and no radiated power is derived from it; the parabola is the general-relativistic prediction as published, and the points are that prediction scaled by the published ratio of observed to predicted decay, 0.997 ± 0.002 (Hulse & Taylor 1975; Weisberg & Huckins 2016), which is how the agreement is quoted. The error bars are drawn: at the last point the bar is 1.25 px tall against a dot 9 pixels across, so they are invisible, and their invisibility is the result. The campaign has run 45 years.

An orbit measured to be shrinking

A 7.75-hour orbital period that shortens by 68 nanoseconds each turn is beyond any single measurement, and unmissable after fifty thousand of them, because the shift accumulates as the square of the elapsed time. Forty-five years of pulse arrival times have made it 87.5 seconds, which is why the rate is a measured quantity rather than an inferred one.

9 figures · Relativistic orbits $The Hill radius against the mass ratio, and where four planets keep their outermost moon. The Hill radius R_H = a(μ/3)^⅓ and the exact L₁ distance from the same quintic, both in units of the planet's own semi-major axis, against the mass ratio μ = m/(M+m) from 10⁻⁸ to 0.12. The two are the same line to the width of the stroke across the whole planetary range — the cube-root formula runs 0.33% short of the exact L₁ distance at the Earth's mass ratio and 10.0% short at the right-hand edge — so the departure between them is drawn on its own percentage scale at the right, which is the only place it can be seen. Below the curves are the numerically determined stability limits, 0.5 R_H for a prograde satellite and 0.7 R_H for a retrograde one, and the outermost known satellite of each of four planets, each at a Hill fraction computed from that planet's own mass and orbit: Moon 0.257, Sinope 0.451, Phoebe 0.198, Neso 0.424. Of the four, three are retrograde — Sinope, Phoebe, Neso — and they reach 0.451 of their planet's Hill radius against 0.257 for the one prograde satellite: 1.8 times as far out, in a diagram whose two stability limits stand in the ratio 0.7 to 0.5. No equilibrium argument predicts that, and it is what the drawing makes plain.$

The region a planet may keep a moon in

A satellite is not held by orbiting but by orbiting inside the radius at which the Sun's tidal field would take it away. That radius is the same distance the innermost Lagrange point sits at, asked a different question — and the boundary the sky actually respects is smaller, by a factor that depends on which way the moon is going round.

10 figures · Hill sphere

Forty-two minutes, from anywhere to anywhere. Left, the gravitational field inside the Earth: a straight line for a uniform sphere, because the enclosed mass grows as r³ and the field as r, and the PREM curve for the real one, which is nearly flat through the whole mantle at about 9.94 m/s² because the dense core is already all below. Right, what falls through it. A body dropped down a diametric tunnel through a uniform Earth executes simple harmonic motion with ω = √(g/R), reaching the far side in 42.2 minutes; a body dropped down a chord at 0.6 of the radius feels only the component along the tunnel, which is the same ω times a smaller distance, so it arrives in the same time from a shorter trip. The period does not contain the length of the tunnel, its direction, or where the body starts. And 2π/ω is also the period of a circular orbit grazing the surface, 84.3 minutes — the tunnel and the orbit are one ellipse, seen twice. The real Earth is not uniform and gets there in 38.2 minutes instead, which is the honest number and is 9% quicker.

The tunnel that takes the same time from anywhere

Inside a uniform sphere the field grows in proportion to the distance from the centre, which is Hooke's law. A body dropped down any straight tunnel arrives in the same time — and that time is the period of an orbit skimming the surface.

8 figures · Shell theorem

Energy error over 240 revolutions, at one step size. The relative error in total energy against revolution number, for three integrators run on the same Kepler orbit at e = 0.5 with the same step of 200 per revolution. The exact energy is a constant, so every curve here is the method rather than the problem. Euler climbs steadily: its energy at the end is 106.5% wrong, and the orbit it draws has spiralled outwards. Runge–Kutta 4 begins 4.7e+2 times more accurate than leapfrog and ends at 2.22e-4, having grown by a factor of 10 across the run: the error is SECULAR. Leapfrog oscillates inside a band and stays there — worst error 2.62e-3, and the second half of the run is no worse than the first, which is measured here rather than claimed. That is the property that decides whether a five-billion-year integration means anything, and it is not accuracy: a symplectic method is the exact solution of a Hamiltonian a step-size away from the intended one, so its energy cannot wander, while a more accurate non-symplectic method has no such constraint and eventually wanders further.

Wrong about where, and right about how much

Runge–Kutta is the more accurate method and loses energy steadily; leapfrog is cruder and its energy error never leaves a band. Over five billion years only one of those properties survives — and neither method knows where the planet is.

8 figures · Numerical integration

The system that gets hotter as it loses energy

The virial theorem makes a self-gravitating body's total energy equal to minus its kinetic energy. Radiating heat away therefore raises the temperature, there is no equilibrium to settle into, and every star and every cluster is running away from one.

8 figures · Virial theorem

A separation of 10⁻¹⁰ becomes 0.0285. The distance between two copies of Burrau's problem, started 10⁻¹⁰ apart in one coordinate of one body and then stepped in lockstep — one loop, one step size, both states advanced by it — for 60 time units. The vertical axis is logarithmic, so the straight stretch is exponential growth, and its slope is the Lyapunov exponent: 0.3413 per time unit, fitted over the 604 samples lying between ten times the starting offset and a tenth of the system's own size. That is an e-folding every 2.93 time units, so the separation multiplies by ten every 6.75 and by 7.8·10⁸ across the whole run. It has not saturated within the drawn interval: the growth rate over the last tenth of the run is still 0.6002 per time unit, 176% of the fitted exponent, and the separation has reached only 0.90% of the system's own size. The exponent is what sets a prediction horizon, and this run is drawn short of it on purpose: the straight line is the measurement, and the flat part that follows is arithmetic about how far apart two bounded systems can get.

A prediction with an expiry date

The inner solar system's Lyapunov time is about five million years, so a centimetre of error becomes an orbit in a hundred million. The ephemeris dies while the system survives — because the elements stay bounded when the phase does not, and only one of those is what stability means.

8 figures · Chaos

The contours a planetary encounter cannot cross. Contours of the Tisserand parameter with respect to Jupiter in the plane of semi-major axis and eccentricity, drawn for a coplanar orbit. An encounter with Jupiter moves a comet along one of these curves and never across one, because T is what the encounter conserves. The heavy contours are at T = 3 and T = 2, and they are the boundaries the comet families are defined by: T > 3 means no encounter is possible at all, since v∞²/v_J² = 3 − T and a negative squared speed is not a trajectory; 2 < T < 3 is the Jupiter family; below 2 the approach speed exceeds Jupiter's own orbital speed and the orbits are the nearly isotropic ones. The shaded boundary is the crossing condition — an orbit whose pericentre is outside Jupiter's, or whose apocentre is inside it, never meets the planet whatever its T. The five comets are placed at their JPL elements and labelled with the T the literature quotes; all five with inclination sit off the coplanar contours by exactly the cos i in the definition, which is why 1P/Halley's is negative — a retrograde orbit meets Jupiter at nearly twice Jupiter's speed.

The number that survives the encounter

A comet that passes Jupiter comes away with every orbital element changed. One combination of them is not changed, and it is enough to recognise the comet afterwards, to sort the comet families, and to say where a spacecraft can and cannot go.

9 figures · Tisserand parameter

An average that arrives, and a snapshot that never does. 2T/|U| against time for Burrau's problem — three masses released from rest, and never repeating, in two forms: the instantaneous ratio, and the ratio of the running time-averages. The instantaneous one runs between 0.03 and 1.97 — the system is a long way from equilibrium at almost every moment, because the bodies are alternately falling together and flying apart. The averaged one settles: after 50 time units it is 1.0142, against the exact 1 the virial theorem requires of any bound system. This system is not periodic and is not even permanently bound — Burrau's problem ejects its lightest body — so the average is drifting rather than converged, and that is the theorem's own condition made visible: it is a statement about bound systems and says nothing about anything that is leaving. What the figure cannot show is the error: a cluster observed once gives 2T/|U| with a scatter of this size, and what makes its mass believable is not the measurement but the relaxation.

An average that weighs what cannot be watched

A cluster's mass can be had from its speeds alone — no orbit followed, no period observed, no distance to any single star. The theorem that allows it is an average over time, which is exactly what a photograph is not.

8 figures · Virial theorem $A thousandth of the field, and all of the precession. Left: the Earth's figure against a sphere of the same equatorial radius, with the flattening drawn 28× its true value. The real difference between the equatorial and polar radii is 21.4 km on 6378 — 1 part in 298 — which at this size would be 0.5 pixels and invisible, so the drawing is a schematic and the number is here instead. Right: the two components of the J₂ perturbation at the surface, each as a fraction of the monopole μ/r², both differentiated from the potential rather than quoted. The radial one strengthens the inward pull by 1.62×10⁻³ over the equator, where the extra mass is, and weakens it by 3.25×10⁻³ over the poles, vanishing at ±35.26° where P₂ does. The transverse one is zero at the equator and at both poles and peaks at ±45°, at 1.62×10⁻³ — and that is the component that does the work. It pulls an inclined orbit back towards the equatorial plane, which is a torque about the line of nodes, and a torque applied to something already turning moves it sideways rather than back. Averaged over an orbit the pair leave a, e and i untouched and turn the whole plane instead, which is why a term a thousandth of the field is the largest single perturbation on almost every satellite ever flown.$

The Earth's shape, read off a satellite's node

The Earth is a thousandth of a part from being a sphere, and that thousandth turns every satellite's orbital plane. Vanguard 1 measured it in 1959 — and one retrograde inclination turns the plane at exactly the rate the Sun moves, which is a perturbation used as a design constraint rather than corrected for.

8 figures · Oblateness

The wake, computed from the streamlines that make it. Left: 26 streamlines past a point mass, in the mass's own frame, each integrated from far upstream with the same speed and a different impact parameter, and mirrored about the axis. Nothing is drawn to converge — every track is the hyperbola its own impact parameter gives it, and they cross downstream because an attraction focuses. Right: the density that focusing produces at 3.2 focusing radii behind the mass, as (b/y)(db/dy) — the Jacobian of the map from starting radius to arrival radius — which peaks at 15.53 times the background at 0.03 radii off the axis. The overdensity is behind the mass, and that is the entire mechanism: the wake pulls backwards on the body that made it. What the figure cannot show is the steady state, because it has no time in it: a real wake is continuously replenished, and the drag is the sum over an infinite train of these encounters, which is where Chandrasekhar's logarithm comes from.

A drag with nothing to drag against

A massive body moving through a sea of light ones raises an overdensity behind itself and is pulled back by it. The force does not depend on the masses of the background bodies at all — and it is strongest at one particular speed.

9 figures · Dynamical friction

The last curve across the cylinder, before and after it breaks. Surfaces of section for the standard map at K = 0.5, 0.9716, 1.4, each 22 trajectories iterated 41 times from a column of starting values. Nothing here is placed: every dot is an iterate. At K = 0.5 the islands are separate and the space between them is filled with curves that run all the way round in θ — a trajectory cannot get from one island to the next, and one started beside the hyperbolic fixed point wanders 2.91 in p — 0.46 of a cylinder — and stops. At K = 0.9716, which is Greene's threshold to four figures, the last of those curves is on the point of going and the same trajectory still only reaches 4.84, which is 0.77 of a cylinder. At K = 1.4 it covers 6.1 cylinders: the barrier is gone and there is nothing left to stop it. That transition is what a chaotic zone in the asteroid belt is, drawn without any asteroids.

Where the chaos comes from

A Lyapunov time says how long a prediction lasts. It does not say what destroyed it. The mechanism is two resonances whose libration widths overlap, and the transition can be watched happening on a surface of section as one number is turned up.

8 figures · Chaos

255 km sees degree 154; 35786 km sees degree 6. The same field spectrum, multiplied by the upward continuation factor (R/r)ˡ⁺¹ for orbits at 255 km, 450 km, 800 km, 35786 km. A harmonic of degree ℓ has ℓ bumps around the planet, so it falls off with height like a wave of that wavelength — fast, and faster the finer it is. The horizontal rule is a measurement floor; where each curve crosses it is the highest degree that orbit can feel at all, and the answer is 154 at 255 km, 101 at 450 km, 65 at 800 km, 6 at 35786 km, or 130 km, 198 km, 308 km, 3340 km of horizontal resolution on the ground. Nothing in that arithmetic is an instrument. It is why GOCE flew at 255 kilometres with drag compensation rather than at a comfortable altitude with a better gradiometer, why the Moon's mascons were not seen until something orbited low over them, and why a geostationary satellite's ephemeris needs a field with four terms in it.

The field a satellite is allowed to feel

Past the flattening, a planet's gravity is a sum of harmonics whose sizes follow a rule with no physics in it. How much of that sum a spacecraft can measure is decided by its altitude and by nothing else — which is why one mission flew at 255 kilometres and had to push itself along.

8 figures · Oblateness

A feeding zone, and the spacing it forces

The radius that decides what a planet may keep also decides what it could reach while it was growing — and measuring the gaps between planets in that unit turns a distribution spanning three orders of magnitude in astronomical units into a band a factor of ten wide, with a floor that is partly a theorem and partly a fit.

9 figures · Hill sphere

Every sequence · Every named object · Search