The Balancing Act — Hydrostatic Equilibrium

Complete lesson

Concept Throughline

After completing this reading, you should be able to:

Explain
why a star needs a pressure gradient (not just pressure) to resist gravitational collapse, and identify thermal pressure as the primary source.
Distinguish
thermal gas pressure, radiation pressure, and degeneracy pressure as different sources of support inside stars.
Derive
and interpret the equation of hydrostatic equilibrium as a force balance between gravity and the pressure gradient.
Apply
the virial theorem to relate a star’s total thermal energy to its gravitational potential energy.
Estimate
a star’s central temperature from its mass and radius, and explain why gravity alone predicts a core temperature of order $10^7\,\mathrm{K}$ .
Explain
the negative-heat-capacity paradox — why a self-gravitating system gets hotter as it loses energy.
Connect
hydrostatic equilibrium to the dynamical timescale, and explain why departures from balance are corrected almost instantly.

Concept Throughline

Gravity never stops pulling inward. For a star to survive, a force must balance gravity at every radius — not uniformly, but with a strength that increases toward the center, where the weight of the overlying material is greatest. That force comes from a pressure gradient, and the condition that pressure exactly balances gravity at every point inside the star is called hydrostatic equilibrium.

Pressure gradient

The rate at which pressure changes with position, $dP/dr$ . A uniform pressure exerts no net force; only a pressure gradient produces one. Inside a star the gradient points outward (pressure falls with radius) and supports the weight of the overlying gas.

Hydrostatic equilibrium

The condition in which the outward pressure-gradient force exactly balances the inward pull of gravity at every radius, so the gas has no net radial acceleration. It is the foundational force-balance equation of stellar structure.

This force balance is the foundation of stellar structure theory. Combined with the virial theorem and the ideal-gas picture, it lets us estimate a solar core temperature of order $10^7\,\mathrm{K}$ , or about $15\,\mathrm{MK}$ . No nuclear-reaction physics is needed to estimate this temperature scale — only gravity, force balance, energy balance, and the thermal behavior of gas.

Horizontal five-step flowchart for hydrostatic equilibrium with rounded boxes labeled Gravity, Force Balance, Pressure Scale, Gas Pressure, and Temperature Scale, each containing the key equation and connected by arrows. — Figure 1Stellar structure is a reasoning chain: gravity sets the inward pull, hydrostatic equilibrium sets the required pressure gradient, the pressure scale implies a central pressure, and the ideal-gas picture turns that into a core temperature scale.ASTR 201 (generated)

Why the Sun Does Not Collapse

Part 1: Why Doesn’t the Sun Collapse?

A deceptively simple question

In Reading 1 you calculated the Sun’s dynamical timescale: $\tau_{\text{dyn}} \sim 50\,\mathrm{min}$ . If gravity were the only force acting, the Sun would collapse in less than an hour. Yet the Sun has been shining for $4.6\,\mathrm{Gyr}$ — about $4 \times 10^{13}$ dynamical times. Somehow it maintains an almost perfect balance between inward pull and outward support.

Dynamical timescale

The characteristic time for a star to respond mechanically to a force imbalance — roughly the free-fall time, $\tau_{\text{dyn}} \sim 1/\sqrt{G\bar\rho}$ . For the Sun it is about 50 minutes; departures from hydrostatic equilibrium are corrected on this timescale.

What force opposes gravity? The answer is pressure — specifically, the pressure of the hot gas inside the star. But saying “pressure holds the star up” is not quite right, and understanding why requires some care. If the balance fails, even briefly, the gas cannot remain static: some layer must accelerate inward or outward on roughly the dynamical timescale.

Pressure vs. pressure gradient

Two-panel didactic diagram comparing a gas shell with equal pressure on top and bottom to a gas shell with larger pressure below than above, alongside downward gravity arrows and net-force labels. — Figure 2Uniform pressure produces no net support. Only when the lower side pushes harder than the upper side does a shell feel a net outward pressure force.ASTR 201 (generated)

Imagine a thin shell of gas at some radius $r$ inside a star. This shell feels gravity pulling it inward and pressure pushing it from both sides — inward from the gas above and outward from the gas below.

If the pressure were the same everywhere inside the star, the inward and outward pressure forces on the shell would exactly cancel, and gravity would win unopposed. For pressure to resist gravity, the pressure below the shell must exceed the pressure above it. What matters physically is the difference in pressure across the shell:

$F_{P,\text{net}} = \left[P(r) - P(r+dr)\right]A.$

If $P(r+dr) = P(r)$ , then $F_{P,\text{net}} = 0$ no matter how large the common pressure is. A star is supported only when the lower side pushes harder than the upper side. The physical point is that it is the pressure gradient, $dP/dr$ , not the pressure $P$ by itself, that provides the outward force.

Think first

At the bottom of a swimming pool, the water pressure is higher than at the surface. Why? Commit to a guess before reading on, then state the analogy to a star.

What makes the pressure rise with depth in the pool, and what plays the role of 'depth' inside a star?

Hold this prediction — the answer is just below.

Quick check

Check your prediction: what makes the pressure rise with depth in the pool, and what plays the role of “depth” inside a star?

Quick check

Predict the motion of a gas shell if the pressure were extremely large but exactly uniform everywhere inside the star. Be explicit about the pressure forces and gravity.

Types of pressure in stars

