Stellar Interior Structure — How Stars Support Their Own Weight

1. Hydrostatic Equilibrium

Consider a thin spherical shell of stellar material at radius r, with thickness dr, density ρ(r), and area 4πr². Its weight (gravitational force) is:

where M(r) = ∫₀ʳ 4πr'² ρ(r') dr' is the mass enclosed within radius r. The outward pressure force on this shell is the difference in pressure across it:

Setting the net force to zero (static equilibrium) gives the fundamental equation of hydrostatic equilibrium:

This single equation is the cornerstone of stellar structure theory. It says that pressure must decrease outward — the interior of a star is under enormous pressure precisely because it must support the weight of everything above it. At the Sun's centre, the pressure is approximately 2.5 × 10¹⁶ Pa — 250 billion times atmospheric pressure — and the temperature is about 1.5 × 10⁷ K.

The Virial Theorem

A powerful consequence of hydrostatic equilibrium is the virial theorem. Integrating the hydrostatic equation over the entire star yields a relationship between the total thermal energy E_th and the gravitational potential energy E_grav:

The virial theorem has a counterintuitive consequence: as a proto-star contracts and loses energy (E_total decreases), half the released gravitational energy heats the interior, and the other half is radiated away. Stars become hotter as they lose energy — they have a negative heat capacity.

2. The Four Equations of Stellar Structure

A complete description of a spherically symmetric star in thermal equilibrium requires four coupled first-order ordinary differential equations, using radius r as the independent variable:

Here L(r) is the luminosity (energy flux through a sphere of radius r), ε(ρ,T) is the nuclear energy generation rate per unit mass, κ is the opacity (cross-section per unit mass for photon absorption), a is the radiation constant, c is the speed of light, and γ is the adiabatic index. The energy transport equation takes one of two forms depending on whether the local conditions favour radiative or convective transport (discussed in Section 6).

To close this system, we need an equation of state P(ρ, T) and prescriptions for κ(ρ, T) and ε(ρ, T) from atomic physics and nuclear reaction rates. The boundary conditions are P = 0, L = 0, ρ = 0 at r = R (stellar surface), and M = 0, L = 0 at r = 0 (centre). Given the total mass M_star, integrating these equations inward from the surface (shooting method) or outward from the centre determines the entire stellar structure.

3. The Lane-Emden Equation and Polytropes

Before computers, astronomers needed analytic solutions to the stellar structure equations. A powerful simplification is the polytropic equation of state:

where K is a constant and n is the polytropic index. Substituting into the hydrostatic equation combined with Poisson's equation (∇²Φ = 4πGρ) and introducing dimensionless variables:

yields the Lane-Emden equation:

The stellar surface is defined by the first zero of θ: θ(ξ₁) = 0. Analytic solutions exist for three special cases:

• n = 0 (incompressible fluid): θ = 1 − ξ²/6, ξ₁ = √6 ≈ 2.449. Uniform density throughout.
• n = 1: θ = sin(ξ)/ξ, ξ₁ = π. The density profile is a sinc function.
• n = 5: θ = (1 + ξ²/3)^(−1/2), ξ₁ → ∞. The star has infinite radius but finite mass.

Other values of n require numerical integration. The physically most important polytropes are:

• n = 3/2 (γ = 5/3): ideal monatomic gas, adiabatic; describes fully convective stars and non-relativistic white dwarfs.
• n = 3 (γ = 4/3): Eddington's "standard model" for massive stars where radiation pressure dominates; also describes relativistic white dwarfs and marks the Chandrasekhar limit.
• n = 3: The Chandrasekhar mass M_Ch = 5.83 M_☉/μ_e² follows from the condition that the n = 3 polytrope has finite mass but zero pressure at the surface.

4. Nuclear Energy Generation

Stars are powered by nuclear fusion in their cores. The energy generation rate ε (erg per gram per second) depends sensitively on temperature and density because fusion requires nuclei to tunnel through the Coulomb barrier. For the pp chain (the dominant process in the Sun):

The steep temperature dependence (T⁴ for the pp chain, T¹⁸ for the CNO cycle in more massive stars) is the key to stellar self-regulation. If the core overheats slightly, nuclear burning accelerates, the core expands, cools, and burning slows back down — a negative feedback loop that keeps the star stable for billions of years.

The CNO Cycle

In stars more massive than about 1.3 M_☉ (with core temperatures above ~2 × 10⁷ K), the CNO cycle dominates. Carbon, nitrogen, and oxygen act as catalysts, producing the same net reaction (4H → He) but with ε_CNO ∝ ρ X X_CNO T¹⁸. The much steeper temperature dependence means CNO stars have very compact, hot cores and tend to be convective in their centres (because the steep temperature gradient drives convection).

