How pressure holds up a star

Choose a small radial column

Let $r$ be distance from the centre, increasing outward, and let $m(r)$ be the mass inside that radius. Use metres, kilograms and pascals. The local density is $\rho$ in kilograms per cubic metre; pressure is $P$ . Spherical symmetry lets us describe the structure with one radial coordinate.

Consider a narrow radial column with cross-sectional area $A$ and thickness $dr$ . The gas just below it pushes outward slightly harder than the gas just above it pushes inward. The pressure force is minus the outward pressure change times area. The derivative $dP/dr$ measures that outward change per metre.

The column’s mass is density times volume. Its gravitational acceleration is $Gm/r^2$ inward, where $G$ is the gravitational constant. Equate the outward pressure force to that inward weight and cancel the column’s area and thickness.

The minus sign records the direction: pressure decreases as radius increases. It does not mean pressure itself is negative. This equation states mechanical balance at one position; it does not determine the temperature, composition or energy source by itself.

-\frac{dP}{dr}\,A\,dr=\rho\,A\,dr\,\frac{Gm}{r^2}\quad\Longrightarrow\quad\frac{dP}{dr}=-\rho\frac{Gm}{r^2}

The small column has outward-positive radial coordinate. The expression assumes spherical Newtonian gravity and negligible radial acceleration.

Sources: Paxton et al., MESA IV, arXiv:1710.08424v1, energy equations and convective boundaries

Count mass as well as radius

A thin spherical shell occupies volume $4\pi r^2dr$ . Multiplying by density gives its mass. This connects the outward radius coordinate with the enclosed-mass coordinate.

Divide the pressure equation by the mass equation. Density cancels, and the weight of the overlying material appears as a pressure change per kilogram of enclosed mass. A saved MESA profile supplies the enclosed masses, radii and pressures needed to check this relation with finite differences.

The massive-star examples in the atlas rotate. Their radial coordinate is volume-equivalent: a sphere of that radius has the same enclosed volume as the averaged rotating surface. The source includes a rotational factor in the pressure equation. Its saved pressure differences agree much better with the corrected equation than with the spherical expression alone.

\frac{dm}{dr}=4\pi r^2\rho,\qquad\frac{dP}{dm}=-\frac{Gm}{4\pi r^4}

These are the spherical, nonrotating relations. The atlas’s rotating family uses the recorded factor fp_rot in its pressure equation; the original polytrope mode uses the spherical limit.

Sources: MESA 25.12.1, history, profile and restart output documentation · Keszthelyi et al., arXiv:2209.06350v1, rotating stellar grids and physical prescriptions · Keszthelyi (2022), fixed reproduction package, Zenodo 7069766, CC BY 4.0

Solve one sphere all the way to the centre

For a first exact calculation, choose constant density. Let the sphere have total mass $M$ and radius $R$ , with negligible external pressure. This is a mathematical reference, not a thermal model of a real star. The enclosed mass then increases as the cube of radius.

Insert that mass into hydrostatic balance. The pressure gradient is proportional to minus r, so integrate from a chosen radius out to the surface, where the pressure is zero. The resulting pressure is largest at the centre and falls smoothly to zero at the edge.

For a sphere with solar mass and radius, use approximately $M=1.99\times10^{30}\,\mathrm{kg}$ and $R=6.96\times10^8\,\mathrm{m}$ . The constant-density central pressure is about $1.34\times10^{14}\,\mathrm{Pa}$ . This number belongs to the uniform sphere; an actual stellar model has a strongly varying density.

The formula is also a good scaling test. At fixed mass, doubling the radius reduces central pressure by a factor of sixteen. The gas occupies a larger volume, and gravity at a given radius fraction is weaker. Both changes reduce the pressure gradient needed for support.

\rho=\frac{3M}{4\pi R^3},\quad m(r)=M\frac{r^3}{R^3},\quad\frac{dP}{dr}=-\frac{4\pi G\rho^2}{3}r

Constant density fixes both the mass profile and the linear radial pressure gradient.

P(r)=\frac{2\pi G\rho^2}{3}(R^2-r^2),\qquad P_c=\frac{3GM^2}{8\pi R^4}

The boundary condition is

P(R)=0

. The pressure unit is force per area; the right-hand side has SI units

\mathrm{kg\,m^{-1}\,s^{-2}}

Sources: Paxton et al., MESA IV, arXiv:1710.08424v1, energy equations and convective boundaries

Mechanical balance can adjust quickly

A star need not wait for a large fraction of its fuel to burn before it reacts to a pressure imbalance. Gravity and pressure can rearrange the gas on a much shorter mechanical time. Estimate that time by asking how long gravity takes to move material over a distance comparable to the radius.

The dynamical time is roughly $t_{\mathrm{dyn}}\sim\sqrt{R^3/(GM)}$ . With solar mass and radius it is about sixteen hundred seconds, or twenty-seven minutes. This is an order-of-magnitude adjustment time, not the period of every stellar oscillation.

A slowly evolving star can therefore move through a long sequence of nearly balanced structures. Each snapshot has local pressure support even though its temperature and composition change gradually. The word equilibrium here refers to the fast mechanical balance.

Sources: Paxton et al., MESA IV, arXiv:1710.08424v1, energy equations and convective boundaries

Losing heat sets a different clock

A star radiates energy into space. If gravity supplied that radiation through contraction, the available energy scale would be gravitational binding energy. Divide that scale by luminosity to obtain the Kelvin–Helmholtz time. Luminosity is energy per second, so this ratio has units of seconds.

With solar mass, radius and luminosity, the estimate is about $3\times10^7\,\mathrm{yr}$ . It is vastly longer than the mechanical adjustment time. A star can consequently remain close to hydrostatic balance throughout substantial thermal evolution.

To follow that thermal evolution we need an equation for energy, an equation for heat transport and a description of the material. A polytropic relation between pressure and density closes a useful static model, but does not supply those missing evolutionary equations. It cannot tell us how much hydrogen has burnt or assign a stellar age.

t_{\mathrm{KH}}\sim\frac{GM^2}{RL}

M

R

and

L

denote stellar mass, radius and luminosity. This thermal timescale omits a structure-dependent numerical factor. A nuclear lifetime instead depends on the available fuel.

Sources: Paxton et al., MESA IV, arXiv:1710.08424v1, energy equations and convective boundaries · Jermyn et al., MESA VI, arXiv:2208.03651v2, Section V energy consistency

The source of pressure matters too

For an ordinary hot gas, temperature and density help determine pressure. Radiation also contributes in a luminous star. In a white dwarf, densely packed electrons contribute degeneracy pressure, associated with the quantum states the electrons occupy. Its mechanical support can remain strong while the remnant loses heat.

This separates two questions that are easy to mix together: what supports the gas, and what supplies the emitted energy? Hydrostatic balance answers the first through a pressure gradient. The next chapter follows how energy crosses the layers of the star.

Sources: Jermyn et al., MESA VI, arXiv:2208.03651v2, Section V energy consistency · Kumar and Townsley (2025), fixed white-dwarf cooling reproduction package, Zenodo 17765239, CC BY 4.0

How pressure holds up a star

Choose a small radial column

Count mass as well as radius

Solve one sphere all the way to the centre

Mechanical balance can adjust quickly

Losing heat sets a different clock

The source of pressure matters too

Try the reasoning

Check the uniform-sphere pressure scaling

Primary sources