Stars

Two stars only a Fourier transform can tell apart

A giant burning hydrogen in a shell and one burning helium in its core sit at the same luminosity, the same temperature and the same mean density. Every scaling relation returns the same mass and radius for both. The gravity-mode period spacing is fifty seconds for one and three hundred for the other.

Assumes Asteroseismology, Stellar evolution and Degeneracy.

The rung below this one took two numbers off a power spectrum — the spacing between consecutive radial modes and the frequency at which the oscillation power peaks — and inverted them for a mass and a radius. The inversion is exact, it is calibrated on the Sun, and it has been applied to some twenty thousand red giants.

It has one blind spot, and it is a large one. Those two numbers describe the envelope: the large separation is an inverse sound-crossing time and therefore a mean density, and the peak frequency is an acoustic cutoff and therefore a surface gravity. Neither says anything about what is happening at the centre.

That matters because two entirely different stars share an envelope. A star that has exhausted hydrogen in its core and is burning it in a shell, and a star that has ignited helium in a core that is now convective, can sit at the same luminosity, the same temperature, the same radius and the same mean density. Every measurement in this collection returns the same answer for both.

Two stars at the same point of every other diagram, 4.1 times apart in one. Above: the gravity-mode period spacing against the large frequency separation, for 90 shell-burning giants and 55 core-burning ones. These are the same stars in every other measurement. They have the same luminosity, the same temperature, the same colour, the same surface gravity and — inside the band drawn — the same Δν, which means the same mean density; the scaling relations of the rung below return the same mass and the same radius for both. What separates them is a quantity that comes from nowhere near the surface: ΔΠ₁ is set by the buoyancy frequency integrated across the core, and a core that has ignited helium is expanded and convective, so its integral is smaller and its period spacing larger. The two sequences do not touch — 61 seconds against 248 — and the gap sorts a catalogue of tens of thousands of giants into stars burning hydrogen in a shell and stars burning helium in a core, by a Fourier transform of a light curve. Below: what that looks like in the spectrum itself, drawn over two radial orders. The number of mixed ℓ = 1 modes between consecutive radial modes is Δν/(ΔΠ₁ν²), so the shell burner has about 12 of them and the core burner about 3: the star with the larger period spacing has the sparser spectrum. The picture cannot show what the same modes are also used for and cannot settle — the splitting of each mixed mode gives the rotation rate of the core separately from the envelope, and the cores come out spinning some ten times faster than the surface and a hundred times slower than any model of angular-momentum transport predicts.
Fig. 1 The measurement that separates them. Above: the gravity-mode period spacing against the large separation, for ninety shell-burning giants and fifty-five core-burning ones. Inside the band where both exist the two sequences do not touch — about eighty seconds against three hundred — and the gap sorts a catalogue of tens of thousands into two evolutionary states by a Fourier transform of a light curve. Below: what that looks like in the spectrum itself. The number of mixed modes between consecutive radial ones is Δν/(ΔΠ1ν2)\Delta\nu/(\Delta\Pi_1\nu^2), so the star with the larger period spacing has the sparser spectrum.

Two cavities in one star

A star supports two families of standing wave and they live in different places.

Pressure modes are sound: the restoring force is compression, the wave speed is the sound speed, and the modes are trapped between the surface and the depth at which the wave refracts back — deeper for lower angular degree. They are nearly equally spaced in frequency, which is why the spectrum looks like a comb.

Gravity modes are buoyancy waves: the restoring force is the density stratification, and the frequency scale is the buoyancy frequency NN. They propagate only where NN exceeds the wave frequency, which in a giant means only in the core, and they are nearly equally spaced in period rather than in frequency.

Between the two cavities is a region where neither wave propagates. In a main-sequence star that barrier is thick and the two families are entirely separate: the Sun’s gravity modes are confined to its interior, have never been convincingly detected, and would be the single most valuable measurement in the subject if they were.

In a red giant the barrier is thin. The core has contracted to something the size of the Earth while the envelope has swollen to the size of Mercury’s orbit, so the buoyancy frequency in the core is enormous — it reaches into the range where the envelope’s pressure modes live. The two cavities couple, and the resulting modes are neither one thing nor the other.

