Astronomy & Celestial Mechanics
the sky, computed
Concept

Convective envelope — where it appears

The outer region of a cool star, where energy is carried by overturning gas rather than by radiation. It is what generates a magnetic field, what dissipates a tide, and what disappears above about 1.3 solar masses.

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

A lever 12 radii long, and the spin it removes. A magnetised stellar wind, drawn with the Alfvén surface at 12 stellar radii — the schematic distance at which the wind's inertia finally beats the field. Inside it the gas is forced to turn with the star, so every gram that leaves carries the specific angular momentum of the radius at which it broke free rather than of the surface it came from, and the lever squares: J̇ = (2/3) Ṁ Ω r_A². Beyond the surface the streamlines curve backwards, because angular momentum conservation makes the azimuthal speed fall as 1/r while the radial speed does not. With a moment of inertia coefficient of 0.073 and a mass loss of 2.3·10⁻¹⁴ solar masses a year, the star loses a fraction 2.3·10⁻¹⁴ of its mass and a fraction 3·10⁻¹¹ of its angular momentum in the same year — a ratio of 1,315, which is (2/3)(r_A/R)²/k² and nothing else. The e-folding time for the spin is 3.3·10¹⁰ years against 4.3·10¹³ years for the mass. Nothing here is to scale in one respect that matters: the wind's density falls by more than ten orders of magnitude across the drawn region, so the streamlines are drawn as though the flow were visible when almost none of it is.

A wind that takes no mass and all the spin

The Sun loses about a ten-thousandth of itself to its own wind over its whole life. It loses most of its rotation to the same wind, and the whole of that asymmetry is one geometric factor — the gas is forced to keep turning with the star until it is a dozen or more radii out from the surface it left.

stars · Magnetic braking A period, a colour, and an age. Rotation period against colour for stars of 125, 625, 1000, 2500, 4570 million years, under the empirical relation P = t^0.5189 × 0.7725(B−V − 0.4)^0.601. The isochrones do not cross and are separated at every colour by exactly the age ratio raised to 0.5189, which is what allows a single measured period to be inverted for an age once the colour is known. The Sun, at B−V = 0.653 and 4570 million years, is placed by the relation at 26.8 days against the 25.4 days it is observed to have. Three clusters are marked at a common colour to show the spacing directly. The dashed boundary at the left is where the Rossby number — the period divided by the convective turnover time — passes 2 on the oldest isochrone, at B−V = 1.35: past that point the braking weakens and the relation is known to over-predict the age, which is the one place a rotation period stops being a clock. What the figure cannot show is the scatter, which is a few days at fixed colour and age and is the real error bar on any single star.

The clock that starts by forgetting

An ordinary star tells nothing about its age. It sits on the main sequence for billions of years at almost fixed brightness and colour, and the one property that changes monotonically is how fast it turns — but only because the braking law destroys the initial condition first, and only until it stops.

stars · Magnetic braking A period below which every orbit is round. Orbital eccentricity against period for binaries in four clusters of 0.125, 0.625, 6, 4 billion years, with the eccentricities drawn from one seeded distribution and then damped by exp(−age/τ), where τ rises as the 5.333 power of the period. Each cluster shows the same thing: below a boundary period nothing survives eccentric, above it the original distribution is untouched, and there is almost nothing in between because the timescale is so steep. The boundary is a clock. It moves as the three-sixteenths power of the age, which the figure checks against the drawn curves, and the calibration puts it at 6.5 days at 125 million years, 8.8 at 625 million and 12.5 at four billion — against measured cut-offs near 7.2, 8.5 and 12.5 days in the Pleiades, the Hyades and M67. The boundaries are also read back off the plotted points rather than trusted, and required to move outward with age. This is the cleanest measurement of tidal dissipation in ordinary stars that exists, and its cleanliness comes from the ages: a cluster's age is read off its main-sequence turn-off and owes nothing whatever to the tide being measured.

A cut-off period that is an age

Plot eccentricity against orbital period for the binary stars of one cluster and the picture has a wall in it. Below a certain period every orbit is circular; above it, the original spread survives untouched. The wall moves outward as the cluster ages, and where it stands is a measurement of how stars dissipate a tide.

