Stellar evolution: main sequence to red giant, HR diagram tracks

Intuition Beginner

Stars are not eternal. They are born, they live, and they die. A star spends most of its life on the main sequence, fusing hydrogen into helium in its core. The Sun has done this for 4.6 billion years and will continue for about 5 billion more. When the core hydrogen runs out, the star swells into a red giant. The Sun will grow so vast that it engulfs Mercury and Venus, and possibly Earth.

Massive stars, those weighing more than about eight times the Sun, go further. After the red giant phase they fuse heavier and heavier elements in concentric layers — carbon, neon, oxygen, silicon — until iron forms at the centre. Iron cannot be fused for energy. With no fuel left, the core collapses in a fraction of a second and the outer layers rebound outward in a supernova, an explosion that can briefly outshine an entire galaxy.

The Hertzsprung-Russell diagram plots every star by its luminosity against its surface temperature. As a star ages, it traces a path across this diagram, recording each life stage: main sequence, red giant, and beyond. By reading the tracks of star clusters — groups of stars born together — astronomers measure how old the cluster is from where its stars turn off the main sequence.

Visual Beginner

A star's life traces a path across the Hertzsprung-Russell diagram. The diagram below labels the main phases a solar-mass star passes through, from birth on the main sequence through the red giant and asymptotic giant branches to its final white-dwarf cooling track.

Phase	Core state	Surface change
Main sequence	Hydrogen fuses to helium	Stable, slowly brightening
Subgiant	Core hydrogen exhausted; shell burning begins	Envelope expands, cools
Red giant branch	Degenerate helium core; hydrogen shell	Huge, red, luminous
Helium flash	Helium ignites degenerately	Brief; little surface change
Horizontal branch	Helium fuses to carbon	Smaller, hotter, stable
Asymptotic giant branch	Helium and hydrogen shells	Large, pulsating, heavy mass loss
White dwarf	No fusion; radiates stored heat	Small, hot, fading

Massive stars skip the degenerate helium flash and the horizontal branch. Instead they burn successively heavier fuels until iron forms, then collapse. Their tracks sweep across the top of the diagram to the blue supergiant region before ending in a supernova.

Worked example Beginner

A star's lifetime depends on how much fuel it has and how fast it burns it. Fuel scales with mass. The burning rate scales with luminosity, and luminosity rises steeply with mass — roughly as mass to the 3.5 power for main-sequence stars. So lifetime scales as mass divided by luminosity, giving mass to the power minus 2.5. A star twice the Sun's mass lives about one-sixth as long; one half its mass lives six times longer.

Take the Sun's 10-billion-year main-sequence life as the benchmark. A star of 5 solar masses lives about 180 million years — fast, bright, and short-lived. A star of 0.5 solar masses lives roughly 56 billion years, longer than the present age of the universe. Every small red dwarf ever formed is still on the main sequence today. This single scaling explains why massive stars are rare: they simply do not last.

Check your understanding Beginner

Formal definition Intermediate+

A star is a self-gravitating nuclear engine that passes through a sequence of equilibrium configurations as its core composition changes. At each stage the four structure equations of 28.02.02 pending hold, but the composition profile evolves on the nuclear timescale, and the star drifts across the Hertzsprung-Russell (HR) diagram. This section traces that drift from the pre-main-sequence contraction through the asymptotic giant branch and binary mass transfer.

Pre-main-sequence contraction: the Hayashi and Henyey tracks

A protostar forms from a collapsing molecular cloud and contracts toward hydrostatic equilibrium along the Hayashi track, a nearly vertical path on the cool side of the HR diagram ^{[Hayashi 1961]}. The Hayashi track is a forbidden-region boundary: no star in hydrostatic equilibrium can lie to its right (cooler at fixed luminosity), because the envelope would become convectively unstable. Fully convective protostars descend the track at nearly constant temperature as they shrink and brighten, with deuterium burning providing a brief additional energy source above roughly $0.013M_{\odot }$.

As the core heats, a radiative zone develops and the contraction path bends onto the Henyey track, moving the star leftward and down toward the main sequence. The star settles onto the zero-age main sequence (ZAMS), the locus of stars that have just begun core hydrogen burning, when hydrostatic and thermal equilibrium are both established.

Main-sequence lifetime and composition changes

On the main sequence the star fuses hydrogen to helium. As the core helium abundance grows, the mean molecular weight $\mu$ rises and (through the virial relation of 28.02.02 pending) the core temperature climbs, so the nuclear rate and luminosity increase slowly even at fixed mass. A solar-mass star brightens by about 40% over its main-sequence life. The lifetime is set by the available fuel divided by the consumption rate:

