Astronomy & Celestial Mechanics
the sky, computed
Generator

The expansion-history generator

Four expansion histories that agree exactly today
Four expansion histories that agree exactly today. The scale factor against time, with the present at zero and every model normalised to a = 1 and an expansion rate of 67.36 km/s/Mpc there. That normalisation is the figure: four universes that are indistinguishable from a measurement made now, separated entirely by what is behind and ahead of them. Each curve is integrated from da/dt = aH₀E(a) rather than from its own closed form, so the four are compared through one routine. The age each implies is where its curve meets zero: empty (Ω = 0) 14.52 Gyr, matter only (Ω = 1) 9.68 Gyr, ΛCDM (Planck 2018) 13.80 Gyr, closed (Ω = 2) 8.29 Gyr. The empty universe's 14.52 Gyr is the Hubble time 1/H₀ exactly, which is what makes it the natural thing to measure an acceleration against. The closed model turns over and is stopped at its turnaround.

The scale factor against time, with the present at zero and every model normalised to a = 1 and an expansion rate of 67.36 km/s/Mpc there. That normalisation is the figure: four universes that are indistinguishable from a measurement made now, separated entirely by what is behind and ahead of them. Each curve is integrated from da/dt = aH₀E(a) rather than from its own closed form, so the four are compared through one routine. The age each implies is where its curve meets zero: empty (Ω = 0) 14.52 Gyr, matter only (Ω = 1) 9.68 Gyr, ΛCDM (Planck 2018) 13.80 Gyr, closed (Ω = 2) 8.29 Gyr. The empty universe's 14.52 Gyr is the Hubble time 1/H₀ exactly, which is what makes it the natural thing to measure an acceleration against. The closed model turns over and is stopped at its turnaround.

15 essays call expansion-history. The drawing above is what it returns with no arguments at all; every call below passes it something, because a placement that passes nothing draws whichever member of the family the generator happens to default to rather than the one its essay argues about.

Where it is called

Every figure listed here is the same construction drawn at different numbers, so a correction to one is a correction to all of them.

The same recession law from two different galaxies. Twenty galaxies at fixed comoving positions after the whole picture has been multiplied by 1.34, with each galaxy's displacement drawn from where it was to where it is. Every arrow's tail sits at the old separation and its tip at the new one, so the arrow is exactly 0.34 times its tail's distance from the highlighted galaxy — proportional to separation for one reason and no other: a uniform scaling moves everything in proportion to its distance from whatever point the scaling is measured about. The right-hand panel measures from a different galaxy and gets the identical law with the identical constant. That is the content of a linear velocity–distance relation. It is the signature of an expansion with no centre, and the observation that every galaxy recedes is therefore not evidence that this one is the centre — it is evidence that none of them is. Cosmology

The shift that is not a Doppler shift

Every galaxy beyond the Local Group has its lines shifted to the red, by an amount proportional to its distance. Read as a velocity that looks like a confession that everything is fleeing from here; read as a change of scale it says the opposite, because a uniform expansion produces exactly the same law measured from any galaxy in it.

 Four expansion histories that agree exactly today. The scale factor against time, with the present at zero and every model normalised to a = 1 and an expansion rate of 67.36 km/s/Mpc there. That normalisation is the figure: four universes that are indistinguishable from a measurement made now, separated entirely by what is behind and ahead of them. Each curve is integrated from da/dt = aH₀E(a) rather than from its own closed form, so the four are compared through one routine. The age each implies is where its curve meets zero: empty (Ω = 0) 14.52 Gyr, matter only (Ω = 1) 9.68 Gyr, ΛCDM (Planck 2018) 13.80 Gyr, closed (Ω = 2) 8.29 Gyr. The empty universe's 14.52 Gyr is the Hubble time 1/H₀ exactly, which is what makes it the natural thing to measure an acceleration against. The closed model turns over and is stopped at its turnaround. Cosmology

One number that is an age, a size and a density