To say that “pressure supports a star” raises a deeper question: what physical process creates that pressure? In ordinary stars like the Sun, the answer begins with the random thermal motions of particles. In hotter, more massive stars, trapped radiation also contributes. Much later in the course, when we discuss white dwarfs and stellar remnants, we will meet a third kind of support that comes from quantum mechanics rather than temperature.

Thermal gas pressure: momentum from particle collisions

A gas is made of particles moving in many directions. When those particles strike a surface, they bounce and transfer momentum. Pressure is the rate of momentum transfer per unit area. Faster particles hit harder; more particles per cubic centimeter hit more often. That is why heating a gas or compressing it both raise its pressure. For a thermal gas, this is summarized by the ideal-gas pressure law.

In ordinary stars like the Sun, this thermal gas pressure is the main source of support against gravity: at fixed density a hotter gas pushes harder, and at fixed temperature a denser gas pushes harder.

Symbol	Meaning
$P_{\text{gas}}$	thermal gas pressure, in $\mathrm{dyn\,cm^{-2}}$
$n$	number density, in particles per $\mathrm{cm^3}$
$k_B$	Boltzmann’s constant
$T$	temperature, in $\mathrm{K}$
$\rho$	mass density, in $\mathrm{g\,cm^{-3}}$
$m_p$	proton mass
$\mu$	mean molecular weight
$\mu m_p$	mean mass per particle

Deep Dive: Deep dive — what gas pressure really is

Gas pressure is a macroscopic description of an enormous number of microscopic collisions. A gas contains a vast ensemble of particles moving in many directions; every time they strike a surface, they transfer momentum, and pressure is the cumulative rate of that transfer per unit area.

Temperature is not “how much heat a single particle contains.” It is a macroscopic description of the typical kinetic-energy distribution of the whole ensemble. A hotter gas has particles that typically move faster, so collisions transfer more momentum. A denser gas has more particles per cubic centimeter, so collisions happen more often. Both effects raise the pressure.

Read term by term, $P_{\text{gas}} = \rho k_B T/(\mu m_p)$ says: $n$ counts how many particles can collide per unit volume; $T$ sets the typical kinetic-energy scale; $\rho$ tracks how much matter is packed in; and $\mu m_p$ converts mass density to number density. Gas pressure rises when there are more particles per unit volume or when those particles have larger typical speeds.

Mean molecular weight

The average mass per particle in a gas, in units of the proton mass ( $\mu m_p$ is the mean mass per particle). For fully ionized solar-composition gas $\mu \approx 0.6$ ; a smaller $\mu$ means more particles per gram and therefore more pressure at fixed density and temperature.

For a fully ionized mixture with hydrogen, helium, and metal mass fractions $X$ , $Y$ , $Z$ (which satisfy $X + Y + Z = 1$ ), the mean molecular weight is given below.

Hydrogen contributes one proton and one electron when fully ionized, so it gives many particles per gram; helium contributes one nucleus and two electrons packed into four nucleons, so fewer; heavier elements give fewer still. For ionized solar-composition gas ( $X \approx 0.70$ , $Y \approx 0.28$ , $Z \approx 0.02$ ) this gives $1/\mu \approx 1.62$ , so $\mu \approx 0.6$ — the value we use throughout this reading.

Multiple choice

If $\mu$ decreases while $\rho$ and $T$ stay fixed, does the gas pressure increase or decrease?

IncreaseDecrease

Temperature describes a distribution, not a single speed

Temperature does not assign one “correct speed” to every particle. It describes a distribution of particle speeds and kinetic energies. Some particles move slower than the average, some faster, and the whole distribution shifts as the gas heats up.

Log-scale plot of normalized Maxwell-Boltzmann speed distributions for protons at 10^4 K, 10^6 K, and 1.5 times 10^7 K. The bottom axis shows speed in centimeters per second, and the top axis shows the equivalent kinetic temperature scale. Hotter curves are broader and shifted toward higher speeds. — Figure 3Temperature describes a distribution of particle speeds, not a single speed. Higher temperatures broaden the Maxwell-Boltzmann distribution and shift it toward faster particles; the top axis translates speed into an equivalent single-particle kinetic temperature.ASTR 201 (generated)

The upper axis of that plot translates speed into the single-particle kinetic-temperature scale $T_{\mathrm{eq}} = m_p v^2/(2k_B)$ . It is a way to compare speeds and thermal-energy scales, not a claim that all protons at temperature $T$ move at one speed.

Radiation pressure: momentum from light

Photons carry momentum even though they have no rest mass. If photons are trapped inside a star and scatter repeatedly from matter, they exert a pressure on the gas — called radiation pressure. For a radiation field in local thermal equilibrium, the radiation energy density is $u_{\text{rad}} = a T^4$ , and for an isotropic field the associated pressure is one third of that.

Radiation pressure

The pressure exerted by a trapped, isotropic photon field, $P_{\text{rad}} = \tfrac{1}{3}aT^4$ . Because it scales as $T^4$ while gas pressure scales as $T$ , radiation pressure grows far more rapidly with temperature and dominates only in hot, massive stars.

The important physics is the temperature dependence: $P_{\text{gas}} \propto T$ while $P_{\text{rad}} \propto T^4$ at fixed density. Radiation pressure therefore grows much more rapidly with temperature than thermal gas pressure does. For a star like the Sun, thermal gas pressure dominates; in very hot, very massive stars, radiation pressure becomes much more important.

We can make that trend explicit with a scaling argument. The ratio of the two pressures scales as

$\frac{P_{\text{rad}}}{P_{\text{gas}}} \propto \frac{T^4}{\rho T} = \frac{T^3}{\rho}.$