$${\tau }_{\mathrm{MS}}\propto \frac{M}{L}\propto \frac{M}{M^{3.5}}=M^{-2.5},$$

using the empirical mass-luminosity relation $L\propto M^{3.5}$. Normalising to the Sun gives ${\tau }_{\mathrm{MS}}\approx 10^{10}(M/M_{\odot }{)}^{-2.5}$ years.

Main-sequence turn-off and cluster ages

In a coeval star cluster the most massive stars exhaust their core hydrogen first and peel off the main sequence, while lower-mass stars remain. The main-sequence turn-off (MSTO) marks the present-day luminosity of the most massive stars still burning core hydrogen, and hence fixes the cluster age through ${\tau }_{\mathrm{cluster}}\approx 10^{10}(M_{\mathrm{TO}}/M_{\odot }{)}^{-2.5}$ years. This method dates open clusters (the Pleiades at $\sim 10^8$ yr, the Hyades at $\sim 6\times 10^8$ yr) and globular clusters ($\sim 10^{10}$ yr), and historically provided the first reliable lower bound on the age of the universe.

The Schönberg-Chandrasekhar limit and the subgiant branch

When core hydrogen is exhausted, the star develops an isothermal helium core surrounded by a hydrogen-burning shell. An isothermal core cannot support the weight of the overlying envelope indefinitely. The Schönberg-Chandrasekhar limit caps the fraction of the total mass that such an inert core can hold in equilibrium ^{[Schönberg & Chandrasekhar 1942]}:

$$\frac{M_{\mathrm{core}}}{M_{\star }}\lesssim 0.37{\left(\frac{{\mu }_{\mathrm{env}}}{{\mu }_{\mathrm{core}}}\right)}^2,$$

where ${\mu }_{\mathrm{env}}$ and ${\mu }_{\mathrm{core}}$ are the mean molecular weights of the envelope and core. For typical compositions ${\mu }_{\mathrm{env}}\approx 0.6$ and ${\mu }_{\mathrm{core}}\approx 1.3$, giving $M_{\mathrm{core}}/M_{\star }\approx 0.1$. Once the core grows past this fraction the envelope contracts and the core heats, the hydrogen shell ignites, and the star expands and cools onto the subgiant branch.

The red giant branch, the helium flash, and the horizontal branch

On the red giant branch (RGB) the luminosity is generated almost entirely by a thin hydrogen-burning shell around a growing degenerate helium core. Because the core is electron-degenerate, its pressure is independent of temperature, so it cannot expand and cool as it gains mass; instead it grows denser and hotter until helium ignites. The star climbs the RGB at nearly constant temperature but rapidly increasing luminosity, with the envelope expanding to tens or hundreds of solar radii. A deep convective envelope dredges up processed material to the surface (the first dredge-up), lowering the surface hydrogen abundance and raising the helium and nitrogen abundances.

For stars below about $2M_{\odot }$ the helium core is degenerate when helium ignites, producing the helium flash: a thermonuclear runaway whose peak luminosity briefly exceeds that of the whole galaxy, though the energy is absorbed by the envelope and leaves no visible trace. The flash lifts the degeneracy, the core expands and cools, and the star settles onto the horizontal branch burning helium to carbon in a convective core. Horizontal-branch stars lie at roughly constant luminosity across a range of temperatures; in globular clusters they trace a horizontal band whose blue end hosts the pulsating RR Lyrae variables, standard candles used to measure distances within the Local Group. Stars above $2M_{\odot }$ ignite helium non-degenerately and quietly, then occupy the analogous red clump at the bright end of the red giants.

The asymptotic giant branch and thermal pulses

Once core helium is exhausted, a carbon-oxygen core forms surrounded by a helium-burning shell and, further out, a hydrogen-burning shell. The star ascends the asymptotic giant branch (AGB), a track that runs close to but slightly hotter than the earlier RGB. The two shells alternate dominance, producing cyclic thermal pulses: the helium shell ignites in a runaway, drives a convective episode that dredges carbon and s-process material to the surface (the third dredge-up), and then subsides while the hydrogen shell rebuilds the helium layer for the next pulse. In stars above about $4\text{–}5M_{\odot }$ the convective envelope reaches hot enough temperatures at its base for hot bottom burning, processing dredged-up carbon into nitrogen via the CNO cycle. The s-process (slow neutron capture) nucleosynthesis during AGB thermal pulses is the dominant site for half the elements heavier than iron, seeded by neutrons released in a ${}^{13}\mathrm{C}$ pocket.

