Astronomy & Celestial Mechanics
the sky, computed
The collection

Every essay — page 12

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 221–240 of 514.

Galaxies

Where the unknown stops being a number and becomes a profile — and most of it is not light.

A delay of 81 days, and a sheet nobody can see that moves H₀ to 82.4. Above: the arrival-time surface of a lensed source, along the line through the lens. The curve is the Fermat potential in days — the geometric cost of taking a longer path, minus the gravitational cost of climbing out of the potential — and the images sit at its stationary points, at -1.20″ and 2.04″, which for an isothermal sphere is β ± θ_E. Fermat's principle is doing all of the work here: light does not take the shortest path or the quickest one, it takes every stationary one, and the number of images is the number of stationary points. The vertical distance between the two is 81 days, and it is measurable — the source is a quasar, quasars vary, and the same wiggle appears in one image and then the other. That single number carries an absolute distance: the delay is D_Δt/c times a dimensionless function of the lens model, and D_Δt goes as 1/H₀, so a monitoring campaign gives the Hubble constant with no rung of any ladder beneath it. Below: the two light curves, shifted by exactly that delay. The dashed curve is the second image with the delay removed, and the agreement is the measurement. What the picture also shows is the reason the answer keeps moving. The second arrival-time curve is the same lens with a uniform sheet of convergence added and the source moved to compensate: every image sits at the same place, every flux ratio is the same, every image shape is the same, and the delay is λ = 0.85 times as long. A lens model fitted to positions alone cannot see the sheet, and inferring H₀ from the same delay under it gives 82.4 instead of 70 — a 15 per cent shift with no observable attached. Breaking it needs a mass measured some other way: the velocity dispersion of the deflector, or a count of everything else along the line of sight.

A distance measured with a stopwatch

The images of a lensed quasar sit at the stationary points of an arrival-time surface, and the height between two of them is a delay in days. That delay is proportional to a distance, and the distance is proportional to one over the Hubble constant — so a flickering quasar gives H₀ with no ladder under it.

8 figures · Strong lensing Two densities in thermal balance at one pressure, 121 times apart. Thermal equilibrium for interstellar gas, drawn as pressure against density with both axes logarithmic. Every point on the curve is a temperature between 40 and 9000 K at which cooling exactly balances the 2·10⁻²⁶ erg per second per hydrogen nucleus that grain photoelectrons deliver. The curve is not monotonic: it rises to 5007 K cm⁻³, falls to 1597, and rises again, so a horizontal line anywhere between those two crosses it three times. At the 3000 K cm⁻³ of the local medium the three crossings are a warm phase at 0.47 cm⁻³ and 6354 K, a cold phase at 57 cm⁻³ and 52 K, and one in between drawn dashed because it cannot survive — a parcel there that is squeezed cools and keeps contracting, and one that expands heats and keeps expanding. The two survivors differ by a factor of 121 in density and by exactly the same factor in temperature — necessarily the same, because their product is the pressure both are held at — and yet they press on one another equally, which is why they can share the same volume of the disc indefinitely rather than mixing.

Two temperatures, and nothing in between

The gas between the stars is not one gas at one temperature. Thermal balance has three solutions at any ordinary interstellar pressure, the middle one cannot survive being nudged, and the two that can differ by a factor of over a hundred in both density and temperature — which is why a map of the interstellar medium looks like clouds rather than like fog.

8 figures · Interstellar medium The critical mass falls as the cloud contracts, with a slope of −0.50. The Jeans mass — the least mass of molecular gas that its own pressure cannot support — against number density, for gas at 10, 20, 50 K, both axes logarithmic. The curves are straight because the mass goes as T^(3/2) ρ^(−1/2), and the slope measured off the drawn 10 K curve is −0.500 against the −0.5 that exponent requires. The consequence is the one that matters and it is a matter of sign: a cloud collapsing at constant temperature moves to the right along one of these lines, so the mass it takes to be unstable keeps falling, and sub-regions that were individually stable when the collapse began become individually unstable during it. A giant molecular cloud at 100 particles per cubic centimetre has a Jeans mass of 99 solar masses; a dense core at 10⁵ has one of 1.70. That is why a cloud of ten thousand solar masses makes a cluster rather than a star, and why the question about star formation is not what makes gas collapse but what stops the fragmenting.