Later in this reading we estimate $T_c \sim M/R$ and $\rho \sim M/R^3$ . Substituting those hydrostatic scaling estimates gives

$\frac{P_{\text{rad}}}{P_{\text{gas}}} \propto \frac{(M/R)^3}{M/R^3} \propto M^2.$

This scaling uses simplified hydrostatic estimates and ignores detailed structure and composition. It shows a trend, not exact stellar interiors.

Log-log plot of normalized gas pressure and radiation pressure versus temperature. The gas-pressure curve rises linearly with temperature, while the radiation-pressure curve rises much more steeply with a fourth-power dependence. — Figure 4Normalized pressure scalings at fixed density: thermal gas pressure rises linearly with T, while radiation pressure rises as T^4. The goal is to compare slopes, not to build a full stellar model.ASTR 201 (generated)

Deep Dive: Deep dive — photons as a fluid: pressure, diffusion, and why stars look like blackbodies

So far we have treated radiation pressure as a formula, $P_{\rm rad} = \tfrac{1}{3} a T^4$ . Where does it come from physically? Three ideas connect: photons carry momentum, photons interact repeatedly with matter inside stars, and those interactions drive the radiation field toward thermal equilibrium.

1. Photons carry momentum → radiation can push. Even though photons have no rest mass, they carry energy and momentum, $p = E/c$ . When absorbed or scattered, they transfer momentum to matter. Pressure is momentum transfer per unit area per unit time, so photons hitting matter from all directions exert a real pressure.

2. Inside a star, photons do not free-stream — they diffuse. In space photons travel freely; inside a star they are constantly absorbed and re-emitted or scattered. The mean free path is tiny compared with the stellar radius, so photons execute a random walk — radiative diffusion. Two consequences follow: the radiation field becomes nearly isotropic, and energy transport is slow and governed by gradients rather than free streaming.

3. Why the factor of one third? The radiation energy density is $u_{\rm rad} = a T^4$ , but pressure is the momentum flux along one direction. Averaging an isotropic field over all directions in three dimensions contributes only one third of the energy density to any axis, giving $P_{\rm rad} = \tfrac{1}{3} u_{\rm rad}$ .

4. Thermal equilibrium → blackbody radiation. Because photons are repeatedly absorbed, re-emitted, and scattered, the radiation field relaxes to local thermal equilibrium (LTE). In LTE the spectrum becomes a blackbody spectrum and the energy density depends only on temperature. This is why stars behave approximately like blackbodies: deep inside, radiation and matter are tightly coupled; the surface layers then emit something close to that blackbody radiation.

5. Pressure and energy transport are linked. The same photons that push outward (pressure) also carry energy outward (luminosity). That is why stellar structure and stellar luminosity are deeply connected.

Synthesis — why radiation pressure matters. Radiation pressure is not a separate “extra” effect. It is the natural consequence of hot matter emitting photons, photons interacting with matter, and the system reaching thermal equilibrium. In massive stars, where temperatures are high, this radiation field becomes strong enough to help support the star against gravity.

Radiative diffusion

The slow, random-walk transport of radiation through a stellar interior, where the photon mean free path is tiny compared with the stellar radius. Repeated absorption, re-emission, and scattering make the radiation field nearly isotropic and drive it toward a local blackbody spectrum.

Degeneracy pressure

A quantum-mechanical pressure arising from the Pauli exclusion principle, independent of temperature. It supports white dwarfs (electron degeneracy) and neutron stars (neutron degeneracy), and dominates at high density and low temperature.

Problem

Suppose the density stays the same but the temperature doubles. How does the thermal gas pressure change? How does the radiation pressure change? Which source becomes relatively more important in a hotter star?

Hydrostatic Equilibrium

Part 2: The Equation of Hydrostatic Equilibrium

Setting up the force balance

Generated shell-force diagram with a rectangular gas shell labeled by density rho, area A, and thickness dr, plus arrows and equations for the inner pressure force, outer pressure force, and inward gravitational force. — Figure 5The hydrostatic equation comes directly from shell bookkeeping: pressure on the inner face pushes outward, pressure on the outer face pushes inward, and gravity pulls the shell inward.ASTR 201 (generated)

Consider a thin shell of gas at radius $r$ inside a star, with thickness $dr$ , cross-sectional area $A$ , and density $\rho(r)$ . Its mass is $dm = \rho(r)A\,dr$ . If the inward and outward forces do not cancel, the shell accelerates — which is exactly why this equation matters: hydrostatic equilibrium is the condition for a star to remain nearly static instead of beginning a rapid global readjustment.

Two forces act on the shell in the radial direction. The pressure on the inner face pushes outward, $F_{\text{in}} = P(r)A$ ; the pressure on the outer face pushes inward, $F_{\text{out}} = P(r+dr)A$ . The net pressure force is

$F_P = P(r)A - P(r+dr)A.$

Using the first-order Taylor expansion $P(r+dr) \approx P(r) + \frac{dP}{dr}\,dr$ , this becomes

$F_P \approx -A\frac{dP}{dr}\,dr.$

Because pressure decreases outward, $dP/dr < 0$ , so the pressure-gradient force points outward. Gravity pulls the shell inward with $F_g = -g(r)\,dm = -\frac{G M(r)\rho(r)}{r^2}A\,dr$ , using $g(r) = G M(r)/r^2$ . For hydrostatic equilibrium the net force vanishes, $F_P + F_g = 0$ :

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