seven radial orders of a red giant, at 0.8623 μHz apart. The p-mode spectrum of a red giant — 1.1 solar masses in 30 solar radii — from the asymptotic relation with its second-order term, drawn as seven radial orders of ℓ = 0, 1 and 2 under a Gaussian envelope centred on ν_max = 4.43 μHz. Two numbers are marked and they do very different work. The large separation, 0.8623 μHz, is the spacing between consecutive ℓ = 0 modes and fixes the mean density. The small separation, 0.108 μHz, is 9.9 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 45.9 m s⁻¹ at the peak, 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.
Fig. 2 What that does to a spectrum. A giant’s oscillation power is concentrated at low frequency, the modes are of large amplitude and short lifetime, and the ℓ = 1 ridge — normally a single peak per radial order — is replaced by a cluster. Each member of that cluster is a mixed mode: a pressure wave in the envelope, coupled through the barrier to a gravity wave in the core, with a different share of its energy in each region.
The same comb folded at 0.8623 μHz — three ridges and their curvature. Frequency against frequency modulo Δν, for seven radial orders of a red giant. Folding at 0.8623 μHz stacks the orders into three near-vertical ridges, one for each degree, and that is what makes Δν a fact about the data rather than a fitted parameter: get it wrong and the ridges lean. What is left over is the curvature — the ℓ = 0 ridge wanders 7.7 per cent of Δν across the drawn range, which is the departure from the asymptotic relation and the part of the spectrum that knows about the star's outer layers. The small separation is here too and here it is visible: 0.108 μHz between the ℓ = 0 and ℓ = 2 ridges, 70 pixels on this axis against 70 — the same quantity that is three pixels wide in the unfolded spectrum. A constant offset of 1 μHz has been subtracted before folding so that no ridge wraps round the edge.
Fig. 3 The same data folded at the large separation, which is how the structure is actually read. Radial modes stack into a straight vertical ridge; quadrupole modes sit just beside them; and the dipole modes, instead of forming a third ridge, spread across the diagram. The spread is the coupling. A mode with most of its energy in the envelope stays near where a pure pressure mode would be; one with most of its energy in the core is displaced far from it, and the pattern of displacements is the thing that carries the core’s structure.

Why the spacing is a core measurement

The asymptotic period spacing of pure gravity modes is

ΔΠ1  =  2π22(Nrdr)1,\Delta\Pi_1 \;=\; \frac{2\pi^2}{\sqrt{2}}\left(\int \frac{N}{r}\,dr\right)^{-1},

with the integral taken across the region where gravity waves propagate — that is, across the core and nowhere else. Nothing about the envelope appears in it.

The buoyancy frequency is large where the density gradient is steep and zero wherever the gas is convective, because a convective region is by definition one in which a displaced parcel keeps going. That single fact is the whole diagnostic:

A shell-burning giant has a radiative, degenerate helium core. The gradient is steep, NN is large through the whole core, the integral is large, and the period spacing is small — 40 to 90 seconds.

A core-helium-burning giant has a convective core. Helium burning is violently temperature-sensitive, so the core convects, and throughout that convective region NN is zero and contributes nothing to the integral. The integral is much smaller, and the period spacing is 200 to 400 seconds.

The two states differ by a factor of three to five in a quantity that has nothing to do with the surface, and they are indistinguishable in everything that does.

A furnace with the thermostat taken out. Left, the two pressures in a red giant's helium core at 10⁶ g/cm³. The ideal-gas pressure rises with temperature, as the whole regulation of a star depends on it doing; the degenerate electron pressure is a horizontal line, because the exclusion principle does not know the temperature. At the ignition point near 10⁸ K the degenerate term is 5.1 times the gas term, so a temperature rise there raises the total pressure by almost nothing and the core does not expand. Right, what that costs. Both curves start at the same temperature with the same triple-alpha rate, whose logarithmic slope is 41 — the famous T⁴⁰, differentiated here rather than quoted. The regulated core settles; the degenerate one, with nowhere to put the extra energy, goes vertical. The flash reaches 10¹⁰ solar luminosities for a few seconds and not one photon of it is seen: every erg goes into lifting the degeneracy, and the star's visible response is to become fainter and settle on the horizontal branch.
Fig. 4 The event that separates the two states. In a star below about two solar masses the helium core becomes degenerate before it is hot enough to burn, and degenerate matter does not expand when heated — so when ignition comes it is a runaway, releasing in minutes a luminosity comparable with a whole galaxy’s, all of it absorbed by the overlying envelope, which is the same reason a core-collapse shock is invisible from outside. Nothing about the flash reaches the surface. The star’s outside changes almost imperceptibly, and the only evidence that it has happened is that the core is now convective, which is exactly what the period spacing measures.

