Astronomy & Celestial Mechanics
the sky, computed
Depth

Series — page 9

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.
Io's measured heat needs k₂/Q = 0.016, and Enceladus's needs 0.011. Tidal surface heat flux against orbital eccentricity, from Ė = (21/2)(k₂/Q)GM_p²R⁵ne²/a⁶ evaluated at each satellite's own orbit, drawn at a common k₂/Q of 0.015. Every line has slope 2 because the dissipation is quadratic in e and nothing else on this axis varies. The filled marks are each body at its actual eccentricity; the two ringed ones are the bodies with a measured surface heat flux, and they are the only points here that are observations. Solving each of those for k₂/Q gives 0.016 for Io and 0.011 for Enceladus — within a factor of 1.5 of one another, for a warm silicate body and a 500-kilometre ball of ice, which ought to be a coincidence and is instead the sharpest problem in the subject: nothing about Enceladus's ice can plausibly dissipate at 0.011, and the number is what the heat requires all the same. The dashed line is the Earth's measured surface heat flux, 0.087 W/m², which Io exceeds by a factor of 28 — the most volcanically active body in the solar system is the fourth largest moon of the fifth planet, and the reason is entirely in the orbit.

Tidal heating

  1. 1 A moon heated by not being allowed to relax
  2. 2 One heat flow, and two viscosities
  3. 3 A heat flow that depends on a number nobody can compute
  4. 4 A migration timed by the planet and not the moon
4 essays · gravitation Tully–Fisher, from two assumptions and no fitting. 44 model galaxies spanning two and a half decades in luminosity, each built from a constant disc surface brightness of 380 L☉ per square parsec with 0.13 dex of scatter and a stellar mass-to-light ratio of 1.4 with 0.1 dex. Nothing about a luminosity–speed relation is put in: the scale length follows from the surface brightness, the mass from the light, and the speed from v² = GM/R. The construction makes L ∝ v⁴ exactly, so the true slope is −10.0 magnitudes per decade of speed. Least squares of magnitude on log speed returns -8.45, and the reverse regression -9.15: the scatter here is scatter in the speed at fixed luminosity, and error in the abscissa flattens a fitted slope, so neither fit returns the value the family was built with and the distance between them is the size of the effect. Residual scatter about the drawn line is 0.54 magnitudes. Observed slopes run from about −7.5 in blue light to −10 in the near infrared, where the mass-to-light ratio is steadiest — which is the same statement as the second assumption above.

Tully–Fisher

  1. 1 A line width that is a distance
  2. 2 The speed a line width stands in for
  3. 3 A disc that turns slower than its mass requires
  4. 4 A stirring supernovae cannot pay for
4 essays · galaxies 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
  3. 3 A few thousand clouds outrun a hundred billion stars
  4. 4 Torques that add because the orbits hold still
4 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
  3. 3 Four numbers that carry a state without an orbit
  4. 4 A check that cannot see the clock
4 essays · orbits 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
  3. 3 The one radius where the orbits stop mattering
  4. 4 Motions across the sky measure the orbits first
4 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
  3. 3 An equilibrium reached without a single encounter
  4. 4 A relaxation that stops before it is finished
4 essays · gravitation A coherent one-per-cent distortion, invisible on every galaxy in the picture. 150 background galaxies behind a lens of Einstein radius 14″, each drawn at its own ellipticity: an intrinsic shape with a dispersion of 0.3 per component, plus the reduced shear the lens adds. The strongest shear on any galaxy here is 0.035, one part in 8 of the intrinsic scatter, so no object in this field is measurably distorted and the tangential alignment cannot be seen by eye at all. Averaged over these 150, the mean tangential ellipticity is 0.0088 ± 0.0245 against the 0.0130 the lens model predicts — consistent with the lens and equally consistent with nothing, because 150 galaxies buy a precision of 0.024 and the signal is 0.013. Detecting it at five sigma takes about 13,275 of them, which is not a picture anybody can draw. The cross component — every shape rotated by 45°, which gravitational lensing cannot produce — averages −0.0016 ± 0.0245, consistent with nothing, and that null is what separates a mass from a badly figured optic. The signal is not in any galaxy. It is in the sum, and the whole design of a lensing survey follows from that.

Weak lensing

  1. 1 A one-per-cent distortion, and a million galaxies to see it
  2. 2 Two parameters that lensing measures as one
  3. 3 A measurement of dark matter that needs a model of galaxies
  4. 4 A lens counted instead of measured