Dividing through by $A\,dr$ and substituting $g(r)$ gives the equation of hydrostatic equilibrium.

Read as a sentence: at every radius inside the star, the pressure must decrease outward at exactly the rate needed to support the weight of the overlying gas. Since $\rho > 0$ and $g > 0$ , the right-hand side is negative, so $dP/dr$ must also be negative — pressure is highest at the center and falls to nearly zero at the surface. This is a local force-balance law; solving for the full structure $P(r)$ , $\rho(r)$ , $M(r)$ needs additional equations (the subject of Reading 5).

Quick check

Classify each statement as a local law or a global scaling estimate, then say what physical question each one answers:

$dP/dr = -\rho g$
$P_c \sim G M^2/R^4$

Observable

The Sun holds its size for billions of years

The Sun keeps nearly the same size for many billions of years, far longer than its $\sim 50$ -minute dynamical timescale.

Model

Force balance on a thin shell

Treat the interior as a stack of thin shells and require the pressure-gradient force to balance gravity on each one.

Inference

Each shell satisfies dP/dr = -ρg

Each shell must satisfy $dP/dr = -\rho g$ . If that balance is violated, the shell accelerates inward or outward instead of staying in place.

Quick check

Hydrostatic equilibrium says $dP/dr < 0$ everywhere inside a star. What would happen if $dP/dr > 0$ at some radius (pressure increasing outward)? What if $dP/dr = 0$ (uniform pressure)?

Estimating the Central Pressure

Part 3: Estimating the Central Pressure

Hydrostatic equilibrium tells us the slope of the pressure profile, but not yet the pressure scale itself. In astronomy we often know a star’s mass and radius before we know its internal structure. A scaling argument lets us infer interior conditions directly from those observable quantities.

Reading the Math: the central pressure

Hydrostatic equilibrium gives the slope of the pressure profile, not the pressure itself. To turn that slope into a number we use the move that drives this whole module — Reading the Math: approximate the derivative, extract the scaling, then name the assumption you just made. We will run it again for the core temperature, for fusion, for radiation transport, and for the entire main sequence. Here is its first full pass.

① Approximate the derivative. We do not know $P(r)$ in detail, but we know its two endpoints: at the center ( $r = 0$ ) the pressure is $P_c$ ; at the surface ( $r = R$ ) it has dropped to essentially nothing, $P_{\text{surf}} \approx 0$ . Approximate the derivative by the average slope between those endpoints:

$\frac{dP}{dr} \;\approx\; \frac{\Delta P}{\Delta r} \;=\; \frac{P_{\text{surf}} - P_c}{R - 0} \;\approx\; \frac{0 - P_c}{R} \;=\; -\frac{P_c}{R}.$

The minus sign is not bookkeeping — it is the physics: pressure decreases outward. This single replacement, a derivative turned into a ratio of global scales, is the engine of every scaling in Module 3.

$Single-panel plot of normalized pressure versus fractional radius for a smooth toy stellar profile, with labels marking central pressure, near-zero surface pressure, a vertical delta-P arrow, a horizontal delta-r arrow, and a dashed secant representing the scale estimate P_c over R.$

Figure 6The pressure profile is smooth across a scale of order the stellar radius, so the gradient scale is a total drop of order P_c across a distance of order R. That is why dP/dr ~ P_c/R is a sensible scaling estimate.ASTR 201 (generated)

② Extract the scaling. Put that approximate gradient into hydrostatic equilibrium, $dP/dr = -\rho g$ , and replace the two remaining local quantities by their global scales — the mean density $\rho \sim M/R^3$ and the surface gravity $g \sim GM/R^2$ :

$-\frac{P_c}{R} \;\sim\; -\rho g \;\sim\; -\frac{M}{R^3}\cdot\frac{GM}{R^2} \;=\; -\frac{GM^2}{R^5}.$

The minus signs match on both sides — the approximation respects the physics — so cancel them and multiply through by $R$ to read off the central-pressure scaling.

Stronger gravity demands higher internal pressure, and smaller radii make that demand rise sharply: $P_c \propto R^{-4}$ . That is why compact stars require enormous internal pressure even when their total mass is not especially large.

③ Name the assumption. Every arrow above hid an approximation. Collect them — this is the audit you return to whenever a scaling disagrees with a real star, and the list Reading 5 will stress-test one row at a time:

We assumed	by replacing	What it costs
pressure vanishes at the surface	$P_{\text{surf}} \to 0$ in $\Delta P$	negligible: $P_{\text{surf}} \ll P_c$
the star has one density	$\rho(r) \to \bar\rho \sim M/R^3$	real stars are centrally concentrated → underestimates $P_c$ (here by $\sim 20\times$ )
one length scale	$dr \to R$	yields a scale, not the true profile $P(r)$

The exponents survive all three approximations; only the coefficient suffers. That is the deal a scaling makes — right exponents, approximate coefficient — and it is exactly why the worked estimate below lands a factor of $\sim 20$ low yet still nails the $P_c \propto M^2/R^4$ dependence.

Two-panel schematic. Left panel shows a thin shell inside a star with an outward pressure force on the inner face, a smaller inward pressure force on the outer face, and inward gravity. Right panel shows a whole star labeled with mass M, radius R, and a high central pressure scale. — Figure 7Conceptual comparison between a local differential law (dP/dr = -rho g) and the global scaling estimate it motivates (P_c ~ GM^2/R^4).ASTR 201 (generated)

Observable

Stellar masses and radii are measurable

Stellar masses and radii can be measured, even when central conditions cannot.

