Astronomy & Celestial Mechanics
the sky, computed
Concept

Opacity — where it appears

How strongly matter blocks radiation per unit mass, which decides how steep a temperature gradient is needed to carry a given luminosity. It rises steeply where a species is partly ionised, and those narrow windows are the only places in a star where a pulsation can be driven.

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

Opacity against temperature, and the three things that supply it. The Rosseland mean opacity of a gas of composition X = 0.7, Y = 0.28, Z = 0.02, on logarithmic axes, at 10⁻⁷ g/cm³ and 10⁻⁶ g/cm³. The three faint curves are the separate processes at the first density — electron scattering, the Kramers bound-free and free-free term, and the negative hydrogen ion — and the solid curve is their sum. Every one of them is multiplied by the fraction of hydrogen the Saha equation says is ionised at that temperature and density, or by one minus it for H⁻, which is the only reason the low-temperature end is a picture of a star rather than of a formula outside its range: ungated, Kramers alone gives 10,343 cm²/g at 5,800 K, against the 0.40 drawn here. The peak sits at 15,400 K, where hydrogen is 78% ionised — that bump is not a detail, it is the engine of a Cepheid — and the flat floor at high temperature is electron scattering, which is the one term with no temperature in it at all.

The surface that is a depth

A star has no surface. What looks like one is the level at which the optical depth reaches about two-thirds — and the edge is sharp only because the opacity climbs so steeply that the transition takes a ten-thousandth of the radius.

starlight · Opacity Two means of one opacity, 75 times apart, and only the smaller one is in the equation. Above: a synthetic opacity across the frequencies that carry a star's flux, drawn against x = hν/kT, with a continuum falling as ν⁻³ and a forest of 6 lines per unit x on top of it. The shaded curve is the Rosseland weighting function, ∂B_ν/∂T, which peaks at x = 3.83 and is what decides which frequencies matter. The two horizontal lines are the two ways of averaging. The Planck mean is an ordinary average and lands high, among the lines, because that is where most of the opacity is. The Rosseland mean is a harmonic average — it averages 1/κ rather than κ, because what carries the flux out of a star is transparency and transparencies add — and it lands 75 times lower, close to the continuum, because a harmonic mean is dominated by the smallest values in it. In other words the opacity that appears in the equation of radiative transport is a measurement of the gaps between the lines. Below: what that means for a table. The Planck mean rises as the first power of the line density, slope 0.76 as drawn — every line added is another contribution to an ordinary average. The Rosseland mean does almost nothing at first, slope 0.13, and then turns up sharply, slope 0.99, once the lines are close enough to blanket the windows. That is why adding several million atomic transitions to an opacity table in the early 1990s changed nothing for decades and then changed stellar structure: the new lines were not the first lines, they were the ones that finally closed the gaps.

A mean dominated by the gaps

The opacity in the equation of radiative transport is not an average of the opacity. It is a harmonic average weighted by the temperature derivative of the Planck function, which makes it a measurement of the transparent windows between the lines rather than of the lines — and that single fact decides what adding a million spectral lines to a table does.

starlight · Opacity The main sequence against its own ceiling. Luminosity against mass, in solar units, with the Eddington limit drawn on the same axes. The limit is exactly linear in mass. The main sequence is not — until it is: above about 55 solar masses the relation flattens to a straight proportionality as well, and the two lines then run parallel at a fixed ratio of 83 per cent. They never cross on this range, which is the more interesting outcome: the most massive stars do not approach their limit gradually, they are built at a fixed fraction of it. The Sun, by contrast, radiates 0.0026 per cent of its own ceiling.

A brightness that would blow the star apart

Light pushes. For a star of a given mass there is a luminosity at which the outward push on its own outer layers equals the inward pull of gravity, and it depends on nothing but the mass and the opacity — not on the star's structure, its composition or its age.

stars · Eddington limit Which gradient is steeper, and where. The temperature gradient radiation would need in order to carry the whole luminosity, against the gradient a rising blob of gas follows, through an n = 3 polytrope with a Kramers opacity and an energy generation going as ρT^4. Wherever the first exceeds the second the layer is unstable and convects. This model puts the outer boundary at 89 per cent of the radius; helioseismology measures the base of the Sun's convection zone at 71.3 per cent, and the gap between the two numbers is the price of a polytrope with one opacity law rather than a solar model. The curve also rises above the adiabatic value inside 4 per cent of the radius, which is a convective core — it appears when the energy generation is concentrated enough, and it is the structure a star burning by the CNO cycle has.

The part of a star that boils

Energy leaves a star's interior either by radiation or by bulk motion, and which one happens is settled by comparing two temperature gradients. The comparison decides that the Sun's outer third churns and its core does not, and that a massive star does exactly the opposite.