The slope of the velocity–distance relation has units of inverse time, so it can be read as an age, multiplied by c to give a length, or squared to give a density. All three readings are natural, all three are quoted, and not one of them is the quantity it appears to be.

 Distance modulus against redshift, for three universes. The distance modulus μ = 5 log₁₀(D_L/10 pc) against redshift for three universes, all with H₀ = 67.36 km/s/Mpc, with 60 model supernovae drawn from the ΛCDM curve with 0.15 magnitudes of scatter. The point of the figure is how little difference there is: across two decades of redshift the three curves stay within a few tenths of a magnitude, and at z = 0.5 the accelerating and decelerating cases differ by 0.387 mag. A cosmology is not read off this plot. It is read off the residual, which is the next figure. Cosmology

An expansion that was supposed to be slowing

Gravity is attractive, so an expanding universe full of matter must be decelerating, and the only question was by how much. Two teams set out to measure the deceleration and both found a quarter of a magnitude of extra faintness at redshift half — which is the wrong sign.

 The same universe, four times, as a fraction of itself. Each bar is the fractional contribution of the four components to the total density at one epoch, computed from the Planck 2018 parameters by scaling each component from today: radiation as (1+z)⁴, both kinds of matter as (1+z)³, and Λ as a constant. The familiar figure — five per cent baryons, twenty-six dark matter, sixty-nine dark energy — is the top bar and only the top bar. At recombination the same universe is three-quarters dark matter and Λ is one part in ten million; before matter–radiation equality it is mostly radiation. A pie chart of the contents of the universe is therefore a statement about a moment, and the moment is the one it happens to be drawn in. Cosmology

A budget whose familiar part is five per cent

Five per cent ordinary matter, twenty-six per cent dark matter, sixty-nine per cent dark energy. The figures are quoted everywhere and each one comes from a different measurement, the denominator they are fractions of is itself built out of the expansion rate, and the whole chart is a statement about one instant that was a different chart at every earlier time.

 Three horizons, and a light cone that bulges. Cosmic time upwards, comoving distance sideways, for the Planck 2018 cosmology. Galaxies sit still in these coordinates, so their worldlines are the vertical grey lines: comoving distance is defined to take the expansion out. The solid inner curve is the past light cone — the set of events whose light reaches here and now — and it reaches out to the particle horizon at 46.1 billion light years, which is the diagram's central number and the one that sounds impossible. The universe is 13.80 billion years old and light has travelled for 13.80 billion years, and yet the material that emitted the oldest light is now 46 billion light years away. Nothing has outrun light: the comoving distance covered is ∫c dt/a, and dividing by a scale factor that was small early on makes the integral three times ct. The dashed curve is the Hubble sphere, where the recession speed equals c, at 14.5 Gly today — and the light cone lies outside it for most of its length, which is exactly why a galaxy can be observed while receding faster than light. The outer dot-dash curve is the event horizon, at 16.7 Gly: a signal sent from here today never reaches anything beyond it. Cosmology

A horizon three times larger than the age allows

The universe is 13.8 billion years old and light travels one light year a year, so the observable universe should be 13.8 billion light years across. It is 46.1, nothing has outrun light, and the discrepancy is a piece of arithmetic rather than a paradox.

 w = −0.9 is 34 millimagnitudes, and one supernova scatters by 120. Above: how much the distance modulus moves when the dark energy is not a constant. Each curve is a universe with the same Ωₘ = 0.315 and a different equation of state w, drawn as a difference from w = −1 in magnitudes. At redshift a half, w = −0.9 is worth 34 millimagnitudes — the shaded band is the 0.12-magnitude intrinsic scatter of a single standardised type Ia supernova, and the signal is a fifth of it. Nothing about one object can see this; the measurement is the mean of 1500, whose error on the mean is 3.1 millimagnitudes, and even that only works because the shape of the curve in redshift is different from every systematic anybody has thought of. Below: why the supernovae are not enough on their own. Each locus is the set of (Ωₘ, w) that a measurement cannot tell apart from the fiducial model — computed, not sketched: the supernova curve is the ridge of the same sum of squares a fit would minimise over 0.02–1 in redshift, and the acoustic-scale curve is the exact set of models with the same comoving distance to last scattering, which is what fixes the angle the microwave background's first peak subtends. They cross at 36 degrees. Neither is a measurement of w and the pair is, which is why the constraint on the equation of state is a picture of two loci crossing rather than a number read off a curve — and why −1.03 ± 0.03 is a statement about how well they cross rather than about how well anything was measured. Cosmology

