Hydrostatic equilibrium is the fundamental condition that describes the mechanical stability of a star, where the inward pull of gravity is precisely balanced by an outward pressure gradient. This balance ensures that the star remains static and does not undergo rapid expansion or collapse. Mathematically, for a spherically symmetric star, this is expressed by the hydrostatic equilibrium equation: $$\frac{dP}{dr} = -G \frac{M_r\rho}{r^2} = -\rho g$$ where $P$ is pressure, $r$ is the radial distance from the center, $M_r$ is the interior mass within radius $r$, $\rho$ is the local density, and $g$ is the local acceleration of gravity. This equation dictates that pressure must increase toward the center of the star to support the weight of the overlying layers.
Pressure Equation of State
A pressure equation of state relates the internal pressure of stellar material to its density, temperature, and chemical composition. The general pressure integral derived from the momentum transferred by particles (massive or massless) is: $$P = \frac{1}{3} \int_0^\infty p n_p v_p dp$$ where $p$ is momentum, $v$ is velocity, and $n_p$ is the number density of particles per momentum interval. In most stellar interiors, the total pressure is the sum of gas pressure ($P_g$) and radiation pressure ($P_{rad}$): $$P_t = \frac{\rho kT}{\mu m_H} + \frac{1}{3} a T^4$$ where $k$ is Boltzmann’s constant, $\mu$ is the mean molecular weight, $m_H$ is the mass of a hydrogen atom, and $a$ is the radiation constant.
Ideal Gas Law
The ideal gas law applies to non-degenerate, non-relativistic stellar gases where particles interact only through perfectly elastic collisions. Derived using the Maxwell-Boltzmann velocity distribution, it is expressed as: $$P_g = n k T = \frac{\rho kT}{\mu m_H}$$ where $n$ is the particle number density. The mean molecular weight ($\mu$) depends on the ionization state and composition of the gas. For a completely ionized gas, it is approximated by: $$\frac{1}{\mu_i} \approx 2X + \frac{3}{4}Y + \frac{1}{2}Z$$ where $X, Y, \text{ and } Z$ are the mass fractions of hydrogen, helium, and metals, respectively.
Fermi-Dirac and Bose-Einstein Statistics
When classical assumptions fail, quantum statistics are required to describe particle distributions:
Fermi-Dirac Statistics: Apply to fermions (e.g., electrons, protons, neutrons), which obey the Pauli exclusion principle. In extremely dense environments like white dwarfs, electron degeneracy pressure becomes dominant.
Non-relativistic degeneracy: $P \approx \frac{(3\pi^2)^{2/3}}{5} \frac{\hbar^2}{m_e} \left( \frac{Z}{A} \frac{\rho}{m_H} \right)^{5/3}$.
Relativistic degeneracy: $P \approx \frac{(3\pi^2)^{1/3}}{4} \hbar c \left( \frac{Z}{A} \frac{\rho}{m_H} \right)^{4/3}$.
Bose-Einstein Statistics: Apply to bosons (e.g., photons), which do not obey the exclusion principle. This leads to the formulation of radiation pressure: $$P_{rad} = \frac{1}{3} a T^4$$ which can dominate in very massive, hot stars.
Stars are powered by two primary energy sources: gravitation and nuclear fusion.
Gravitational Energy: Released during contraction, with the total potential energy approximated as $U_g \approx -\frac{3}{5} \frac{GM^2}{R}$. According to the virial theorem, half is radiated and half heats the interior.
Nuclear Energy: Principally the fusion of light nuclei into heavier ones, releasing energy based on $E=mc^2$. The nuclear energy generation rate ($\epsilon$) is often modeled as a power law:
$$\epsilon = \epsilon_0 X_i X_x \rho^\alpha T^\beta$$
where $\beta$ varies significantly by reaction type (e.g., $\beta \approx 4$ for the pp chain and $\beta \approx 20$ for the CNO cycle).
Timescales:
Three primary timescales characterize stellar life:
1. Dynamic (Free-fall) Timescale: The time for a cloud to collapse under gravity if pressure is removed: $t_{ff} = \sqrt{\frac{3\pi}{32} \frac{1}{G \rho_0}}$.
2. Thermal (Kelvin-Helmholtz) Timescale: The time to radiate away a star’s total gravitational energy: $t_{KH} \approx \frac{3}{10} \frac{GM^2}{RL}$.
3. Nuclear Timescale: The time to exhaust nuclear fuel: $t_{nuc} \approx \frac{E_{nuc}}{L}$. For the Sun, $t_{nuc} \approx 10^{10}$ years.
Quantum Tunneling:
Classical physics suggests stellar cores are too cool for nuclei to overcome the Coulomb barrier through thermal motion alone. Quantum mechanical tunneling allows particles to penetrate this barrier even when their kinetic energy is insufficient. The probability of tunneling is exponential: $$P_{tunnel} \propto e^{-b E^{-1/2}}$$ where $b$ is a constant related to particle charges and masses. This probability combined with the Maxwell-Boltzmann distribution defines the Gamow peak, the narrow energy range where most nuclear reactions occur.
Nucleosynthesis:
Nucleosynthesis is the sequence of nuclear reactions that transform elements:
Hydrogen Burning: Occurs via the proton-proton (pp) chains (dominant in low-mass stars) or the CNO cycle (dominant in massive stars).
Helium Burning: Occurs via the triple-alpha process ($3\alpha \to ^{12}C$) at temperatures $\sim 10^8$ K.
Advanced Burning: Successive stages (carbon, oxygen, neon, and silicon burning) produce elements up to the iron peak. Elements heavier than iron are produced via the s-process (slow neutron capture) or r-process (rapid neutron capture).
Energy is transported from the core to the surface via radiation, convection, or conduction. Radiative Transport: Driven by the radiation pressure gradient: $$\frac{dT}{dr} = -\frac{3}{4ac} \frac{\kappa \rho}{T^3} \frac{L_r}{4\pi r^2}$$ where $\kappa$ is the opacity.
Convection: Occurs when the temperature gradient is superadiabatic. The Schwarzschild criterion for convection (neglecting composition changes) is: $$\left| \frac{dT}{dr} \right|_{act} > \left| \frac{dT}{dr} \right|_{ad}$$ where the adiabatic gradient for an ideal gas is $\frac{dT}{dr}|_{ad} = -\left(1 - \frac{1}{\gamma}\right) \frac{\mu m_H}{k} \frac{GM_r}{r^2}$.
Thermodynamics: The first law ($dU = dQ - dW$) relates internal energy changes to heat and work. For an adiabatic process, $P V^\gamma = \text{constant}$.
Stellar models are constructed by numerically solving the four fundamental differential equations (hydrostatic equilibrium, mass conservation, energy generation, and energy transport).
Numerical Modeling: The star is divided into concentric shells (zones), and differential equations are replaced by difference equations. Solutions require matching boundary conditions at the center ($M_r \to 0, L_r \to 0$) and surface ($P \to 0, T \to 0$).
Vogt-Russell Theorem: States that a star’s mass and composition uniquely determine its structure and evolution.
Polytropes: Simplified analytical models where pressure is a power of density ($P = K\rho^{(n+1)/n}$); they are solved using the Lane-Emden equation:
$$\frac{1}{\xi^2} \frac{d}{d\xi} \left( \xi^2 \frac{dD_n}{d\xi} \right) = -D_n^n$$
where $\xi$ is a dimensionless radius and $D_n$ is a dimensionless density function.