stars · Energy transport The period spans 492× and the pulsation constant 2.1×. The pulsation constant Q = P√(ρ̄/ρ̄_⊙) for 7 radial pulsators, against their periods. A radial pulsation is a standing sound wave across a star and the sound speed in a self-gravitating body is set by its own gravity, so the period should go as (Gρ̄)^−1/2 and Q should be the same for every star pulsating in the same mode. Across a factor of 492 in period — from 1.3 hours to 27 days — Q varies by 2.07, with a mean of 0.0393 days. That is the whole reason a period–luminosity relation can exist: a period is a density, a density is a mass and a radius, and a radius with a temperature is a luminosity. The masses here are the weak link and the caption should say so — for the Cepheids they come from evolutionary models rather than from a dynamical measurement, and the drift of Q with period is where that assumption is showing.

The valve that has to sit at the right depth

A period–luminosity relation only exists because a period is a density in disguise. Which stars have a period at all is decided by whether a partial ionisation zone sits deep enough to have mass behind it and shallow enough that convection has not taken it over.

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

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.

exoplanets · Exoplanet atmospheres A better measurement that made the model worse: 0.9 per cent in the sound speed. The fractional difference between the Sun's sound speed as its own oscillations measure it and as a structural model predicts it, against fractional radius. Zero would be agreement. The lower curve is the model built on the solar abundances used until the mid-2000s, and it hugs the axis: a part in a thousand across most of the interior, which was for a long time the best-tested piece of stellar physics anybody had. The upper curve is the same model with the abundances re-measured using three-dimensional atmospheres and without assuming local thermodynamic equilibrium — better measurements by every methodological standard, which lowered carbon, nitrogen and oxygen by around thirty per cent. The disagreement grows to 0.9 per cent, and it is not spread through the star: it peaks at 0.683 of the radius, just beneath the base of the convection zone at 0.713. The same substitution moves the model's own convection-zone base from 0.715 to 0.729, against a seismic value known to about a thousandth. What is being tested here is not really the abundances but what converts a composition into a structure, which is the opacity: the metals whose abundances fell are exactly the ones whose bound–free absorption dominates at those temperatures, and an opacity larger by some fifteen per cent near that boundary would restore the agreement. Laboratory measurements of iron at those conditions have since come in high by about that much, which is a satisfying result to have arrived at by way of a discrepancy in the sound speed of the Sun. The curves are published inversions and model differences rather than anything computed here; what the figure adds is where they peak and by how much.

A better measurement that made the model worse

The Sun's composition was re-measured with better atmospheres and better physics, and the carbon, nitrogen and oxygen abundances fell by about thirty per cent. The improved model then disagreed with the Sun's own oscillations by ten times as much as the model it replaced, and it still does.

starlight · Opacity A spot is dark because it is squeezed. Pressure against depth below the quiet photosphere's optical surface, drawn as a ratio to the pressure there. The rising curve is the surrounding gas, which grows exponentially with a 140-kilometre scale height because that is what hydrostatic equilibrium in an ideal gas produces. The flat pair of bands is the spot's own budget: a magnetic pressure of 3.58·10⁵ dyn/cm² from a 3000-gauss field, which is 25 per cent of the total, plus the gas pressure left over. Horizontal balance requires the two columns to reach the same total at the same geometric level, and the level at which they do is 348 kilometres below the quiet surface — so the spot's own optical surface sits in a hollow. Measured Wilson depressions, obtained from the foreshortening of a spot near the limb, are four to six hundred kilometres. Nothing about the darkness was assumed: a field strength read off a Zeeman splitting fixes the magnetic share, the share fixes the depression, and the depression fixes the temperature the deeper layer must have to carry the reduced flux.

The darkness a field pays for

A sunspot is not cool because something is missing. It is cool because a three-thousand-gauss field supplies part of the pressure that holds the column up, the gas therefore supplies less, and the level at which that gas becomes opaque sits some hundreds of kilometres deeper than the surface around it.

stars · Plasma beta A flux that measures a temperature to 0.12 per cent. Neutrino flux against central temperature, in units of the standard model's, for the three main solar channels. The exponents are not arbitrary: each reflects how far up the Gamow peak the reaction has to reach, so the channel with the largest Coulomb barrier is the steepest. The boron-8 flux goes as roughly the twenty-fourth power, which means a measurement good to 3 per cent constrains the Sun's central temperature to 0.12 per cent — better than any other technique by an order of magnitude. The same steepness is why the flux is useless as a check on anything else: a stellar model whose central temperature is uncertain at the half-per-cent level predicts this flux to within a factor, and the disagreement between two model families is far larger than the measurement.