Mass loss, planetary nebulae, and the initial-final mass relation

AGB stars shed mass at enormous rates through dust-driven winds. The Reimers formula parameterises the RGB and early-AGB mass-loss rate,

$$\dot{M}=-4\times 10^{-13}\eta \frac{(L/L_{\odot })(R/R_{\odot })}{M/M_{\odot }}{M}_{\odot }{\mathrm{yr}}^{-1},$$

with $\eta \sim 0.4$ fitted to observed RGB stars. On the AGB the rate rises to a superwind of $10^{-5}$ to $10^{-4}{M}_{\odot }{\mathrm{yr}}^{-1}$ that strips the envelope in under a million years. The ejected shell, ionised by the exposed hot core, glows as a planetary nebula, while the core, below the Chandrasekhar limit, cools into a white dwarf. The relation between a star's initial mass and the resulting white-dwarf mass is the initial-final mass relation, which shows that stars up to about $8M_{\odot }$ lose enough material to leave remnants between roughly $0.5$ and $1.1M_{\odot }$.

Binary evolution: Roche lobe overflow, common envelope, and Type Ia supernovae

More than half of all stars live in binary systems, and mass transfer reshapes their evolution. As a star expands to fill its Roche lobe — the teardrop-shaped region of gravitational dominance around each component — matter spills through the inner Lagrange point onto the companion ^{[Paczynski 1971]}. Stable Roche lobe overflow can rejuvenate the accretor with fresh hydrogen, producing Algol-type systems and, when the donor is massive, Wolf-Rayet binaries.

When mass transfer is unstable, the companion and the donor's core spiral together inside a shared gaseous common envelope, which is ejected by the orbital energy released. The outcome is a close binary of compact remnants — white dwarfs, neutron stars, or black holes — that cannot form through single-star evolution. A white dwarf accreting from a non-degenerate companion, or merging with another white dwarf, can be driven to the Chandrasekhar limit and ignite carbon degenerately throughout its volume, producing a Type Ia supernova. Because these explosions occur at a nearly fixed mass, Type Ia supernovae are nearly uniform in peak luminosity and serve as the standard candles that revealed the accelerating expansion of the universe.

Key derivation: main-sequence lifetime scaling and the Schönberg-Chandrasekhar limit Intermediate+

Two quantitative results govern the gross structure of post-main-sequence evolution: how long a star stays on the main sequence, and what triggers its departure. The first is a scaling argument; the second is a hydrostatic support argument. Together they set the clock and the trigger for red-giant expansion.

The main-sequence lifetime scaling

A star has a fuel reservoir proportional to its mass, $E_{\mathrm{fuel}}\propto M$, and consumes it at a rate proportional to its luminosity, $L$. Its main-sequence lifetime is therefore

$${\tau }_{\mathrm{MS}}\sim \frac{E_{\mathrm{fuel}}}{L}\propto \frac{M}{L}.$$

The mass-luminosity relation, derived from the structure equations and confirmed empirically across the main sequence, is approximately $L\propto M^{3.5}$ for stars between $0.5$ and $50M_{\odot }$. (The exponent is steeper for the most massive stars, where radiation pressure dominates, and flatter for the lowest masses, where convection changes the interior structure.) Substituting,

$${\tau }_{\mathrm{MS}}\propto \frac{M}{M^{3.5}}=M^{-2.5}.$$

Normalising to the Sun, whose main-sequence lifetime is $\sim 10^{10}$ years, gives

$${\tau }_{\mathrm{MS}}\approx 10^{10}{\left(\frac{M}{M_{\odot }}\right)}^{-2.5}\mathrm{yr}.$$

A $10M_{\odot }$ star lives $\sim 3\times 10^7$ years; a $0.5M_{\odot }$ star lives $\sim 6\times 10^{10}$ years. Because the lowest-mass stars outlive the present age of the universe, no M dwarf has yet left the main sequence, and the full red-giant evolution of such stars is inferred from models rather than observed.

The Schönberg-Chandrasekhar isothermal-core limit

When core hydrogen burning ceases, the helium core is nearly isothermal (no energy source within it) and is wrapped in a hydrogen-burning shell that adds mass to it. The question is whether the core can remain in hydrostatic equilibrium as it grows. Consider the core of mass $M_c$ and radius $R_c$, surrounded by an envelope of mass $M_{\mathrm{env}}\approx M_{\star }$. The pressure at the core-envelope interface must balance the weight of the envelope, $P_c\sim G{M}_{\mathrm{env}}{M}_c/(4\pi R_c^4)$. Within the isothermal core the ideal-gas pressure is $P_c\sim {\rho }_{c}k{T}_c/({\mu }_{\mathrm{core}}{m}_{\mathrm{H}})$, and the virial relation fixes the core temperature at $T_c\sim G{M}_c{\mu }_{\mathrm{core}}{m}_{\mathrm{H}}/(k{R}_c)$.

