Astronomy & Celestial Mechanics
the sky, computed
Depth

Series — page 8

A field says what an essay is about. A series follows one idea essay by essay — from the question that introduces it to the one that assumes all the others.
Refraction against apparent altitude. How far the air lifts an object, in arcminutes, against where the object appears to be. The lift is 34.5′ at the horizon — larger than the Sun's own diameter — and falls below one arcminute above 45°. At the horizon it is also the least reliable number in positional astronomy, because it depends on the temperature profile of the air the sight line crosses.

Refraction

  1. 1 The Sun sets before it sets
  2. 2 The atmosphere is a prism as well as a lens
2 essays · sky 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.

Relativistic orbits

  1. 1 Forty-three arcseconds, after everything else
  2. 2 An orbit measured to be shrinking
2 essays · gravitation Catching a target 40° ahead. A phasing manoeuvre. Dropping into an orbit 6% lower shortens the period to 0.9553 of the target's, so the chaser gains 16.1° each lap and closes 40° in 3 revolutions. Speeding up would have lost ground instead.

Rendezvous

  1. 1 Catching up by slowing down, which cost Gemini 4 its fuel
  2. 2 A burn that moves the wrong way
2 essays · spaceflight The rotation curve of NGC 3198, decomposed. Circular speed against radius for a three-component model of NGC 3198: a Hernquist bulge of 1.0×10⁹ M☉, an exponential disc of 2.20×10¹⁰ M☉ with a scale length of 2.6 kpc, and a pseudo-isothermal halo whose asymptotic speed is not chosen but solved for — it is whatever brings the total to the measured 150 km/s at 30 kpc, and comes out at 171 km/s. The components add in quadrature because accelerations add. The disc alone peaks at 5.7 kpc and falls away; the total does not.

Rotation curves

  1. 1 A rotation curve that refuses to fall
  2. 2 The same curve, two galaxies
2 essays · galaxies Earth's eccentricity is a sum of 8 sinusoids. The eccentricity of Earth over 800 thousand years, from the Laplace–Lagrange solution for all eight planets — the secular matrix built from the JPL masses and semi-major axes, symmetrised, and diagonalised by Jacobi rotations. It runs between 0.0035 and 0.0436, and it has no period, because it is a sum of 8 incommensurable frequencies. The two largest contributions to this planet are the modes at 3.73 and 7.33 arcseconds per year, drawn as the flat lines: those are constants, and everything moving in the figure is their beat. The check is the quadratic form ½ΣΛe², which a symmetric secular matrix conserves exactly and which drifts by 6.7e-16 across the whole interval — computed from the same curves the figure draws and from nothing the eigenvalues were fitted to. The true angular momentum deficit, Σ Λ(1 − √(1−e²)), drifts by 7.1e-4, and that difference is not an error either: the two agree only to fourth order in e, and Mercury at 0.206 supplies almost all of the gap.

Secular theory

  1. 1 No planet has an eccentricity of its own
  2. 2 The bound that holds only in the linear theory
2 essays · orbits Why a shell pulls a point inside it not at all. A double cone from a point inside a uniform shell. In the narrow-cone limit the far patch is 2.80 times further away and 2.80 times wider — the same number, checked to a part in a thousand before this figure is drawn. Its mass is greater by the square of that ratio and its pull weaker by the same square, so the two cancel in every direction.

Shell theorem

  1. 1 A sphere pulls exactly like a point, and the proof is a pair of cones
  2. 2 The tunnel that takes the same time from anywhere
2 essays · gravitation 4 cycles of wings, and the polarity reverses at every boundary. Sunspot latitude against date, one mark per spot group, over 4 cycles of 11 years. The pattern is the reason the plot is called a butterfly diagram, and both of its features are laws with names. Spörer's law is the downward slope: spots emerge near ±28° at the start of a cycle and near ±7° at the end, so each wing narrows towards the equator and never crosses it. Hale's law is what the two mark shapes say: the leading spot of a pair has one magnetic polarity in the north and the other in the south, and both reverse when a new cycle starts — so a diagram that repeats every 11 years in appearance repeats only every 22 in magnetism. The wings overlap: the first high-latitude spots of a cycle appear about 1.6 years before the last low-latitude spots of the one before, which is why counting spots gives a cycle length slightly different from measuring one between polarity reversals. The dot density in time is the sunspot number itself, drawn from the standard skewed fitting function — the rise to maximum takes about four years and the decline about seven, in every cycle ever recorded.

Solar cycle

  1. 1 A magnetic clock read off a butterfly
  2. 2 One step of memory kept at the poles