A flux that is a thermometer to a tenth of a per cent

The boron-8 neutrino flux from the Sun's core rises as roughly the twenty-fourth power of the central temperature. That makes it the sharpest thermometer in astrophysics and simultaneously the most fragile prediction — a model uncertain in its central temperature by half a per cent predicts the flux to within a factor.

stars · Solar neutrinos An exponent between 2.3 and 4, not a number. The local logarithmic slope of the mass–luminosity relation against mass — the exponent that "L ∝ M^3.5" stands in for — read off the same fit to eclipsing binaries that the relation itself is drawn from. It is a step function with two corners, and the steps are 2.3 below 0.43 M☉, 4 0.43 to 2 M☉, 3.5 2 to 55 M☉. The horizontal lines are what homology predicts: a star in radiative equilibrium has L ∝ μ⁴M³/κ̄, which gives an exponent of 3 when the opacity is electron scattering and a constant, 5.5 when it is Kramers and depends on density and temperature, and 1 where radiation pressure holds the star up and the luminosity is pinned to the Eddington value. The measured steps sit between those predictions rather than on them, which is what a real star does — no star is homologous, the low-mass ones are convective throughout and the heavy ones have convective cores, and the fitted exponents are what is left after all of that. The useful statement is not that the exponent is 3.5 but that it is between 2 and 5 everywhere and 4 in the middle, and that the reason it moves is a change of which opacity is carrying the energy out.

An exponent that is a slope, not a law

The phrase “L goes as M to the three and a half” stands in for a curve with three straight pieces and two corners. The exponent is 2.3 below half a solar mass, 4 in the middle and falls towards 1 at the top — and each change of exponent is a change in which opacity is carrying the energy out.

stars · The mass–luminosity relation Where titanium oxide comes apart, and what the pressure does to it. The fraction of titanium oxide dissociated into Ti + O, against temperature, for a gas of equal parts at 10, 300, 6000 N/m² of total pressure. The equilibrium is Saha's with a free atom in place of the free electron and the pair's reduced mass in place of the electron's — the two atoms' partial pressures multiplied together, divided by the molecule's, equal an equilibrium constant Kₚ(T) — which for equal abundances gives a dissociated fraction of √(Kₚ/(Kₚ+P)) exactly. A molecule is one particle becoming two, so pressure suppresses the reaction exactly as it suppresses ionisation — the half-dissociation point moves from 3145 K at 10 N/m² to 4164 K at 6000. The bond energy is 6.87 eV against hydrogen's ionisation potential of 13.6, which is why this happens at a third of the temperature. Rotational and vibrational partition functions are included, classically and harmonically; the atoms keep their ground-term weights.

A spectrum with no continuum left in it

Below about four thousand kelvin a photosphere stops being a gas of atoms. Titanium oxide, water and carbon monoxide take over, the same equilibrium settles how much of each survives with a bond energy where an ionisation potential used to be — and the bands are so crowded that the level every line depth is measured against is nowhere on the plate.

starlight · Ionisation Hayashi lines: a luminosity that moves and a temperature that does not. The effective temperature a fully convective star must have, against its luminosity, for 0.5, 1, 2 solar masses, drawn as a Hertzsprung–Russell diagram with temperature falling to the right. The temperature follows from matching the star's convective adiabat to a photosphere whose opacity, from the negative hydrogen ion, rises as the ninth power of temperature. Over a factor of thirty thousand in luminosity each line moves by only 363 K, 398 K, 436 K respectively: the lines are nearly vertical, at 3650 K, 4000 K, 4383 K for one solar luminosity. A contracting young star, or a swelling giant, moves up or down its line almost without changing colour. To the right of the lines no star in hydrostatic equilibrium can sit.

A wall built by a steep opacity

A star that is convective all the way to its surface cannot choose its surface temperature. The interior fixes a relation between pressure and temperature, the photosphere fixes another through its opacity, and the negative hydrogen ion's opacity rises so steeply with temperature that the two meet at nearly the same temperature whatever the star's size. Young stars contracting and old stars swelling slide up and down a line in the diagram almost without changing colour — and nothing in equilibrium can sit to its right.

starlight · Opacity The edge of the Sun, computed. The brightness of the Sun at 500 nanometres across its limb, as a fraction of the brightness at the centre of the disc, against height above the level where the photosphere's radial optical depth is one, from the radiative transfer along each line of sight through a grey model atmosphere. Inside the disc the brightness falls gently with limb darkening; at the edge it falls to nothing over about two hundred kilometres. The steepest point of the fall — the inflection point, which is what a solar radius conventionally means — is 335 km above the radial τ = 1 level, because a ray grazing the limb passes through 270 times more atmosphere than a vertical one and becomes opaque that much higher up, at about 336 km. The brightness falls to one per cent 555 km above the τ = 1 level and to a thousandth at 694.

