Astronomy & Celestial Mechanics
the sky, computed
Concept

Power spectrum — where it appears

How much variance sits at each frequency or scale, whether in the comb of a star's oscillation modes or in the clustering of cosmic matter. It discards all phase information, so two very different-looking fields can share one — which is why the microwave background's non-Gaussianity is measured separately.

Named by 11 essays across 4 fields — each of them below, with the objects they name alongside it.

14 radial orders of the Sun, at 135.1 μHz apart. The p-mode spectrum of the Sun — 1 solar mass in 1 solar radius — from the asymptotic relation with its second-order term, drawn as 14 radial orders of ℓ = 0, 1 and 2 under a Gaussian envelope centred on ν_max = 3,090 μHz. Two numbers are marked and they do very different work. The large separation, 135.1 μHz, is the spacing between consecutive ℓ = 0 modes and fixes the mean density. The small separation, 9.00 μHz, is 2.7 pixels on this axis — it fixes the age, and it is why the échelle diagram exists rather than being a convenience. The vertical axis is the measurement: each mode moves the surface by about 20.0 cm s⁻¹ at the peak, and brightens it by a few parts per million, which is why this was impossible before a decade-long velocity series. Each mode is one line: its true width is set by its lifetime and is far below a pixel here.

The interior read from a comb of frequencies

A star's surface moves by about twenty centimetres a second, in thousands of overlapping sound modes at once. Two numbers off that spectrum give a mass and a radius with almost no stellar model in the chain, and a third gives an age.

stars · Asteroseismology 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.

cosmology · Large-scale structure The damping tail measures a thickness: ℓ_D = 1400 for a shell 80 deep in redshift. The suppression of small-scale structure in the microwave background, against multipole, on logarithmic axes. Every other feature of the power spectrum measures the last-scattering surface as a surface — its distance, the sound horizon written on it, the ruler it provides. This one measures how thick it is. Recombination takes time: the ionised fraction falls over a range of redshift rather than all at once, and while it is falling the photons are still scattering, so each one random-walks. A photon taking N steps of a mean free path λ diffuses √N λ, which is much further than a single step and much less than the whole interval, and any temperature fluctuation smaller than that distance is mixed away before it can be frozen in. What is left is a Gaussian cut-off, drawn here for three shell thicknesses. A thicker shell means more steps and a longer walk, so it damps at a smaller multipole: Δz = 40 gives ℓ_D = 1980, Δz = 80 gives ℓ_D = 1400, Δz = 160 gives ℓ_D = 990. The observed cut-off is near ℓ = 1400, corresponding to a diffusion length of about 0.03 proper megaparsecs at the time — a scale the reader should compare with the sound horizon, some hundred and fifty comoving megaparsecs, which is what the peaks measure. The tail is therefore a genuine probe of the inside of the transition rather than of its position, and because the damping depends on the free-electron density it is also one of the cleanest constraints on anything that changes it: extra relativistic species, a varying fine-structure constant, or an early energy injection all move ℓ_D while leaving the peak positions nearly alone. The normalisation here is set once by the observed value, so what the figure claims is the scaling with thickness and the shape of the cut-off, not the absolute number.

A blur that measures a depth

Every other feature of the microwave background measures the last-scattering surface as a surface — its distance, the ruler written on it, the geometry between. The damping of its small-scale structure measures how thick it is, because a photon random-walking through a finite transition smears away anything smaller than its walk.

cosmology · Microwave background One stream after three different populations of dark subhaloes. The density of stars along the middle sixteen kiloparsecs of a thin stream, relative to its mean, after four billion years of encounters with dark subhaloes drawn from a mass function falling as the −1.9 power of the mass between 10⁵ and 10⁸ solar masses. Top: the abundance cold dark matter predicts — in this realisation 26 encounters, 4 of them above a million solar masses — which leave 32 per cent of the stream's length more than 30 per cent below the mean density. Middle: 0.3 of that abundance, 8 encounters and 7 per cent of the length in gaps. Bottom: none at all, where the only variation is the counting noise of the 60,000 model stars. No single dip in the top trace is a clean gap: they overlap, and a large one and several small ones are superposed at every point. What differs between the three is the statistics of the whole trace, not any feature in it.

A spectrum that counts what no single gap can

