Hydrostatic Equilibrium in Stars

Hydrostatic equilibrium is the single most-used idea in the theory of stars, and it wears three faces at once. Mechanically it is a force balance: at every radius the outward push of the pressure gradient exactly cancels the inward pull of gravity, so no mass shell accelerates. Structurally it is the first of the four equations of stellar structure: the momentum equation with its inertial term deleted, coupling to mass continuity, an equation of state, and energy transport to close the problem that produces a whole star. And dynamically it is a statement about timescales: it holds not because a star is frozen but because the free-fall time on which gravity would collapse it — about half an hour for the Sun — is fourteen orders of magnitude shorter than the nuclear time on which it evolves, so any mechanical imbalance is erased almost the instant it appears. Left to gravity alone a self-gravitating gas ball would fall in on itself in minutes to hours. That the Sun has instead burned steadily for four and a half billion years means some agent almost perfectly opposes gravity everywhere inside it, and that agent is the gradient of gas-plus-radiation pressure. This article takes the bare force balance, writes it in both the radius and mass coordinates that stellar-structure codes actually integrate, and shows the two great consequences that follow from it alone with no equation of state at all: the virial theorem, recalled by name in aq-astro-virial-theorem-stars, and a rigorous lower bound on the central pressure. The explicit run of temperature and the role of radiation pressure are developed in the sibling nodes and cited here, not re-derived.

Key takeaways

Hydrostatic equilibrium balances the inward weight of each mass shell against the outward force of the pressure gradient; when they cancel, the star neither collapses nor expands.
The governing equation is dP/dr = −G Mᵣ ρ / r², the first of the four coupled equations of stellar structure.
It is exact only because the dynamical (free-fall) time of the Sun, ≈ 1.6×10³ s ≈ 27 min, is ~2×10¹⁴ times shorter than its ~10¹⁰ yr nuclear lifetime — any imbalance is corrected almost instantly.
Combining it with mass continuity gives an order-of-magnitude central pressure Pᶜ ∼ G M² / R⁴ ≈ 10¹⁵ Pa for the Sun; detailed solar models give 2.5×10¹⁶ Pa.
Multiplying by volume and integrating yields the virial theorem, 3∫(P/ρ)dMᵣ = −E_grav, linking a star's thermal and gravitational energy.
With no equation of state at all it forces a model-independent bound Pᶜ > G M² / (8π R⁴) ≈ 4.5×10¹³ Pa for the Sun.
Key idea: at every radius in a static star the outward pressure-gradient force per unit volume exactly cancels the inward gravitational force per unit volume, dP/dr = −G Mᵣ ρ / r².

The core insight

At every radius in a static star the outward pressure-gradient force per unit volume exactly cancels the inward gravitational force per unit volume: $\dfrac{\mathrm{d}P}{\mathrm{d}r} = -\dfrac{G M_r \rho}{r^2}$ .

Where this lives today

Hydrostatic equilibrium is not a museum piece — it is the load-bearing assumption of essentially every stellar model computed today, and its precision is now directly testable. Helioseismology measures the Sun's acoustic oscillation frequencies to parts in $10^4$ using instruments such as SOHO/MDI, GONG, and SDO/HMI, and inverting them reconstructs the interior sound speed $c_s^2=\Gamma_1 P/\rho$ as a function of radius; the agreement with a hydrostatic Standard Solar Model is better than half a percent through most of the interior, a stunning confirmation that the Sun really is in the balance equation (1) predicts. Asteroseismology with Kepler, TESS, and PLATO extends the same test to thousands of other stars. Codes like MESA integrate the Lagrangian form (2) alongside continuity, energy generation, and transport to evolve stars from birth to death. At the extreme end, the relativistic generalisation — the Tolman–Oppenheimer–Volkoff equation — governs neutron stars, where NICER's X-ray timing and LIGO/Virgo's measurement of tidal deformability in GW170817 are pinning down the dense-matter equation of state through the equilibrium structure it produces. And the failure of hydrostatic equilibrium is its own frontier: three-dimensional supernova simulations must abandon the static approximation entirely, following the full momentum equation as a massive star's core collapses in less than a second and rebounds — the one regime where $\rho\,\ddot r$ is the whole story.

