Stellar structure: hydrostatic equilibrium, nuclear burning, the PP chain and CNO cycle

Intuition Beginner

A star is a giant ball of hot gas held together by gravity. The enormous weight of its outer layers pushes inward, while the energy from nuclear fusion in the core pushes outward. These two forces balance perfectly. This balance is called hydrostatic equilibrium, and it keeps a star stable for billions of years. If the core cools even slightly, gravity squeezes it harder, heating it and speeding up fusion; if it overheats, the star expands and cools. This self-correcting thermostat is what holds a star steady.

In the Sun's core, at about 15 million degrees and crushing pressure, hydrogen nuclei smash together and fuse into helium through the proton-proton chain, or PP chain. Every second the Sun consumes roughly 600 million tons of hydrogen. About 0.7% of that mass becomes pure energy, following Einstein's relation E = mc². Produced deep in the core, this energy takes many thousands of years to fight its way outward, finally escaping at the surface as sunlight.

Stars heavier than about 1.3 solar masses burn hotter. Their cores reach temperatures where a different fusion pathway takes over, the CNO cycle. Here carbon, nitrogen, and oxygen act as catalysts, shuttling protons through a chain of reactions that still turns four hydrogens into one helium. The CNO cycle is far more sensitive to temperature than the PP chain, so it dominates in massive stars and barely runs in the Sun. A star's mass decides not only how brightly it shines but which engine lights its fire.

Visual Beginner

A star's interior splits into concentric zones, each moving energy outward by a different mechanism. The diagram shows the balance of forces and the runaway temperature sensitivity that separates the two fusion engines.

Zone	What happens	Energy transport
Core	Hydrogen fuses to helium; energy is released	Radiation
Radiative zone	Photons random-walk outward, absorbed and re-emitted	Radiation
Convective zone	Hot gas rises, cools, and sinks in rolling cells	Convection
Photosphere	The visible surface; photons escape to space	Radiation

The Sun has a radiative zone and a convective envelope. Red dwarfs are fully convective, while high-mass stars are convective only in a small core. This zoning is set by how steeply temperature rises inward and by how opaque the gas is to radiation.

Worked example Beginner

The Sun radiates about 3.8 × 10²⁶ watts. That energy comes from mass turned into energy, so the Sun loses mass at a rate fixed by E = mc². Dividing the luminosity by c² gives a mass loss of roughly 4.3 billion kilograms every second. This is the mass actually converted to energy; the hydrogen consumed is far larger, because only 0.7% of the fused mass becomes energy and the rest becomes helium.

Only the hydrogen in the hot core can fuse, about 10% of the Sun's initial 70% hydrogen supply. The convertible mass is 0.7% of that reservoir, near 10²⁷ kilograms. Dividing by the loss rate of 4.3 billion kilograms per second gives roughly 2 × 10¹⁷ seconds, or about 7 billion years. This is close to the Sun's expected main-sequence lifetime of around 10 billion years. The rough agreement is a triumph of the simple idea that a star shines by converting mass into energy.

Check your understanding Beginner

Formal definition Intermediate+

A star in hydrostatic equilibrium is described by four coupled ordinary differential equations. Treat the star as a spherically symmetric fluid with pressure $P(r)$, density $\rho (r)$, temperature $T(r)$, enclosed mass $M(r)$, and local luminosity $L(r)$, all functions of the radial coordinate $r$.

Hydrostatic equilibrium balances the inward pull of gravity against the outward pressure gradient:

$$\frac{dP}{dr}=-\frac{GM(r)\rho (r)}{r^2}.$$

Mass continuity states how enclosed mass grows with radius:

$$\frac{dM}{dr}=4\pi r^2\rho (r).$$

Energy generation accumulates the luminosity produced by nuclear reactions and gravitational contraction:

$$\frac{dL}{dr}=4\pi r^2\rho (r)\epsilon (r),$$

where $\epsilon$ is the nuclear energy generation rate per unit mass.

Energy transport sets the temperature gradient. In a radiative region,

$$\frac{dT}{dr}=-\frac{3\kappa (r)\rho (r)L(r)}{16\pi acT(r{)}^3{r}^2},$$