orbits · Circularisation The Sun's internal rotation: differential above 0.693 R, rigid below, and the 15 nHz between two modes. Angular velocity against fractional radius, in nanohertz, from the inversion of hundreds of thousands of measured rotational splittings. Two features had to be discovered rather than deduced. The convection zone rotates differentially in latitude — 473 nanohertz at the equator, 330 at sixty degrees — and that latitude dependence persists all the way down through it, in surfaces that are very nearly radial rather than the cylinders a rotating convective fluid was expected to produce. And below about 0.693 of the radius the latitude dependence stops: the radiative interior turns as a single rigid body at 430 nanohertz, which is a remarkable thing for a fluid with no strength to do, and requires something — most likely a weak internal magnetic field — to be enforcing it. The two regimes are joined by a shear layer a few per cent of the radius thick, the tachocline, and it is where the solar magnetic cycle is generally thought to be generated, because it is the only place in the Sun with the shear a dynamo of that strength needs. The horizontal marks are what individual modes would report: the rotation averaged over the cavity each one occupies, weighted by the time the wave spends at each radius. A mode of degree 100 is trapped near the surface and reports 473 nanohertz; one of degree 1 passes through the deep interior and reports 458. The 15-nanohertz difference between them is the measurement, and splittings are measured to about a nanohertz — which is why the profile can be resolved at all. The kernels used here have the right support and the right sense; a real inversion uses computed eigenfunctions, and its resolution below 0.2 R is poor for the reason the previous figure gives.

A shear layer that should have spread

The Sun's convection zone turns differentially — its equator laps its poles about once every three months — and the radiative interior below turns as one rigid piece. Between them is a transition four per cent of the radius thick. Nothing in hydrodynamics keeps a velocity discontinuity that thin for four and a half billion years.

stars · Internal rotation The same orbit, realigned by one star and not by the other. The time an equilibrium tide takes to bring a planet's orbit into the plane of its star's equator, against orbital separation in stellar radii, for a planet of 1e-3 stellar masses. Both axes are logarithmic. The two curves differ only in how efficiently the star dissipates the tide, by a factor of 10⁴ — the contrast between a star with a convective envelope, where turbulence turns the tidal flow into heat, and one hotter than about 6,250 K, which has almost none. The lower curve is calibrated so that a Jupiter at 8 stellar radii realigns a cool star's orbit in 1 billion years, which is what the aligned systems require; the tidal quality factor of a star is not known from first principles and this is the honest way to say so. Everything else follows from the sixth power of the separation, which is measured off the drawn curve as 6.000. The two curves cross a Hubble time at 12.4 and 2.7 stellar radii, and their ratio is the sixth root of the dissipation contrast. Hot Jupiters sit between those two numbers. So the same arrival distribution of orbital tilts is erased around cool stars and preserved around hot ones, and a survey that finds cool hosts aligned and hot hosts scattered has measured the filter rather than the arrivals.

A misalignment only cool stars forget

A third of hot Jupiters orbit at a large angle to their star's equator, and some go round backwards. Sort the same planets by the temperature of their host and the picture changes — below about 6,250 kelvin almost all are aligned, and above it almost none are. The boundary is not about the planets.

exoplanets · Spin–orbit alignment A hundredfold expansion is a four-order slowdown. Equatorial rotation speed against radius for three stars leaving the main sequence at 2, 10, 100 km/s, on the single assumption that nothing exerts a torque. Both axes are logarithmic. The specific angular momentum is held fixed along each track and checked at four radii rather than asserted, so the speed falls as the reciprocal of the radius and then faster, as the moment of inertia coefficient slides from 0.073 on the main sequence to 0.02 in a centrally condensed giant envelope. A star at 100 km/s crosses the 8 km/s line — the usual boundary for calling a giant a rapid rotator — at 45.6 solar radii, and everything larger is slower. That is why the observed giants are almost all under two kilometres a second, and why the one or two per cent that are not cannot be explained by anything the star did on its own: the angular momentum has to have been delivered, by a swallowed companion or a merger. The figure assumes no mass loss, which for the largest radii drawn is the weakest of its assumptions.

A surface that slowed because the star grew

A star leaving the main sequence expands by a factor of a hundred. Conserve its angular momentum and its equator slows by four orders of magnitude, which is why every red giant is a slow rotator. One or two per cent are not, and a giant turning at fifty kilometres a second has to have been given the angular momentum by something else.

stars · Stellar evolution The same equation, run away at q = 0.8 and settled at q = 0.5. The donor's overfill of its own Roche lobe against time, integrated for 3 mass ratios at a donor adiabatic response of −0.33 — a star with a deep convective envelope, which expands as it loses mass. The overfill sets the transfer rate and the transfer rate changes the overfill, and the sign of that feedback is the stability criterion: where the lobe shrinks faster than the star does, ζ_L > ζ_ad, the overfill grows and the growth is exponential. At q = 0.8 the lobe responds at 0.03 and the overfill runs away; At q = 1.4 the lobe responds at 1.32 and the overfill runs away. At q = 0.5 it responds at -0.62 and the transfer throttles itself back. The critical ratio for this donor is 0.63 — below equal masses, which is the result the whole subject turns on: a giant transferring to a lighter companion is already unstable before the mass ratio has reversed, and what follows is not accretion but a common envelope. Nothing here is a fitted rate. The vertical scale is logarithmic and the runaway is a straight line on it, which is what an exponential is.

Three clocks and a runaway