The number that would say whether it is a constant

Whether dark energy is a cosmological constant is the question of whether w is exactly −1, and w = −0.9 changes a distance modulus by thirty-four millimagnitudes at redshift a half — a fifth of the scatter of a single supernova, along a degeneracy only the acoustic scale can cut across.

 A lumpy universe makes more of both, and one of them was already too much. The deuterium–lithium plane, with the curve a homogeneous universe traces as its baryon density varies and the points a two-zone universe reaches at a fixed mean density of η₁₀ = 6.13. The dense zone occupies 15 per cent of the volume, and the contrast between the zones runs from 1 — which is the homogeneous case — to 100. Every mixture lies up and to the right of its own homogeneous point, and that is not a modelling choice: an abundance is measured per baryon, so what a telescope averages is η times the abundance — and for both deuterium and lithium that product is a convex function of the density, whose average therefore exceeds its value at the average. The excess is large. At a contrast of 100 the mixture gives D/H = 1.13e-4 against 2.51e-5 smooth, and ⁷Li/H = 1.77e-8 against 4.70e-10, a factor of 37.7. The extra deuterium is exactly what the proposal was for: it lets the mean baryon density be raised while the observed D/H is still matched, which in the 1980s was the one way to make the baryons account for all the matter that dynamics required. The extra lithium is what it costs. The measured abundance was already a factor of 2.9 below the homogeneous prediction, and every step toward the lumpy universe that fixes the density makes that discrepancy worse — which is why the answer to "the baryons are lumpy" turned out to be that the extra matter is not baryons. Cosmology

The universe that was lumpy at one second

If the baryons were unevenly spread when the network fired, each region ran its own nucleosynthesis and what is observed is an average. For a decade that was the one way to make ordinary matter account for all the matter — and the reason it fails is a theorem about convex curves.

 Whether an event horizon exists at all, against one number. The comoving event horizon today — the distance a signal sent now will ever cover — against the equation of state of the dark energy, for a flat universe with the measured matter density. The curve runs away at w = −1/3 and does not exist above it: that is where the expansion stops accelerating, and in a universe that does not accelerate the integral ∫da/a²E diverges and every galaxy is eventually reachable, however far away. Below −1/3 the horizon is finite and shrinks as w falls, because a more negative equation of state makes the dark energy density grow with time rather than stay constant. At the cosmological constant's w = −1 the horizon is 16.7 billion light years against a particle horizon of 46.1, so 4.7 per cent of the volume now observable is still reachable. The band is the measured −1.03 ± 0.03. What the figure is for is the asymmetry in what the measurement still allows: two sigma toward zero puts the horizon at 17.5 Gly and two sigma the other way at 14.7, and the shape of the curve means that the closer the true value sits to −1/3 the more violently the answer moves. The reachable fraction is not a robust number in the way the particle horizon is. Cosmology

Whether there is a horizon at all

An event horizon exists precisely when the expansion accelerates, and its size is not a smooth function of how much. The integral that defines it runs away as the equation of state approaches minus a third, so two sigma either way on a measured number are two very different futures.

 How much a galaxy's redshift changes in 10 years, and which way. The change in a galaxy's apparent recession velocity over 10 years of the observer's time, c ż/(1 + z), against the galaxy's redshift, for three universes with the same present expansion rate of 67.36 km/s/Mpc. The drift is (1 + z) H₀ − H(z): positive if the expansion rate at the galaxy's epoch was less than (1 + z) times today's, which is to say if the expansion has been accelerating since. In the empty universe the expansion rate is exactly (1 + z) H₀ at every epoch and nothing drifts at all. In the matter-only universe the expansion has only ever slowed, and every redshift falls: −8.6 cm/s over 10 years at z = 1 and −25.5 at z = 4. In ΛCDM the drift is positive nearby, largest at z = 0.63 where it reaches 2.51 cm/s, changes sign at z = 1.91, and is −5.5 cm/s at z = 4. The whole signal is a few centimetres per second in a decade, against the thirty kilometres per second of the Earth's own orbital motion that has to be removed from every spectrum first. Cosmology

