Lecture 7: After the Main Sequence — Low-Mass Evolution and White Dwarfs
When hydrogen runs out, why does a star swell before it shrinks?
Learning Objectives
After completing this reading, you should be able to:
- Infer from star-cluster HR diagrams why low-mass stars evolve onto the red giant branch, horizontal branch, and white dwarf cooling sequence.
- Derive from the virial theorem why an inert core that loses energy contracts and heats up.
- Explain causally why shell hydrogen burning drives red giant expansion instead of simple fading.
- Calculate a red giant radius from the Stefan-Boltzmann law using explicit dimensionless ratios.
- Compare the energy released per unit mass in hydrogen burning and helium burning, and use that comparison together with luminosity arguments to explain why helium-burning phases are shorter.
- Explain, using the temperature sensitivity of the triple-alpha reaction and the weak temperature dependence of degeneracy pressure, why helium ignition in a degenerate core runs away.
- Use the non-relativistic degeneracy-pressure scaling and the white-dwarf mass-radius relation to infer why more massive white dwarfs are smaller.
Concept Throughline
The HR diagram tells us what stars do after the main sequence. This reading explains why. When core hydrogen runs out, the core contracts and heats, shell burning turns on, the envelope expands into a red giant, helium ignites, the star sheds its envelope, and a white dwarf remains. Every stage is an inference from observables tied to a physical model.
At each stage, the key reasoning task is: what new interior structure or energy source is required to explain the observed position of the star on the HR diagram?
Track A (Core, ~35 min): Read the opening observable section and then Parts 1–6 in order. Skip boxes marked Enrichment.
Track B (Full, ~45 min): Read everything, including the enrichment boxes connecting helium flashes, Type Ia supernovae, and the white-dwarf endpoint.
Both tracks cover all core learning objectives.
We will use four different mathematical relationships, and they do not mean the same thing:
- \(=\) means an exact statement within the model being used.
- \(\approx\) means a numerical approximation after values are inserted.
- \(\sim\) means an order-of-magnitude estimate.
- \(\propto\) means proportionality or scaling, with constants suppressed.
When we move from one of these to another, we will say so explicitly.
This short Hubble visualization is a good way to begin because it reminds us that the HR diagram starts from measurements of brightness and color. Globular clusters are especially useful because the stars lie at nearly the same distance and many of them are old enough that post-main-sequence populations are easy to identify.
Constructing the Hertzsprung-Russell Diagram for Globular Star Cluster Omega Centauri. A Hubble Space Telescope color image of the core of the globular star cluster Omega Centauri is used to construct a Hertzsprung-Russell diagram of the stellar populations in the cluster. When stars are sorted by brightness and color they can be used to create a graph that astronomers use to trace stellar evolution.
Credit: NASA, ESA, and J. Anderson, R. van der Marel, G. Bacon, and M. Estacion (STScI)
Observable: What Do We See When Stars Leave the Main Sequence?
Observable
Gaia’s Hertzsprung-Russell diagram. Millions of stars trace the main sequence, giant branch, and white dwarf sequence. What to notice: a coeval stellar population will lose its highest-mass main-sequence stars first, so the turnoff point acts like a clock. (Credit: ESA/Gaia)
Evolutionary tracks off the main sequence for stars of different masses. A solar-mass star ascends the Red Giant Branch (RGB), undergoes the helium flash, moves to the Horizontal Branch (HB), climbs the Asymptotic Giant Branch (AGB), and ends as a white dwarf. More massive stars (5, 10 M☉) take wider loops through the supergiant region. The labeled phases (RGB, HB, AGB) correspond to the physical stages described in this reading.
Astronomers do not watch a single Sun-like star for billions of years. Instead, we infer stellar evolution from populations. In old stellar systems we observe:
- a main-sequence turnoff, where the most massive stars that can still burn core hydrogen define an age scale,
- a heavily populated red giant branch (RGB), where stars are cool but luminous,
- a horizontal branch (HB), where stars are less luminous than the RGB tip but hotter,
- planetary nebulae, glowing shells of gas around hot compact central stars,
- and a white dwarf cooling sequence, where hot but faint remnants cool to lower luminosities over time.
These are not random points scattered across the HR diagram. They are organized structures.
Because stars in a cluster are observed at nearly the same distance and formed at roughly the same time, differences in their HR-diagram positions can be interpreted primarily as differences in stellar mass and evolutionary stage.
Model
The model idea is simple before it becomes detailed: a star changes position on the HR diagram when its luminosity, surface temperature, and radius change. Those surface properties change because the interior energy source, pressure support, and energy-transport mechanism change.
So the problem is not merely to label post-main-sequence populations. The problem is to explain what interior physics sends a low-mass star from:
\[ \text{main sequence} \rightarrow \text{subgiant} \rightarrow \text{RGB} \rightarrow \text{HB} \rightarrow \text{AGB} \rightarrow \text{planetary nebula} \rightarrow \text{white dwarf}. \]
In older, metal-poor populations this core-helium-burning phase often appears as a horizontal branch; in more metal-rich populations it may appear instead as a red clump. In both cases, the key inference is that a new central energy source has turned on.
Inference
The observations imply that low-mass stars do not simply fade when hydrogen runs out in the core. They undergo a sequence of structural reorganizations. The driving question for this reading is:
Why do stars follow this path instead of simply fading?
Part 1: Core Hydrogen Exhaustion Starts the Evolutionary Track
Observable
In old star clusters, we do not see stars vanish from the main sequence and disappear. We see them peel away onto the subgiant branch and then ascend the red giant branch. That means the end of core hydrogen burning does not shut the star off immediately. It launches a new phase.
Model
On the main sequence, a low-mass star is in:
- hydrostatic equilibrium: pressure gradients balance gravity,
- thermal equilibrium: luminosity leaving the surface is replenished by nuclear burning in the core.
Once core hydrogen is exhausted, the core can no longer replenish the energy it loses. The correct starting point is the stellar virial theorem:
\[ 2 K_\text{th} + U_\text{grav} = 0 \tag{1}\]
Virial theorem (stellar)
What it predicts
The relationship between total thermal energy \(K_\text{th}\) and gravitational potential energy \(U_\text{grav}\) in a star.
What it depends on
\(K_\text{th} = -U_\text{grav}/2\). Total energy \(E = U_\text{grav}/2 < 0\) (bound).
What it's saying
In hydrostatic equilibrium, thermal energy and gravitational potential energy are linked by the virial theorem. During contraction, a decrease in gravitational potential energy is split between increasing the gas thermal energy and energy that can ultimately be radiated away. Losing total energy makes a self-gravitating core hotter.
Assumptions
- Star is in hydrostatic equilibrium
- Ideal gas (non-relativistic particles)
- Time-averaged or in quasi-static equilibrium
See: the equation
For this section, use the virial theorem as a quasi-static, ideal-gas argument for the contracting core. We are not claiming that every detail of the post-main-sequence interior can be read off from this one equation alone.
Let us define the symbols:
- \(K_{\text{th}}\) is the thermal kinetic energy of the gas, in
erg. - \(U_{\text{grav}}\) is the gravitational potential energy, in
erg, and it is negative.
The total energy of the core is
\[ E = K_{\text{th}} + U_{\text{grav}}. \]
Now use the virial theorem step by step:
\[ 2K_{\text{th}} + U_{\text{grav}} = 0 \]
\[ U_{\text{grav}} = -2K_{\text{th}} \]
Substitute into the total-energy expression:
\[ E = K_{\text{th}} + U_{\text{grav}} = K_{\text{th}} - 2K_{\text{th}} = -K_{\text{th}}. \]
So we also have
\[ E = -K_{\text{th}} = \frac{U_{\text{grav}}}{2}. \]
This line is the key to the whole section. For a self-gravitating system, losing total energy makes the bound state deeper.
If the core radiates energy away, then
\[ \Delta E < 0. \]
Because \(E = U_{\text{grav}}/2\), making \(E\) more negative means
\[ \Delta U_{\text{grav}} < 0. \]
More negative gravitational potential energy means the potential well deepens, which means the core must contract.
Because \(E = -K_{\text{th}}\), making \(E\) more negative means
\[ \Delta K_{\text{th}} > 0. \]
So contraction raises the thermal kinetic energy. For an ideal gas,
\[ K_{\text{th}} = \frac{3}{2}N k_B T, \]
where
- \(N\) is the number of particles,
- \(k_B\) is Boltzmann’s constant in
erg K^{-1}, - \(T\) is temperature in
K.
If \(N\) is fixed for the contracting core, then higher \(K_{\text{th}}\) means higher \(T\).
This is the negative-heat-capacity behavior of self-gravitating objects:
- the core loses energy,
- the core contracts,
- the core gets hotter.
The core response is therefore the opposite of what students expect from ordinary objects like stovetops or coffee mugs.
After central hydrogen is exhausted, the core is mostly helium ash left behind by earlier fusion. Because the temperature is still too low for helium fusion, this helium core is inert: it contributes mass and gravity, but not new nuclear power.
At the outer edge of that contracting core, a hydrogen-rich layer that was previously too cool for fusion gets compressed and heated. The shell that ignites is the layer just outside the inert helium core because it is the first hydrogen-rich region compressed strongly by the contracting core. The center itself cannot resume hydrogen burning because its hydrogen fuel has already been exhausted.
Inference
The turnoff and subgiant observations imply that core hydrogen exhaustion does not remove gravity from the problem. It lets gravity regain control. The virial theorem then predicts the next stage: contraction of the inert helium core and heating of the shell around it.
A stellar core loses energy, so \(\Delta E < 0\). Using
\[ E = -K_{\text{th}} = \frac{U_{\text{grav}}}{2}, \]
decide whether the core temperature rises or falls. Write your answer as a chain of algebraic statements and a one-sentence physical interpretation.
Start with
\[ \Delta E < 0. \]
Because
\[ E = \frac{U_{\text{grav}}}{2}, \]
we have
\[ \Delta U_{\text{grav}} = 2\,\Delta E < 0. \]
So the gravitational potential energy becomes more negative, which means the core contracts.
Because
\[ E = -K_{\text{th}}, \]
we also have
\[ \Delta K_{\text{th}} = -\Delta E > 0. \]
Finally, since
\[ K_{\text{th}} = \frac{3}{2}N k_B T, \]
with \(N\) fixed, larger \(K_{\text{th}}\) means larger \(T\).
Inference: a self-gravitating core that loses energy contracts and heats up. That is negative heat capacity in action.
Part 2: Shell Burning Drives the Red Giant Branch
Observable
Hydrogen shell burning on the Red Giant Branch. When core hydrogen is exhausted, a shell of hydrogen burns around the inert helium core. The helium ash accumulates, growing the core mass. Meanwhile, the enormous energy output from the shell drives the outer envelope to expand and cool — the star becomes a red giant. The expanded, cool, hydrogen-rich outer envelope is hundreds of times larger than the tiny burning region at its heart.
Evolutionary tracks off the main sequence for stars of different masses. A solar-mass star ascends the Red Giant Branch (RGB), undergoes the helium flash, moves to the Horizontal Branch (HB), climbs the Asymptotic Giant Branch (AGB), and ends as a white dwarf. More massive stars (5, 10 M☉) take wider loops through the supergiant region. The labeled phases (RGB, HB, AGB) correspond to the physical stages described in this reading.
On the HR diagram, the star does not move down and to the right as a fading ember. It climbs the red giant branch: luminosity rises strongly while surface temperature falls to roughly \(3{,}500\text{–}4{,}500~\text{K}\). Observationally, that means the star becomes both larger and brighter even though core hydrogen burning has ended.
From
\[ L = 4\pi R^2 \sigma T_{\text{eff}}^4, \]
we can write the schematic radius scaling
\[ R \propto L^{1/2} T_{\text{eff}}^{-2}. \]
So if a star becomes much more luminous while its surface temperature drops, its radius must increase dramatically. That is the defining signature of a red giant.
Model
The model has two linked pieces: why shell burning is powerful, and why that power inflates the envelope.
Why can shell hydrogen burning exceed the old core luminosity?
The luminosity from any burning region is
\[ L = \int \varepsilon \, dm, \]
where
- \(\varepsilon\) is the energy generation rate per unit mass, in
erg g^{-1} s^{-1}, - \(dm\) is the mass element, in
g.
For pedagogical purposes, we use the rough local scaling
\[ \varepsilon_{\text{pp}} \propto \rho X^2 T^4, \]
where
- \(\rho\) is density,
- \(X\) is the hydrogen mass fraction,
- \(T\) is temperature.
The exact dependence varies with conditions, but the main point is that compression and heating of the shell make the burning rate rise sharply.
When the helium core contracts, the hydrogen-rich layer just outside it is compressed. Both \(\rho\) and \(T\) in that shell rise. Because the reaction rate depends strongly on temperature, even a modest rise in \(T\) increases \(\varepsilon_{\text{pp}}\) substantially. The shell is geometrically thin, but it is located exactly where the temperature and density have become large enough for intense burning. The result is a large luminosity emerging from a small region.
Why does the envelope expand instead of just glowing harder at the same size?
Here is the causal chain, step by step:
- Core contraction raises the pressure at the core-envelope boundary.
- That compresses and heats the hydrogen shell.
- The shell luminosity rises.
- That enhanced outward energy flux must be transported through the base of the envelope.
- The envelope opacity is high, so purely radiative transport struggles to carry the new luminosity.
- The radiative temperature gradient steepens.
- Convection turns on when the gradient exceeds the stability threshold.
- The envelope expands until the star finds a new equilibrium with a larger radius, lower surface temperature, and high luminosity.
So the correct physical story is not “convection makes giants big” in isolation. It is:
\[ \text{core contraction} \rightarrow \text{strong shell burning} \rightarrow \text{energy injection at envelope base} \rightarrow \text{high-opacity transport bottleneck} \rightarrow \text{convection} \rightarrow \text{envelope expansion}. \]
That is why red giant structure is an envelope response to an interior burning-shell problem.
As the envelope expands, the radiating surface becomes much larger. A larger surface area allows the star to remain very luminous while cooling at the photosphere, moving it up and to the right onto the red giant branch.
Radius example from the Stefan-Boltzmann law
We can now quantify what “giant” means. Use the Stefan-Boltzmann law:
\[ L = 4\pi R^2 \sigma T^4 \tag{2}\]
Stefan-Boltzmann law
What it predicts
Given \(R\) and \(T\), it predicts the luminosity \(L\).
What it depends on
Scales as \(L \propto R^2 T^4\).
What it's saying
Luminosity depends on surface area (\(R^2\)) and temperature (\(T^4\)). Double the temperature, get 16× the luminosity.
Assumptions
- Blackbody radiation
- Spherical, uniformly radiating surface
- Effective surface temperature
See: the equation
Take the ratio of the giant to the Sun:
\[ \frac{L}{L_\odot} = \frac{4\pi R^2 \sigma T_{\text{eff}}^4}{4\pi R_\odot^2 \sigma T_{\text{eff},\odot}^4} = \left(\frac{R}{R_\odot}\right)^2 \left(\frac{T_{\text{eff}}}{T_{\text{eff},\odot}}\right)^4. \]
The luminosity ratio \(L/L_\odot\) is dimensionless, so the right-hand side must also be dimensionless, as required for \(\left(R/R_\odot\right)^2\).
Now solve for the radius ratio:
\[ \left(\frac{R}{R_\odot}\right)^2 = \frac{L/L_\odot}{\left(T_{\text{eff}}/T_{\text{eff},\odot}\right)^4} \]
\[ \frac{R}{R_\odot} = \left[ \frac{L/L_\odot}{\left(T_{\text{eff}}/T_{\text{eff},\odot}\right)^4} \right]^{1/2}. \]
For a red giant with
\[ L = 1{,}000\,L_\odot, \qquad T_{\text{eff}} = 3{,}500~\text{K}, \qquad T_{\text{eff},\odot} = 5{,}800~\text{K}, \]
we have
\[ \frac{T_{\text{eff}}}{T_{\text{eff},\odot}} = \frac{3{,}500~\text{K}}{5{,}800~\text{K}} = 0.603. \]
Notice that the kelvin units cancel in the ratio:
\[ \frac{\text{K}}{\text{K}} = 1. \]
Now raise that ratio to the fourth power:
\[ \left(\frac{T_{\text{eff}}}{T_{\text{eff},\odot}}\right)^4 = (0.603)^4 \approx 0.133. \]
Substitute into the radius equation:
\[ \left(\frac{R}{R_\odot}\right)^2 = \frac{1{,}000}{0.133} \approx 7.52 \times 10^3. \]
Take the square root:
\[ \frac{R}{R_\odot} = \sqrt{7.52 \times 10^3} \approx 86.7. \]
Since \(86^2 = 7396\) and \(87^2 = 7569\), the square root is indeed about \(86.7\).
So the star’s radius is
\[ \boxed{\frac{R}{R_\odot} \approx 87}. \]
Inference
An RGB star that is luminous but cool must have a very large radius. The red giant branch therefore implies not simple dimming, but large-scale envelope expansion driven by shell burning and transport physics.
Use the Stefan-Boltzmann law to determine whether a star with
\[ L = 300\,L_\odot \qquad \text{and} \qquad T_{\text{eff}} = 4{,}000~\text{K} \]
is likely to be main-sequence-sized or giant-sized. Show the radius-ratio calculation and interpret the result physically.
Use
\[ \frac{R}{R_\odot} = \left[ \frac{L/L_\odot}{\left(T_{\text{eff}}/T_{\text{eff},\odot}\right)^4} \right]^{1/2}. \]
First compute the temperature ratio:
\[ \frac{T_{\text{eff}}}{T_{\text{eff},\odot}} = \frac{4{,}000~\text{K}}{5{,}800~\text{K}} = 0.690. \]
Then
\[ (0.690)^4 \approx 0.227. \]
So
\[ \left(\frac{R}{R_\odot}\right)^2 = \frac{300}{0.227} \approx 1.32 \times 10^3 \]
and therefore
\[ \frac{R}{R_\odot} \approx \sqrt{1.32 \times 10^3} \approx 36. \]
So
\[ \boxed{\frac{R}{R_\odot} \approx 36}. \]
Inference: a star with this luminosity and temperature is much too large to be an ordinary main-sequence star. It is a giant.
Part 3: Helium Ignition, the Triple-Alpha Process, and the Helium Flash
Observable
In old clusters, the red giant branch does not continue forever. Stars reach the RGB tip, then populate the horizontal branch or, in some populations, the red clump. That tells us a new central energy source must turn on. The question is what that source is, and why it behaves so differently from hydrogen burning.
Model
Why helium ignition requires much higher temperatures
Helium nuclei have charge \(Z = 2\), so the Coulomb barrier is much larger than it was for proton-proton burning. The helium core therefore has to contract to much higher temperature:
\[ T \sim 10^8~\text{K} \]
before helium burning can proceed efficiently.
But there is another problem. There are no stable nuclei at mass number \(A=5\) or \(A=8\). So nature cannot climb smoothly upward from helium by single-particle additions. That forces the triple-alpha pathway. The equation below is the net reaction:
\[ 3\,{^4\text{He}} \rightarrow {^{12}\text{C}} + \gamma \tag{3}\]
Triple-alpha process
What it predicts
Net triple-alpha reaction: three helium-4 nuclei fuse to form one carbon-12 nucleus plus a gamma-ray photon.
What it depends on
Requires \(T \gtrsim 10^8\) K and \(\rho \gtrsim 10^4\) g/cm³. Rate scales as \(\varepsilon_{3\alpha} \propto T^{40}\).
What it's saying
The only way to bridge the mass-5 and mass-8 gaps in nuclear stability. The physical pathway proceeds through unstable \(^8\)Be and then alpha capture into an excited state of \(^{12}\)C (the Hoyle state at 7.65 MeV) — without which carbon would be exceedingly rare and carbon-based life impossible.
Assumptions
- Occurs in red giant cores and horizontal branch stars
- Two-step process: \(2\alpha \rightarrow\ ^8\)Be (unstable) + \(\alpha \rightarrow\ ^{12}\)C
- Competes with alpha capture: \(^{12}\)C + \(\alpha \rightarrow\ ^{16}\)O
See: the equation
The triple-alpha process — how stars build carbon from helium. Two ⁴He nuclei fuse to form unstable ⁸Be, which normally decays in ~10⁻¹⁶ s. But at the extreme densities and temperatures of helium burning (~10⁸ K), a third ⁴He captures onto ⁸Be before it decays, producing stable ¹²C and releasing gamma rays. Fred Hoyle predicted the nuclear resonance state that makes this reaction possible — from the mere existence of carbon in the universe.
The reaction proceeds in two steps:
\[ {^4\text{He}} + {^4\text{He}} \rightleftharpoons {^8\text{Be}} \]
followed by
\[ {^8\text{Be}} + {^4\text{He}} \rightarrow {^{12}\text{C}}^* \rightarrow {^{12}\text{C}} + \gamma. \]
The Hoyle state in \(^{12}\text{C}\) makes this chain fast enough to matter astrophysically.
Energy per unit mass: hydrogen burning versus helium burning
To understand why helium-burning phases are so short, compare the energy released per unit mass of fuel.
For hydrogen burning, four protons eventually become one helium nucleus and release about
\[ E_{\text{H}} \approx 26.7~\text{MeV} \]
from about
\[ 4~\text{u} \]
of fuel. So the energy per unit mass is
\[ \frac{E_{\text{H}}}{m_{\text{H fuel}}} = \frac{26.7~\text{MeV}}{4~\text{u}} = 6.675~\text{MeV/u}. \]
For helium burning, three helium nuclei produce one carbon nucleus and release about
\[ E_{\text{He}} \approx 7.3~\text{MeV} \]
from
\[ 12~\text{u} \]
of helium fuel. So
\[ \frac{E_{\text{He}}}{m_{\text{He fuel}}} = \frac{7.3~\text{MeV}}{12~\text{u}} = 0.608~\text{MeV/u}. \]
Now take the ratio:
\[ \frac{\left(E_{\text{He}}/m_{\text{He fuel}}\right)} {\left(E_{\text{H}}/m_{\text{H fuel}}\right)} = \frac{0.608~\text{MeV/u}}{6.675~\text{MeV/u}} \approx 0.091. \]
So helium burning releases only about
\[ \boxed{9\%} \]
as much energy per unit mass of fuel as hydrogen burning.
That is one important reason helium-burning phases are shorter than hydrogen-burning phases. The total lifetime also depends on how much fuel is available and on the luminosity of the star during that phase.
A rough lifetime estimate follows
\[ t \sim \frac{\text{available nuclear energy}}{L}. \]
So even if a phase has some fuel available, its duration is reduced if the star is very luminous during that stage.
Why the triple-alpha process is unstable in a degenerate core
The triple-alpha energy generation rate has an extremely steep temperature dependence. For pedagogical purposes, we use the schematic local scaling
\[ \varepsilon_{3\alpha} \propto T^{40}. \]
Now perturb the temperature slightly. Let
\[ T' = T(1+\delta), \]
where \(\delta\) is a small fractional perturbation.
Then the new energy generation rate satisfies
\[ \frac{\varepsilon'_{3\alpha}}{\varepsilon_{3\alpha}} = \left(\frac{T'}{T}\right)^{40} = (1+\delta)^{40}. \]
For small \(\delta\), the binomial approximation gives
\[ (1+\delta)^{40} \approx 1 + 40\delta. \]
So a \(1\%\) temperature rise, with \(\delta = 0.01\), gives
\[ \frac{\varepsilon'_{3\alpha}}{\varepsilon_{3\alpha}} \approx 1 + 40(0.01) = 1.4. \]
The exact value is
\[ (1.01)^{40} \approx 1.49, \]
so a \(1\%\) temperature rise increases the energy generation rate by roughly \(40\text{–}50\%\).
That is already an extremely sensitive response in an ordinary gas.
Now add degeneracy. In a non-relativistic degenerate electron gas,
\[ P_\text{deg} \sim \frac{\hbar^2}{m_e}\,n_e^{5/3} \tag{4}\]
Electron degeneracy pressure (non-relativistic)
What it predicts
Given electron number density \(n_e\), it predicts the pressure from quantum confinement (Pauli exclusion).
What it depends on
Scales as \(P \propto n_e^{5/3}\) (non-relativistic). No temperature dependence.
What it's saying
Compressing fermions forces them into higher momentum states (Pauli exclusion + Heisenberg uncertainty). The resulting pressure exists even at \(T = 0\) — fundamentally different from thermal pressure. This is what holds up white dwarfs.
Assumptions
- Non-relativistic electrons (\(p_F \ll m_e c\))
- Fully degenerate gas (\(T \ll T_F\))
- Relativistic regime: \(P \propto n_e^{4/3}\) (softer — leads to Chandrasekhar limit)
See: the equation
The important physical point is that
\[ P_{\text{deg}} \propto \rho^{5/3} \]
with no explicit temperature dependence. This is a schematic non-relativistic degeneracy scaling, not a full equation of state with all constants restored.
That means:
- in an ordinary thermal gas, heating raises pressure and the core expands,
- in a degenerate core, the pressure depends mainly on density rather than temperature, so heating does not trigger the usual rapid pressure increase,
- so the core does not expand enough to cool itself,
- so the temperature keeps rising,
- and because \(\varepsilon_{3\alpha} \propto T^{40}\), the burning runs away.
That is the helium flash.
Why the helium flash is not seen directly
The helium flash is violent in the core but not at the surface. The released energy is absorbed by the surrounding stellar interior and goes into:
- lifting the electron degeneracy,
- expanding the core slightly,
- restoring an ordinary thermal thermostat.
The released energy is absorbed deep in the stellar interior and goes mainly into lifting degeneracy and re-expanding the core, rather than emerging promptly from the surface as a bright optical outburst.
So the flash is a deep interior restructuring event, not a giant optical outburst at the photosphere.
Inference
The horizontal branch implies a new core energy source: helium burning. The helium flash implies that helium ignites in a degenerate core, where pressure is largely temperature-independent and the normal stellar thermostat is broken.
If the core temperature in a helium-burning region rises by \(2\%\), estimate the change in \(\varepsilon_{3\alpha}\) using
\[ \varepsilon_{3\alpha} \propto T^{40}. \]
Then explain why the same perturbation is much more dangerous in a degenerate core than in an ordinary thermal gas.
Let
\[ T' = T(1+0.02) = 1.02T. \]
Then
\[ \frac{\varepsilon'_{3\alpha}}{\varepsilon_{3\alpha}} = \left(\frac{T'}{T}\right)^{40} = (1.02)^{40}. \]
Using the small-perturbation approximation,
\[ (1.02)^{40} \approx 1 + 40(0.02) = 1.8. \]
The exact value is
\[ (1.02)^{40} \approx 2.21. \]
So a \(2\%\) temperature rise roughly doubles the energy generation rate.
In an ordinary gas, that stronger burning raises pressure, the core expands, and the temperature drops back. In a degenerate core, \(P_{\text{deg}}\) depends mainly on density, not temperature, so expansion does not happen quickly enough. The extra burning therefore drives still more heating: runaway.
Part 4: The AGB and Planetary Nebula Phase
Observable
Internal structure of an AGB star (close-up of the core region for a ~1 M☉ star). The inert carbon-oxygen core sits at the center, surrounded by a helium-burning shell, a helium layer, and an outer hydrogen-burning shell. The entire core structure fits inside a region smaller than the Earth, while the stellar envelope extends to ~1–1.5 AU. Not to scale.
The Ring Nebula as a planetary-nebula example. This is not a star-forming nebula. It is an envelope ejected by a low-mass star near the end of its life. The glowing ring is gas photoionized by the hot exposed remnant core at the center. Observations like this are the reason we infer that low-mass stars end by envelope loss, not by core collapse. (Credit: NASA, ESA, CSA, and STScI)
After the horizontal branch, low-mass stars return to a cool, luminous giant phase: the asymptotic giant branch (AGB). We also observe planetary nebulae, which are glowing shells of gas around hot compact central stars. That means the late giant phase must somehow end by ejecting the envelope and exposing the core.
Model
The AGB resembles the RGB in that the star again has a compact inert core, surrounding shell burning, and an enormous convective envelope. The difference is that two burning shells are now present instead of one.
Once core helium is exhausted, the star has an inert carbon-oxygen core. The evolution is again driven by core contraction:
\[ \text{C/O core contracts} \rightarrow \text{He shell burning turns on} \rightarrow \text{H shell burning continues farther out} \rightarrow \text{envelope re-expands}. \]
Now there are two burning shells instead of one:
- an inner helium-burning shell,
- an outer hydrogen-burning shell.
This shell structure is thermally unstable. Because helium burning occurs in a thin shell with strong temperature sensitivity, the helium shell can become thermally unstable and ignite in pulses, causing:
- episodic luminosity spikes in the shell,
- convective dredge-up that mixes carbon-rich material upward,
- strong mass loss from the outer envelope.
As mass loss strips away the envelope, the hot compact core is exposed. That bare core emits ultraviolet photons that ionize the expelled gas. The planetary nebula therefore glows because it is photoionized by the hot remnant core.
Because the ejected gas expands and its density drops rapidly, the nebula is only a brief phase compared with the much longer-lived white dwarf that remains.
Despite the name, a planetary nebula has nothing to do with planets; the term is historical.
So the central-star luminosity source and the nebular glow must be distinguished carefully:
- the central star is a hot exposed post-AGB core,
- the nebula glows because that core ionizes the ejected gas.
Inference
Planetary nebulae imply that low-mass stars do not end by core collapse. They end by envelope ejection. The observed nebula plus hot central star together point to an AGB star that lost its envelope and exposed a compact remnant core.
Explain why a planetary nebula is short-lived even though the white dwarf remnant persists for billions of years. Your answer must identify the observable, the physical mechanism making the gas glow, and the reason the nebula disappears long before the white dwarf does.
The observable is a glowing shell of ejected gas around a hot compact central star.
The gas glows because ultraviolet photons from the exposed post-AGB core ionize the ejected material. Recombination and line emission then produce the visible nebula.
The nebula is short-lived because the ejected gas expands and thins out. As it spreads, its density drops and the surface brightness falls rapidly. After roughly
\[ \sim 10^4~\text{yr}, \]
the shell has dispersed enough that it is no longer bright as a planetary nebula.
The white dwarf remnant persists much longer because it is a compact degenerate object with a large reservoir of thermal energy and a very long cooling time.
Inference: the nebula is a transient envelope phenomenon, not the long-term remnant itself.
Part 5: White Dwarfs Are Degenerate Stellar Remnants
Observable
White dwarf mass-radius relation. More massive white dwarfs are smaller — the counter-intuitive result of R ∝ M⁻¹/³ from degeneracy pressure. The full relativistic curve (solid) plunges to R = 0 at the Chandrasekhar limit (1.44 M☉), where electron degeneracy can no longer support the star. Beyond this mass, gravity wins. (Credit: ASTR 201 (generated))
In HR diagrams, white dwarfs are hot but faint. That means they must have small radii. Spectra and binary masses show that they can contain roughly half a solar mass or more inside an Earth-sized volume. Observationally, that is an extraordinary combination: high mass, tiny radius, and no active fusion.
Model
A white dwarf is supported by electron degeneracy pressure, not by ordinary thermal gas pressure.
Degeneracy pressure comes from quantum state packing, not from thermal agitation.
The key scaling for a non-relativistic degenerate electron gas is
\[ P_{\text{deg}} \propto \rho^{5/3}. \]
Assumptions:
- the electrons are degenerate,
- they are still non-relativistic,
- the pressure is dominated by the electron gas.
Because this pressure comes from quantum state filling, it is largely independent of temperature. That is why a white dwarf can cool without losing its pressure support.
Electron degeneracy pressure supports the white dwarf mechanically against gravity, but it is not a continuing energy source. The white dwarf shines because it is still hot and slowly cooling.
Why more massive white dwarfs are smaller
We can combine the degeneracy-pressure scaling with a characteristic gravitational pressure scaling.
For a star of mass \(M\) and radius \(R\),
\[ \rho \sim \frac{M}{R^3}. \]
Insert that into the degeneracy-pressure scaling:
\[ P_{\text{deg}} \propto \rho^{5/3} \propto \left(\frac{M}{R^3}\right)^{5/3} = \frac{M^{5/3}}{R^5}. \]
A characteristic gravitational pressure scale is
\[ P_{\text{grav}} \sim \frac{GM^2}{R^4}. \]
Here \(P_{\text{grav}} \sim GM^2/R^4\) is a characteristic self-gravitational pressure scale, not an exact local pressure formula.
Equilibrium requires these to scale together:
\[ \frac{M^{5/3}}{R^5} \propto \frac{GM^2}{R^4}. \]
Multiply both sides by \(R^5\):
\[ M^{5/3} \propto G M^2 R. \]
Now divide both sides by \(G M^2\):
\[ R \propto \frac{M^{5/3}}{M^2} = M^{-1/3}. \]
So the white-dwarf mass-radius relation is
\[ R_\text{WD} \propto M^{-1/3} \tag{5}\]
White dwarf mass-radius relation
What it predicts
Given a white dwarf's mass, it predicts the radius (inversely).
What it depends on
Scales as \(R \propto M^{-1/3}\). More massive white dwarfs are smaller.
What it's saying
The opposite of ordinary objects: adding mass shrinks the star because stronger gravity compresses the degenerate electrons harder. This trend holds for non-relativistic white dwarfs and breaks down near the Chandrasekhar limit.
Assumptions
- Non-relativistic regime (well below Chandrasekhar limit)
- Zero-temperature (degenerate) equation of state
- Breaks down as \(M \rightarrow M_\text{Ch}\) (relativistic corrections)
See: the equation
The algebra matches the physical expectation: adding mass strengthens gravity, so equilibrium requires a higher density. In a degenerate object, achieving higher pressure means packing electrons into a smaller volume, so the radius decreases.
More massive white dwarfs are smaller because stronger gravity demands higher density, and higher density means larger degeneracy pressure.
Density example with units
Take a white dwarf with
\[ M = 0.6\,M_\odot \qquad \text{and} \qquad R = R_\oplus = 6.4 \times 10^8~\text{cm}. \]
Use
\[ \bar{\rho} = \frac{M}{\frac{4}{3}\pi R^3}. \]
Convert the mass. Using \(M_\odot \approx 1.99 \times 10^{33}~\text{g}\), which we round here to \(2.0 \times 10^{33}~\text{g}\) for simple arithmetic,
\[ M = 0.6 \left(2.0 \times 10^{33}~\text{g}\right) = 1.2 \times 10^{33}~\text{g}. \]
Now compute the volume:
\[ V = \frac{4}{3}\pi R^3 = \frac{4}{3}\pi \left(6.4 \times 10^8~\text{cm}\right)^3 \]
\[ V \approx 1.10 \times 10^{27}~\text{cm}^3. \]
Now divide:
\[ \bar{\rho} = \frac{1.2 \times 10^{33}~\text{g}}{1.10 \times 10^{27}~\text{cm}^3} \approx 1.09 \times 10^6~\text{g}\,\text{cm}^{-3}. \]
So the average density is
\[ \boxed{\bar{\rho} \approx 1.1 \times 10^6~\text{g}\,\text{cm}^{-3}}. \]
The units are correct because
\[ \frac{\text{g}}{\text{cm}^3} = \text{g}\,\text{cm}^{-3}. \]
Inference
An Earth-sized, hot, faint stellar remnant with no fusion source points to a white dwarf supported by electron degeneracy pressure. The inverse mass-radius relation is a signature of quantum pressure, not ordinary gas support.
A white dwarf has more mass than Earth but a radius of order Earth’s radius. Use the scaling
\[ R_{\text{WD}} \propto M^{-1/3} \]
to explain why increasing the mass makes the remnant smaller rather than larger. Your answer must mention both gravity and degeneracy pressure.
Increasing the white dwarf mass strengthens gravity, so equilibrium requires a larger inward pressure scale.
In a degenerate electron gas, the way to get more pressure is to increase the density. Higher density means electrons are confined to a smaller volume, which increases their momentum spread and therefore the pressure.
That is why the star responds by becoming smaller and denser rather than larger. The scaling
\[ R_{\text{WD}} \propto M^{-1/3} \]
captures that balance mathematically.
Inference: an inverse mass-radius relation is evidence for quantum support, not ordinary thermal support.
If a carbon-oxygen white dwarf in a binary system accretes enough mass to approach the Chandrasekhar limit, the density and temperature in its interior rise until carbon ignites. Because the star is degenerate, the pressure does not respond normally to the temperature increase. The thermostat is broken. Carbon burning runs away, and the white dwarf is disrupted in a Type Ia supernova.
This is the same basic logic as the helium flash, but at much higher density and with much more destructive consequences.
Part 6: Putting the Full Low-Mass Path Together
Observable
The combined HR-diagram evidence from clusters and field stars shows a path:
- main-sequence turnoff,
- subgiant branch,
- red giant branch,
- horizontal branch or red clump,
- AGB stars,
- planetary nebulae,
- white dwarf cooling sequence.
These are the observable signposts of low-mass stellar evolution.
Model
The physical chain is now complete:
- Core hydrogen is exhausted.
- The helium core contracts and heats by the virial theorem.
- Hydrogen shell burning ignites around the inert core.
- The envelope expands because shell luminosity and transport bottlenecks reorganize the outer star.
- Helium ignites when the core reaches \(T \sim 10^8~\text{K}\).
- A helium flash occurs if the core is degenerate.
- Core helium burning creates carbon, and some carbon captures helium to form oxygen.
- The carbon-oxygen core contracts after helium is exhausted.
- Double-shell burning on the AGB drives instability and mass loss.
- The envelope is ejected as a planetary nebula.
- A white dwarf remains, supported by electron degeneracy pressure.
The overall timing is also part of the model:
| Phase | Typical duration for a Sun-like star | Dominant energy source |
|---|---|---|
| Main sequence | \(\sim 10~\text{Gyr}\) | Core H burning |
| RGB ascent | \(\sim 10^8\text{–}10^9~\text{yr}\) | H shell burning |
| Core He burning | \(\sim 10^8~\text{yr}\) | Core He burning + H shell |
| AGB | \(\sim 10^7~\text{yr}\) | He shell + H shell |
| Planetary nebula | \(\sim 10^4~\text{yr}\) | Ionized ejecta around hot core |
| White dwarf cooling | \(\gtrsim 10^{10}~\text{yr}\) | No fusion; cooling only |
These timescales reflect the basic rule
\[ t \sim \frac{\text{available energy}}{L}, \]
so long phases have large fuel reservoirs and modest luminosities, while short phases are either high-luminosity burning stages or brief transient transitions.
Inference
Low-mass stars do not end by catastrophic collapse. They end by structural reorganization and envelope loss. The observed HR-diagram path therefore implies a sequence of new burning configurations followed by a degenerate remnant.
When we observe:
- a cool luminous red giant,
- a horizontal-branch helium-burning star,
- a planetary nebula with a hot compact central star,
- or an Earth-sized white dwarf cooling with no fusion,
the physical model lets us infer that these are not unrelated classes of objects, but different stages of the same low-mass evolutionary pathway, all triggered by the exhaustion of core hydrogen.
Reference Tables
Key Equations and What They Mean Here
| Equation | Role in this reading | Core idea |
|---|---|---|
| \(2K_{\text{th}} + U_{\text{grav}} = 0\) | Core response after H exhaustion | Energy loss leads to contraction and heating |
| \(L = 4\pi R^2 \sigma T_{\text{eff}}^4\) | Giant radius inference | Luminous + cool implies large radius |
| \(3\,{^4\text{He}} \rightarrow {^{12}\text{C}} + \gamma\) | Core helium burning | Carbon is built from helium |
| \(\varepsilon_{3\alpha} \propto T^{40}\) | Helium flash stability | Small temperature changes strongly amplify burning |
| \(P_{\text{deg}} \propto \rho^{5/3}\) | White dwarf support | Degeneracy pressure is density-controlled, not temperature-controlled |
| \(R_{\text{WD}} \propto M^{-1/3}\) | White dwarf structure | More massive white dwarfs are smaller |
Symbol Legend
| Symbol | Meaning | Units |
|---|---|---|
| \(K_{\text{th}}\) | Thermal kinetic energy | erg |
| \(U_{\text{grav}}\) | Gravitational potential energy | erg |
| \(E\) | Total energy | erg |
| \(\varepsilon\) | Energy generation rate per unit mass | erg g^{-1} s^{-1} |
| \(T_{\text{eff}}\) | Effective surface temperature | K |
| \(\rho\) | Density | g cm^{-3} |
| \(P_{\text{deg}}\) | Degeneracy pressure | dyn cm^{-2} |
| RGB | Red giant branch | — |
| HB | Horizontal branch | — |
| AGB | Asymptotic giant branch | — |
This reading ends with electron degeneracy holding the line against gravity. But the white-dwarf support law has a limit. In Reading 8, we ask what happens when the electrons become relativistic and the non-relativistic \(P_{\text{deg}} \propto \rho^{5/3}\) scaling is no longer enough. That leads to the Chandrasekhar mass.