The cloud that cannot hold itself up

A cold cloud collapses when gravity beats pressure, and the mass at which that happens falls as the square root of the density — so the collapse makes the condition for collapse easier, over and over, until the gas can no longer get rid of the heat. That is why a cloud of ten thousand solar masses makes a cluster and not a star.

9 figures · Molecular clouds One star in 380 makes essentially all the ionising light. Three cumulative fractions against stellar mass, for a broken power-law initial mass function with slopes 1.3 and 2.3 breaking at 0.5 solar masses. Each curve says what share of one quantity is produced by stars heavier than the mass on the axis, and the three do not resemble one another. Only 0.41 per cent of the hydrogen-ionising photons come from stars below 15 solar masses, because the ionising output of a star climbs by five orders of magnitude between eight and twenty. One star in 380 is above that mass, and between them those stars hold 14 per cent of the mass. Those two numbers are the leverage in every star-formation rate quoted from an Hα line. What is measured is the light of a handful of very massive stars; what is reported is the mass of a whole population; and the number in between is an integral over a part of the mass function that no extragalactic observation reaches. The medians are marked but should be read with care, and the reason is visible in the curves: the mass-weighted median at 1.27 solar masses is a property of the population, while the light-weighted one at 58 is a property of where the plot stops — halving the upper mass limit moves it to 36. An integrand that rises with mass has its median wherever the axis ends.

A birth rate measured from light nothing young emitted

A star-formation rate is quoted in solar masses per year, and nothing in the measurement counts a star or weighs anything. What is measured is a luminosity produced by about one star in four hundred, and the conversion to a mass is an integral over a part of the mass function that no extragalactic observation reaches.

9 figures · Star formation

Cosmology

One object, seen once, from inside — and every number in it the output of a model.

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.

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.

8 figures · Expansion 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.

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.

8 figures · Hubble constant H₀: nine determinations in two families. Published determinations of H₀, each with its quoted one-sigma interval, sorted into two families — measured locally, calibrated by a ladder, against inferred from z ≈ 1100 through a model. The shaded band behind each family is that family's inverse-variance weighted mean: 72.66 ± 0.75 across 5 of them, against 67.40 ± 0.41 across 4. The difference is 5.26 ± 0.85 km/s/Mpc, which is 6.2 standard deviations, computed here from the quoted errors alone. That number is an upper bound on the significance rather than the significance: the determinations within each family share calibrations, samples and in two cases the same supernovae, so they are not independent, and a correlated pair combines to something wider than the formula used here gives. What the figure does establish is that the split is not one discrepant measurement against a consensus — it is two internally consistent groups, and the grouping is by method rather than by result.

The same constant, measured twice, five sigma apart

The distance ladder gives an expansion rate of about 73 kilometres per second per megaparsec. The microwave background gives 67.4. Both quote errors near one per cent, both have been rebuilt from scratch by rival teams, and the gap between them has grown as the measurements have improved.

8 figures · Hubble constant Four distances to the same galaxy. Four quantities all called "the distance", against redshift, at a Hubble constant of 67.36 km/s/Mpc, a matter density parameter of 0.3153 and a dark-energy parameter of 0.6847. They agree below z ≈ 0.1 and then part company completely. Comoving distance is the separation now, and it is what a map of the universe is drawn in. Luminosity distance is what a brightness gives, and it is larger by (1+z) because the photons arrive both redshifted and spread out in time. Light-travel distance is the age difference times c, and it is bounded by the age of the universe. Angular-diameter distance is what an angle gives, and it is the odd one: it rises, turns over at z = 1.59 where it reaches 1.79 Gpc, and falls thereafter. Past that redshift a galaxy of fixed size looks bigger the further away it is, because the universe it is being seen across was smaller when the light left it. At z = 10 the four differ by a factor of 121 between the largest and the smallest, so a sentence quoting a cosmological distance without saying which one has not given a number.