Balancing these two estimates for $P_c$ and solving for the limiting mass fraction yields the Schönberg-Chandrasekhar limit,

$$q_{\mathrm{SC}}\equiv \frac{M_c}{M_{\star }}\lesssim 0.37{\left(\frac{{\mu }_{\mathrm{env}}}{{\mu }_{\mathrm{core}}}\right)}^{2}f(n),$$

where $f(n)$ is a slowly varying factor of order unity depending on the polytropic structure of the envelope, and the mean-molecular-weight ratio appears because a helium-rich core (${\mu }_{\mathrm{core}}\approx 1.3$) is denser than a hydrogen-rich envelope (${\mu }_{\mathrm{env}}\approx 0.6$) at the same pressure. For a Population I star this gives $q_{\mathrm{SC}}\approx 0.10\text{–}0.12$. Once the shell has fed the core past roughly ten percent of the stellar mass, the core can no longer support the envelope, contracts on a Kelvin-Helmholtz timescale, and the envelope expands: the star leaves the main sequence. This is the physical trigger of the subgiant and red-giant phases, and it operates in every star massive enough to build an inert helium core — including the Sun, which will cross this threshold in roughly five billion years.

Exercises Intermediate+

Advanced results Master

Stellar evolution codes: MESA, KEPLER, TYCHO

Modern stellar evolution is computational. The structure equations of 28.02.02 pending are integrated forward in time as the composition profile evolves through nuclear burning and convective mixing. Three codes dominate the historical and current landscape. KEPLER, developed by Stan Woosley and collaborators, pioneered the modelling of massive stars through advanced burning and explosive nucleosynthesis, coupling a large nuclear network directly to the hydrostatic structure. TYCHO (Young and Arnett) emphasised opacity and equation-of-state sensitivity with an explicit hydrodynamics solver. The current workhorse is MESA (Modules for Experiments in Stellar Astrophysics; Paxton et al. 2011, 2013, 2015, 2018, 2019), an open-source, modular code with adaptive mesh refinement, tabulated OPAL and OP opacities, flexible equation-of-state and nuclear-network modules, and robust controls for handling phase transitions like the helium flash and thermal pulses. MESA's public release has made reproducible stellar evolution calculations accessible to the whole community, and its inlists (configuration files) underpin most published evolutionary tracks and isochrones.

Convective overshooting, semiconvection, and mixing

The boundary of a convective zone is not a sharp wall. Convective overshooting carries mixed material past the formal Schwarzschild boundary into the formally stable layer, because convective elements retain momentum as they reach the boundary. Overshooting enlarges the core, supplies fresh fuel, and lengthens the main-sequence lifetime, shifting the main-sequence turn-off and altering the resulting stellar masses inferred from cluster ages. The extent of overshooting is usually parameterised as a fraction ${\alpha }_{\mathrm{ov}}$ of the pressure scale height and calibrated against observed cluster morphology and eclipsing binaries.

In massive stars with a convective core and a composition gradient outside it, a different instability arises. Semiconvection occurs in a region that is stable to the Schwarzschild criterion but unstable to the Ledoux criterion (which accounts for composition gradients). The layer is marginally stable and mixes slowly, producing a partially mixed zone that strongly affects the size of the core entering helium burning and the subsequent evolutionary path. The treatment of semiconvection (diffusive versus overshoot-like) is one of the largest sources of uncertainty in massive-star evolution models.

Nuclear reaction networks and the s-process

The nuclear history of a star is governed by a reaction network whose size varies with the evolutionary phase. Main-sequence hydrogen burning needs only the PP chain or CNO cycle (a dozen isotopes). Helium burning adds the triple-alpha process and the alpha chain through carbon and oxygen. During AGB thermal pulses, the network must include the s-process (slow neutron capture), which builds roughly half the elements heavier than iron along the valley of beta stability. The dominant neutron source in low-mass AGB stars is the ${}^{13}\mathrm{C}(\alpha ,n){}^{16}\mathrm{O}$ reaction operating in a thin ${}^{13}\mathrm{C}$ pocket that forms at the base of the helium intershell during third dredge-up. In more massive AGB stars the ${}^{22}\mathrm{Ne}(\alpha ,n){}^{25}\mathrm{Mg}$ reaction dominates at higher temperatures. The resulting neutron exposure, interwoven with the thermal-pulse cycle and the dredge-up of processed material to the surface, sets the predicted s-process abundance distribution, which compares well with observed abundances in AGB stars and in the isotopic makeup of presolar silicon-carbide grains found in meteorites.