4 essays · galaxies Mestel cooling for white dwarfs of 0.4 to 1 solar masses. Luminosity against age for white dwarfs of 0.4, 0.6, 0.8, 1 solar masses, both axes logarithmic. Each line is the whole of the star's later life: there is no nuclear source, so what is radiated is the thermal energy the ions had when the star was made, leaking out through a thin envelope whose opacity is Kramers'. That one sentence gives L ∝ t^(−7/5), and the lines are straight here because a power law is straight on these axes — the slope is the exponent and can be measured off the picture. Two features are worth more than the numbers. The first is how flat the sequence becomes: a factor of ten in luminosity costs a factor of 5.2 in age, so a white dwarf spends most of its existence faint, and any population of them piles up at the bottom. The second is that the heavier tracks lie above the lighter ones. A heavier white dwarf is smaller, so it has more ions to cool and less surface to lose them through, and at a given age it is the brighter object. The picture does not show what happens at the faint end, where the ions crystallise and the heat capacity stops being constant; that is a separate figure, and it is where this law stops being enough.

White dwarf cooling

  1. 1 A clock with no fuel in it
  2. 2 A clock that stops while its interior freezes
  3. 3 A colour that runs backwards at the end of cooling
  4. 4 A star that distils itself as it freezes
4 essays · stars The field a collapse would arrive with, and the one it has. Field strength against hydrogen density, both logarithmic, over the eight decades between diffuse gas and a protostellar core. The steeper line is what perfect flux freezing demands: a sphere collapsing conserves both mass and flux, so B goes as R⁻² while ρ goes as R⁻³, and therefore B goes as the two-thirds power of the density exactly — an exponent with no free parameter in it. The shallower locus is what Zeeman measurements find: a flat branch at about 10 µG up to 300 per cubic centimetre, where the density is rising and the field is not, and a rise as the 0.65 power of the density above it. By the density of a core the two differ by a factor of 1, and a star built at the frozen-flux value would carry a field four orders of magnitude beyond anything measured on one. The flat branch is the important half: it says the gas is moving along field lines without dragging them, which is what gravity does to a cloud that is still magnetically supported, and it locates where the freezing has to break.

Flux freezing

  1. 1 A field that would have arrived ten thousand times too strong
  2. 2 The same flux in three sizes
  3. 3 A field that waits for the core to freeze
3 essays · stars Two populations, and the valley between them. 420 model galaxies, 42 per cent of them drawn from a tight red sequence tilted by -0.08 magnitudes of colour per magnitude of brightness, the rest from a broad blue cloud. Only 3.3 per cent land in the band of width 0.20 magnitudes midway between the two sequences, against 17 per cent in the same band laid on the blue cloud and 27 per cent on the red sequence. The bimodality is the strongest statement the galaxy census makes: a galaxy is usually forming stars or has stopped, rarely in between, and since colour tracks the age of the newest stars, the emptiness of the middle is a statement about a timescale — whatever ends star formation does it in a few hundred million years, not a few billion.

Galaxy populations

  1. 1 Two colours, and almost nothing between
  2. 2 Two ways to stop that do not interact
  3. 3 A crossing time that is the filter's
3 essays · galaxies A halo is born with a spin of a few hundredths. The distribution of the dimensionless spin parameter λ = J|E|^½ ÷ (G M⁵ᐟ²) across dark matter haloes, drawn as a lognormal of median 0.035 and logarithmic width 0.5, with the disc scale length each λ implies printed along the lower axis. The distribution is required to integrate to one and to peak at 0.0272, which is the median times e raised to minus sigma squared, and is the signature of a lognormal rather than of a bell curve drawn to look like one. λ is small because a halo is supported by random motion rather than by rotation: a value of 0.035 means the halo turns at about three and a half per cent of the rate it would need to hold itself up centrifugally. It is also nearly independent of halo mass, which is what points at a common origin. Mapping it to a disc through R_d = λ R₂₀₀/√2 with all of the specific angular momentum retained, a 10¹²-solar-mass halo of radius 206 kiloparsecs gives 5.1 kiloparsecs at the median, against the 2.6 kiloparsecs the Milky Way's disc actually has. The gap is not a failure of the estimate; it is the measurement that the baryons arrived with less spin per unit mass than the halo they arrived in.

Galaxy spin

  1. 1 A disc the size its halo was born with
  2. 2 Two spins that cancel, and an orbit that leaves
  3. 3 A halo with too much slow gas for a disc
3 essays · galaxies Inside a star that is holding itself up (polytrope n = 3). Temperature, density and pressure through a star, as fractions of their central values, against fractional radius. All three come from one numerical integration of the Lane–Emden equation at index 3, which is the hydrostatic balance written for a gas whose pressure is a power of its density. The inner half of the radius holds 90% of the mass, and the outer half is nearly weightless — which is why the load, and so the temperature, is concentrated where the burning is.

Hydrostatic equilibrium

  1. 1 A star is held up by its own weight
  2. 2 A floor under the centre that assumes nothing
  3. 3 Mechanics sets a star's temperature before fusion does