What was actually measured

Four years of nearly uninterrupted photometry from the Kepler mission, on some sixteen thousand red giants, at a precision of tens of parts per million per measurement.

What is measured is a light curve. What is extracted is a power spectrum, and from it: the frequency of maximum power, the large separation, the individual mode frequencies where the signal-to-noise allows, and — for the ℓ = 1 modes — the spacing of the periods after the pressure-mode contribution has been removed.

That last step is the difficult one. The observed periods of mixed modes are not equally spaced, because each is shifted by its coupling to the envelope; the underlying gravity-mode spacing has to be recovered by fitting the coupling as well. The standard method stretches the period axis by a function of the mode’s inertia and looks for equal spacing in the stretched variable, which converts the problem into a search for a straight line.

The result, first published in 2011, was immediate and unambiguous: the giants fall into two clean sequences with almost nothing between them.

Around 55 to 65 seconds for the shell burners at the relevant Δν\Delta\nu, and around 250 to 300 seconds for the clump. A handful of stars sit in between and they are the ones caught in the brief transition, which is itself a measurement of how long that transition takes.

The plane the two relations make, and how badly it is conditioned in mass. The large separation against the frequency of maximum power, both logarithmic, for four stars whose masses, radii and temperatures are stated and whose frequencies are computed from them — the Sun at 3,090 and 135 μHz; a subgiant at 671 and 39.0 μHz; a red-clump star at 34 and 4.06 μHz; a red giant at 4.43 and 0.862 μHz. Over the plane are lines of constant radius, at 1, 3, 10, 30 solar radii, and of constant mass, at 0.8 and 2. The two families are not at right angles and that is the finding: the constant-mass lines run at a slope of about 0.77 and the constant-radius lines at 0.5, so a factor of 2.5 in mass moves a star only 0.11 decades across the plane. Mass enters the inversion at the quarter power of the observables and radius at the first, which is why an asteroseismic radius is good to a few per cent and an asteroseismic mass to something nearer ten. Every point here was put on the plane by the relations and then read back off it: solving the two equations for mass and radius returns the stated values to 7.8e-16.
Fig. 5 The relations the rung below this one uses, and what they cannot see. Every star on this diagram has a mass and a radius from two numbers, and the clump and the giant branch overlap in exactly the region where the ambiguity lies. The scaling relations are not wrong about these stars; they are silent about the one thing that distinguishes them, which is a good description of what a rung is for.

The result nobody wanted

The same modes carry a second measurement, and it produced a discrepancy that is still open.

Rotation splits every non-radial mode into a multiplet, and the size of the splitting is a weighted average of the rotation rate over the region the mode samples. For a mixed mode that weighting is the useful part: a mode with most of its energy in the core reports the core’s rotation, and one with most in the envelope reports the envelope’s. Both are present in the same spectrum.

The measurement is unambiguous. Red-giant cores rotate about ten times faster than their envelopes — which sounds like a lot until it is compared with what should happen. A core that contracts by a factor of ten in radius while conserving its angular momentum should spin up by a factor of a hundred; the envelope expands and slows. The predicted contrast is a factor of a thousand or more.

The observed cores are spinning a hundred times too slowly. Something is transporting angular momentum out of the core efficiently, throughout the star’s life, and no proposed mechanism — magnetic torques, internal gravity waves, hydrodynamic instabilities — produces the right amount at the right stage without producing the wrong amount at some other. It is one of the cleanest quantitative failures in stellar physics, and it exists because a Fourier transform of a light curve can weigh the rotation of a region the size of the Earth inside a star a hundred parsecs away.

What the sorting bought

The separation is not an end in itself, and three things followed from it within a few years.

A distance scale. Core-helium-burning stars all arrive with nearly the same core mass, so they all sit at nearly the same luminosity — which makes the clump a standard candle, and a bright one, visible across the Galaxy in the infrared. Its usefulness had always been limited by contamination: any sample of clump stars selected by colour and magnitude contains giant-branch stars at the same place, and they are not standard anything. Seismic sorting removes them one by one, and the resulting candle is one of the cleaner rungs available inside the Galaxy.

A mass loss measurement. A star loses mass on the giant branch, at rates no theory predicts from first principles. Comparing the seismic masses of clump stars in a cluster with those of branch stars in the same cluster measures how much was lost between the two states — the same quantity that the initial–final mass relation reaches from the other end, and the two now agree at about 0.2 solar masses for a solar-mass star.