2 essays · stars Seven sources over twelve decades of flux, and one threshold below the one that matters. The Sun's neutrino spectrum at the Earth, with the continua drawn per unit energy and the two monoenergetic lines as spikes at their own energies. The pp reaction supplies 91 per cent of all of them and its endpoint is at 0.4233 MeV. The vertical lines are experimental thresholds, and they are the figure's argument: only gallium sits below that endpoint. Chlorine, which produced the deficit and held it for twenty years, could not see a single pp neutrino — it counted ⁷Be and ⁸B, which are a rare branch of a rare branch, together under a per cent of the total — and the water detectors that followed were higher still. So the discrepancy that eventually turned out to be a property of the neutrino was measured, for two decades, using the least representative one per cent of the flux available. The total drawn here is 6.54·10¹⁰ cm⁻² s⁻¹, and it is checkable without any stellar model at all: every completed chain turns four protons into helium, releases 26.73 MeV of which 0.59 leaves as neutrinos, and emits two neutrinos — so the solar constant of 1361 W m⁻² fixes the number at 6.5·10¹⁰ cm⁻² s⁻¹, within 0.6 per cent of the sum of the model's own branches. The dominant flux is a consequence of the Sun shining and of nothing else.

Solar neutrinos

  1. 1 The only thing that leaves the centre
  2. 2 A flux that is a thermometer to a tenth of a per cent
2 essays · stars A stream that is not the orbit it came from. 456 stars released in pairs from the two saddles of a 10⁵ solar-mass cluster over 4.0 billion years, integrated in a halo whose circular speed is 220 kilometres a second, drawn with the progenitor's own orbit. The cluster runs between 10 and 25 kiloparsecs and the orbit is the thin closed-looking curve; the stars are everything else. The point of the figure is the discrepancy. Stars leaving through the inner saddle are on slightly smaller orbits, turning round at a median of 24.77 kiloparsecs rather than the progenitor's 25.00, and therefore running ahead; stars leaving through the outer saddle reach 25.22 and fall behind. The whole spread is 1.8 per cent of the apocentre, which is the number worth carrying: an offset far too small to see in this drawing builds the entire stream, because it acts for four billion years. The two arms are therefore not merely displaced along the orbit, they are on different orbits, and the track a survey measures is a family of them rather than any single one. Fitting a Galactic potential by demanding that a stream lie along an orbit is wrong by exactly this much, and the size of the error grows with the mass of the progenitor, because the mass is what sets the distance between the two doors.

Stellar streams

  1. 1 A stream is not the orbit it came from
  2. 2 A hole that says mass times time
2 essays · galaxies A wind that writes its own terminal speed on the blue edge, at 2017 km/s. A P Cygni profile, computed by integrating a spherical wind on a grid rather than drawn. The wind accelerates outward as one minus the inverse radius to the power 0.8, reaching 2000 kilometres a second, and its density falls as the inverse square of the radius over its own speed; the horizontal axis is the Doppler shift in units of that terminal speed, with blue to the left. The column directly in front of the stellar disc is moving toward the observer at every speed from nearly nothing near the surface up to the terminal value far out, so it takes light out of the beam across that whole range and the absorption trough runs from the rest wavelength to the blue edge at -1.01 of the terminal speed. That edge is the measurement: no model of the star, its distance or its mass enters it, only the geometry of a column seen end-on. The emission comes from everywhere else, where the wind moves across the line of sight and can only add photons; it is Infinity times stronger to the red than to the blue, because the part of the shell receding directly away is hidden behind the star and its blueshifted counterpart is not. Two features, one shell, and between them a velocity and a mass-loss rate.

Stellar winds

  1. 1 A speed read off an edge
  2. 2 The flare that arrives from somewhere else
2 essays · starlight Three slopes: 2.5, −0.75 and −1.25, and only the middle one is about the electrons. A radio source's spectrum across five decades of frequency, computed for an electron population with index p = 2.5 in a uniform field. Nothing in this curve is a temperature, because a power law has no scale and therefore nothing a thermometer could read. The straight section between the two bends has slope −0.75, measured here off the drawn curve, and the electron index follows from it and from nothing else: α = (p − 1)/2, so a flux ratio between two frequencies is a measurement of the energy distribution of particles in a place no detector will ever visit. The two bends are the other two measurements. Below 80 MHz the source is opaque to its own radiation and rises as ν^2.5 — a slope fixed at 5/2 by the geometry alone, whatever the electrons are doing — and the frequency at which that happens gives an angular size for a source nothing has resolved, because the turnover is where the brightness temperature meets the electrons' own. Above 12 GHz the spectrum steepens by 0.50: the electrons that radiate at high frequency lose their energy fastest, so the top of the distribution has already emptied, and the frequency of the break is a clock. What no part of this curve gives is the magnetic field. Only the field and the particle density together enter the emission, and every field strength ever quoted for a radio source comes from assuming the two share the energy equally.

Synchrotron radiation

  1. 1 A spectrum with no temperature in it
  2. 2 A minimum that was mistaken for a principle