The physical picture

Picture a column of bricks: each brick carries the weight of every brick above it, so the compressive load rises steadily toward the bottom. A star is the same picture wrapped into a sphere and made of gas. Each thin layer must hold up the weight of everything above it, and it does so by being at higher pressure than the layer just outside it. The support therefore comes not from the pressure itself but from its difference across a layer — the gradient — pointing outward because pressure rises inward. Add a constant to the pressure everywhere and nothing changes, exactly as buoyancy in a fluid comes from pressure differences and not from absolute pressure. The first figure isolates this balance on a single shell and reads off the equation it implies.

Figure

A spherical shell of gas at radius r feels the enclosed mass Mᵣ pulling it inward with weight g ρ dr dA while the higher pressure on its inner face pushes it outward by −(dP/dr) dr dA. Equating the two per unit volume gives dP/dr = −G Mᵣ ρ / r².

Why is the star not violently oscillating between these two forces? Because it is an exquisite thermostat operating far faster than it evolves. If the pressure momentarily sags, gravity wins for an instant, the core contracts and heats, the pressure climbs, and balance is restored — all within the dynamical time. That time is set by how fast gravity alone could collapse the star, t_dyn ∼ (G ρ̄)^{−1/2} ∼ √(R³/GM), and for the Sun it is only about half an hour. The nuclear time on which the star actually changes — burning its hydrogen — is ~10¹⁰ years. The second figure places these two clocks on one logarithmic axis; their fourteen-decade separation is the entire justification for setting the acceleration to zero.

Figure

The two clocks of a star on a logarithmic time axis: the free-fall time t_dyn ≈ 1.6×10³ s (≈ 27 min) on which gravity would collapse the Sun, versus the nuclear time t_nuc ≈ 3×10¹⁷ s (≈ 10¹⁰ yr) on which it evolves. Their ratio ≈ 2×10¹⁴ is why the inertial term is dropped.

The core equations

The starting point is the force balance on a shell, written with radius as the independent variable. It says that the pressure must drop as you move outward, and by exactly the amount needed to carry the local weight.

\frac{\mathrm{d}P}{\mathrm{d}r} = -\frac{G M_r \rho}{r^2} \qquad (1)

where: $P$ — total pressure (gas + radiation) at radius $r$ ; $r$ — radial coordinate from the centre; $G$ — Newton's constant, $6.674\times10^{-11}\ \mathrm{N\,m^2\,kg^{-2}}$ ; $M_r$ — mass enclosed within $r$ ; $\rho$ — local density. The minus sign is the whole physics: pressure falls outward ( $\mathrm{d}P/\mathrm{d}r<0$ ), so the interior is always hotter and denser than the surface. Stellar-structure codes rarely use radius as the coordinate — mass is conserved as a star evolves while its radius is not, so the natural independent variable is the enclosed mass $M_r$ . Dividing (1) by the continuity relation converts it to the Lagrangian form.

\frac{\mathrm{d}P}{\mathrm{d}M_r} = -\frac{G M_r}{4\pi r^4} \qquad (2)

where: $M_r$ — enclosed mass, now the independent variable; $r=r(M_r)$ — radius of the shell enclosing $M_r$ . The density has cancelled: this is the form actually integrated in codes such as MESA, paired with $\mathrm{d}r/\mathrm{d}M_r = 1/(4\pi r^2\rho)$ from aq-astro-mass-continuity-equation. The factor multiplying $\rho$ in equation (1) is just the local gravitational acceleration produced by the enclosed mass — a consequence of Newton's shell theorem, so only the mass interior to $r$ contributes.

g(r) = \frac{G M_r}{r^2} \qquad (3)

where: $g$ — inward gravitational acceleration at radius $r$ ; the exterior mass exerts no net force by Gauss's law. With this, equation (1) reads compactly as $\mathrm{d}P/\mathrm{d}r=-\rho g$ . Even before solving anything, dimensional analysis of (1) fixes the scale of the central pressure. Replacing derivatives by whole-star differences and the density by the mean density gives an estimate accurate to an order-unity factor.