Model

Hydrostatic equilibrium plus a mean-density approximation

Apply the local force-balance law with the interior replaced by global scales, $\rho \sim M/R^3$ and $g \sim GM/R^2$ .

Inference

The required central pressure scale is P_c ~ GM²/R⁴

The required central-pressure scale is $P_c \sim G M^2/R^4$ . Massive or compact stars therefore need much larger central pressure.

Three-panel generated toy-model plot showing normalized enclosed mass, gravitational acceleration, and pressure as functions of fractional radius inside a uniform-density star. — Figure 8Even a uniform-density toy star teaches the right qualitative lesson: enclosed mass grows outward, gravity rises roughly linearly inside, and pressure must peak at the center.ASTR 201 (generated)

Generated heatmap of logarithmic central pressure in dynes per square centimeter as a function of stellar mass and radius in solar units, with contour labels and a marked Sun point. — Figure 9Central pressure depends strongly on compactness. Holding mass fixed while shrinking the radius drives the required support sharply upward because P_c ~ GM^2/R^4.ASTR 201 (generated)

Numeric answer

A star has the same mass as the Sun but half the Sun’s radius. Using $P_c \sim G M^2/R^4$ , by what factor does the required central pressure increase?

Your answer

Worked example: the Sun’s central pressure

Worked Example 1The Sun's central pressure

Problem

Estimate the Sun’s central pressure from the scaling $P_c \sim G M_\odot^2/R_\odot^4$ , using $G = 6.674\times10^{-8}\,\mathrm{cm^3\,g^{-1}\,s^{-2}}$ , $M_\odot = 2.0\times10^{33}\,\mathrm{g}$ , and $R_\odot = 7.0\times10^{10}\,\mathrm{cm}$ . Express the answer in $\mathrm{dyn\,cm^{-2}}$ and in atmospheres.

StepEvaluate the powers

$(2.0\times10^{33}\,\mathrm{g})^2 = 4.0\times10^{66}\,\mathrm{g^2}$ and $(7.0\times10^{10}\,\mathrm{cm})^4 \approx 2.4\times10^{43}\,\mathrm{cm^4}$ .

StepCombine the coefficients

$P_c \sim \dfrac{(6.674\times10^{-8})(4.0\times10^{66})}{2.4\times10^{43}} \approx 1.1\times10^{16}\,\mathrm{g\,cm^{-1}\,s^{-2}}$ .

Dimensional check

$[G M^2/R^4] = (\mathrm{cm^3\,g^{-1}\,s^{-2}})(\mathrm{g^2})(\mathrm{cm^{-4}}) = \mathrm{g\,cm^{-1}\,s^{-2}}$ . Since $1\,\mathrm{dyn} = 1\,\mathrm{g\,cm\,s^{-2}}$ , this is $\mathrm{dyn\,cm^{-2}}$ — a pressure. ✓

Result

$P_c \sim 1.1\times10^{16}\,\mathrm{dyn\,cm^{-2}} \approx 1.1\times10^{10}\,\mathrm{atm}$ — about ten billion atmospheres (dividing by $1.013\times10^6\,\mathrm{dyn\,cm^{-2}\,atm^{-1}}$ ). Detailed solar models give $\approx 2.5\times10^{17}\,\mathrm{dyn\,cm^{-2}}$ , so this estimate is low by about a factor of 20 — expected, because the Sun is centrally concentrated rather than uniform in density.

Problem

Verify that $G M^2/R^4$ has units of pressure in CGS by tracking the units of each factor.

The Virial Theorem

Part 4: The Virial Theorem for Stars

Force balance tells us what support a star needs at each radius, but not how gravity and thermal energy are linked as the star contracts, radiates, and evolves. For that we need an energy argument. The virial theorem is the key tool: it connects the star’s thermal energy to its gravitational binding energy and explains why contraction heats a self-gravitating object instead of cooling it.

Virial theorem

For a bound, self-gravitating system in equilibrium, $2K_{\text{th}} + U_{\text{grav}} = 0$ , so the total thermal energy is fixed at half the magnitude of the (negative) gravitational potential energy. It implies that gravitational contraction heats a star.

Energy balance in self-gravitating systems

Hydrostatic equilibrium is a force-balance statement; the virial theorem is an energy-balance statement relating the star’s total thermal energy to its gravitational potential energy. For a bound star, the gravitational potential energy has a characteristic scale.

The minus sign matters: a bound self-gravitating object has less energy than the same mass dispersed to infinite separation. Contraction (smaller $R$ ) makes $|U_{\text{grav}}|$ larger.

For a star in hydrostatic equilibrium, the virial theorem links this to the thermal energy.

Rearranging, $K_{\text{th}} = -\tfrac{1}{2}U_{\text{grav}}$ , and the total energy $E_{\text{tot}} = K_{\text{th}} + U_{\text{grav}} = \tfrac{1}{2}U_{\text{grav}} = -K_{\text{th}} < 0$ — the star is bound. In quasi-static contraction, roughly half the released gravitational energy increases the thermal energy of the gas, while roughly half must be radiated away.

Numeric answer

Using $U_{\text{grav}} \sim -G M^2/R$ , if a star contracts to one-third its radius at fixed mass, by what factor does $|U_{\text{grav}}|$ change?

Your answer

The negative heat capacity paradox

