4.4 Equations of stellar structure

Fundamental Equations of Stellar Structure

Stars hold together because of a precise balance between competing physical processes: gravity pulling inward, pressure pushing outward, energy being generated and transported, and mass distributed throughout. The four equations of stellar structure capture this balance mathematically. Together with boundary conditions and an equation of state, they let you build a complete model of a star's interior from core to surface.

The Four Structure Equations

Each equation is a first-order ordinary differential equation in radius $r$, describing how a physical quantity changes as you move outward through the star.

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

$\frac{dP}{dr} = -\frac{G M(r) \rho(r)}{r^2}$

• $P$ is pressure, $G$ the gravitational constant, $M(r)$ the mass enclosed within radius $r$, and $\rho(r)$ the local density.
• The negative sign tells you pressure decreases as you move outward. If this balance breaks, the star either collapses or expands.

Mass Conservation relates how enclosed mass grows with radius:

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

This is just saying that the mass in a thin shell of thickness $dr$ equals its volume ($4\pi r^2 dr$) times its density. It connects the density profile to the total mass distribution.

Energy Generation describes how luminosity builds up through the star:

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

• $L(r)$ is the luminosity passing through radius $r$, and $\epsilon(r)$ is the energy generation rate per unit mass (from nuclear reactions).
• In the core where fusion occurs, $\epsilon$ is large and luminosity climbs steeply. In the envelope, $\epsilon \approx 0$ and $L(r)$ stays roughly constant.

Energy Transport (Radiative Transfer) governs the temperature gradient when energy moves outward by radiation:

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

• $T$ is temperature, $\kappa$ the opacity (how effectively material blocks radiation), $a$ the radiation density constant, and $c$ the speed of light.
• High opacity or high luminosity produces a steeper temperature gradient. If the gradient becomes too steep, convection takes over as the dominant transport mechanism, and you replace this equation with the adiabatic gradient condition instead.

These four equations contain seven unknowns: $P$, $T$, $\rho$, $M$, $L$, $\kappa$, and $\epsilon$. To close the system, you need three additional relations: an equation of state (linking $P$, $\rho$, and $T$), an opacity law ($\kappa$ as a function of $\rho$, $T$, and composition), and nuclear reaction rates ($\epsilon$ as a function of $\rho$, $T$, and composition).

Boundary Conditions

You can't solve differential equations without boundary conditions. Stellar structure requires them at both the center and the surface.

At the center ($r = 0$):

• $M(0) = 0$ — no mass is enclosed at the very center
• $L(0) = 0$ — no luminosity is generated interior to a point

At the surface ($r = R$):

• $P(R) \approx 0$ — pressure drops to (approximately) zero at the photosphere
• $T(R) = T_{\text{eff}}$ — the surface temperature matches the effective temperature

In practice, the surface conditions aren't exactly zero; they connect to an atmospheric model. But for most stellar structure calculations, these approximations work well.

The boundary conditions, combined with initial parameters like total mass and chemical composition (often parameterized by hydrogen fraction $X$, helium fraction $Y$, and metallicity $Z$), uniquely determine the solution. This is why stars of the same mass and composition have the same structure: the Vogt-Russell theorem.

Modeling Stellar Interiors

Solving these equations analytically is only possible for highly idealized cases (like polytropic models). Real stellar models are built numerically:

1. Divide the star into many concentric shells.
2. Assign initial guesses for $P$, $T$, $\rho$, $M$, and $L$ at each shell.
3. Integrate the structure equations inward from the surface and outward from the center.
4. Adjust parameters until the two solutions match smoothly at a fitting point (typically somewhere in the interior). This is the Henyey method, or a variant of it.

The resulting model predicts observable properties: effective temperature, luminosity, radius, and surface composition. By computing a sequence of models at successive time steps (updating composition as nuclear burning proceeds), you generate evolutionary tracks on the Hertzsprung-Russell diagram, showing how a star's temperature and luminosity change over its lifetime.

Internal structure profiles from these models reveal how temperature, density, pressure, and composition vary from core to surface. For example, a solar model shows core temperatures near $1.5 \times 10^7$ K and core densities around $150 \, \text{g/cm}^3$, dropping by orders of magnitude toward the surface.

Limitations of Stellar Structure Models

These models are powerful but rest on several simplifying assumptions. Knowing where they break down matters for interpreting results.

• Spherical symmetry is assumed throughout. Rotation and magnetic fields break this symmetry. For slowly rotating stars like the Sun, the error is small. For rapidly rotating stars (like many O and B stars), centrifugal distortion becomes significant and 1D models lose accuracy.
• Mixing length theory (MLT) is used to model convection. It introduces a free parameter, the mixing length $\alpha$, typically calibrated to the Sun ($\alpha \approx 1.5\text{–}2.0$). MLT is a rough approximation of an inherently 3D, turbulent process, and it handles convective boundaries poorly.
• Local thermodynamic equilibrium (LTE) assumes radiation and matter are in equilibrium locally. This holds well in deep interiors but breaks down in stellar atmospheres where photons can escape freely (the mean free path becomes large).
• Mass loss is often neglected or treated simply. For massive stars (Wolf-Rayet stars can lose several solar masses over their lifetimes) and evolved giants with strong stellar winds, this is a serious omission.
• Static (time-independent) structure is assumed at each timestep. The star adjusts instantaneously to changes. This fails during rapid evolutionary phases like thermal pulses on the asymptotic giant branch or the final stages before core collapse.
• 1D modeling cannot capture 3D effects like turbulence, convective overshooting, and rotational mixing. Modern 3D hydrodynamic simulations address some of these issues but are computationally expensive and limited to short timescales.
• Nuclear reaction rate uncertainties, particularly for reactions like ${}^{12}\text{C}(\alpha, \gamma){}^{16}\text{O}$, propagate into predictions of stellar lifetimes, nucleosynthesis yields, and the final fate of massive stars.
• Opacity tables and equations of state are computed from atomic physics and have their own uncertainties. These are especially challenging in extreme conditions like white dwarf interiors (electron degeneracy) or partially ionized stellar envelopes.

Despite these limitations, the four structure equations remain the foundation of stellar astrophysics. Most of what we know about stellar evolution, from the main sequence through to compact remnants, comes from solving these equations under various assumptions and refining the input physics.