Four distances to the same galaxy

Inside the Local Group the word "distance" has one meaning. Past a redshift of about a tenth it has four, they disagree by factors of a hundred by the time the light is old, and one of them stops increasing and starts coming back.

9 figures · Distance ladder 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.

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.

9 figures · Dark energy A blackbody at 2.7255 kelvin, filling the sky. The Planck function at 2.7255 K in the units the measurement is reported in, with the peak marked where Wien's law in frequency puts it: x = hν/kT = 2.8214, so ν = 160.2 GHz and the intensity there is 384 MJy per steradian. The points are drawn at the twenty-one frequencies across the FIRAS band, displaced from the curve by a Gaussian of 50 parts per million of the peak, which is the root-mean-square deviation the instrument actually reported. At the scale of this plot that displacement is a fifth of a pixel and the points sit on the line — which is the entire finding. Nothing else in astronomy is a blackbody to a part in twenty thousand: a stellar spectrum is a blackbody with absorption lines cut into it and a continuum that is the wrong shape at both ends. A thermal spectrum this exact requires that the radiation was once in equilibrium with matter, which requires that the universe was once opaque, which requires that it was once hot and dense.

The most perfect blackbody ever measured

The whole sky glows at 2.7255 kelvin, and its spectrum matches the Planck curve to fifty parts in a million. Nothing else in astronomy is thermal to a part in twenty thousand, and a spectrum that exact is not a description of the radiation — it is a constraint on the history of everything that could have disturbed it.

8 figures · Microwave background The acoustic peaks, and where the geometry says they should be. The temperature angular power spectrum of the microwave background. The drawn curve is a monotone interpolation through the published positions and heights of the six peaks and five troughs of the Planck 2018 TT measurement — it is a representation of data, and nothing between two extrema is claimed. The marks along the top are not: they are computed from this collection's own cosmology as ℓₐ(m − 0.267), where ℓₐ = π × 13866 / 144.43 = 301.6 — π times the comoving distance to last scattering divided by the sound horizon there — is the angle the sound horizon subtends at last scattering turned into a multipole. The two agree to 2.8 per cent at worst across six peaks, which is the whole of what makes this a measurement of geometry: a wave of known physical wavelength, seen at a known distance, is a protractor. The first peak at ℓ = 220 corresponds to about 0.82 degrees on the sky — roughly twice the width of the full Moon, which is the largest hot and cold patch the sky has.

A standing wave frozen at one instant

The temperature of the microwave background varies across the sky by one part in a hundred thousand, and the sizes of the patches are not random. There is a preferred angular scale near one degree, and it is a sound wave that stopped ringing four hundred thousand years after the beginning.

9 figures · Microwave background Four abundances, one free parameter. The abundances big-bang nucleosynthesis predicts, against the one number it is free to choose: η₁₀, the ratio of baryons to photons in units of 10⁻¹⁰. Four curves spanning nine decades, from a helium mass fraction of about a quarter down to a lithium abundance of one atom in ten billion, and they are not four independent predictions — they all come out of the same reaction network run at the same density. The horizontal bands are what is measured in the sky, each at its published one sigma. The measurement that matters is deuterium, because its curve is the steep one: inverting the drawn curve at D/H = 2.527e-5 gives η₁₀ = 6.11, and the ends of the observed interval give 6.06 to 6.15. The vertical band is what the microwave background gives, 6.13 ± 0.04, from the height of the second acoustic peak relative to the first. Those two agree to 0.4 per cent, and they have nothing whatever in common: one is a nuclear-reaction network run in the first three minutes and read off a quasar absorption line, the other is a fluid oscillation at four hundred thousand years read off a sky map. Lithium is the exception and it is not a small one — the network predicts 4.70e-10 at the microwave background's density and the halo stars show 1.60e-10, a factor of 2.9 too much, which is unresolved.