P_c \sim \frac{G M^2}{R^4} \qquad (4)

where: $M$ — total stellar mass; $R$ — total radius; $P_c$ — order-of-magnitude central pressure. For the Sun this is $\approx 1.1\times10^{15}\ \mathrm{Pa}$ ; the true value is higher because real stars are centrally concentrated, not uniform. The timescale that certifies the static approximation is the free-fall or dynamical time — how long gravity would take to collapse the star with its pressure switched off. It depends only on the mean density.

t_{\rm dyn} \sim \sqrt{\frac{R^3}{G M}} \sim \frac{1}{\sqrt{G\bar\rho}} \qquad (5)

where: $t_{\rm dyn}$ — dynamical (free-fall) time; $\bar\rho=3M/4\pi R^3$ — mean density. The two forms differ by the order-unity factor $\sqrt{4\pi/3}\approx2$ ; for the Sun, $\sqrt{R^3/GM}\approx1.6\times10^3\ \mathrm{s}$ . Near the surface, gravity and enclosed mass are nearly constant, so equation (1) reduces to the plane-parallel barometric law with an exponential pressure profile — the same law that governs planetary atmospheres, here used to model stellar photospheres.

\frac{\mathrm{d}P}{\mathrm{d}r} = -\rho g \;\Longrightarrow\; P(z)=P_0\,e^{-z/H},\quad H=\frac{k_B T}{\mu m_H g} \qquad (6)

where: $z$ — depth below a reference level; $H$ — pressure scale height (the e-folding length); $\mu$ — mean molecular weight in units of $m_H$ ; $k_B$ — Boltzmann's constant. A larger $g$ or heavier gas gives a thinner atmosphere.

Derivation AHydrostatic equilibrium as the static limit of the momentum equation

The goal is to show that hydrostatic equilibrium is not an independent postulate but the radial momentum equation of a self-gravitating gas with the acceleration term deleted, and to make explicit when that deletion is legitimate. We keep Newton's second law for a gas element with its inertial term intact, estimate the size of that term against gravity for a real star, and set it to zero only after confirming it is negligible — so the reader sees the approximation being made rather than assumed.

Step 1. Apply Newton's second law per unit volume to a gas element at radius $r$ : mass times acceleration equals the pressure-gradient force plus gravity, giving $\rho\,\mathrm{d}^2r/\mathrm{d}t^2 = -\mathrm{d}P/\mathrm{d}r - G M_r\rho/r^2$ . Why: this is the exact radial Euler equation for a spherical gas; every term is a force per unit volume, and nothing has yet been approximated.

Step 2. Estimate the inertial term. If the star changes on a timescale $t$ , then $|\rho\,\ddot r| \sim \rho R/t^2$ , while the gravity term is $\sim \rho G M/R^2 = \rho g$ . Their ratio is $(R^2/gt^2)\sim(t_{\rm dyn}/t)^2$ . Why: this recasts "is the acceleration negligible?" as "is the star evolving slower than it can fall?", turning a vague approximation into a sharp timescale comparison.

Step 3. Insert the numbers for the Sun: $t_{\rm dyn}\approx1.6\times10^3$ s while the nuclear evolution time is $t\sim10^{10}$ yr $\approx3\times10^{17}$ s, so $(t_{\rm dyn}/t)^2\sim3\times10^{-29}$ . The inertial term is utterly negligible. Why: the ratio is so extreme that hydrostatic equilibrium is not merely a good approximation but effectively exact for a main-sequence star.

Step 4. Drop the acceleration to obtain $0=-\mathrm{d}P/\mathrm{d}r - G M_r\rho/r^2$ , i.e. equation (1). Why: this identifies hydrostatic equilibrium precisely as the $t_{\rm dyn}\ll t_{\rm evol}$ limit, and tells us it fails — the full equation must return — during collapse, explosion, or pulsation, when the star is forced to move on its own dynamical time.