A star radiates energy from its surface, so its total energy $E_{\text{tot}}$ becomes more negative. Because $E_{\text{tot}} = \tfrac{1}{2}U_{\text{grav}}$ , this means $U_{\text{grav}}$ also becomes more negative — the star contracts into a more tightly bound state. But the virial theorem also says $K_{\text{th}} = -\tfrac{1}{2}U_{\text{grav}}$ , so if $U_{\text{grav}}$ becomes more negative, $K_{\text{th}}$ becomes larger. More thermal energy means higher typical particle speeds and a higher temperature.

So the star gets hotter as it loses energy. This is negative heat capacity, a defining feature of self-gravitating systems.

Negative heat capacity

The property of a self-gravitating system whereby losing total energy raises its temperature: radiation drives contraction, contraction deepens the gravitational well, and the virial theorem converts that into a higher thermal energy. It is how gravitational contraction heats a protostar toward fusion.

Two-panel generated figure with a virial energy ledger bar chart for gravitational, thermal, and total energy on the left and a causal chain from radiating energy to core heating on the right. — Figure 10The virial theorem fixes the bookkeeping: K_th = -(1/2) U_grav and E_tot < 0. That is why a star can lose energy, contract, and still get hotter.ASTR 201 (generated)

Observable

Protostars shrink and heat instead of free-falling

Protostars radiate energy while gradually shrinking, instead of collapsing in free fall.

Model

A bound self-gravitating gas obeys the virial theorem

A bound self-gravitating gas obeys $2K_{\text{th}} + U_{\text{grav}} = 0$ .

Inference

Losing energy makes the core hotter

As radiation removes total energy, the star contracts, $U_{\text{grav}}$ becomes more negative, and the thermal energy increases. The core gets hotter as the star loses energy.

Think first — energy logic

A protostar loses energy by radiation. Before using any equations, reason physically: does it expand or contract? Does the temperature rise or fall?

Predict the direction of both changes (size and temperature), using physical reasoning rather than formulas.

Then open the reasoning below to check.

Quick check

Check your prediction: as a protostar radiates energy away, does it expand or contract, and does its core temperature rise or fall? Reason physically.

Estimating the Core Temperature

Part 5: Estimating the Core Temperature

Pressure is the macroscopic requirement. To finish answering our guiding question, we connect that required pressure to microscopic particle motion — and therefore to temperature. We run Reading the Math a second time, with one honest difference: there is no new derivative to approximate here. The derivative work was already spent in Part 3 getting $P_c$ ; now we only need two pressure estimates to agree.

① No new derivative — equate the two pressure scales. Part 3 gave the hydrostatic requirement $P_c \sim G M^2/R^4$ . The ideal-gas law gives the thermal pressure the core actually supplies, $P_c \sim \rho\,k_B T_c/(\mu m_p)$ . A star in balance must satisfy both, so set them equal — using the same mean density $\rho \sim M/R^3$ :

$\frac{M}{R^3}\cdot\frac{k_B T_c}{\mu m_p} \;\sim\; \frac{G M^2}{R^4}.$

② Extract the scaling. Cancel one factor of $M$ and multiply by $R^3$ , leaving $k_B T_c/(\mu m_p) \sim G M/R$ . Rearranging for $T_c$ gives the core-temperature scaling.

The left-hand side is the thermal-energy scale per particle; the right-hand side is the gravitational-energy scale per particle. Hydrostatic support requires these to be comparable — so gravity alone predicts a stellar core temperature of order $10^7\,\mathrm{K}$ , no nuclear physics needed.

③ Name the assumption. This pass inherits every row of the Part 3 audit — it is built on top of $P_c \sim G M^2/R^4$ — and adds one of its own:

We assumed	by replacing	What it costs
gas pressure dominates	$P \to P_{\text{gas}}$ , dropping $P_{\text{rad}} = \tfrac{1}{3}aT^4$	fine for the Sun; in massive stars $P_{\text{rad}}/P_{\text{gas}} \propto M^2$ grows and lowers the required $T_c$

The audit only ever grows: each Reading-the-Math pass stacks its assumptions on the ones before. Reading 5 collects the full stack and tests where it finally breaks.

Generated log-log plot of core temperature scale in megakelvin versus stellar mass, showing a hydrostatic scaling curve based on T proportional to M over R, a Sun point, and a dashed comparison for the unrealistic fixed-radius case. — Figure 11Hydrostatic equilibrium sets a fusion-scale temperature for every main-sequence star, but the core temperature rises only modestly with mass because more massive stars are also larger.ASTR 201 (generated)

Worked example: the Sun’s core temperature

Worked Example 2The Sun's core temperature

Problem

Estimate the Sun’s core temperature from $T_c \sim \mu G M_\odot m_p/(k_B R_\odot)$ , using $\mu = 0.6$ , $G = 6.674\times10^{-8}\,\mathrm{cm^3\,g^{-1}\,s^{-2}}$ , $M_\odot = 2.0\times10^{33}\,\mathrm{g}$ , $m_p = 1.67\times10^{-24}\,\mathrm{g}$ , $k_B = 1.38\times10^{-16}\,\mathrm{erg\,K^{-1}}$ , and $R_\odot = 7.0\times10^{10}\,\mathrm{cm}$ .

StepEvaluate the numerator

$(0.6)(6.674\times10^{-8})(2.0\times10^{33})(1.67\times10^{-24}) \approx 1.34\times10^{2}\,\mathrm{cm^3\,s^{-2}}$ .

StepEvaluate the denominator

$(1.38\times10^{-16})(7.0\times10^{10}) \approx 9.66\times10^{-6}\,\mathrm{erg\,K^{-1}\,cm}$ .