Four abundances and one free parameter

In the first three minutes the universe ran a nuclear reaction network with exactly one adjustable number in it. That number predicts four abundances spanning nine orders of magnitude, three of them are observed and match, and the fourth is wrong by a factor of three and has been for twenty-five years.

8 figures · Nucleosynthesis 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.

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.

8 figures · Density parameters A bump at a hundred megaparsecs. The two-point correlation function of galaxies, multiplied by the square of the separation so that the interesting part is not buried under the power law. The smooth dashed curve is the broad-band clustering — the ordinary fact that galaxies are near other galaxies, which carries no cosmological information and is treated as a nuisance term in the real analysis. The bump on top of it is the whole measurement, and its position here is not fitted: it is the comoving sound horizon at the drag epoch, integrated in this file from ∫c_s da/a²H with c_s = c/√(3(1+R)), which comes out at 146.9 Mpc — 99.0 in the h⁻¹ Mpc a survey works in. The excess says that a galaxy is very slightly more likely to have a companion at that separation than at 90 or 115, by about one part in two hundred, and the reason is that a pressure wave in the photon–baryon fluid ran outward from every overdensity for four hundred thousand years and stopped where it was when the photons let go. The points are drawn with the error a survey of a million galaxies achieves; a single pair of galaxies at 100 Mpc means nothing, which is why the measurement waited for the surveys.

The same ruler measured twice, ten billion years apart

The sound wave that produced the acoustic peaks in the microwave background also left a faint excess in how galaxies are spaced, at a separation of about a hundred megaparsecs. It is the only cosmological distance indicator whose length is set by physics rather than by a chain of calibrations.

8 figures · Baryon acoustic oscillations Every shell contributes the same. Concentric shells of equal thickness around an observer, with stars scattered uniformly in volume. The number of stars in a shell grows as its radius squared and the flux from each falls as the radius squared, so the two cancel exactly and every shell delivers the same total light — the drawn counts are 6, 15, 31, 55, 88, in proportion to r³ − r₀³, and the drawn sizes fall as 1/r. That cancellation is the paradox, and it is why no amount of dust helps: dust absorbs the light and then re-radiates it, and in a universe old enough for the sum to converge it would come to the same temperature as the stars. The sum diverges linearly with radius, so something has to stop it — and the only two candidates are that the shells eventually overlap, or that there are no shells beyond a certain distance because there has not been time for their light to arrive.

Why the sky is dark

In an infinite universe of stars every sight line ends on a stellar surface, so the whole sky should be as bright as the Sun's disc. It is not, and the usual answer — that the expansion redshifts the light away — accounts for a factor of six out of a hundred trillion.

8 figures · Olbers's paradox 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.

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.

9 figures · Horizons One length, leaving and returning. The comoving Hubble radius c/aH against the scale factor, both logarithmic. Everything to the right of the kink is exact and is the same cosmology as every other figure here: after inflation the comoving Hubble radius grows, as a in the radiation era and as √a in the matter era, then turns over once Λ takes hold. Everything to the left is a schematic — the energy scale of inflation is unmeasured, so neither the height of the plateau nor the 62 e-folds drawn is a number anybody has — but the shape is not negotiable: in any accelerated expansion H is nearly constant, so c/aH falls as 1/a. The consequence is the mechanism. The solid horizontal line is a fixed comoving length of 209 Mpc: it starts inside the Hubble radius, is carried outside during inflation, sits frozen while nothing can act across it, and re-enters at z = 1090, which is last scattering — it is the scale the first acoustic peak is made of. The dashed line above it is the comoving size of the whole observable universe, 14148 Mpc, and the figure's quiet second finding is that it is still outside the Hubble radius: the very largest angular scales in the microwave background have never re-entered, which is why they show the primordial spectrum with no acoustic processing on it at all. One crossing does two jobs: it makes the sky uniform, because everything now visible was once in causal contact, and it makes it not quite uniform, because a quantum fluctuation stretched beyond the horizon has nothing left that can smooth it out.

