Modern Physics · Advanced Topics

Astrophysics & Stellar Structure

Stars are self-gravitating balls of plasma in hydrostatic equilibrium, powered by nuclear fusion. Understanding stellar structure reveals the origin of chemical elements, the life cycle of stars, and the endpoints — white dwarfs, neutron stars, black holes.

PrerequisitesGravitation(Ch.7)Nuclearphysics(Ch.NP)Statisticalmechanics(Ch.S)GeneralrGravitation (Ch. 7) \cdot Nuclear physics (Ch. NP) \cdot Statistical mechanics (Ch. S) \cdot General relativity (Ch. GR)
Learning Goals
  • Derive the hydrostatic equilibrium equation and use it to estimate central pressures and temperatures in stars.
  • Apply the virial theorem to compare the Kelvin-Helmholtz timescale with the nuclear burning timescale of the Sun.
  • ExplainthemassluminosityrelationLM3.5andcomputemainsequencelifetimesasafuExplain the mass-luminosity relation L \propto M^3.5 and compute main-sequence lifetimes as a function of stellar mass.
  • Describe stellar endpoints (white dwarfs, neutron stars, black holes) and the mass thresholds that determine them.
  • Use the Chandrasekhar limit and the distance ladder (parallax, Cepheids, Type Ia SNe) to connect stellar physics to cosmological measurements.

AS.1 Hydrostatic Equilibrium

A star's structure is determined by the balance between gravity (inward) and pressure gradient (outward). For a spherical shell at radius r containing mass M(r):

dP/dr=GM(r)ρ/r2(hydrostaticequilibrium)dP/dr = -G M(r) \rho / r^{2} \qquad (hydrostatic equilibrium)(AS.1)
dM/dr=4\pir2ρ(masscontinuity)dM/dr = 4\pir^{2} \rho \qquad (mass continuity)(AS.2)

These two equations govern the pressure and density profiles. The virial theoremfor a self-gravitating system states: 2K + U_grav = 0, so E_total = −K = U_grav/2. Contracting a star releases gravitational energy, half radiated away, half heating the gas.

The Kelvin-Helmholtz timescale (gravitational contraction time):

τKH=GM2/(RL)1.5×107yr(Sun)(KelvinHelmholtztimescale)\tau_KH = GM^{2}/(RL) \approx 1.5\times10^{7} yr (Sun) \qquad (Kelvin-Helmholtz timescale)(AS.3)

The Sun's actual age is 4.6×10⁹ yr — 300× longer. This proved that nuclear burning, not gravity, powers the Sun (Thomson/Kelvin's classical estimate was wrong).

AS.2 Energy Transport

Heat flows outward through a star by three mechanisms:

Radiative transport: photons random-walk outward. The opacity κ (m²/kg) determines how far they travel. The radiative temperature gradient:

dT/drrad=(3κρL)/(16\piacr2T3)(radiativegradient)dT/dr|_rad = -(3\kappa\rho L)/(16\piacr^{2}T^{3}) \qquad (radiative gradient)(AS.4)

where a = 4σ/c is the radiation constant and L(r) is the luminosity at r.

Convection: when |dT/dr|_rad exceeds the adiabatic gradient, the gas is convectively unstable (Schwarzschild criterion). Convection is efficient — nearly adiabatic in stellar interiors. The outer ∼30% of the Sun (by radius) is convective; the core is radiative.

AS.3 Nuclear Burning in Stars

The proton-proton (pp) chain dominates in solar-mass stars (T_core ≈ 1.5×10⁷ K):

pp I: 4 p → ⁴He + 2e⁺ + 2νe + 26.7 MeV. Rate ∝ T⁴ (steep T dependence). CNO cycle dominates for M > 1.5 M☉ (T_core > 1.7×10⁷ K): ε ∝ T¹⁶ to T²⁰.

εppρT4εCNOρT17(nuclearenergygenerationrates)\varepsilon_pp \propto \rho T^{4} \qquad \varepsilon_CNO \propto \rho T^{17} \qquad (nuclear energy generation rates)(AS.5)

The steep temperature dependence is a thermostat: if T rises, ε rises, pressure rises, star expands, T falls — negative feedback keeps stars stable. Main sequence lifetime:

τMS0.1×fH×Mc2/L1010yr×(M/M)/(L/L)(mainsequencelifetime)\tau_MS \approx 0.1 \times f_{H} \times M c^{2}/L \approx 10^{10} yr \times (M/M\odot)/(L/L\odot) \qquad (main sequence lifetime)(AS.6)

where f_H ≈ 10% of the mass is burned (core hydrogen fraction). Since L ∝ M³ to M⁴, massive stars (10 M☉) live only ∼30 Myr; low-mass stars (0.1 M☉) live > 100 Gyr.

Example AS.1The Mass-Luminosity Relation

DeriveLM3fromstellarstructureequationsandexplainthemassluminosityrelationDerive L \propto M^{3} from stellar structure equations and explain the mass-luminosity relation.

Pressure scale:Central pressure from hydrostatic equilibrium: P_{c} ∼ GM^{2}/R^{4}. For ideal gas: P = \rhok_BT/(\mum_H),soTcGM\mumH/(kBRH), so T_{c} ∼ GM\mum_H/(k_{B} R.
Opacity:Electronscatteringopacityκ0.2(1+X)m2/kgconst.Theradiativeluminosity:L=(4\picElectron scattering opacity \kappa \approx 0.2(1+X) m^{2}/kg \approx const. The radiative luminosity: L = (4\picGM)/(κ)×βwhereβ=Pgas/Ptotal(EddingtonluminositylimitGM)/(\kappa) \times \beta where \beta = P_{gas}/P_{total} (Eddington luminosity limit.
Eddington limit:LEdd=4\picGM/κ3.2×104(M/M)L.StarsmusthaveL<LEddorradiationpressureblowsL_{Edd} = 4\picGM/\kappa \approx 3.2\times10^{4} (M/M\odot) L\odot. Stars must have L < L_{Edd} or radiation pressure blows the envelope away.
Mass-luminosity:Combiningpressureandopacityarguments:LM3forintermediatemassstars.ObservationaCombining pressure and opacity arguments: L \propto M^{3} for intermediate-mass stars. Observationally:LL(M/M)3.5(empiricalfit).A10Mstar:L103.5L3000Llly: L \approx L\odot(M/M\odot)^3.5 (empirical fit). A 10M\odot star: L \approx 10^3.5 L\odot \approx 3000 L\odot.
Main sequence:SinceLM3.5andτM/LM2.5:a10Mstarlives316×shorterthantheSun:101Since L \propto M^3.5 and \tau \propto M/L \propto M^{-2.5}: a 10M\odot star lives ∼316\times shorter than the Sun: 10^10/31630Myr.ConsistentwiththeageoftheOrionOBassociation0/316 \approx 30 Myr. Consistent with the age of the Orion OB association.

AS.4 Stellar Evolution and Death

After main sequence, the hydrogen shell burns, and the core contracts as the envelope expands → Red Giant. Helium ignites (triple-alpha process: 3 ⁴He → ¹²C + 7.27 MeV), then carbon, oxygen in more massive stars.

Stellar endpoints by mass:

M < 8 M☉: planetary nebula + white dwarf (supported by electron degeneracy pressure). Chandrasekhar limit: M_Ch = 1.44 M☉. Above this, no stable white dwarf exists.

8 M☉ < M < ~20 M☉: core-collapse supernova → neutron star(supported by neutron degeneracy + repulsive nuclear force). Density ∼ nuclear: ρ ≈ 5×10¹⁷ kg/m³. Tolman-Oppenheimer-Volkoff limit ≈ 2–3 M☉.

M > ~20 M☉: collapse to black hole. Event horizon at r_s = 2GM/c². For 10 M☉ BH: r_s ≈ 30 km.

Theorem AS.1Chandrasekhar Limit
White dwarfs are supported by electron degeneracy pressure. For non-relativistic electrons: P_ρ5/3.Themassradiusrelation:RM1/3moremassiveWDsaresmaller.AsMinc\propto \rho^{5/3}. The mass-radius relation: R \propto M^{-1/3} — more massive WDs are smaller. As M increases,electronsbecomerelativistic(Pdegρ(4/3reases, electrons become relativistic (P_{deg} \propto \rho^(4/3:MCh=(5.87/\mue2)(\hbarc/G)3/2mH21.44MM_{Ch} = (5.87/\mue^{2}) (\hbarc/G)^{3/2} m_{H}^{-2} \approx 1.44 M\odotwhere\mueisthemeanmolecularweightperelectron(2forC/O).AtM=MCh,theradiusgwhere \mue is the mean molecular weight per electron (\approx2 for C/O). At M = M_{Ch}, the radius goes to zero — no stable configuration. Type Ia supernovae occur when a white dwarf accretes to this limit and undergoes thermonuclear explosion — they are standard candles (discovered dark energy, Nobel 2011).

AS.5 Cosmological Distance Ladder

Measuring cosmic distances requires a hierarchy of methods:

Parallax (up to ∼1 kpc, Gaia satellite): d = 1/p arcsec (parsec).Cepheids: period-luminosity relation (Leavitt 1912): log L = a log P + b, P = 1–100 days. Used to ∼100 Mpc (HST Key Project, H₀ calibration).Type Ia supernovae: standard candles (M_V ≈ −19.3), usable to z ∼ 2. Used to discover accelerating expansion (dark energy).

The Hubble constant H₀ = v_rec/d (recession velocity per distance) measures the expansion rate. Current tension: H₀ = 73.0 ± 1.0 km/s/Mpc (local, Cepheids/SNe Ia) vs H₀ = 67.4 ± 0.5 km/s/Mpc (CMB/Planck). 5σ discrepancy — possibly new physics.

Definition AS.1Common Traps
  • Luminosity and brightness are different: observed flux falls with distance squared.
  • Magnitude is logarithmic: smaller magnitudes are brighter.
  • Hydrostatic balance is dynamic equilibrium: gravity is balanced by pressure gradients, not absent.
  • Stellar lifetimes depend strongly on mass: massive stars burn brighter and die sooner.
Exercises — AS.1–AS.5 Astrophysics
1.
Calculate the number of pp-chain reactions per second in the Sun, the mass converted to energy per year, and the neutrino flux at Earth.
reactions/s
Straightforward
2.
EstimatetheradiusandmeandensityofawhitedwarfwithmassM=1.2M(composedofC/Estimate the radius and mean density of a white dwarf with mass M = 1.2 M\odot (composed of C/O,\mue=2).HowdoesthiscomparetotheChandrasekharlimitO, \mue = 2). How does this compare to the Chandrasekhar limit?
kg/m³
Intermediate
3.Calculate the surface gravity,escapevelocity,andgravitationalredshiftatthesurfaceofaneutronstar(M=gravity, escape velocity, and gravitational redshift at the surface of a neutron star (M =stimate its rotation period.
Intermediate
4.Quantify the Hubble tension between local (Cepheid/SNeIa) and CMB measurements. What are the main proposed resolutions — systematic or new physics?
Challenging
Key Takeaways
  • Hydrostaticequilibrium:dP/dr=GMρ/r2.Starsliveonthermalnuclearenergy,notgravitHydrostatic equilibrium: dP/dr = -GM\rho/r^{2}. Stars live on thermal nuclear energy, not gravity (KH timescale too short).
  • Massluminosity:LM3.5.Mainsequencelifetime:τM2.5.MassivestarsdieyoungMass-luminosity: L \propto M^3.5. Main sequence lifetime: \tau \propto M^{-2.5}. Massive stars die young.
  • Nuclearreactions:ppchain(T4)forsolartype;CNOcycle(T17)dominatesabove1.5MNuclear reactions: pp chain (T^{4}) for solar-type; CNO cycle (T^{17}) dominates above 1.5 M\odot.
  • Stellar endpoints: WD (M < 8M☉, Chandrasekhar 1.44M☉), NS (8–20M☉, TOV ~2.5M☉), BH (>20M☉).
  • WDsupportedbyelectrondegeneracy(Pρ5/3ρ4/3nearMCh).RM(1/3WD supported by electron degeneracy (P \propto \rho^5/3 \to \rho^4/3 near M_{Ch}). R \propto M^(-1/3.
  • HubbleconstantH0:5σtensionbetweenlocal(73)andCMB(67.4)measurementsunresolvedHubble constant H_{0}: 5\sigma tension between local (73) and CMB (67.4) measurements — unresolved.