Dimensional check

$\dfrac{\mathrm{cm^3\,g^{-1}\,s^{-2}}\cdot\mathrm{g}\cdot\mathrm{g}}{\mathrm{erg\,K^{-1}\,cm}} = \dfrac{\mathrm{g\,cm^2\,s^{-2}}}{\mathrm{erg}}\,\mathrm{K} = \mathrm{K}$ , because $1\,\mathrm{erg} = 1\,\mathrm{g\,cm^2\,s^{-2}}$ . ✓

Result

$T_c \sim (1.34\times10^{2})/(9.66\times10^{-6}) \approx 1.4\times10^{7}\,\mathrm{K} = 14\,\mathrm{MK}$ . Detailed solar models give $\approx 1.5\times10^{7}\,\mathrm{K} = 15\,\mathrm{MK}$ , so this stripped-down estimate is already in the correct solar ballpark — a scaling success, not an exact stellar-structure solution.

Observable

The Sun's measured mass and radius

The Sun’s mass ( $M_\odot = 2.0\times10^{33}\,\mathrm{g}$ ) and radius ( $R_\odot = 7.0\times10^{10}\,\mathrm{cm}$ ) — both measured from binary orbits and angular size plus distance.

Model

Hydrostatic equilibrium plus the ideal-gas picture

Combine hydrostatic equilibrium, gas-pressure dominance, the ideal-gas law, and the mean-density scaling $\rho \sim M/R^3$ .

Inference

A core temperature of order 10⁷ K

$T_c \sim 10^7\,\mathrm{K}$ — about 15 million kelvin, hot enough for nuclear fusion. The temperature needed for fusion is set by gravity.

What this temperature means

A core temperature of $T_c \approx 1.5\times10^7\,\mathrm{K}$ corresponds to an average thermal energy per particle of $\tfrac{3}{2}k_B T_c \approx 3.1\times10^{-9}\,\mathrm{erg} \approx 1.9\,\mathrm{keV}$ , or about $2\,\mathrm{keV}$ (using $1\,\mathrm{keV} = 1.602\times10^{-9}\,\mathrm{erg}$ ). At these temperatures atoms are fully ionized, so the solar core is a plasma of free electrons, protons, and helium nuclei.

But this is still not enough thermal energy to overcome the proton–proton Coulomb barrier classically. That barrier is of order $1\,\mathrm{MeV} = 10^3\,\mathrm{keV}$ , hundreds of times larger than the typical thermal energy. So classical thermal motion alone should not allow fusion — the missing ingredient is quantum tunneling, the subject of Reading 3.

Problem

Using $T_c \propto M/R$ and the main-sequence mass–radius relation $R \propto M^{0.8}$ , how does core temperature scale with mass? Is a $10\,M_\odot$ star’s core hotter or cooler than the Sun’s?

Step	Physics	What it determines
Gravity	$g(r)$	inward pull
Hydrostatic equilibrium	$dP/dr$	pressure gradient
Pressure scaling	$P_c$	central pressure
Ideal gas	$P \leftrightarrow T$	temperature
Virial theorem	$K,\ U$	energy balance

Synthesis and Reference

Synthesis: The Stellar Reasoning Ladder

Step back and look at the chain we have built. Each rung uses one physical idea to infer the next:

1. Gravity sets the local inward pull:

$g(r) = \frac{G M(r)}{r^2}.$

2. Force balance on each shell gives hydrostatic equilibrium:

$\frac{dP}{dr} = -\rho g.$

3. That force balance implies a central pressure scale:

$P_c \sim \frac{G M^2}{R^4}.$

4. The ideal-gas relation turns a pressure requirement into a temperature requirement:

$T_c \sim \frac{\mu G M m_p}{k_B R}.$

5. The virial theorem explains why contraction heats the core:

$2K_{\text{th}} + U_{\text{grav}} = 0.$

6. Gravity therefore drives stars toward fusion temperatures of order $10^7\,\mathrm{K}$ .

Starting from gravity alone, and adding only force balance, energy balance, and a thermal-gas model, we find that stars naturally develop hot cores. Fusion is not an arbitrary extra ingredient pasted onto stars later — gravity itself drives the star to the temperature scale where fusion becomes possible.

Quick check

This reading handed you four tools. Without re-deriving anything, name which one answers each question — hydrostatic equilibrium, the virial theorem, the ideal-gas law, or the dynamical timescale:

How quickly does a star restore balance after a small squeeze?
Why does a contracting protostar get hotter as it radiates energy away?
What central pressure must a star of a given mass and radius sustain?
What core temperature does that pressure imply?

Reference Tables

Key results from hydrostatic equilibrium

Quantity	Scaling / Formula	Sun value
Central pressure	$P_c \sim \frac{G M^2}{R^4}$	$\sim 10^{16}\,\mathrm{dyn\,cm^{-2}}$
Core temperature	$T_c \sim \frac{\mu G M m_p}{k_B R}$	$\sim 1.5\times10^{7}\,\mathrm{K}$
Core thermal energy/particle	$\frac{3}{2}k_B T_c$	$\sim 2\,\mathrm{keV}$
Virial relation	$2K_{\text{th}} + U_{\text{grav}} = 0$	—

Symbol legend