with $\kappa$ the opacity, $a$ the radiation constant, and $c$ the speed of light. When the radiative gradient exceeds the adiabatic gradient, the gas becomes convectively unstable and energy is carried by bulk fluid motion instead.

These four equations close with an equation of state $P=P(\rho ,T,\text{composition})$, an opacity law $\kappa =\kappa (\rho ,T,\text{composition})$, and a nuclear generation law $\epsilon =\epsilon (\rho ,T,\text{composition})$. Together with boundary conditions, they determine the run of every physical quantity from centre to surface.

Boundary conditions and equations of state

At the centre, $r=0$ forces $M(0)=0$ and $L(0)=0$, with $P$ and $T$ taking their central (finite) values. At the surface, one imposes $R$ as the radius where $P\to 0$ and $T\to T_{\mathrm{eff}}$, with $M(R)=M_{\star }$ and $L(R)=L_{\star }$. The problem is a two-point boundary value problem, solved numerically by relaxation or shooting methods.

The pressure receives up to three contributions. For an ideal ionised gas, $P_{\mathrm{gas}}=\rho kT/(\mu m_{\mathrm{H}})$ with $\mu$ the mean molecular weight. Radiation contributes $P_{\mathrm{rad}}=\frac{1}{3}a{T}^4$. At very high density and low temperature, electrons fill all low quantum states and exert degeneracy pressure independent of temperature, $P_{\mathrm{deg}}\propto {\rho }^{5/3}$ (non-relativistic) or ${\rho }^{4/3}$ (relativistic). Degeneracy pressure dominates in white dwarfs and in the cores of low-mass giants, and its temperature-independence is what permits thermonuclear runaway.

Opacity and the nuclear energy generation rates

Opacity measures how strongly the stellar gas blocks photons. The leading contributions are electron (Thomson) scattering ${\kappa }_{\mathrm{es}}\approx 0.2(1+X){\mathrm{m}}^2\mathrm{k}{\mathrm{g}}^{-1}$, important in hot massive stars, and free-free and bound-free absorption (Kramers opacity), ${\kappa }_{\mathrm{ff}},{\kappa }_{\mathrm{bf}}\propto \rho T^{-3.5}$. Detailed opacity tables (OPAL, OP) tabulate $\kappa$ across all relevant conditions.

Nuclear burning converts mass to energy with extreme temperature sensitivity. The proton-proton chain rate scales as

$${\epsilon }_{\mathrm{pp}}\propto \rho X^2{T}^4,$$

while the CNO cycle scales as

$${\epsilon }_{\mathrm{CNO}}\propto \rho X{X}_{\mathrm{CNO}}{T}^{17},$$

where $X$ is the hydrogen mass fraction and $X_{\mathrm{CNO}}$ the combined carbon-nitrogen-oxygen abundance. The steep CNO temperature dependence concentrates the energy production of a massive star into a tiny central volume, which drives a steep temperature gradient and a convective core.

The PP chain and the CNO cycle

The net PP chain reaction fuses four protons into helium-4:

$$4{}^1\mathrm{H}⟶{}^4\mathrm{He}+2e^{+}+2{\nu }_e+\gamma ,$$

releasing about $26.7\mathrm{MeV}$, of which roughly $2\%$ escapes in neutrinos. The CNO cycle has the same net reaction but routes the protons through proton captures and beta decays on carbon, nitrogen, and oxygen nuclei. In the Sun the PP chain dominates; in stars above about $1.3M_{\odot }$ the CNO cycle takes over.

The Standard Solar Model and the neutrino problem

The Standard Solar Model (Bahcall and collaborators) integrates the four structure equations calibrated to the Sun's present mass, radius, luminosity, age, and composition. For decades Ray Davis's Homestake experiment detected only about one third of the electron neutrinos the model predicted. This solar neutrino problem was resolved in 2001-2002 by the Sudbury Neutrino Observatory, which showed that the missing neutrinos had oscillated into muon and tau flavours in transit ^{[Bahcall et al. 2001]}. The model was correct; the neutrinos have mass and change identity.

Key result: the virial theorem and the Kelvin-Helmholtz timescale Intermediate+