Whether mass transfer between two stars is stable is a comparison of two logarithmic derivatives. How fast it runs is a separate question with three possible answers fourteen orders of magnitude apart — and the answer decides whether the companion accretes, is buried, or is swallowed.

stars · Mass transfer A factor of 2.9, and the density that would remove it. The predicted ⁷Li abundance against the baryon density, with the halo-star measurement as a horizontal band and the microwave background's density as a vertical one. At η₁₀ = 6.13 the network gives 4.70e-10 and the stars show 1.60e-10, a factor of 2.94 — 0.47 in the logarithm, which is the number every proposed resolution has to produce. The obvious cure is drawn as the second marker: ⁷Li rises as η², so the density that would reproduce the observation is η₁₀ = 3.58, a third below the measured one. That is where the cure fails, and it fails on a different element. Deuterium falls as η^(−1.6), so at 3.58 the predicted D/H is 5.946e-5 against the 2.527e-5 measured in quasar absorbers — 2.35 times too much, and 114 times deuterium's own error bar. Helium, meanwhile, moves to 0.2445 and stays inside its measurement, because its curve is flat. The three light elements do not fail together, and that is what rules out a single wrong parameter: whatever is wrong is wrong about mass seven specifically.

A factor of three, and the flatness that prices every cure

The oldest stars in the Galaxy show a third of the lithium the first three minutes should have left. Three kinds of resolution have been offered, and the datum that rules on all three is not how much lithium is missing but how uniformly it is missing.

cosmology · Nucleosynthesis A bound on the baryon density from an abundance nobody can extrapolate. Primordial deuterium against the baryon density, with the upper bounds that the solar system's own deuterium and helium-3 place on it. The inequality is D_p ≤ D_obs + ³He_obs/g₃, and it holds for any star-formation history whatever: a deuteron entering a star becomes ³He, so the pair can only be moved from one member to the other and then destroyed, never increased. With pre-solar values of D/H = 2.0e-5 and ³He/H = 1.5e-5, the bounds are 3.50e-5 at g₃ = 1, 5.00e-5 at g₃ = 0.5, 8.00e-5 at g₃ = 0.25 — and because deuterium falls with density, each upper bound on the abundance is a lower bound on η: η₁₀ > 4.98, η₁₀ > 3.99, η₁₀ > 2.97. The weakest of them, at a survival fraction of 0.25, still excludes 52 per cent of the density range below the answer. That was the state of the measurement for most of the 1980s, and the shape of it is the thing worth carrying: an abundance too poorly understood to be extrapolated at all still constrained the quantity, because its direction of change under processing was known even though its magnitude was not. The microwave background's 6.13 sits above every bound drawn, which is not a coincidence and is not evidence — a bound that excluded the answer would have been an error, and a bound that admits it is only a bound.

An abundance with no direction to correct in

Every primordial abundance is measured today and extrapolated backwards, and the extrapolation works because processing moves each species one way. Helium-3 is made by small stars and destroyed by large ones, so the sign of its correction is not merely uncertain — it is unknown.

cosmology · Nucleosynthesis Two damping times, one crossing, and the slope that separates them. Circularisation timescale against orbital period for the two tidal mechanisms, both logarithmic, normalised to 1.217 Gyr at 10 days. On logarithmic axes a power law is a straight line and the index is its slope, so the figure's content is that the two lines have different slopes and one crossing. The equilibrium tide gives 5.33, the bulge raised on a convective envelope being dragged ahead by a viscosity that is turbulent convection itself; the dynamical tide gives 7, gravity waves launched at the convective boundary carrying angular momentum to wherever they break. The horizontal lines are the ages of populations a boundary can be read in: where each curve crosses one is the wall that population shows — at 0.125 Gyr, 6.5 days against 7.2; at 0.625 Gyr, 8.8 days against 9.1; at 4 Gyr, 12.5 days against 11.9; at 10 Gyr, 14.8 days against 13.5. One cluster measures one number and both theories have a free normalisation, so one cluster cannot choose. The separation across every age available is a factor of 1.11. The discriminator that does not depend on the normalisation is mass: the equilibrium tide needs a convective envelope, and above about 1.3 solar masses there is not one, so the two predict different behaviour on either side of a boundary the theory names in advance. That is a measurement about where the wall stops behaving, not about where it is — and it is the reason the samples had to grow from tens of binaries per cluster to hundreds.

Two damping times, one crossing, and the slope that separates them

The equilibrium tide gives a damping time going as the sixteen-thirds power of the period and the dynamical tide as the seventh. One cluster measures one number and both theories have a free normalisation, so one cluster cannot choose.

orbits · Circularisation

Named alongside it

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

Angular momentumEquilibrium tideMagnetic brakingMain sequence turn-offTidal quality factorAngular momentum transportAsteroseismologyBaryon-to-photon ratioBig bang nucleosynthesisCluster ageDeuterium abundanceDynamical tide

All concepts