Go deeper · Derivation B · 1 minThe virial theorem from the mass-form hydrostatic equation

Derivation BThe virial theorem from the mass-form hydrostatic equation

The goal is to extract a global energy relation from the local force balance by multiplying the Lagrangian hydrostatic equation by volume and integrating over the whole star. This is the deepest single consequence of hydrostatic equilibrium: with no equation of state assumed, it ties the volume-averaged pressure to the gravitational binding energy, and for an ideal gas it becomes the statement that a self-gravitating star has negative heat capacity. We integrate by parts, kill the surface term, and recognise the gravitational-energy integral.

Step 1. Start from the mass form, $\mathrm{d}P/\mathrm{d}M_r = -G M_r/(4\pi r^4)$ , and multiply both sides by the enclosed volume $V=\tfrac{4}{3}\pi r^3$ , then integrate over the star from $M_r=0$ to $M$ . Why: multiplying by volume before integrating is the trick that turns a pressure gradient into a spatially averaged pressure, which is what a global theorem must involve.

Step 2. Integrate the left side by parts: $\int_0^M V\,(\mathrm{d}P/\mathrm{d}M_r)\,\mathrm{d}M_r = [VP]_0^M - \int_0^M P\,(\mathrm{d}V/\mathrm{d}M_r)\,\mathrm{d}M_r$ . The boundary term $[VP]_0^M$ vanishes because $r=0$ at the centre and $P\approx0$ at the surface. Why: the vanishing surface term is exactly the statement that the star is self-bound with negligible external pressure — the boundary condition doing real physical work.

Step 3. Use $\mathrm{d}V/\mathrm{d}M_r = (\mathrm{d}/\mathrm{d}M_r)(\tfrac43\pi r^3) = 4\pi r^2\,\mathrm{d}r/\mathrm{d}M_r = 1/\rho$ , so the remaining left-hand integral is $-\int_0^M (P/\rho)\,\mathrm{d}M_r$ , while the right side is $\tfrac13\int_0^M (G M_r/r)\,\mathrm{d}M_r = -\tfrac13 E_{\rm grav}$ . Why: recognising the right-hand integral as one-third of the gravitational potential energy is what closes the loop between mechanics and energetics.

Step 4. Equate and rearrange to $3\int_0^M (P/\rho)\,\mathrm{d}M_r = -E_{\rm grav}$ . For an ideal monatomic gas $P/\rho = \tfrac23 u$ per unit mass, so the left side is $2U$ , giving the classic $2U + E_{\rm grav}=0$ . Why: this is why a star that loses energy contracts and heats — the hallmark of gravitational thermodynamics, developed fully in aq-astro-virial-theorem-stars.

Go deeper · Derivation C · 1 minA model-independent lower bound on the central pressure

Derivation CA model-independent lower bound on the central pressure

The goal is to prove, using only hydrostatic equilibrium and the fact that radius increases outward, that the central pressure of any star in equilibrium must exceed a definite value set by its mass and radius alone — with no equation of state, no assumption about composition, and no numerical model. The trick is to add to the pressure a cleverly chosen function of the enclosed mass and radius so that the sum is guaranteed to decrease monotonically from centre to surface, after which the bound follows by inspection.

Step 1. Consider the combination $Q(M_r)=P + G M_r^2/(8\pi r^4)$ and differentiate with respect to $M_r$ , using the mass-form hydrostatic equation for $\mathrm{d}P/\mathrm{d}M_r$ . The added term's $\mathrm{d}/\mathrm{d}M_r$ contributes $+G M_r/(4\pi r^4)$ , which exactly cancels $\mathrm{d}P/\mathrm{d}M_r$ . Why: the second term was reverse-engineered precisely to annihilate the pressure gradient, leaving a residue whose sign we can read off.

Step 2. The surviving residue is $\mathrm{d}Q/\mathrm{d}M_r = -G M_r^2/(2\pi r^5)\,(\mathrm{d}r/\mathrm{d}M_r) < 0$ , because $\mathrm{d}r/\mathrm{d}M_r>0$ — radius grows as you enclose more mass. So $Q$ strictly decreases outward. Why: monotonicity is the entire content of the bound; once $Q$ is known to fall from centre to surface, its central value must be the largest.

