Astronomy & Celestial Mechanics
the sky, computed
Concept

Red noise — where it appears

Noise whose power rises towards low frequencies, so successive measurements are correlated rather than independent. It makes a false-alarm probability computed for white noise far too optimistic, and stellar activity, pulsar timing noise and instrumental drift all have this character.

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

Nothing visible in any pulsar, and a quadrupole in the angle between them. Above: 4 millisecond pulsars' timing residuals over 15 years, at the few hundred nanoseconds a good one reaches. Each wanders, and none of them shows anything a reader could call a signal; a gravitational-wave background of amplitude 2.4·10⁻¹⁵ at one cycle per year contributes a common part to all of them that is smaller than each pulsar's own red noise. Below: the correlation between pairs, against the angle on the sky between them. 2211 pairs out of 67 pulsars, binned into 15 angles, against three curves with no free parameters between them. A quadrupolar background gives the Hellings–Downs shape — positive for nearby pulsars, negative near 83°, and back up to exactly half its zero-separation value at 180° because a background looks the same in opposite directions. An error in the observatory clock would give a flat line, because it shifts every pulsar identically. An error in the solar-system ephemeris would give a cosine, because it moves the barycentre in one direction. The drawn points prefer the quadrupole over the flat line by Δχ² = 358. That is the detection: not a waveform, not an event, not a moment — a shape in an angle, accumulated over fifteen years, on data taken for another purpose entirely.

A detector the size of the galaxy

At a nanohertz no instrument can be built, so the clocks already in the sky are used instead. The signal is in no single pulsar's data — it is in the correlation between pairs as a function of the angle between them, and that curve has no free parameters at all.

gravitation · Gravitational waves A true peak at 0.3103 and a false one at 0.6897, from the same data. The Lomb–Scargle periodogram of 162 simulated observations taken over 419 nights from one site, of a star carrying a 4.2-unit sinusoid at 0.3103 cycles per day under 3 units of Gaussian noise per point. The injected signal is recovered at 0.3103 cycles per day, within the 0.0024 resolution element the baseline allows. The second peak, at 0.6899, is 87 per cent as tall and corresponds to nothing: it is 1 − f, the signal reflected in the one-cycle-per-day spike of the sampling. Neither peak is more real than the other in this picture — deciding between them needs a second site at a different longitude, or a run long enough for the seasonal window to separate them. The dashed line is the power a pure-noise series would exceed once in 1000 trials, computed from 586 independent frequencies rather than from the 4691 grid points searched; using the grid count instead would put the line 2.08 higher and reject a real detection.

A period found in the gaps

Astronomical time series are sampled when the sky is dark and clear and the target is up, which is a schedule with a spectrum of its own. That spectrum is convolved with the real one, so a single sinusoid produces several peaks — and the tallest is not always the true one.

starlight · Periodograms Three ways for a timing model to be wrong, and three shapes that say which. Timing residuals over 6 years for the Crab pulsar, one curve per kind of error in the model, in microseconds. A position error of 1.2 mas leaves a sinusoid of period exactly one year — measured off the drawn curve as 1.000 — with amplitude (a/c)·δθ·cos β = 2.7 μs, because the error is being projected onto a baseline that is the Earth's own orbit and nothing about the pulsar. An unmodelled proper motion of 0.9 mas/yr leaves the same sinusoid with an envelope growing linearly: twice as large at 6 years as at 3. An error of one part in 10⁹ in Ṗ leaves a parabola, the second integral of a frequency drift, whose second derivative is constant to 4e-12 across the span — and that fractional error is deliberately minute, because anything larger produces a residual thousands of times the other two and draws them as flat lines. The shapes do not resemble each other, which is the whole reason a pulsar is an instrument rather than a clock: fitting them simultaneously delivers a position, a proper motion and — from the annual curvature term, not drawn here — a parallax, all from the arrival times of pulses and no image of anything. What is left when every known shape has been removed is the science: glitches, red noise, and the correlated residual between pairs of pulsars that a timing array exists to find.

A clock read against a model of everything in between

A pulsar delivers arrival times and nothing else. Everything else is a model, and what is measured is the difference between the model and the arrivals — a residual whose shape says which term is wrong, and whose annual sinusoid is a position measured from pulses rather than from an image.

stars · Pulsars A peak worth 10.8 in a narrow search is worth nothing in a wide one. The probability that noise alone produces a peak at least as tall as a given power, for searches over four different numbers of independent frequencies. A single frequency examined in isolation gives a one-per-cent chance at a power of 4.6; searching fifty thousand frequencies for the same one-per-cent chance requires 15.4. The threshold rises as the logarithm of the width of the search, which is why the penalty is survivable — but it is a penalty, it is often not applied, and the number of independent frequencies in an unevenly sampled time series is not the number of frequencies on the grid. Overestimating that count is conservative and underestimating it is not, which is the one asymmetry worth remembering.

The tallest peak in nothing at all

A periodogram of pure noise has peaks in it, and the tallest is not small. How tall it has to be before it means something depends on how many frequencies were searched and on what the noise actually is — and astronomical noise is almost never the white noise the standard formula assumes.

starlight · Periodograms

Named alongside it

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

False-alarm probabilityIndependent frequenciesA pulsar timing arrayTiming residualAliasBarycentric correctionBootstrapBox least squaresCharacteristic strainChirp massCorrelationDetection threshold

All concepts