The Sun is three hundred kilometres larger seen edge-on

The radius of the Sun measured from images and eclipses and the radius measured from its oscillations disagree by about three hundred kilometres, and neither is wrong. A line of sight grazing the limb passes through hundreds of times more atmosphere than a vertical one and becomes opaque higher up — so the edge a telescope sees sits above the surface the waves feel, and every other definition of the edge, from the detector's threshold to the colour of the light, moves it again.

sky · Phases and eclipses Brighter by mass, fainter by temperature. Main sequences from 0.5 to 1.3 solar masses in luminosity against effective temperature, hotter to the left, for solar composition, [Fe/H] = −1.5, from homology applied to the solar-composition relation. The dots join stars of the same mass: lowering the metals moves each one up and to the left, so a 0.8 solar-mass star goes from 0.41 to 0.66 solar luminosities and from 5049 to 5879 K. The main sequence is steep, so a move of that size to the left carries a star past solar stars much brighter than itself: at a fixed temperature of 5800 K the metal-poor sequences lie 0.59 magnitudes below the solar one. A metal-poor star is brighter than a solar one of the same mass and fainter than a solar one of the same temperature — which is why such stars were called subdwarfs.

Brighter for its mass and fainter for its colour

Two stars of the same mass are almost the same star only if they are made of the same material. Take most of the metals out and the opacity falls, the radiation leaks out faster and the star settles brighter and much hotter — so far hotter that, at a given colour, it sits below the ordinary main sequence. The relation between mass and light has composition inside it, and every mass, age and distance read off it inherits that.

stars · The mass–luminosity relation The gradient radiation would need, as a young star shrinks. The radiative temperature gradient at the centre of a contracting, fully convective pre-main-sequence star — the gradient radiation alone would need to carry the luminosity — against the star's radius, which shrinks to the right, for 0.3, 0.5, 1, 2 solar masses, logarithmic. The centre stays convective while the gradient exceeds the adiabatic 0.4 (dashed). As the star shrinks its centre gets hotter and denser, Kramers opacity, ∝ ρT⁻³·⁵, falls, and the gradient falls as the radius to the power 2.5; where it crosses 0.4 a radiative core appears. 0.3 M☉: not before the main sequence; 0.5 M☉: not before the main sequence; 1 M☉: at 1.94 R☉; 2 M☉: at 7.62 R☉. Dots: each mass's main-sequence radius, where contraction stops. A Sun-like star leaves its convective line at 1.9 solar radii; the lightest never do.

Where a shrinking star stops boiling at its centre

A young star descends the Hertzsprung–Russell diagram at almost fixed temperature, convective from centre to surface, until something at its centre changes and it turns sharply left toward the main sequence. The something is the opacity. As the star shrinks its centre grows hotter and denser, the Kramers opacity there falls, and radiation becomes able to carry the luminosity without the steep gradient that drives convection. A Sun-like star crosses that line at about twice its present radius, two million years old; a star of a third of the Sun's mass never crosses it at all.

starlight · Opacity An envelope that boils away in a million years, or keeps for longer than the universe. The time a 5-Earth-mass planet takes to lose an envelope of 5.0% of its mass through an isothermal wind driven by its own heat and the star's light, against the escape parameter Λ at the envelope's outer edge — the ratio of the gravitational binding of a molecule there to its thermal energy — for equilibrium temperatures of 500, 800, 1200 K, on a logarithmic time axis. The rate is exponential in Λ, so the time grows by a factor of e for every unit: each step of ten in Λ is more than four decades of time. The envelope goes in under a million years below Λ ≈ 19, 18, 18 and lasts longer than a billion years above Λ ≈ 26, 26, 26. The window between a million and a billion years is a band only about seven units of Λ wide, so a planet is either boiling or effectively sealed.

An envelope that boils down before the starlight strips it

The usual story of how small planets lose their hydrogen starts with the star's ultraviolet light and runs for a hundred million years. It skips the first million. When the gas disc that fed a young planet disappears, an envelope that was swollen and held in by the disc's pressure is suddenly exposed, and it pours off as a wind powered by ordinary heat — until it has shrunk to the size at which the same number near thirty that governs every escaping atmosphere shuts the flow off.

exoplanets · Atmospheric escape

Named alongside it

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

ConvectionHelioseismologyKramers opacityOptical depthMean molecular weightMetallicityPhotosphereRosseland meanEddington limitElectron scatteringHayashi trackHomology

All concepts