3 essays · stars The Sun's internal rotation: differential above 0.693 R, rigid below, and the 15 nHz between two modes. Angular velocity against fractional radius, in nanohertz, from the inversion of hundreds of thousands of measured rotational splittings. Two features had to be discovered rather than deduced. The convection zone rotates differentially in latitude — 473 nanohertz at the equator, 330 at sixty degrees — and that latitude dependence persists all the way down through it, in surfaces that are very nearly radial rather than the cylinders a rotating convective fluid was expected to produce. And below about 0.693 of the radius the latitude dependence stops: the radiative interior turns as a single rigid body at 430 nanohertz, which is a remarkable thing for a fluid with no strength to do, and requires something — most likely a weak internal magnetic field — to be enforcing it. The two regimes are joined by a shear layer a few per cent of the radius thick, the tachocline, and it is where the solar magnetic cycle is generally thought to be generated, because it is the only place in the Sun with the shear a dynamo of that strength needs. The horizontal marks are what individual modes would report: the rotation averaged over the cavity each one occupies, weighted by the time the wave spends at each radius. A mode of degree 100 is trapped near the surface and reports 473 nanohertz; one of degree 1 passes through the deep interior and reports 458. The 15-nanohertz difference between them is the measurement, and splittings are measured to about a nanohertz — which is why the profile can be resolved at all. The kernels used here have the right support and the right sense; a real inversion uses computed eigenfunctions, and its resolution below 0.2 R is poor for the reason the previous figure gives.

Internal rotation

  1. 1 A fluid that turns as one piece
  2. 2 A shear layer that should have spread
  3. 3 A dynamo that did not need the shear layer
3 essays · stars Measured moment-of-inertia factors, from the Sun's 0.07 to the Moon's 0.3931. Eleven bodies whose interiors have never been sampled, arranged by the one interior quantity that has been measured for all of them. C/MR² is 2/5 for a uniform sphere and falls as mass is concentrated toward the centre, and the values here span from 0.07 to 0.3931. That spread is the content. The Moon at 0.3931 is barely differentiated — whatever iron core it has is a few per cent of its radius, which is why the Moon is the one large body in the inner solar system without a magnetic field of its own. Mercury at 0.346 is nearly as low as the Earth despite being an eighth of its mass, and for a body that small the only way to get there is an iron core filling most of the radius. The Sun at 0.07 is off the scale of anything a two-layer model describes; a star is not a planet with a bigger core but a body whose density falls by five orders of magnitude between centre and surface. The faint curves behind are the two-layer relation at a few density contrasts, drawn to show what kind of interior each value is consistent with — and the horizontal placement of each body on them is an illustration rather than a result, since one factor never fixes one core. Every number here was obtained by watching the body turn: a precession rate, a libration amplitude, a gravity field sampled on a flyby, or in the Sun's case the frequencies of its own oscillations.

Moment of inertia

  1. 1 Whether the heavy material sank
  2. 2 The part of a shape the spin cannot explain
  3. 3 A mountain's gravity says what holds it up
3 essays · gravitation 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.

Oblateness

  1. 1 The Earth's shape, read off a satellite's node
  2. 2 The field a satellite is allowed to feel
  3. 3 A core weighed by something that never went in
3 essays · gravitation A spot is dark because it is squeezed. Pressure against depth below the quiet photosphere's optical surface, drawn as a ratio to the pressure there. The rising curve is the surrounding gas, which grows exponentially with a 140-kilometre scale height because that is what hydrostatic equilibrium in an ideal gas produces. The flat pair of bands is the spot's own budget: a magnetic pressure of 3.58·10⁵ dyn/cm² from a 3000-gauss field, which is 25 per cent of the total, plus the gas pressure left over. Horizontal balance requires the two columns to reach the same total at the same geometric level, and the level at which they do is 348 kilometres below the quiet surface — so the spot's own optical surface sits in a hollow. Measured Wilson depressions, obtained from the foreshortening of a spot near the limb, are four to six hundred kilometres. Nothing about the darkness was assumed: a field strength read off a Zeeman splitting fixes the magnetic share, the share fixes the depression, and the depression fixes the temperature the deeper layer must have to carry the reduced flux.

Plasma beta

  1. 1 The darkness a field pays for
  2. 2 A corona held by its field and opened by its gas
  3. 3 A wind that carries the ledger of its loops
3 essays · stars 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
  3. 3 The part of a delay no orbit can imitate
3 essays · gravitation 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
  3. 3 A disc weighed by how far its stars climb
3 essays · galaxies 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
  3. 3 A Sun read through a shield that also changed
3 essays · stars How long a star lasts, against its mass. Main-sequence lifetime against mass, on logarithmic axes. Fuel grows in proportion to mass and consumption grows as its three-and-a-half power, so the lifetime falls steeply — a star of thirty solar masses lives for a few million years.

Stellar lifetimes

  1. 1 Why the biggest stars die first, and take the galaxy with them
  2. 2 A lifetime that is still a forecast
  3. 3 A lifetime that stops getting shorter
3 essays · stars

All essays