The virial theorem links a star's internal thermal energy to its gravitational binding energy and yields the timescale on which a star would fade if its fuel were cut off. It is the cleanest single result of stellar structure theory, and historically it is what proved that gravity alone cannot power the Sun.

Start from hydrostatic equilibrium, $\frac{dP}{dr}=-GM(r)\rho /r^2$, and multiply both sides by $4\pi r^3$, then integrate from the centre to the surface radius $R$. Using $dM=4\pi r^2\rho dr$ on the right-hand side gives

$${\int }_0^{R}4\pi r^3\frac{dP}{dr}dr=-{\int }_0^M\frac{Gm}{r(m)}dm=\Omega ,$$

where $\Omega =-{\int }_0^{M}Gm/rdm<0$ is the total gravitational potential energy. The left-hand side is integrated by parts; the boundary term $[4\pi r^{3}P{]}_0^R$ vanishes because $P$ is zero at the surface and $r^3$ is zero at the centre, leaving

$$-3{\int }_0^{R}4\pi r^{2}Pdr=\Omega ,$$

that is, $3\int PdV=-\Omega$. For a monatomic ideal gas the internal thermal energy is $K=\frac{3}{2}\int PdV$, so

$$2K+\Omega =0.$$

This is the virial theorem for a self-gravitating gas sphere. Two consequences follow immediately. First, the mean temperature satisfies $\bar{T}\sim GM\mu m_{\mathrm{H}}/(kR)$: gravity sets the core temperature. Second, the total energy $E=K+\Omega =\Omega /2$ is negative, so as a star radiates energy away and $E$ grows more negative, $K$ must increase, the star must heat up as it contracts. A contracting star does not cool; it heats.

The Kelvin-Helmholtz timescale estimates how long the Sun could shine on gravitational contraction alone:

$${\tau }_{\mathrm{KH}}\sim \frac{|\Omega |}{L}\sim \frac{G{M}^2}{RL}\approx 3\times 10^7\mathrm{years}$$

for the Sun. Against the geological evidence for an Earth far older than this, nineteenth-century physics (Kelvin, Helmholtz) could not reconcile the Sun's brightness with its age. Only the discovery of nuclear energy and $E=m{c}^2$ supplied a power source lasting the required billions of years, giving the nuclear timescale ${\tau }_{\mathrm{nuc}}\sim 0.007M{c}^2/L\approx 10^{10}$ years. The three stellar timescales, dynamical ($\sim 1$ hour), thermal ($\sim 10^7$ years), and nuclear ($\sim 10^{10}$ years), are separated by many orders of magnitude, which is precisely why stars can be treated as passing through a sequence of hydrostatic equilibria.

Exercises Intermediate+

Advanced results Master

The proton-proton chain branches

Hydrogen burning through the PP chain is not a single reaction but a network of branches distinguished by how the intermediate helium-3 is consumed. In the dominant ppI branch, two helium-3 nuclei combine: ${}^3\mathrm{He}+{}^3\mathrm{He}\to {}^4\mathrm{He}+2{}^1\mathrm{H}$. The ppII and ppIII branches instead capture ${}^3\mathrm{He}$ on an alpha particle to form beryllium-7, which then either captures an electron (ppII, yielding lithium-7 that is broken up by a proton) or captures a proton (ppIII, yielding boron-8 that beta-decays). The ppIII branch produces the high-energy boron-8 neutrinos, the principal target of the chlorine and water Cherenkov detectors. Two side branches complete the picture: pep ($p+e^{-}+p\to {}^2\mathrm{H}+\nu$, producing a monoenergetic 1.44 MeV neutrino) and hep (${}^3\mathrm{He}+p\to {}^4\mathrm{He}+e^{+}+\nu$, producing the highest-energy solar neutrinos at up to 18.8 MeV but at a negligible rate). The branching ratios between ppI/ppII/ppIII are sensitive thermometers of the solar core and are reproduced to within experimental uncertainty by the Standard Solar Model.

The CNO bi-cycle and explosive variants