Step 3. Evaluate the endpoints. At the centre $M_r^2/r^4\to0$ , so $Q(0)=P_c$ . At the surface $Q(R)=P_{\rm surf}+G M^2/(8\pi R^4)\approx G M^2/(8\pi R^4)$ . Since $Q$ decreases, $Q(0)>Q(R)$ , giving $P_c > G M^2/(8\pi R^4)$ . Why: the inequality follows purely from mechanics — any object held up by pressure against its own gravity is forced to this minimum central pressure.

Step 4. Insert solar values: $P_c > 4.5\times10^{13}\ \mathrm{Pa}$ . The true solar central pressure, $2.5\times10^{16}\ \mathrm{Pa}$ , comfortably respects the bound and sits far above it because the Sun is centrally concentrated. Why: the bound is a floor, not an estimate; that the real value exceeds it by ~500× quantifies how far from a uniform sphere a real star is.

Limiting cases

Hydrostatic equilibrium is one limit of the full momentum equation, and its extensions map out where it applies and where it must be replaced.

Case	Result	Meaning
Full dynamical equation	ρ d²r/dt² = −dP/dr − G Mᵣ ρ / r²	Setting the inertial term to zero recovers HSE. Valid whenever accelerations are tiny beside g ∼ 10²–10⁵ m s⁻².
Free-fall (no pressure)	t_dyn ∼ (G ρ̄)^{−1/2} ≈ 27 min (Sun)	If pressure vanished the star would collapse this fast; its brevity beside t_nuc ∼ 10¹⁰ yr certifies the static limit.
Radiation-pressure limit	P → P_rad = ⅓ a T⁴	In very massive stars radiation carries much of the weight; near the Eddington limit outward radiation force approaches gravity and equilibrium becomes marginal.
Polytrope (Lane–Emden)	(1/ξ²) d/dξ(ξ² dθ/dξ) = −θⁿ	With P = K ρ^{1+1/n}, HSE plus Poisson's equation collapse to one dimensionless ODE, yielding analytic stellar models.
Relativistic (TOV)	dP/dr = −G(ρ+P/c²)(Mᵣ+4πr³P/c²) / [r²(1−2GMᵣ/rc²)]	For neutron stars, general relativity corrects HSE; pressure itself gravitates, and no equilibrium exists above ~2–3 M⊙.

Common trapIt is the pressure gradient, not the pressure, that holds a star up

A common slip is to imagine that "high pressure pushes the star outward against gravity." A uniform pressure pushes equally on every face of a fluid element and produces exactly zero net force — you could bathe a star in an ocean of constant pressure and change nothing. Support comes only from the difference in pressure across a layer, the gradient $\mathrm{d}P/\mathrm{d}r$ , higher on the inner face than the outer. A second, subtler trap: hydrostatic equilibrium does not mean the star is static in every sense. It evolves continuously on the nuclear timescale; the equation only asserts that mechanical accelerations are negligible. It genuinely fails — and the full time-dependent momentum equation must return — during gravitational collapse, supernova explosion, or stellar pulsation, whenever $\rho\,\mathrm{d}^2r/\mathrm{d}t^2$ grows to rival gravity and the pressure gradient.

Worked examples

1 · The central pressure of the Sun, three ways (fully worked)

Estimate the Sun's central pressure from the order-of-magnitude formula, bound it rigorously from below, and compare both with the Standard Solar Model. Use M⊙ = 1.99×10³⁰ kg, R⊙ = 6.96×10⁸ m, G = 6.674×10⁻¹¹ SI.

Step 1. The dimensional estimate from equation (4): $P_c \sim G M^2/R^4$ . Numerically $G M^2 = 6.674\times10^{-11}\times(1.99\times10^{30})^2 = 2.64\times10^{50}$ , and $R^4 = (6.96\times10^8)^4 = 2.35\times10^{35}$ , so $P_c \sim 1.1\times10^{15}\ \mathrm{Pa}$ . The estimate gives Pᶜ ∼ 10¹⁵ Pa — about ten billion times atmospheric pressure.