And a galactic archaeology. A giant’s mass is its age, because the mass fixes how long it took to get there. Sorting giants by evolutionary state and weighing each one turns a photometric survey into a catalogue of ages across the disc — which is how the thick disc was shown to be uniformly old and the thin disc to have a wide age spread, a distinction that colours and abundances alone had never settled.

Where the model stops

The coupling is a model. Extracting ΔΠ1\Delta\Pi_1 requires a description of how strongly the two cavities communicate, and the standard treatment is asymptotic — valid when the barrier is thick and the modes are many. For the most evolved giants it is not, and the spacings quoted there carry a systematic nobody has fully characterised.

A convective core has a boundary that is not sharp. Overshooting mixes material beyond the formally convective region, and where the boundary is placed changes the integral that sets the period spacing. The measured spacings of clump stars are in fact used to constrain overshooting rather than predicted from it, which makes them an input to stellar models rather than a test of them in that particular respect.

And the sample is what the instrument saw. A four-year data set resolves period spacings down to a limit set by its own length, and the shortest-spacing stars — the most evolved shell burners — are exactly the ones whose modes are hardest to resolve. The clean separation drawn at the top of this essay is partly a statement about what four years of photometry can measure.

Why these stars ring at all

Nothing in the account so far says where the oscillations come from, and the answer decides which stars can be studied this way and with what instrument.

The modes are not driven by a valve of the kind that makes a Cepheid pulsate. They are damped — every one of them would die away if left alone — and they are kept going by being continuously excited by the turbulent convection beneath the surface. Convective cells rising and overturning are a broadband source of acoustic noise, and the star’s resonances pick out the frequencies they can support.

Three consequences follow, and each shapes the observation.

The amplitudes are tiny and the modes are many. A stochastic source excites everything within its bandwidth, so hundreds of modes are present at once, each at an amplitude of a few parts per million in brightness for a solar-type star and up to hundreds of parts per million for a giant. That is a photometric precision no ground-based observation achieves through the atmosphere.

The frequencies are smeared. A damped mode driven by noise has a finite linewidth, set by its damping rate — so a mode is not a spike but a Lorentzian a fraction of a microhertz wide, and its measured frequency has an uncertainty that no length of observation reduces below that.

And the amplitude scales with the star. Larger, more luminous stars have more vigorous convection and slower oscillations, so a giant’s modes are hours long and comparatively easy to detect while a dwarf’s are minutes long and hard. That is why the giants were the first population measured in bulk, and why the technique’s reach across the Hertzsprung–Russell diagram is set by the excitation rather than by the physics being measured.

The star’s inclination, from the amplitudes

There is one more quantity the spectrum gives up, and it is geometric rather than structural.

A mode of angular degree \ell splits into 2+12\ell+1 components, and how much power appears in each depends on the angle between the star’s rotation axis and the line of sight. A star seen pole-on shows only the components whose patterns are visible from the pole; one seen equator-on shows a different distribution; and the ratio of the component heights is a function of the inclination alone.

So fitting the relative amplitudes within a split multiplet measures the inclination of the star’s rotation axis — a quantity that is otherwise almost unobtainable, since a spectrum gives only the product of the equatorial speed with the sine of that angle.

The use is immediate for planetary systems. A transiting planet’s orbit is nearly edge-on by construction, so comparing the orbit’s inclination with the star’s tells whether the system is aligned — and misalignment is evidence about how the planet arrived. The seismic route works where the usual spectroscopic method fails, on slow rotators and on systems with long-period planets that transit rarely.

A pattern of amplitudes in a power spectrum therefore answers a question about the geometry of a planetary system, which is a fair example of how much a resolved set of modes contains once the modes are identified.

The measurement is not free of assumptions — it requires the modes to be excited equally in all the components of a multiplet, which is expected and is not directly checked — but it needs no spectroscopy and no transit, and it applies to any star whose modes are resolved.

There is a further use of the same geometry that has nothing to do with planets. For a star in a cluster, the distribution of measured inclinations across many members says whether the cluster’s stars formed with their spins aligned or randomly oriented — a question about how the cloud that made them fragmented, answered by a statistic over power spectra rather than by any image. The samples large enough to attempt it exist for two clusters, and the answers so far disagree with each other, which is where such measurements usually begin.

The obstacle is not the method but the sample: measuring an inclination this way needs modes resolved well enough to separate the components of a multiplet, which needs a long time series on a bright star, and a cluster’s members are neither bright nor few.

