Astronomy & Celestial Mechanics
the sky, computed
Generator

The blackbody generator

Blackbody curves at 3000, 5800, 10000 K
Blackbody curves at 3000, 5800, 10000 K. Thermal emission against wavelength, each curve scaled to its own peak so the shift can be seen on one plot. The peak moves to shorter wavelengths as the temperature rises, which is why colour is a thermometer.

Thermal emission against wavelength, each curve scaled to its own peak so the shift can be seen on one plot. The peak moves to shorter wavelengths as the temperature rises, which is why colour is a thermometer.

13 essays call blackbody. 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 magnitude scale, plotted. The logarithm of the received light against magnitude from -27 to 32, with the vertical scale left unlabelled because only its slope matters. The relation is a straight line of slope −0.4, which is what makes five magnitudes exactly a hundredfold — and the numbers run backwards, so brighter is smaller. 9 landmarks are marked, from the Sun to the deepest exposures, spanning a factor of 1.2·10²³ in received light. Starlight

A scale that runs backwards, multiplies, and works

Brighter stars have smaller magnitudes, and five steps is a factor of a hundred. A scale invented by eye in the second century BC turned out to be logarithmic, because eyes are.

 Blackbody curves at 3000, 5800, 10000 K. Thermal emission against wavelength, each curve scaled to its own peak so the shift can be seen on one plot. The peak moves to shorter wavelengths as the temperature rises, which is why colour is a thermometer. Starlight

Colour is a thermometer, and it reads across the galaxy

A star's colour gives its surface temperature, from two brightness measurements and no other information. It is the cheapest useful measurement in astronomy.

 The extinction law, for three kinds of dust. How much of a star's light dust removes, against inverse wavelength, normalised to one at the V band. Blue light is extinguished more than red, which is why reddening and extinction are the same measurement — and the hump at 4.6 inverse microns is a feature of the grains themselves, present on almost every sight line and still without an agreed carrier. Larger grains give a flatter law and a larger R_V. Starlight

Dust makes everything look further away

A star behind dust is fainter, so a distance taken from its brightness comes out too large. The same dust also makes it redder, and the reddening is measurable where the dimming is not — which is the only reason the correction can be applied at all.

 The Hertzsprung–Russell diagram. Luminosity against surface temperature, both in solar units and on logarithmic axes, with temperature increasing to the left. The main sequence is computed from the mass–luminosity and mass–radius relations; the dashed diagonals are lines of constant radius. Stars

The diagram that sorted the stars, by plotting two things against each other

Plot brightness against colour for a few thousand stars and they do not scatter. They fall on a narrow band with two islands off it, and explaining that structure is most of stellar astronomy.

 A spectrum belonging to no temperature at all. The disc's summed emission, with the individual annuli drawn faintly beneath it. Each ring is a blackbody at its own temperature, and each is drawn at the area it actually has — the outer rings are cool and enormous, the inner ones hot and small. The sum has three parts and only the two ends belong to a temperature: a Rayleigh–Jeans rise of slope 2 from the outermost ring, a Wien cutoff at the hottest, and between them a stretch of slope 0.316, against the 1/3 that comes out of integrating ν²T(r)r dr with T ∝ r^−3/4. That middle section is the observational signature of a disc: no single blackbody produces it, no photosphere produces it, and its width rather than its peak is what says how far in the disc goes. What the figure cannot show is that a real disc's innermost rings are neither thin nor blackbodies, which is where the model's clean edges stop. Stars

A spectrum that is a stack of temperatures

An accretion disc is not hot. Its inner edge is, its outer edge is not, and the temperature runs continuously between them as a power of radius — so what leaves the disc is the sum of a great many blackbodies at different temperatures, which is a spectrum with no temperature in it and a slope no single body can produce.

 A radius that depends on the colour it is measured in. Transit depth against wavelength for a planet of 1.38 Jupiter radii at 1400 K, whose atmosphere has a scale height of 538 km — computed from H = kT/µg, not assumed. One scale height of extra opacity adds 153 parts per million to a transit of 1.40 per cent, so the whole spectral signal is 862 ppm at its strongest: one part in 16 of the transit that carries it. The features are at real band centres — sodium at 0.589 µm, water at 1.4 µm, carbon dioxide at 4.3 µm — with the rise at the blue end the Rayleigh slope of scattering off the smallest particles. Exoplanets