2 essays · starlight 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.

The three-body problem

  1. 1 Three bodies, and what "no solution" actually means
  2. 2 The curve that says where a body cannot go
2 essays · gravitation A clock that has lost 5.7 hours in 2720 years. ΔT = TT − UT1, the accumulated difference between a uniform time scale and the Earth's own rotation, from 700 BC to 2020, on a logarithmic scale. The points are the published record; the smooth curve is the parabola 31.94·u² fitted to the entries at or before 1000, with u in centuries from 1820. It is a parabola and not a line because the Earth is not merely slow, it is slowing: a day lengthening at a constant rate makes a clock fall behind by the integral of the lag. That coefficient is a length-of-day rate and nothing else — 31.94 = ½ × r × 36525 days per century gives r = 1.75 ms per century, against the 1.78 ms computed from tidal angular momentum with no eclipse anywhere in the derivation. At 700 BC the offset is 5.7 hours, which is 85° of the Earth's rotation, and that is why the ancient measurements are records of the place a total eclipse was seen from rather than of the hour: the hour was never written down accurately enough to matter, and the shadow's track on the ground was.

Timescales

  1. 1 Six kinds of second
  2. 2 A second that depends on where the clock is
2 essays · sky Two bodies at a mass ratio of 3 to 1. Both bodies orbit their common centre of mass, on similar ellipses whose sizes are in inverse proportion to the masses — here 3 to 1, so the heavier body's path is 3 times smaller.

The two-body problem

  1. 1 Neither body is still, and the wobble is how planets are found
  2. 2 Two bodies replaced by one that does not exist
2 essays · gravitation 7 decades of relaxation time, and a Hubble time between them. Crossing time and two-body relaxation time for 5 self-gravitating systems, against the number of bodies in each. The lower marks are the crossing time R/σ, which spans 3 decades; the upper ones are t_relax ≈ 0.1N/ln N times it, which spans 8. The horizontal rule is a Hubble time, and the same expression puts these systems on opposite sides of it. Two of them relax — an open cluster and a globular cluster — so their stars have exchanged energy, their heavy stars have sunk, and their present structure is not the one they were born with. The rest do not: a dwarf spheroidal, an elliptical galaxy, a cluster of galaxies, the elliptical missing it by a factor of 10⁶, which is why a galaxy is modelled as a smooth potential with no stars in it at all. The dependence is on N and hardly at all on anything else, because the Coulomb logarithm makes every decade of impact parameter contribute equally. A cluster that relaxes also evaporates: about a stellar mass in a hundred and forty leaves per relaxation time, so the globular cluster empties in 6.1·10¹¹ years.

Two-body relaxation

  1. 1 How long a system takes to forget
  2. 2 A cluster that boils itself away
2 essays · gravitation 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.

Universal variables

  1. 1 The one solve that does not ask which conic it is
  2. 2 The singularity that is a change of variable
2 essays · orbits Leavitt's law, through the Milky Way calibrators. Mean absolute magnitude against the logarithm of the pulsation period, for classical Cepheids with independently known distances, and RR Lyrae itself for comparison. The line is M = −2.81 log P − 1.43: a tenfold longer period is 2.81 magnitudes brighter, a factor of 13 in luminosity. The vertical axis runs the astronomers' way, with brighter upward.

Variable stars

  1. 1 A star that tells its distance by how slowly it blinks
  2. 2 The valve that has to sit at the right depth
2 essays · stars Mass from a velocity dispersion, across nine decades. The virial estimate M = 5 σ²R_e/G, drawn for half-light radii from 0.2 to 20 kpc, with four systems placed on it: M32 at σ = 75 km/s and R_e = 0.1 kpc, giving 6.5e+8 M☉; the Milky Way's bulge at σ = 105 km/s and R_e = 1 kpc, giving 1.3e+10 M☉; M87 at σ = 320 km/s and R_e = 7 kpc, giving 8.3e+11 M☉; the Coma cluster at σ = 1000 km/s and R_e = 1500 kpc, giving 1.7e+15 M☉. The same expression spans from a dwarf elliptical to a cluster of galaxies — nine decades of mass — because the virial theorem does not care what the moving objects are. In the last case they are whole galaxies, and the estimate exceeds everything visible in the cluster by a factor of about fifty, which is where the problem was first noticed.

Velocity dispersion

  1. 1 A galaxy held up by disorder
  2. 2 A dispersion inflated by orbits nobody resolved
2 essays · galaxies 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.

Virial theorem

  1. 1 An average that weighs what cannot be watched
  2. 2 The system that gets hotter as it loses energy
2 essays · gravitation 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.

Vis-viva

  1. 1 One equation for the speed anywhere, and the eccentricity is not in it
  2. 2 The same burn is worth more when moving fast
2 essays · orbits

All essays