One gap in a stellar stream measures a subhalo's mass multiplied by an unknown time. A stream crossed by dozens of them measures something better — how many subhaloes of each size there are — but only as a statistic of the whole stream's density, where no feature can be pointed to and one stream is one throw of the dice.

galaxies · Stellar streams A velocity error that grows with distance. The uncertainty in one galaxy's peculiar velocity, v = cz − H₀d, against its distance in Mpc/h, when the distance comes from an indicator with a fractional error of 5 per cent, 20 per cent, 25 per cent (type Ia supernovae, 5%; Tully–Fisher, 20%; fundamental plane, 25%). The redshift is exact to a few kilometres a second, so the whole error is the distance error times H₀, and it grows in proportion to the distance. The horizontal line is the rms velocity of a galaxy along one direction in linear theory for the model drawn, 302 km/s. Each galaxy's own velocity is measured to better than its own size only inside 60 Mpc/h for 5 per cent, 15 Mpc/h for 20 per cent, 12 Mpc/h for 25 per cent; beyond that a single velocity is mostly noise, and a velocity survey is useful only as an average over many galaxies.

Velocities that see past the edge of the map

A redshift survey infers the motions of galaxies from the distortions of its own map. A peculiar-velocity survey measures them directly, one galaxy at a time, by subtracting a distance from a redshift — and every one of those measurements is dominated by its own error beyond a few tens of megaparsecs. What makes the survey worth doing is that a velocity answers for structure far larger than the volume it is measured in.

cosmology · Large-scale structure The same tilt, and ice ages of two different lengths. Top: the obliquity anomaly in degrees from the five largest terms of the standard series, over 2 million years of model time; it is dominated by the 41,000-year cycle and has no power at all near 100,000 years. Bottom: ice volume in a threshold model driven by it — the ice grows steadily, and it collapses whenever volume plus a multiple of the obliquity anomaly crosses a threshold, marked by dots. For the first 1 million years the threshold is low and the ice collapses at nearly every obliquity maximum, on average every 43 thousand years. At the vertical line the threshold is raised, and nothing else changes: the ice now survives one or two obliquity maxima and collapses at the second or third, every 108 thousand years on average. The forcing is identical on both sides; the ice sheet counts it differently.

An ice sheet that skips a beat of the tilt

The last eight hundred thousand years of ice ages come roughly every hundred thousand years, and the Earth's orbit has no forcing of any strength at that period. A model in which an ice sheet collapses only when it is both large and warmed by a high-obliquity summer reproduces the rhythm with nothing but the forty-one-thousand-year tilt — the ice sheet grows through one or two maxima and falls on the next — and a slow rise in how large the ice must grow turns the older forty-one-thousand-year world into the newer one without any change in the orbit.

orbits · Orbital averages 40 galaxies detected, and a map of all of them: the summed light traces the structure the catalogue misses. A 160 by 160 Mpc/h slice of a simulated universe at redshift 6, its galaxies clustered on a density field and given luminosities from a Schechter function with a faint-end slope of -1.5. Left: the 40 galaxies bright enough to detect individually, above 1.5 times the characteristic luminosity. Right: the same slice as a map of summed line intensity in 32 by 32 cells, shaded linearly between the 5th and 97th percentiles of the cells, with no galaxy detected at all. The catalogue is a sparse scatter in which the structure is barely visible; the map shows it, and its cells correlate with the underlying density at 0.65. The map's light comes overwhelmingly from galaxies below the detection limit.

A map of galaxies none of which was seen

A survey that lists galaxies one by one stops at the faintest it can detect, and at high redshift most of the light comes from galaxies fainter than that. Map one spectral line in summed intensity instead, pixel by pixel and frequency by frequency, and every galaxy contributes whether it is detected or not. The fluctuations of that map have two parts — clustering, which weighs the total light, and graininess, which weighs how unevenly it is shared — so one power spectrum measures two moments of a luminosity function nobody has counted.