AGB thermal pulse cycle and hot bottom burning

An AGB thermal pulse is a limit cycle. During the quiescent interpulse phase, the hydrogen shell burns outward, depositing helium onto the intershell. When enough helium accumulates, the helium shell ignites in a flash, driving a convective episode that extends from the helium shell to the hydrogen shell. This convection dredges carbon (and s-process isotopes) to the surface — the third dredge-up — converting an oxygen-rich star into a carbon star once the surface C/O ratio exceeds unity. After the pulse, the hydrogen shell reignites and the cycle repeats, with pulse intervals of $10^4$ to $10^5$ years depending on core mass. In stars above roughly $4\text{–}5M_{\odot }$ the convective envelope penetrates deep enough that its base reaches temperatures above $\sim 4\times 10^7\mathrm{K}$, activating hot bottom burning, which processes the freshly dredged-up carbon back into nitrogen through the CNO bi-cycle and prevents the star from becoming a carbon star. Hot bottom burning is the leading site of primary nitrogen production in the universe.

Post-AGB evolution and the Born-Again scenario

After the AGB envelope is ejected, the exposed core moves rapidly across the HR diagram at nearly constant luminosity, ionising the ejected shell into a planetary nebula and then descending onto the white-dwarf cooling track. A small fraction of central stars undergo a late thermal pulse after they have already left the AGB, reigniting helium shell burning and driving the star back onto the AGB — the Born-Again scenario. The prototype is Sakurai's object (V4334 Sagittarii), which in the 1990s re-expanded and cooled on a timescale of years, demonstrating that the pulse cycle can operate after post-AGB departure and that the relevant nuclear and mixing timescales are far shorter than standard models predicted. Such objects are natural laboratories for mixing physics and nucleosynthesis on observable timescales.

Blue stragglers and the initial-final mass relation

Blue stragglers are stars that appear on or above the main-sequence turn-off of a cluster, where they should not exist if all cluster members are coeval. Their formation channels are mass gain: either direct stellar collisions in the dense cores of globular clusters, or Roche-lobe mass transfer in a close binary that rejuvenates the accretor with fresh hydrogen. The collisional channel produces a merged star whose interior is chemically anomalous (no normal core-envelope structure), while the mass-transfer channel can produce blue stragglers with depleted carbon and oxygen at the surface from the processed donor material. Their presence biases simple cluster-age estimates from the main-sequence turn-off.

The initial-final mass relation (IFMR) maps a star's zero-age main-sequence mass to the mass of the white dwarf it leaves behind, after all mass loss on the RGB, AGB, and superwind phases. Empirically determined from white dwarfs in open clusters, the IFMR is roughly linear above $\sim 3M_{\odot }$ but flattens toward a final mass of $\sim 0.5M_{\odot }$ for the lowest-mass progenitors. The relation is a key ingredient in white-dwarf cosmochronology (dating stellar populations from their coolest white dwarfs) and in determining the chemical enrichment budget of galaxies.

Rotation and magnetic fields in evolution

Stellar rotation induces meridional circulation and shear instabilities that mix material beyond the convective boundaries, in competition with composition gradients that resist mixing. Rotational mixing can bring fresh hydrogen into the core during main-sequence burning, extending the lifetime and altering the surface abundances (producing the nitrogen enhancements and helium enrichments observed in fast rotators). It also transports angular momentum inward, with consequences for the spin of the eventual compact remnant. Magnetic fields couple to rotation through the Spruit-Tayler dynamo, which generates a toroidal field from differential rotation and uses it to transport angular momentum far more efficiently than purely hydrodynamic processes. The interplay of rotation, magnetic torques, and mass loss determines the angular-momentum content of stellar remnants and is essential for understanding the spin distribution of neutron stars and black holes detected via gravitational waves.

Massive star evolution: Wolf-Rayet stars, LBVs, and pair-instability

Stars above about $8M_{\odot }$ ignite helium, carbon, neon, oxygen, and silicon in successive convective cores, each stage shorter than the last (for a $25M_{\odot }$ star: hydrogen $\sim 7\mathrm{Myr}$, helium $\sim 0.5\mathrm{Myr}$, carbon $\sim 600\mathrm{yr}$, oxygen $\sim 0.5\mathrm{yr}$, silicon $\sim 1\mathrm{day}$). Intense radiation-driven winds strip the hydrogen envelope. Once the hydrogen envelope is gone, the star appears as a Wolf-Rayet star, spectroscopically dominated by broad emission lines of helium, nitrogen (WN subtype, CNO-processed), carbon, and oxygen (WC/WO subtypes, helium-burning products). The most luminous hydrogen-bearing supergiants are unstable to the luminous blue variable (LBV) phase, characterised by giant eruptions (Eta Carinae being the prototype) that can eject many solar masses in decades.