A radius that depends on the colour it is measured in

Measure a transit in one colour and then another, and the planet is a different size. The difference is a few atmospheric scale heights, which is a few hundred parts per million of an already tiny signal.

 The habitable zone, and the radius inside which a day never ends. The conservative habitable zone — the runaway-greenhouse and maximum-greenhouse limits of Kopparapu's parameterisation — as a band in stellar mass against orbital distance, both logarithmic. For the Sun it runs from 0.99 to 1.71 AU, which is the published result and the check on the arithmetic here. The band moves inward far faster than the mass falls, because luminosity goes as roughly the fourth power of mass: a 0.2 M☉ star's zone is at 0.082–0.158 AU. The dashed line is the distance inside which a planet is tidally locked within 4.5 billion years, computed from τ = ω₀αmQa⁶/(3GM⋆²k₂R³) with an initial ten-hour spin, Q = 100 and k₂ = 0.3. It crosses the inner edge of the zone at 0.67 M☉ — so around every star below that, which is the great majority of stars, a planet in the habitable zone has one hemisphere in permanent daylight. The Earth is outside its own locking radius of 0.53 AU, and TRAPPIST-1e at 0.029 AU is inside its star's by a factor of 8. Exoplanets

One number sets the zone

The habitable zone is a band of stellar flux, so it scales as the square root of luminosity and moves inward far faster than mass falls. For most stars it lies inside the radius at which a planet is tidally locked.

 One whole orbit. The system's total brightness through one orbit: the transit at phase 0, the slow rise and fall of the planet's illuminated hemisphere between, and the secondary eclipse at phase 0.5 where the planet's own light is removed. The transit is 1.05%; the secondary eclipse is 1800 ppm, about 6 times shallower. Exoplanets

The planet is seen when it disappears

Half an orbit after the transit the planet passes behind its star, and the light that vanishes is the planet's own. Subtracting two brightnesses taken hours apart isolates a body nothing has ever resolved.

 The classical law gives every star the same colour. Planck's law and the Rayleigh–Jeans law at 3,000 K, 5,772 K, 10,000 K, both normalised to the 5,772 K Planck peak, on logarithmic axes. The classical law comes from counting standing waves in a cavity — 8πλ⁻⁴ of them per unit volume per unit wavelength — and giving each the kT that equipartition allows. It agrees with Planck's where the modes are crowded and each holds much less than kT, and it runs away where they are not: at 80 nm it exceeds the real spectrum by a factor of 1.1·10¹² while agreeing to within 55.9 per cent at 3000 nm, and the integral under it does not converge at all. That is the ultraviolet catastrophe, and it is the half everybody knows. The quieter half is that in 2ckT/λ⁴ the temperature is an overall factor, so the ratio of the law at two wavelengths is independent of it: the B − V index of a classical star comes out identical at 3,000 K, 5,772 K, 10,000 K — the same -0.968 magnitudes, to the last digit the quadrature carries — while Planck's law spreads the same three stars over 1.47 magnitudes. A classical universe has stars of every brightness and one colour. Colour is a thermometer only because the exponential in the denominator does not cancel, and the quantum of energy that put it there was fitted to this shape before anybody knew what it meant. Starlight

The classical law gives every star one colour

The ultraviolet catastrophe is the famous half. The quieter half is that in 2ckT/λ⁴ the temperature is an overall factor, so the ratio of the law at two wavelengths has no temperature in it — and a classical universe has stars of every brightness and one colour.

 Two laws that are the peak and the area of one curve. Planck curves at 3,000 K, 5,772 K, 9,600 K on logarithmic axes, with each peak marked. Normalised by its own peak, Planck's law is a universal function of x = hc/λkT, and three exponents follow from that alone and are fitted here off the drawn curves rather than quoted: the peak wavelength goes as T^-1.000, which is Wien's displacement law with a constant of 2.897772 mm K obtained by solving 5(1 − e^(−x)) = x for x = 4.965114; the peak HEIGHT goes as T^5.000; and the area goes as T^4.000. The third is the first two multiplied. A peak five powers high on a curve one power narrow encloses four powers of area, so Stefan–Boltzmann is not an independent fact about radiation — it is Wien's law and the height of the peak, taken together. That is also why the two are worth having at once. A colour gives the temperature and a flux gives the luminosity, and L = 4πR²σT⁴ then gives a radius: for the Sun at 5,772 K receiving 1361 W/m² at 1.000 AU, the arithmetic returns 6.957·10⁸ m against a measured 6.957·10⁸. A thermometer alone cannot do that, because a colour is a ratio and a ratio has no size in it; the radius comes from the one law that is an absolute quantity rather than a shape. Starlight