Step 2. The rigorous bound from the third derivation: $P_c > G M^2/(8\pi R^4) = 1.1\times10^{15}/(8\pi) = 4.5\times10^{13}\ \mathrm{Pa}$ . This is a hard floor, Pᶜ > 4.5×10¹³ Pa, true for any equation of state — the estimate in Step 1 already sits a factor ~25 above it, as it must.

Step 3. The Standard Solar Model gives $P_c = 2.5\times10^{16}\ \mathrm{Pa}$ , about 23× the dimensional estimate. The gap is not an error in the physics but a statement about geometry: the naive estimate assumed a uniform sphere, whereas the Sun's density rises from $\sim1.4\times10^3\ \mathrm{kg\,m^{-3}}$ (mean) to $\sim1.5\times10^5\ \mathrm{kg\,m^{-3}}$ at the centre. Central concentration multiplies the pressure by ~20, which is exactly the correction a real solved model supplies over the back-of-envelope scaling.

2 · The Sun's central temperature from equilibrium alone

Combine hydrostatic equilibrium with the ideal-gas law to estimate the central temperature, and show it lands where hydrogen fusion ignites. Take mean molecular weight μ ≈ 0.6 (ionised solar mix), m_H = 1.67×10⁻²⁷ kg, k_B = 1.38×10⁻²³ J K⁻¹.

Step 1. The virial-style estimate sets the thermal energy per particle comparable to the gravitational energy per particle: $k_B T_c \sim G M \mu m_H / R$ . Equivalently, this is what you get by writing $P_c \sim \rho_c k_B T_c/(\mu m_H)$ and using the pressure scale from Example 1. This route needs no detailed model — only equilibrium and the ideal-gas equation of state from aq-astro-equation-of-state-stellar-matter.

Step 2. Insert the numbers: $T_c \sim \dfrac{G M \mu m_H}{R k_B} = \dfrac{6.674\times10^{-11}\times1.99\times10^{30}\times0.6\times1.67\times10^{-27}}{6.96\times10^8\times1.38\times10^{-23}}$ . The numerator is $1.33\times10^{-7}$ and the denominator $9.6\times10^{-15}$ , giving T_c ≈ 1.4×10⁷ K.

Step 3. The Standard Solar Model central temperature is $1.57\times10^7$ K — within 10% of an estimate that used nothing but mechanical balance and the ideal-gas law. A pure equilibrium argument predicts a ~10⁷ K core, precisely the temperature at which the proton–proton chain ignites. Hydrostatic equilibrium alone, without any nuclear input, already tells you the Sun's core must be hot enough to burn hydrogen.

3 · How fast would the Sun collapse without pressure?

Compute the dynamical timescale to confirm the static approximation, using the mean density ρ̄⊙ ≈ 1.4×10³ kg m⁻³.

Step 1. Use the compact form $t_{\rm dyn}=\sqrt{R^3/(G M)}$ : with $R^3=3.37\times10^{26}$ and $G M=1.33\times10^{20}$ , $t_{\rm dyn}=\sqrt{2.54\times10^6}=1.6\times10^3\ \mathrm{s}$ . The Sun would collapse in ≈ 27 minutes if its pressure vanished.

Step 2. The mean-density form gives $1/\sqrt{G\bar\rho}=1/\sqrt{6.674\times10^{-11}\times1.4\times10^3}=3.3\times10^3\ \mathrm{s}\approx54\ \mathrm{min}$ . The two differ only by the order-unity factor $\sqrt{4\pi/3}\approx2.05$ ; both say tens of minutes, and the precise coefficient depends on the collapse geometry.

Step 3. Compare with the nuclear lifetime $t_{\rm nuc}\approx10^{10}\ \mathrm{yr}=3.2\times10^{17}\ \mathrm{s}$ . The ratio is $t_{\rm nuc}/t_{\rm dyn}\approx2\times10^{14}$ . The star re-establishes mechanical balance 10¹⁴ times faster than it evolves, which is why treating it as a sequence of hydrostatic states is not just convenient but essentially exact.