Symbol	Meaning	CGS units
$P$	gas pressure	$\mathrm{dyn\,cm^{-2}}$ ( $= \mathrm{g\,cm^{-1}\,s^{-2}}$ )
$n$	number density	$\mathrm{cm^{-3}}$
$\rho$	mass density	$\mathrm{g\,cm^{-3}}$
$M(r)$	mass enclosed within radius $r$	g
$g(r)$	local gravitational acceleration	$\mathrm{cm\,s^{-2}}$
$u_{\text{rad}}$	radiation energy density	$\mathrm{erg\,cm^{-3}}$
$a$	radiation constant	$\mathrm{erg\,cm^{-3}\,K^{-4}}$
$k_B$	Boltzmann constant	$1.38\times10^{-16}\,\mathrm{erg\,K^{-1}}$
$m_p$	proton mass	$1.67\times10^{-24}\,\mathrm{g}$
$\mu$	mean molecular weight	dimensionless ( $\approx 0.6$ for ionized solar gas)
$K_{\text{th}}$	total thermal energy	erg
$U_{\text{grav}}$	gravitational potential energy	erg (negative)

Conservation laws at work

Conservation law	Where it appears	What it constrains
Momentum conservation	Hydrostatic equilibrium ( $dP/dr = -\rho g$ )	Force balance at every radius — net force on each shell is zero
Energy conservation	Virial theorem ( $2K_{\text{th}} + U_{\text{grav}} = 0$ )	Relationship between thermal and gravitational energy in equilibrium

Summary: Gravity vs. Pressure — Round 1

The most important ideas from this reading:

Hydrostatic equilibrium, $dP/dr = -\rho g$ , is the foundational force-balance equation of stellar structure: the pressure gradient at each radius exactly balances the local weight of the overlying gas.
The virial theorem, $2K_{\text{th}} + U_{\text{grav}} = 0$ , connects thermal to gravitational energy and explains negative heat capacity: as a star loses energy, it contracts and gets hotter.
The core temperature follows from mass and radius, $T_c \sim \mu G M m_p/(k_B R)$ ; for the Sun this gives $T_c \sim 10^7\,\mathrm{K} \approx 15\,\mathrm{MK}$ .
Stars are dynamically stable: departures from hydrostatic equilibrium are corrected on the dynamical timescale — only minutes for the Sun.

Where we stand · Reading 2 of 5

The through-line of Module 3: what holds a star up, and what makes it shine?

✓ Settled. A pressure gradient — not pressure itself — balances gravity at every radius: $dP/dr = -\rho g$ . Running Reading the Math on that balance gives the central pressure $P_c \sim G M^2/R^4$ , and coupling it to the ideal gas law gives a core temperature $T_c \sim \mu G M m_p/(k_B R) \sim 15~\text{MK}$ — gravity alone fixes the core-temperature scale, no nuclear physics required. The virial theorem explains why contraction heats the core toward that scale.

? Still open. That hot core radiates. Without a power source the star cools and contracts on the Kelvin–Helmholtz time, $\sim 30~\text{Myr}$ — far short of the Sun’s $4.6~\text{Gyr}$ . What keeps the core hot?

→ Next. The energy source. Reading 3 — nuclear fusion, and why fusing protons at “only” 15 MK (when the Coulomb barrier demands $\sim 10~\text{GK}$ ) takes all four fundamental forces of nature.

Quick retrieval: why does a star need a pressure gradient, not just high pressure?

Glossary

Degeneracy pressure: A quantum-mechanical pressure arising from the Pauli exclusion principle, independent of temperature. It supports white dwarfs (electron degeneracy) and neutron stars (neutron degeneracy), and dominates at high density and low temperature.
Dynamical timescale: The characteristic time for a star to respond mechanically to a force imbalance — roughly the free-fall time, $\tau_{\text{dyn}} \sim 1/\sqrt{G\bar\rho}$ . For the Sun it is about 50 minutes; departures from hydrostatic equilibrium are corrected on this timescale.
Hydrostatic equilibrium: The condition in which the outward pressure-gradient force exactly balances the inward pull of gravity at every radius, so the gas has no net radial acceleration. It is the foundational force-balance equation of stellar structure.
Mean molecular weight: The average mass per particle in a gas, in units of the proton mass ( $\mu m_p$ is the mean mass per particle). For fully ionized solar-composition gas $\mu \approx 0.6$ ; a smaller $\mu$ means more particles per gram and therefore more pressure at fixed density and temperature.
Negative heat capacity: The property of a self-gravitating system whereby losing total energy raises its temperature: radiation drives contraction, contraction deepens the gravitational well, and the virial theorem converts that into a higher thermal energy. It is how gravitational contraction heats a protostar toward fusion.
Pressure gradient: The rate at which pressure changes with position, $dP/dr$ . A uniform pressure exerts no net force; only a pressure gradient produces one. Inside a star the gradient points outward (pressure falls with radius) and supports the weight of the overlying gas.
Radiation pressure: The pressure exerted by a trapped, isotropic photon field, $P_{\text{rad}} = \tfrac{1}{3}aT^4$ . Because it scales as $T^4$ while gas pressure scales as $T$ , radiation pressure grows far more rapidly with temperature and dominates only in hot, massive stars.
Radiative diffusion: The slow, random-walk transport of radiation through a stellar interior, where the photon mean free path is tiny compared with the stellar radius. Repeated absorption, re-emission, and scattering make the radiation field nearly isotropic and drive it toward a local blackbody spectrum.
Virial theorem: For a bound, self-gravitating system in equilibrium, $2K_{\text{th}} + U_{\text{grav}} = 0$ , so the total thermal energy is fixed at half the magnitude of the (negative) gravitational potential energy. It implies that gravitational contraction heats a star.