Two coincidences with one mechanism

Opposite sides of the microwave sky were never in causal contact and have the same temperature to one part in a hundred thousand. The total density sits on the knife edge of flatness, which is an unstable equilibrium. Two fine-tunings of quite different kinds, and one epoch of accelerated expansion removes both.

10 figures · Inflation A slice 1,000 megaparsecs across. 872 galaxies in a wedge 1000 comoving megaparsecs deep, generated as a Thomas process — Poisson centres at 0.00012 per square megaparsec, 13 galaxies each on average, scattered about their centre with a Gaussian of 22 Mpc. That is not a simulation of how structure formed and does not pretend to be; it is the simplest clustered process whose correlation function is known exactly, which is what lets the next figure check a measurement against something rather than against itself. What it does reproduce is the impression: at this scale the distribution is obviously not uniform, there are groups and there are gaps, and the eye finds patterns in it readily — including several that are not there, because the eye finds patterns in a Poisson field too. The question the next two figures ask is the only one that settles it: at what separation does the clustering stop, and is there a scale above which a box of this universe looks like any other box?

Homogeneous above a hundred megaparsecs

Every cosmological calculation in this field assumes the universe is the same everywhere. Looked at on any scale a person can picture, it plainly is not — it is stars in galaxies in groups in clusters in filaments around voids. The assumption is not a claim about appearance; it is a claim about a statistic, and the statistic has a scale attached.

8 figures · Large-scale structure One ionising photon per atom, and where that happens. The number of photons energetic enough to ionise hydrogen, per baryon, against temperature — with temperature falling to the right, so the universe ages to the right. The quantity is the fraction of a Planck distribution above 13.6 eV, integrated numerically, divided by the 6.13·10⁻¹⁰ baryons there are per photon. It crosses one at 5,844 K, which is z = 2143, and it falls by 19 orders of magnitude across the plot because it is the tail of an exponential. That crossing is the answer to why recombination waits until three thousand kelvin. At hydrogen's own ionisation temperature of 157,803 K there are two billion ionising photons per atom and the gas has no chance; the temperature has to fall by a factor of forty before the supply runs out, and it runs out suddenly because an exponential tail does. The number that sets the factor of forty is not an energy at all — it is the photon-to-baryon ratio, which is to say the entropy per baryon, which is to say a number fixed long before any of this and measured today from the second acoustic peak and from deuterium alike.

The surface the background actually is

Hydrogen ionises at 157,800 kelvin and the universe became transparent at 3,000. The factor of fifty between them is not an error, and it is not about energy — it is about there being two billion photons for every atom, so the far tail of the distribution can keep the gas ionised long after the typical photon has become useless.

9 figures · Microwave background A neutral fraction of 2.4·10⁻⁶ is already opaque. The Gunn–Peterson optical depth against the neutral fraction of the intergalactic medium, at z = 3, 5, 6.3, for Ω_b = 0.0493 and h = 0.674. Note the range of the vertical axis. A fully neutral medium at z = 6.3 gives τ = 4.1·10⁵, which is not absorption but extinction of everything; the medium reaches τ = 1 — the point at which it stops transmitting most of the light — at neutral fractions of 6.1·10⁻⁶ at z = 3, 3.3·10⁻⁶ at z = 5, 2.4·10⁻⁶ at z = 6.3. That is why the argument runs from the flux that survives rather than from the flux that does not. A spectrum showing any transmission at all between Lyman α and Lyman β is a measurement that the medium is ionised to better than one part in 164,727, and no fit to any absorption line is needed to establish it.

A trough that proves the forest survived

A uniform neutral medium at redshift six would absorb Lyman alpha with an optical depth of four hundred thousand. So the existence of any transmitted light in a quasar's spectrum is a measurement — of a neutral fraction below one part in ten thousand — made from a detection rather than from an absorption.

8 figures · Reionisation

Every sequence · Every named object · Search