The CNO cycle is itself a bicycle. The primary (CNO or CN) branch cycles through ${}^{12}\mathrm{C}\to {}^{13}\mathrm{N}\to {}^{13}\mathrm{C}\to {}^{14}\mathrm{N}\to {}^{15}\mathrm{O}\to {}^{15}\mathrm{N}$, where ${}^{15}\mathrm{N}(p,\alpha ){}^{12}\mathrm{C}$ closes the loop. The slowest link is the beta decay of ${}^{15}\mathrm{O}$, which bottlenecks the cycle and concentrates most of the catalytic material in ${}^{14}\mathrm{N}$. A secondary branch leaks at ${}^{15}\mathrm{N}(p,\gamma ){}^{16}\mathrm{O}$, feeding the NO cycle that produces ${}^{16}\mathrm{O}$ and ${}^{17}\mathrm{O}$ and returning to ${}^{14}\mathrm{N}$, while a further arm reaches fluorine and neon-sodium. In the hot, hydrogen-rich envelopes of novae and X-ray burst accretors, the beta-decay bottleneck is bypassed by proton captures faster than the decays: the hot CNO cycle and the rp-process (rapid proton capture) process seed material up to the iron-group and beyond, releasing energy at rates that power the observed outbursts. These networks are now mapped by detailed reaction-rate libraries such as JINA REACLIB.

Neutrino energy loss

At temperatures above about $10^9\mathrm{K}$, reached in advanced burning stages and in stellar collapse, neutrinos carry away energy far faster than photons. Four mechanisms operate: pair neutrinos ($\gamma +\gamma \to e^{+}+e^{-}\to \nu +\bar{\nu }$), photoneutrinos ($\gamma +e^{-}\to e^{-}+\nu +\bar{\nu }$), plasma neutrinos (decay of a plasmon into a neutrino pair), and bremsstrahlung neutrinos (electron-ion scattering radiating a neutrino pair). Because neutrinos escape the star without interaction, they act as a powerful cooling channel that drives the rapid late-stage evolution of massive stars and removes most of the binding energy of a collapsing core in a single burst.

Polytropes and the Lane-Emden equation

A polytropic equation of state $P=K{\rho }^{1+1/n}$ reduces the structure equations to a single dimensionless ordinary differential equation, the Lane-Emden equation:

$$\frac{1}{{\xi }^2}\frac{d}{d\xi }\left({\xi }^2\frac{d\theta }{d\xi }\right)=-{\theta }^n,$$