A redshift that changes while it is watched

A galaxy's redshift is a ratio of two sizes of the universe, and the second one is still growing while the light is being collected. So every redshift drifts, by a few centimetres per second in a decade, and the direction of the drift says whether the expansion has been speeding up since the light left — the one test of that question that needs no distance and no model of any source.

 The expansion rate and the acceleration, and which of them reaches inside an orbit. Two quantities per unit distance through cosmic time, both in units of today's H₀². The square of the expansion rate, H², falls steeply from the big bang and is 1.00 today by definition. The acceleration of the expansion, ä/a, is the sum of a matter-and-radiation part, −Ωₘ/(2a³) − Ωᵣ/a⁴, drawn dashed, and the cosmological constant's part, ΩΛ = 0.685, which is the same at every time. The sum was negative — the expansion decelerating — until the universe was 7.7 Gyr old, at a = 0.614 or redshift 0.63, and is 0.527 today. The equation of motion of anything orbiting inside a bound system carries ä/a and never H: the rate at which distant galaxies recede does not appear in it at all. Of ä/a, the matter part is the mean density of the universe, which inside a galaxy or a planetary system is already counted in the mass that is doing the holding, many million times over. What is left is the constant: a fixed outward acceleration per unit distance, ΩΛ H₀², that does not grow with time and does not care how fast the universe is expanding. Cosmology

An orbit feels the acceleration and never the rate

If space expands, it is natural to ask why the Earth's orbit does not. The answer is not that gravity resists the stretching. It is that the expansion rate never appears in the equation of motion of a bound orbit at all — only the acceleration does, and of that only the cosmological constant's part survives, as a fixed outward push that moves the Earth's orbit once, by twelve picometres, and never again.

 Two galaxies that stopped moving apart, and the mass that did it. The separation of two galaxies on a radial orbit that began together at the big bang, against cosmic time, solved so that after 13.797 Gyr they are 770 kpc apart and approaching at 110 km/s — the present separation and approach speed of the Milky Way and the Andromeda galaxy. The curve is a cycloid, r = A(1 − cos θ) and t = B(θ − sin θ), and only one cycloid passes through that point with that slope. It rose to 1037 kpc, turned round when the universe was 8.5 Gyr old, and on this purely radial orbit the two meet 3.3 Gyr from now. Its period fixes the mass: A³/(GB²) = 4.2·10¹² solar masses. The dashed curve is the same calculation with the cosmological constant's outward push included, integrated rather than solved; to arrive at the same place at the same speed against that push it needs 4.76·10¹² solar masses, 13 per cent more. Far more than the stars of the two galaxies, it is the timing argument's measurement of the Local Group's dark matter. Cosmology

The age of the universe weighs the Local Group

The Andromeda galaxy is approaching the Milky Way, and in an expanding universe that means the two once moved apart, stopped and turned round. One radial orbit passes through their present separation with their present speed after exactly the age of the universe, and its period fixes the mass that turned them — four trillion suns, twenty-five times what their stars can account for.

 The dipole a 370 km/s motion has to put into counts of distant sources. The amplitude of the dipole in the number of sources per unit solid angle expected from the Sun's motion at 369.82 km/s relative to the microwave background, β = 1.234e-3, split into its two parts: aberration, 2β, which crowds sources towards the direction of motion, and the Doppler boost, x(1 + α)β, which brightens sources there and lifts fainter ones above the flux limit. For radio sources, with counts steepening as S^(−1) and spectra falling as ν^(−0.75), the expected dipole is 0.0046; for mid-infrared quasars, with counts steepening as S^(−1.7) and spectra falling as ν^(−1.26), the expected dipole is 0.0072, and the measured one is 0.0155 — 2.16 times larger, which would need a speed of 797 km/s. A dipole of a few parts in a thousand needs a catalogue of more than a million sources to see at all. The disagreement is not with the direction, which lies close to the microwave background's, but with the size, and it is not yet explained: either the samples carry a systematic nobody has found, or the matter and the radiation do not share one rest frame on these scales — in which case the assumption that the universe looks the same from everywhere is wrong in a way it has never been caught being wrong before. Cosmology