Two laws that are one curve read twice

Wien's displacement and Stefan–Boltzmann are the peak and the integral of the same function. The peak is five powers high and one power narrow, so the area is four — and the fourth power that is taught as a separate law is the first two multiplied.

 Three shifts larger than the error bar, and two that cancel. The blackbody B − V index against temperature, with three systematic shifts marked at 5,772 K. All three are computed from the same Planck integrals the relation itself is: a 3,800 K companion contributing 25 per cent of the V light reddens the index by 0.093 magnitudes, because the companion is relatively brighter in the redder band; 0.35 magnitudes of visual extinction at a total-to-selective ratio of 3.1 adds 0.113 directly, since the colour excess is the extinction divided by that ratio; and a metallicity of -1 dex subtracts 0.200, because the metal lines that eat the B band are the ones a metal-poor star is short of. Read as temperatures, the same star comes back at 5331 K, 5243 K and 7021 K against a true 5,772. Two of the three have opposite signs, and that is the difficulty rather than the relief. All three together shift the index by 0.005 magnitudes against 0.405 of combined magnitude — very nearly nothing, because the opposing pair removes almost all of it — and return 5744 K, which is 28 K from the truth by cancellation and not by accuracy. A reddened metal-poor star and an unreddened solar-metallicity one are the same point on this curve, and no amount of photometry in two bands separates them. What does separate them is a third band, or a spectrum — which is where the cheapest measurement in astronomy stops being cheap. Starlight

Three shifts larger than the error bar, and two that cancel

An unresolved companion, a reddening and a metallicity each move a colour index by more than any modern photometer's precision. Two of them move it in opposite directions, so the three together can return the right temperature by cancellation rather than by accuracy.

 How steeply the ratio of all the light to one band rises with temperature. The exponent n in (bolometric flux ÷ band flux) ∝ Tⁿ, for a blackbody, against temperature, for bands at V (0.55 µm), J (1.25 µm), H (1.65 µm), K (2.2 µm). The bolometric flux rises as T⁴; the flux in a band rises as T on the long-wavelength side of the Planck peak and much faster on the short side, so n is the difference. For the Sun's temperature it is 2.33 in K, 1.70 in J and −0.56 in V. An error of one per cent in the modelled band flux becomes an error of 1/n per cent in the temperature: 0.43 per cent, 25 K, in K. In the visible the ratio barely depends on temperature and carries no information about it; in the infrared it does, and the temperature is read with its model error divided. Starlight

A temperature read against the light's long tail

A star's total light depends on its temperature to the fourth power; its light in the infrared depends on it barely more than linearly, because the infrared lies on the long, flat tail of the Planck curve. The ratio of the two cancels the star's size and distance and leaves a steep function of temperature — and the only model it needs is the infrared surface flux, whose errors are divided by two and a half on the way to the answer. That is the infrared flux method, one of the two foundations of the stellar temperature scale, and the other foundation disagrees with it by up to a hundred and fifty kelvin.

 What a hundred kelvin moves, for a star like the Sun. The change in eight quantities derived for a Sun-like star and its planet when every temperature is raised by 100 K — 1.7% — as the uncertainty in the absolute temperature scale allows, against the precision each is quoted to (grey bars). The star's radius from its flux and parallax falls by 3.4%, because radius goes as the inverse square of temperature; the planet's radius follows, and its density rises by +12.8%. The oscillation-based mass rises by +2.6%, the main-sequence mass by +3.5% and the age falls by 9.8%. The iron abundance read from neutral iron lines rises by 0.08 dex, +20.7%. In every case the shift is comparable to or larger than the precision quoted, and in every case it moves all stars of a sample together. Starlight

A hundred kelvin under every radius

The temperatures of cool stars are quoted to tens of kelvin, and the scale they are all measured on is uncertain by about a hundred. That hundred is under almost everything else: a star's radius moves by twice its fractional error, a transiting planet's density by six times, a main-sequence age by ten per cent and an iron abundance by a fifth. And because a scale error moves every star of a sample the same way, no amount of averaging removes it.

The whole library · All essays