At very high masses ($\gtrsim 150M_{\odot }$), the core reaches temperatures where photons pair-produce electron-positron pairs, softening the equation of state and reducing radiation pressure support. The resulting pair-instability can pulsate (pulsational pair-instability supernovae, shedding shells) or, above roughly $260M_{\odot }$, disrupt the entire star in a pair-instability supernova that leaves no remnant and synthesises large amounts of nickel-56. The predicted mass gap in remnant masses between $\sim 50$ and $\sim 120M_{\odot }$ is an active target of gravitational-wave observations.

Electron-capture supernovae, failed supernovae, and silent collapse

At the boundary between white-dwarf progenitors and core-collapse progenitors ($8\text{–}10M_{\odot }$), stars develop degenerate oxygen-neon-magnesium cores. When electron captures on magnesium and neon reduce the electron pressure, the core collapses before silicon burning can stabilise it, producing an electron-capture supernova. These events are fainter than standard iron-core-collapse supernovae and may account for some of the diversity in observed transients. In the ongoing failed supernova programme, surveys of nearby massive stars have identified candidates that brighten briefly and then vanish, consistent with a star whose core collapsed into a black hole without a luminous explosion — a silent collapse. If confirmed at a rate of tens of percent among massive stars, failed supernovae would resolve the long-standing discrepancy between the observed core-collapse supernova rate and the massive-star death rate.

Stellar populations I, II, and III

Stellar evolution is framed against the chemical history of the galaxy through the Population scheme. Population I stars (like the Sun) are young, metal-rich, and found in the galactic disk; their high metal content reflects billions of years of prior nucleosynthesis. Population II stars are old, metal-poor, and found in the halo and globular clusters; they formed early and carry the imprint of only a few generations of prior stars. Population III stars are the hypothesised first generation, formed from primordial Big Bang material (hydrogen, helium, trace lithium) with effectively zero metals. Because metal-free opacity is far lower, Population III stars are predicted to be very massive ($\gtrsim 100M_{\odot }$) and to follow entirely different evolution tracks — no CNO cycle at ignition, no dust-driven mass loss, and predominantly pair-instability or direct-collapse deaths. No Population III star has been directly observed; the search proceeds via deep JWST imaging of the early universe and via the chemical fingerprints (large carbon and oxygen enhancements) imprinted on the most iron-poor Population II stars, which may be the polluted descendants of Population III supernovae.

Nebular diagnostics of evolved stars

The circumstellar material ejected during the AGB, planetary-nebula, and supernova phases carries a record of the star's nucleosynthetic history. Planetary nebulae provide abundance diagnostics through collisionally excited optical and ultraviolet lines, which reveal the products of first, second, and third dredge-ups (helium, nitrogen, carbon, and s-process enrichments) and allow direct tests of AGB nucleosynthesis models. The [O III], [N II], and [S II] forbidden-line ratios classify nebulae by excitation and constrain the central star temperature and luminosity. Supernova remnants show ejecta-dominated interiors rich in oxygen, neon, and magnesium (from massive-star interiors) or iron and nickel (from explosive silicon burning and Type Ia detonation), observable in X-ray lines that map the layered structure predicted by the progenitor's evolution. These diagnostics close the loop between stellar evolution theory and the chemical enrichment of the interstellar medium that seeds the next generation.

Connections Master

Connection to stellar structure 28.02.02 pending

This unit is the direct continuation of 28.02.02 pending. The four structure equations and the equation of state derived there supply the instantaneous physics; stellar evolution is what happens as the composition profile changes on the nuclear timescale. The hydrostatic equilibrium and energy-transport equations fix where a star of given mass and composition sits on the HR diagram, while the nuclear generation laws set the rate at which that position drifts. The virial theorem and degeneracy pressure, both introduced in 28.02.02 pending, are what make the helium flash, the core-mass luminosity relation on the RGB, and the Chandrasekhar-limited ignition of Type Ia supernovae physically possible.

Connection to the Hertzsprung-Russell diagram 28.02.01