Physics of the portrait

Hydrostatic equilibrium is where mechanics became astrophysics: the moment when Newton's balance of forces, applied to a sphere of hot gas, was made to predict the inside of a star no telescope could ever see. No one did more to turn that balance into a working theory than Arthur Eddington.

Origins · 1 minArthur Eddington and how this was discovered

Arthur Eddington

1882–1944 · the architect of stellar structure theory

Arthur Stanley Eddington took the bare force balance that Lane and Emden had written for a gas sphere and built it into a physical theory of the stars. His central insight was that the pressure holding a star up is not only gas pressure but radiation pressure — the momentum of the photons streaming outward — and that in massive stars the two must be balanced together against gravity. From hydrostatic equilibrium plus radiative transport he derived, in the early 1920s, the mass–luminosity relation, predicting a star's brightness from its mass before anyone knew what powered it. He recognised that radiation pressure sets a ceiling on stellar luminosity, the Eddington limit $L_{\rm Edd}=4\pi G M c/\kappa$ , beyond which a star cannot stay in equilibrium and blows its outer layers away as a wind — for the Sun's mass, $L_{\rm Edd}\approx3.8\times10^4\,L_\odot$ . His 1926 treatise The Internal Constitution of the Stars codified the subject and argued, against fierce scepticism, that stars must be powered by subatomic energy released in their hot cores — a decade before Bethe identified the fusion cycles. Eddington also led the 1919 eclipse expedition that confirmed the bending of starlight by the Sun, giving Einstein's general relativity its first observational triumph. The equations a modern stellar-structure code integrates are, in their bones, Eddington's.

History

1870

— Lane models the Sun as a gas sphere held up by pressure against gravity

1907

— Emden's Gaskugeln systematizes polytropic models (the Lane–Emden equation)

1916

— Eddington adds radiation pressure to the hydrostatic support of stars

1926

— Eddington's The Internal Constitution of the Stars codifies stellar structure

1939

— Chandrasekhar's textbook sets the modern mathematical foundations

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., §10.1–10.2 · Kippenhahn, Weigert & Weiss, Stellar Structure and Evolution, 2nd ed., Ch. 2–3 · Prialnik, An Introduction to the Theory of Stellar Structure and Evolution, Ch. 5 · Hansen, Kawaler & Trimble, Stellar Interiors, 2nd ed., Ch. 1 · Chandrasekhar, An Introduction to the Study of Stellar Structure, §3. The key numerical benchmarks for the Sun ( $M_\odot=1.99\times10^{30}$ kg, $R_\odot=6.96\times10^8$ m, $G=6.674\times10^{-11}$ SI) are: mean density $\bar\rho_\odot\approx1.4\times10^3\ \mathrm{kg\,m^{-3}}$ ; the order-of-magnitude central-pressure estimate $P_c\sim GM^2/R^4\approx1.1\times10^{15}$ Pa; the model-independent lower bound $P_c>GM^2/(8\pi R^4)\approx4.5\times10^{13}$ Pa; the Standard Solar Model value $2.5\times10^{16}$ Pa; the dynamical time $t_{\rm dyn}=\sqrt{R^3/GM}\approx1.6\times10^3$ s; and the virial central-temperature estimate $T_c\approx1.4\times10^7$ K (for $\mu=0.6$ ), compared with the SSM value $1.57\times10^7$ K.

Prerequisites. Mass Continuity Equation · Equation of State for Stellar Matter · Newtonian Gravity and Gauss's Law. Leads to. Equations of Stellar Structure · Virial Theorem for Stars · Radiative Temperature Gradient · Radiation Pressure in Stellar Interiors.

Ask the tutor: Want to see how hydrostatic equilibrium joins mass continuity, energy generation, and radiative transport to build a complete stellar model — or how it bends into the TOV equation for neutron stars?

Is a star really "static" if it is in hydrostatic equilibrium?