cosmology · Olbers's paradox [C II] power is round in the plane of k; CO(4–3) mapped at the wrong distance is elongated 5.1 times along k∥. Contours of the clustering power spectrum in the plane of wavenumber across the sky, k⊥, and along the line of sight, k∥, in h/Mpc, at 10², 10²·⁵, 10³, 10³·⁵ (Mpc/h)³ in arbitrary units, for the target line, [C II] from redshift 6.60 (black), and for the CO(4–3) interloper from redshift 0.84 placed at the target's distance (colour). The target's contours are quarter-circles, because its structure has no preferred direction. The interloper's are ellipses, drawn in towards the k∥ axis and elongated along it by the ratio of the two rescalings, α∥/α⊥ = 5.1: a structure that is round where it really is looks magnified across the sky and flattened along the line of sight where the map puts it. That difference of shape, not of amplitude, is what separates the two in a map where both have been summed.

A line from the wrong redshift lands at the wrong distance

A map of ionised carbon from redshift six is built by assuming every photon at a given frequency came from that line. Carbon monoxide from galaxies at redshift one sends photons to the same frequency, and the map places them five billion light-years too far away. Their structure arrives magnified across the sky and squashed along the line of sight — round where it really is, flattened in the map — and that distortion of shape, which no amount of brightness can hide, is how the impostor is found.

cosmology · Olbers's paradox The same ruler, and a wave that arrives early. The acoustic oscillation in the galaxy power spectrum — the ratio to a smooth spectrum, minus one — against wavenumber, for no free-streaming species and for three neutrinos, with the sound horizon held at 99.0 h⁻¹ Mpc in both so that only the phase differs. The two agree on the largest scales and separate as the wavenumber grows: by 0.3 h/Mpc the wave with neutrinos runs 0.21 radians ahead, so its crests sit at smaller wavenumber — near 0.2 h/Mpc by 2.2 thousandths of an h/Mpc. The amplitude 0.05 and the damping of the feature by bulk flows over 7 h⁻¹ Mpc are illustrative; the phase is computed.

A wave pulled ahead by particles that outran it

The neutrinos of the early universe moved at the speed of light while the sound wave that became the acoustic feature moved at little more than half of it, and their gravity dragged the wave forward. The drag survives as a phase in the galaxy power spectrum, a shift that no change in the ruler's length can copy, so galaxy surveys can count the species that streamed freely without asking the microwave background.

cosmology · Baryon acoustic oscillations How deep a hole the Hubble tension needs, and how deep holes are. The typical departure of the matter density from the mean, averaged over a sphere — one standard deviation, σ, from the linear power spectrum with σ₈ = 0.81 — against the sphere's radius, both on logarithmic axes. It falls from 0.55 at 20 Mpc to 0.022 at 300 Mpc. A local expansion rate 8.3 per cent above the global one, the size of the Hubble tension, needs the sphere the supernovae sample to be underdense by 0.47 (dashed), with a growth rate of 0.53. Across the radii the supernovae span (shaded, 100–600 Mpc) that is 4 to 67 standard deviations.

A hole too deep for the sky to have made

If the Milky Way sat inside a vast underdense region, everything around it would be expanding faster than the universe as a whole, and the locally measured Hubble constant would come out high — which is what Cepheids and supernovae find. The idea needs no new physics. It needs a void about half empty across the six hundred megaparsecs the supernovae sample, and the same cosmology that predicts the microwave background prices such a void at dozens of standard deviations.

cosmology · Hubble constant A clump that does not grow until the radiation is outweighed. The amplitude of a dark-matter overdensity inside the horizon against the scale factor in units of its value at matter–radiation equality, both logarithmic, from y = 10⁻⁴ to today near y = 3400. The solution of the growth equation, δ = 1 + 3y/2, stays flat until equality and then rises in proportion to the scale factor; matter domination alone would have grown it in proportion to the scale factor from the start. Over the whole range the true growth is 5.1e+3 against 3.4e+7 — radiation withheld a factor of 6666, which is the interval from y = 10⁻⁴ to about y = 0.67.

A clump that waits for the radiation to thin

Dark matter feels no pressure from light, so nothing should have stopped it clumping in the radiation era. Almost nothing did stop it, and yet a clump inside the horizon at that time barely grows. The radiation does not push on it. It only makes the universe expand too fast for the clump's own gravity to keep up — and the length of that wait is the bend in the spectrum of every structure in the sky.

cosmology · Large-scale structure

Named alongside it

The objects these essays reach for when they reach for this one.

Cosmic varianceGalaxy biasShot noiseCold dark matterGrowth rateLine-intensity mappingPeculiar velocityReionisationSound horizonAlcock paczynski testAnisotropic stressAsteroseismology

All concepts