The HR diagram of 28.02.01 is the canvas on which stellar evolution is drawn. The main sequence, the red giant branch, the horizontal branch, the asymptotic giant branch, and the white-dwarf cooling sequence are all regions of that diagram, and the evolution tracks of this unit are the curves connecting them. The main-sequence turn-off method for dating clusters, introduced qualitatively in 28.02.01, is here given its quantitative basis through the main-sequence lifetime scaling ${\tau }_{\mathrm{MS}}\propto M^{-2.5}$ and the Schönberg-Chandrasekhar trigger. The mass-luminosity relation of the earlier unit becomes the input that fixes the rate at which stars consume their fuel.

Forward connection to stellar remnants 28.02.04 pending

The terminal phases described here set the initial conditions for the stellar remnants of 28.02.04 pending. A star's mass and its mass-loss history decide whether it leaves a white dwarf (initial mass below $\sim 8M_{\odot }$), a neutron star (roughly $8$ to $25M_{\odot }$), or a black hole (above $\sim 25M_{\odot }$, modulo rotation and mass loss). The initial-final mass relation governs the white-dwarf mass spectrum; the iron-core mass and the explosion energy determine whether core collapse succeeds in launching a supernova or fails silently into a black hole. Binary evolution, especially common-envelope ejection and white-dwarf mergers, produces the close compact-object binaries that 28.02.04 pending treats as gravitational-wave sources.

Connection to nuclear physics and nucleosynthesis

Stellar evolution is, at root, a sequence of nuclear burning stages, each governed by cross sections measured (or extrapolated) in the laboratory. The s-process nucleosynthesis during AGB thermal pulses, the alpha chain in massive-star cores, and the explosive nucleosynthesis in supernovae together account for the origin of essentially every element heavier than helium. Reaction-rate uncertainties — particularly for the ${}^{12}\mathrm{C}(\alpha ,\gamma ){}^{16}\mathrm{O}$ rate that sets the carbon-to-oxygen ratio at the end of helium burning — propagate directly into the predicted stellar evolution tracks, the final remnant masses, and the chemical enrichment of galaxies. Underground facilities such as LUNA push direct measurements closer to stellar energies and are steadily reducing these uncertainties.

Connection to cosmology and the first stars 28.04.01

Stellar evolution feeds back into cosmology through chemical enrichment, ionising radiation, and the formation of the first stars. The Population III stars of the early universe, formed from primordial material, followed evolution tracks entirely different from present-day stars because of their near-zero metallicity, and their pair-instability supernovae seeded the interstellar medium with the first heavy elements. The cosmic star-formation history, the reionisation of the intergalactic medium, and the chemical evolution of galaxies are all governed by how stellar populations of different masses and metallicities evolve and die. The Hubble constant tension between cluster ages and the expansion age, historically a major puzzle, was resolved only once stellar evolution models, opacities, and the discovery of cosmic acceleration were all brought into agreement.

Historical and philosophical context Master

The recognition that stars evolve — that they are not fixed lights but objects with histories — was a twentieth-century achievement. The Hertzsprung-Russell diagram, drawn independently by Ejnar Hertzsprung (1911) and Henry Norris Russell (1913), was at first interpreted by Russell as an evolutionary sequence: stars were thought to be born as hot red giants, contract and heat onto the main sequence, then cool down it into red dwarfs. This picture was backwards. Once nuclear energy was identified as the stellar power source (Bethe, 1939), it became clear that the main sequence is not an evolutionary track but the locus of hydrogen-burning stars, and that evolution carries stars off it, not along it.

The quantitative theory of post-main-sequence evolution emerged in the 1950s and 1960s with the arrival of electronic computation. Chushiro Hayashi's 1961 derivation of the forbidden region that bears his name explained why pre-main-sequence and red-giant stars are confined to the cool side of the HR diagram. The Schönberg-Chandrasekhar limit (Schönberg and Chandrasekhar, 1942) had already explained why stars must leave the main sequence when their inert cores reach roughly ten percent of the total mass, but its full implication for red-giant expansion awaited detailed numerical models. Louis Henyey and his collaborators developed the Henyey relaxation method, a stable algorithm for solving the two-point boundary-value problem of stellar structure across time steps, which became the backbone of every subsequent stellar evolution code.

The landmark synthesis was Icko Iben's 1967 review, "Stellar evolution within and off the main sequence" ^{[Iben 1967]}, which presented the first comprehensive set of computed evolution tracks spanning the main sequence, the subgiant branch, the red giant branch, the helium flash, the horizontal branch, and the asymptotic giant branch, tied together by a single consistent numerical framework. Iben's tracks showed how a star's mass determines its path across the HR diagram and the timing of each phase, and they placed the qualitative narrative of stellar lives on a quantitative footing. The same period saw Bohdan Paczynski's work on close-binary evolution and mass transfer ^{[Paczynski 1971]}, which opened the study of how binary interaction reshapes stellar fates — producing cataclysmic variables, X-ray binaries, Type Ia supernovae, and the compact-object mergers now detected by gravitational-wave observatories.