So the measurement exists, works, and has been applied to a handful of individual stars where it settled a question no other technique could reach.

Both applications also share a limitation that is worth stating once: everything here depends on identifying which mode is which, and a spectrum whose modes cannot be labelled is a list of frequencies rather than a measurement.

One more comb covers the second of the two stars the essay is about.

nine radial orders of a red-clump star, at 4.057 μHz apart. The p-mode spectrum of a red-clump star — 1.2 solar masses in 11 solar radii — from the asymptotic relation with its second-order term, drawn as nine radial orders of ℓ = 0, 1 and 2 under a Gaussian envelope centred on ν_max = 34 μHz. Two numbers are marked and they do very different work. The large separation, 4.057 μHz, is the spacing between consecutive ℓ = 0 modes and fixes the mean density. The small separation, 0.467 μHz, is 7.2 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 8.87 m s⁻¹ at the peak, 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.
Fig. 6 The oscillation comb of a red-clump star. It looks like the red-giant comb in every respect a photometric measurement can see, and the period spacing of its mixed modes — which is not visible here at all — is what separates them.

The generalisation

The structure of the argument is one that recurs whenever a system has two regions with very different natural frequencies and a barrier between them.

Two pendulums coupled by a weak spring have normal modes that are neither pendulum’s; where their natural frequencies are far apart the modes are almost pure, and where they cross the modes exchange character over a narrow range. That exchange — an avoided crossing — is exactly what a giant’s dipole modes do as the star evolves and the core’s buoyancy frequency sweeps upwards through the envelope’s acoustic band. Each individual mode is pushed aside as a gravity mode approaches it, hands over its character, and settles back.

The consequence for measurement is the useful part: a mode that is mostly one thing carries a little information about the other, and a spectrum containing many such modes contains the coupling itself. That is why the observable here is not a frequency but a pattern of displacements from where frequencies would otherwise be — the same reason an échelle diagram is drawn at all, and the same reason a small divisor in secular theory is the interesting quantity rather than a nuisance.

The comb and the échelle are the same data twice, and each is worth drawing for a star at a different evolutionary stage, since the whole method is a claim that the pattern moves in a predictable way.

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.
Fig. 7 Fourteen radial orders of the Sun at a large frequency spacing of 135 microhertz. The comb is nearly uniform and the departures from uniformity are the sound-speed structure — the Sun is the calibrating case because its interior is known independently from helioseismology.
The same comb folded at 38.97 μHz — three ridges and their curvature. Frequency against frequency modulo Δν, for twelve radial orders of a subgiant. Folding at 38.97 μHz stacks the orders into three near-vertical ridges, one for each degree, and that is what makes Δν a fact about the data rather than a fitted parameter: get it wrong and the ridges lean. What is left over is the curvature — the ℓ = 0 ridge wanders 8.0 per cent of Δν across the drawn range, which is the departure from the asymptotic relation and the part of the spectrum that knows about the star's outer layers. The small separation is here too and here it is visible: 2.92 μHz between the ℓ = 0 and ℓ = 2 ridges, 45 pixels on this axis against 45 — the same quantity that is three pixels wide in the unfolded spectrum. A constant offset of 4 μHz has been subtracted before folding so that no ridge wraps round the edge.
Fig. 8 And the échelle diagram for a subgiant, folded at thirty-nine microhertz. The ridges bend where a mixed mode couples the envelope to the core, which happens as soon as the star leaves the main sequence — so the diagnostic the essay is about first appears here rather than on the giant branch.

Where this ladder goes next

Later rungs on this anchor: the suppressed dipole modes, a substantial minority of giants whose ℓ = 1 amplitudes are far below the rest and which are thought to have strong magnetic fields trapping the waves in their cores; the direct measurement of those fields from the asymmetry of the splittings, achieved in 2022; the inversion of the rotation profile rather than its two-zone summary; period spacings in subgiants, where the cavities are only beginning to couple and the evolution is fast enough to see; and the use of clump stars as standard candles, which is the practical reason a catalogue wants them identified in the first place.

About the same objects

Not linked from either essay — found by the objects both name.

What links here

Essays that link to this one from their own argument.

The objects this essay names

Each one links to every other essay that touches it.

Angular momentum transportAsteroseismologyBuoyancy frequencyCore-helium burningEvanescent zoneGravity-modeMixed modePeriod spacingPressure modeRed-clumpRed giant branchRotational splitting