5. Opacity and Radiative Transport

Radiation carries energy outward through the stellar interior as a random walk of photons, each absorbed and re-emitted countless times. The mean free path of a photon is ℓ = 1/(κρ), where κ is the Rosseland mean opacity (averaged over frequency, weighted by the temperature derivative of the Planck function):

The main sources of opacity in stellar interiors are:

• Electron scattering (Thomson): κ_es = 0.20(1+X) cm²/g, independent of T. Dominates in hot, highly ionised interiors of massive stars.
• Free-free (bremsstrahlung) absorption: κ_ff ∝ ρ T^(−7/2) (Kramers' opacity). Dominates in hot stars below electron-scattering temperature.
• Bound-free (photoionisation): κ_bf ∝ ρ T^(−7/2). Important in cooler regions where atoms are partially ionised.
• Bound-bound (line) absorption: Dominates in the photosphere and cool outer layers.
• H⁻ opacity: In stellar envelopes below ~8000 K, the negative hydrogen ion is the dominant opacity source.

The radiative diffusion equation for energy flux gives the temperature gradient required to transport luminosity L outward by radiation alone:

High opacity κ means a steep temperature gradient is needed to push the required luminosity through. When this gradient becomes too steep, the fluid becomes convectively unstable.

6. Convective Energy Transport and the Schwarzschild Criterion

Convection sets in when the actual temperature gradient exceeds what a rising fluid blob can sustain adiabatically. The Schwarzschild criterion for convective instability is:

When this condition is met, a fluid element displaced upward (into lower pressure and density) expands faster than the surrounding medium and finds itself denser — it then sinks back, driving a large-scale overturning circulation. Convection is a very efficient energy transporter: the actual temperature gradient in a convective zone stays very close to ∇_ad, far shallower than ∇_rad.

Mixing-Length Theory

The most widely used model of stellar convection is mixing-length theory (MLT), introduced by Ludwig Biermann and formalised by Erika Böhm-Vitense in 1958. It treats convection as blobs of fluid rising and falling a characteristic distance called the mixing length ℓ = α H_P, where H_P = −dr/d(ln P) is the pressure scale height and α ≈ 1.5–2.0 is a dimensionless free parameter calibrated against the Sun. MLT gives the convective flux:

where δ = −(∂ ln ρ/∂ ln T)_P and c_P is the specific heat. Despite its simplicity, MLT reproduces stellar evolution tracks, pulsation periods, and helioseismological sound-speed profiles with impressive accuracy — though modern 3D radiation-hydrodynamics simulations are beginning to reveal its limitations.

7. The Solar Interior: A Case Study

The Sun (mass M_☉ = 1.99 × 10³⁰ kg, radius R_☉ = 6.96 × 10⁸ m, luminosity L_☉ = 3.85 × 10²⁶ W) consists of three main interior regions:

• Core (0–0.25 R_☉): Temperature ~1.5 × 10⁷ K, density ~150 g/cm³. The pp chain converts 600 million tonnes of hydrogen to helium every second. Energy is transported outward by radiation.
• Radiative zone (0.25–0.72 R_☉): Temperature falls from ~7 × 10⁶ K to ~2 × 10⁶ K. A photon's random walk from core to the base of the convection zone takes approximately 170,000 years. The composition is nearly uniform here due to the absence of mixing.
• Convective envelope (0.72–1.0 R_☉): Above the tachocline (the thin shear layer at 0.72 R_☉), opacity rises sharply due to partial ionisation of helium. The Schwarzschild criterion is violated and vigorous convection develops, transporting energy in ~1 month (not 170,000 years). The convective cells visible on the solar surface are called granules (1 Mm across, 10 min lifetime) and supergranules (30 Mm, 1 day).

8. Helioseismology: Seeing Inside the Sun

The Sun oscillates in millions of normal modes — acoustic waves (p-modes) and surface gravity waves (f-modes) trapped in resonant cavities. The frequencies of these modes, observed as Doppler shifts in the solar photosphere by instruments like GOLF and MDI/HMI, encode information about the sound speed c_s = √(γ P/ρ) at every depth:

By inverting the set of observed mode frequencies (solving a Sturm-Liouville eigenvalue problem), helioseismologists can reconstruct the sound-speed profile to better than 0.1% accuracy throughout most of the solar interior. The results confirmed the radiative-convective boundary at exactly 0.713 R_☉, the helium abundance in the convection zone (Y = 0.249), and — crucially — the absence of a rapidly rotating core that some theories had predicted.

The same technique applied to other stars is called asteroseismology, now a major observational programme of the NASA Kepler and TESS space telescopes. Asteroseismology has determined precise stellar ages, masses, radii, and core rotation rates for tens of thousands of stars, revolutionising our understanding of Galactic stellar populations.