The Sun's speed counted in quasars comes out twice too large

The microwave background is warmer in one direction by a part in a thousand, and that dipole is read as the Sun's motion through it at 370 kilometres per second. The same motion must crowd the counts of distant galaxies and quasars towards the same direction by an amount that can be calculated exactly. When a million quasars were counted, the direction agreed and the size came out more than twice too large — as though the Sun were moving at 800 kilometres per second relative to the matter.

 A delay that counts the gas nothing can see. The mean dispersion measure of a radio pulse against the redshift it comes from, for several shares of the baryons residing in diffuse ionised gas between galaxies. A pulse is delayed by free electrons in proportion to the column it crosses, and that column is an integral over the expansion history of a density the baryon budget fixes — so the only unknown in the whole expression is the share itself. At z = 1 the relation gives 1085 pc cm⁻³ if every baryon is out there and 543 at 50 per cent. The measurement runs the other way: a burst with a known host redshift and a measured delay returns the fraction, and the answer came out consistent with nucleosynthesis. A pulse a millisecond long weighs the half of the ordinary matter that no survey could find, and it does so because the thing that delays it is the thing that does not shine. Cosmology

Half the ordinary matter was missing, and a millisecond found it

Nucleosynthesis fixes how many baryons there are to better than a per cent. Every survey of where they are came up about half short for two decades — and what closed the gap was the delay a radio pulse picks up crossing gas too thin and too hot for any telescope to have seen.

 A mass measured by the structure it prevents. The fractional suppression of the small-scale matter power spectrum against the sum of the neutrino masses. The relic density follows from the thermal history with no astrophysics in it — the neutrinos’ share of the critical density, times h², is the mass sum divided by 93.14 eV exactly — and the suppression is about eight times the neutrinos' share of the matter, because a particle moving at a large fraction of the speed of light streams out of an overdensity while it is still growing and takes its own gravity with it. The two mass orderings put a floor under the sum at 58 meV and 100 meV; the cosmological upper bound is near 0.12 eV, which is 7.2 per cent of suppression. The floor and the ceiling are within a factor of two of each other, so this is a measurement about to happen or a model about to break, and the quantity it returns is a sum rather than any individual mass. Cosmology

A mass measured by what it stopped from forming

Neutrinos were relativistic in the early universe and are not now, so they are the one entry in the cosmic budget that changes category. What cosmology measures is not their density but the hole they leave — they stream out of a growing clump and take their gravity with them, and the missing structure bounds a particle mass more tightly than any laboratory has.

 Structure stops forming above about 14 times the observed Λ. The fraction of matter that ever collapses into a bound object, against the cosmological constant in units of the observed one, holding the primordial fluctuation amplitude fixed. Growth of structure stops once Λ dominates the expansion, so the linear growth factor approaches a finite limit rather than rising for ever, and a larger Λ freezes it earlier and smaller. The asymptotic growth factor at the observed Λ is 1.110; at a hundredth of it, 5.152; at a hundred times, 0.239. The collapsed fraction falls from 0.86 to 0.40 to 0.000 across the same range, and passes a tenth of its present value at 14 times the observed constant. The observed value sits a factor of 14 below the largest one that leaves anything at all — which is the whole of the anthropic argument, and that factor is what has to be compared against the 10¹²⁰ by which the naive theoretical estimate misses. Cosmology

A coincidence that is a factor of fourteen

The cosmological constant is famously wrong by a hundred and twenty orders of magnitude, and famously coincidental in sitting near the matter density now. Computing how large it could be and still leave any structure at all turns the second complaint into a number — and the number is fourteen, not a hundred and twenty.

The whole library · All essays