with $\rho ={\rho }_c{\theta }^n$ and $r=\alpha \xi$. The index $n=3$ (Eddington's model) approximates a radiation-pressure-supported massive main-sequence star and yields a mass-radius relation insensitive to central density, giving roughly $R\propto M^{-1/2}$ for massive stars; the index $n=1.5$ describes a fully convective star governed by an adiabatic $\gamma =5/3$ gas law and fits the lower main sequence and the Hayashi track. Analytic solutions exist for $n=0,1,5$; all others are integrated numerically. The Eddington standard model (radiative transport, constant $\kappa /P$ ratio, $n=3$) gives the celebrated approximate mass-luminosity relation $L\propto M^3$, sharpened empirically to $L\propto M^{3.5}$ over much of the main sequence.

Opacity tables, MESA, and helioseismology

Modern stellar modelling replaces Kramers power laws with tabulated opacities from the OPAL (Lawrence Livermore) and OP (Opacity Project) collaborations, which include millions of bound-bound, bound-free, and free-free transitions for arbitrary compositions. The open-source MESA code (Modules for Experiments in Stellar Astrophysics) integrates the full structure equations with adaptive mesh refinement, detailed nuclear networks, and up-to-date opacities and equation-of-state tables, and is now the workhorse of stellar evolution.

The strongest observational test of solar structure comes from helioseismology. The Sun oscillates in millions of acoustic p-modes and gravity g-modes, and the precise frequencies of these oscillations map the sound-speed profile $c_s(r)=\sqrt{\gamma P/\rho }$ from the surface down to the core. The Standard Solar Model reproduces the p-mode spectrum to within a fraction of a percent over most of the solar interior. The residual mismatch, the solar abundance problem, arose when Asplund and collaborators revised the solar metallicity $Z$ downward (around 2009) using three-dimensional hydrodynamic atmosphere models, bringing the model's sound-speed profile into conflict with helioseismology. The tension is unresolved and motivates ongoing work on opacities, diffusion, and additional mixing.

Pre-main-sequence contraction and structural limits

Before hydrogen ignition, a protostar contracts down the Hayashi track, a nearly vertical, fully convective path in the Hertzsprung-Russell diagram set by the requirement that the envelope remain convectively unstable at low surface temperature. As the contracting core heats, radiative cores grow and the star moves toward the main sequence. T Tauri stars are the observational counterparts of this phase. Three structural limits bound the outcomes. Below about $0.08M_{\odot }$ the core never reaches hydrogen-burning temperature and the object becomes a brown dwarf (transient deuterium burning may occur above roughly $0.013M_{\odot }$). The lithium depletion boundary, the mass below which lithium survives because the core never becomes hot enough to burn it, ages clusters independently of stellar models. In an evolved star with an inert helium core surrounded by a hydrogen-burning shell, the Schönberg-Chandrasekhar limit caps the isothermal core at roughly $10\%$ of the total mass; beyond it the core cannot support the envelope and contracts rapidly, triggering the red-giant phase.

Connections Master

Connection to stars and stellar evolution 28.02.01

This unit takes the Hertzsprung-Russell diagram and the mass-luminosity relation of 28.02.01 as its starting point and supplies the interior physics that explains them. Hydrostatic equilibrium plus the ideal-gas equation of state yields, through the virial theorem, the relation between mass, radius, and core temperature that fixes a star's position on the main sequence. The nuclear generation laws derived here are what give the mass-luminosity relation $L\propto M^{3.5}$ its physical content, and they set the main-sequence lifetimes summarised in the companion unit. The forward link is to 28.02.03 pending, where stellar evolution tracks follow directly from how the structure equations respond as the core composition changes.

Connection to planetary interiors 28.01.03 pending

Planets and stars are both self-gravitating fluid bodies obeying a hydrostatic equilibrium equation, and the same machinery, the equation of state, the energy-transport equation, and the distinction between radiative and convective zones, applies in both settings 28.01.03 pending. The differentiating physics is the energy source: stars run nuclear reactions that sustain their luminosity for gigayears, while planets cool by secular contraction and radioactivity. The polytrope formalism developed for stars is routinely adapted to model giant-planet interiors, and the degeneracy pressure that supports white dwarfs also supports the cores of Jupiter and Saturn.

Connection to cosmology and primordial abundances 28.04.01

The CNO cycle depends on the carbon, nitrogen, and oxygen catalysts produced by earlier generations of stars, so the relative importance of the PP chain and the CNO cycle has shifted across cosmic time. The very first stars, formed from primordial Big Bang nucleosynthesis material, were essentially metal-free and could not run the CNO cycle at all, which changes their structure, their lifetimes, and their explosive fates. The primordial hydrogen and helium abundances fixed during the Big Bang 28.04.01 are the raw fuel for all subsequent stellar burning, tying stellar structure directly to the thermal history of the early universe.

Connection to nuclear and particle physics

Stellar interiors are the universe's nuclear physics laboratories, and every reaction rate in the PP chain and CNO cycle is a nuclear cross section extrapolated from laboratory measurements because stellar energies are far below the Coulomb barrier. The solar neutrino problem, resolved by neutrino oscillations, was the first direct evidence that neutrinos have mass and was a decisive extension of the Standard Model of particle physics. The same weak interactions that set the neutrino flux also fix the rates of the pep and hep side branches, so a star is simultaneously a nuclear reactor and a particle-physics experiment.

Historical and philosophical context Master

The problem of the Sun's energy source dominated nineteenth-century physics. Kelvin and Helmholtz argued that the Sun shone by slow gravitational contraction, but the resulting Kelvin-Helmholtz timescale of some tens of millions of years collided with the geological and evolutionary timescales then emerging from the work of Lyell and Darwin. The contradiction was real, and it was unresolved until the twentieth century. Arthur Eddington, in The Internal Constitution of the Stars (1926) ^{[Eddington 1926]}, applied quantum ideas and the virial theorem to argue that stellar cores must reach millions of degrees and that the only energy source adequate to the observed luminosities was the conversion of mass to energy by subatomic processes, though the specific reactions were not yet known.

The missing ingredient was nuclear physics. George Gamow's 1928 theory of quantum tunnelling explained how two positively charged protons could fuse despite their electrical repulsion, and Robert Atkinson and Fritz Houtermans showed that tunnelling opened a viable stellar energy source. The decisive identification came with Hans Bethe, whose 1939 paper "Energy Production in Stars" ^{[Bethe 1939]} laid out both the proton-proton chain and the CNO cycle and matched them to the temperatures and masses of real stars, work recognised by the Nobel Prize in 1967. Bethe's analysis showed that the two engines operate in different stellar mass ranges, explaining at a stroke why the main sequence is a one-parameter family ordered by mass.

The solar neutrino problem then subjected the whole theory to a thirty-year audit. John Bahcall's Standard Solar Model predicted a neutrino flux that Ray Davis's Homestake experiment, running from 1967, consistently failed to match, detecting only about a third of the electron neutrinos expected ^{[Bahcall et al. 2001]}. The resolution by the Sudbury Neutrino Observatory in 2001-2002 transformed the discrepancy from a failure of stellar physics into a discovery about neutrinos themselves: they oscillate among flavours and carry mass. The stellar model survived its longest test intact.

The deeper philosophical lesson of stellar structure is that stars are machines for converting mass into the binding energy of heavier nuclei, and thereby into light. The same equations that describe this conversion also describe the conditions under which carbon, oxygen, and the rest of the periodic table are forged. To understand the interior of a star is to understand where the energy that warms a planet and the atoms that build a body both originate, binding nuclear physics, geology, and cosmology into a single quantitative science.

Bibliography Master

1. Eddington, A. S. (1926). The Internal Constitution of the Stars. Cambridge University Press. The foundational theoretical treatment of hydrostatic equilibrium and stellar interiors.

2. Bethe, H. A. (1939). "Energy Production in Stars." Physical Review 55, 434-456. Identification of the proton-proton chain and the CNO cycle as stellar energy sources.

3. Chandrasekhar, S. (1939). An Introduction to the Study of Stellar Structure. University of Chicago Press. The polytropic and Lane-Emden framework in full generality.

4. Schwarzschild, M. (1958). Structure and Evolution of the Stars. Princeton University Press. The first modern textbook tying structure equations to nuclear burning and evolution.

5. Kippenhahn, R. & Weigert, A. (1990). Stellar Structure and Evolution. Springer. The standard graduate reference for the four structure equations, opacities, and nuclear generation.

6. Clayton, D. D. (1983). Principles of Stellar Evolution and Nucleosynthesis. University of Chicago Press. Nuclear reaction rates and energy generation in detail.

7. Hansen, C. J., Kawaler, S. D. & Trimble, V. (2004). Stellar Interiors (2nd ed.). Springer. Equation of state, opacity, and asteroseismology.

8. Carroll, B. W. & Ostlie, D. A. (2017). An Introduction to Modern Astrophysics (2nd ed.). Cambridge. Chapters 12-16 give the undergraduate treatment of stellar structure and nuclear burning.

9. Bahcall, J. N., Pinsonneault, M. H. & Basu, S. (2001). "Solar Models: Current Epoch and Time Dependences, Neutrinos, and Helioseismological Properties." Astrophysical Journal 555, 990-1012. The Standard Solar Model and its neutrino predictions.

10. Asplund, M., Grevesse, N., Sauval, A. J. & Scott, P. (2009). "The Chemical Composition of the Sun." Annual Review of Astronomy and Astrophysics 47, 481-522. The revised solar abundances and the helioseismology conflict.

11. Paxton, B. et al. (2011). "Modules for Experiments in Stellar Astrophysics (MESA)." Astrophysical Journal Supplement 192, 3. The modern stellar evolution code.

12. Bahcall, J. N. & Davis, R. (2000). "The Solar Neutrino Problem." Scientific American 282, 30-36. The long history of the problem before its SNO resolution.