The computational era, from KEPLER in the 1970s to MESA in the 2010s, has steadily enlarged the domain of stellar evolution that can be modelled from first principles. The helium flash, once a numerical impasse because of the vast disparity between the flash timescale and the evolutionary timescale, is now routinely followed. Thermal-pulse AGB models reproduce the observed carbon-star luminosity functions and s-process abundance patterns. Massive-star models with rotation and magnetic torques reproduce the observed nitrogen enhancements and the spin distribution of young pulsars. Yet central uncertainties remain: the physics of convective boundaries (overshooting, semiconvection), the treatment of mass loss, and the explosion mechanism of core-collapse supernovae are all parameterised rather than derived, and they together set the dominant theoretical error bars on stellar ages, remnant masses, and nucleosynthetic yields.

The philosophical weight of stellar evolution is that it places every heavy atom in a historical context. The carbon in a cell, the oxygen in a breath, the iron in blood, and the iodine in a thyroid were all forged at specific stages of specific stars' lives — most of them during the AGB thermal pulses and supernova explosions traced in this unit. To follow an evolution track across the HR diagram is to follow the origin of the chemical elements of which planets, oceans, and living bodies are made. The starstuff aphorism is, in this precise sense, a statement about stellar evolution: we are what stars, in their dying phases, ejected into the gas from which the Sun and the Earth later condensed.

Bibliography Master

1. Schönberg, M. & Chandrasekhar, S. (1942). "On the Evolution of the Main-Sequence Stars." Astrophysical Journal 96, 161-172. Derivation of the isothermal-core mass limit that triggers red-giant expansion.

2. Hayashi, C. (1961). "Stellar Evolution in Early Phases of Gravitational Contraction." Publications of the Astronomical Society of Japan 13, 450-452. The Hayashi track and the forbidden region on the cool side of the HR diagram.

3. Henyey, L. G., Wilets, L., Böhm, K. H., Lelevier, R. & Levee, R. D. (1959). "A Method for Automatic Computation of Stellar Evolution." Astrophysical Journal 129, 628-636. The Henyey relaxation method that made numerical stellar evolution tractable.

4. Iben, I. (1967). "Stellar Evolution within and off the Main Sequence." Annual Review of Astronomy and Astrophysics 5, 571-626. The first comprehensive computed evolution tracks from the main sequence through the asymptotic giant branch.

5. Paczynski, B. (1971). "Evolutionary Processes in Close Binary Systems." Annual Review of Astronomy and Astrophysics 9, 183-208. The framework for Roche-lobe overflow, common-envelope evolution, and binary mass transfer.

6. Reimers, D. (1975). "Circumstellar Absorption Lines and Mass Loss from Red Giants." Memoires of the Societe Royale des Sciences de Liege 8, 369-382. The Reimers mass-loss formula for red giant branch stars.

7. Kippenhahn, R. & Weigert, A. (1990). Stellar Structure and Evolution. Springer. The standard graduate reference for stellar evolution tracks, the helium flash, and the asymptotic giant branch.

8. Iben, I. & Renzini, A. (1983). "Asymptotic Giant Branch Evolution and Beyond." Annual Review of Astronomy and Astrophysics 21, 271-342. The thermal-pulse cycle, dredge-ups, and the transition to the white-dwarf cooling track.

9. Carroll, B. W. & Ostlie, D. A. (2017). An Introduction to Modern Astrophysics (2nd ed.). Cambridge. Chapter 13 gives the undergraduate treatment of stellar evolution and HR-diagram tracks.

10. Paxton, B. et al. (2011). "Modules for Experiments in Stellar Astrophysics (MESA)." Astrophysical Journal Supplement 192, 3. The open-source stellar evolution code and its application to evolution tracks.

11. Herwig, F. (2005). "Evolution of Asymptotic Giant Branch Stars." Annual Review of Astronomy and Astrophysics 43, 435-479. Modern review of AGB thermal pulses, third dredge-up, and s-process nucleosynthesis.

12. Woosley, S. E., Heger, A. & Weaver, T. A. (2002). "The Evolution and Explosion of Massive Stars." Reviews of Modern Physics 74, 1015-1071. Comprehensive models of massive-star evolution through advanced burning and core collapse.