Not in every sense. Hydrostatic equilibrium says only that the mechanical accelerations are negligible — no mass shell is measurably speeding up or slowing down at any instant. The star is still evolving continuously: it burns hydrogen, its composition changes, and its radius and luminosity drift on the nuclear timescale of billions of years. The right picture is a sequence of hydrostatic snapshots, each in near-perfect force balance, slowly morphing into the next. Only during genuinely fast events — collapse, explosion, pulsation — does the acceleration term stop being negligible and the star cease to be hydrostatic.

Why does only the mass interior to r matter for gravity?

This is Newton's shell theorem, a consequence of the inverse-square law and spherical symmetry. Any spherically symmetric shell of matter exerts zero net gravitational force on a point inside it — the pulls from all directions cancel exactly. So a gas element at radius $r$ feels gravity only from the mass $M_r$ enclosed within its own radius; everything outside is gravitationally invisible to it. That is why the enclosed mass $M_r$ , not the total mass $M$ , appears in the hydrostatic equation, and why $g(r)=GM_r/r^2$ . The same theorem is what lets us treat each shell independently when integrating the structure equations.

How can equilibrium alone predict the Sun's core temperature?

Because hydrostatic equilibrium fixes the pressure needed to hold the star up, and the ideal-gas law converts that pressure into a temperature once you know the density. The chain is purely mechanical: balance sets $P_c\sim GM^2/R^4$ , and $P=\rho k_B T/\mu m_H$ then gives $k_B T_c\sim GM\mu m_H/R\approx1.4\times10^7$ K. No nuclear physics enters. The striking thing is that this "blind" estimate lands within 10% of the true value and squarely at the ignition temperature of hydrogen fusion — telling you that a ball of gas the Sun's mass and size must have a core hot enough to burn, before you know anything about how it burns.

What is the difference between the Eulerian and Lagrangian forms?

They are the same physics in different coordinates. The Eulerian form, $\mathrm{d}P/\mathrm{d}r=-G M_r\rho/r^2$ , uses radius as the independent variable — natural for thinking about a snapshot in space. The Lagrangian form, $\mathrm{d}P/\mathrm{d}M_r=-G M_r/(4\pi r^4)$ , uses the enclosed mass instead, obtained by dividing through by the continuity relation $\mathrm{d}M_r/\mathrm{d}r=4\pi r^2\rho$ . Stellar-evolution codes prefer the mass coordinate because mass is conserved as a star ages while its radius expands and contracts, so the grid points stay attached to the same fluid and the moving outer boundary is trivial to track.

When does radiation pressure matter for the balance?

Gas pressure alone holds up low- and intermediate-mass stars, but in massive stars ( $\gtrsim10\,M_\odot$ ) the radiation pressure $P_{\rm rad}=\tfrac13 aT^4$ becomes a large fraction of the total, because it scales as $T^4$ and their cores are hot. The hydrostatic equation is unchanged in form — you simply use the total pressure $P=P_{\rm gas}+P_{\rm rad}$ — but the physics shifts: if a star's luminosity approaches the Eddington limit $L_{\rm Edd}=4\pi GMc/\kappa$ , the outward radiation force nearly equals gravity, the pressure gradient the gas can supply runs out, and the outer layers are driven off as a wind. This effectively caps how massive a star can be. The full treatment lives in aq-astro-radiation-pressure-stellar-interiors.

Does hydrostatic equilibrium apply to neutron stars and white dwarfs?

Yes, but with important changes. White dwarfs and neutron stars are held up not by thermal gas pressure but by degeneracy pressure — a quantum effect that persists even at zero temperature — yet the balance they satisfy is still hydrostatic equilibrium with that pressure inserted. For neutron stars the gravity is so strong that Newtonian mechanics fails, and the equation must be replaced by its general-relativistic version, the Tolman–Oppenheimer–Volkoff equation, in which pressure itself gravitates and there is a maximum mass (~2–3 $M_\odot$ ) beyond which no equilibrium exists and collapse to a black hole is unavoidable. The Newtonian equation is the weak-gravity limit of the TOV equation.