Virial Theorem for Stars

The virial theorem is not a historical curiosity — it is a live diagnostic in modern astrophysics. In stellar-evolution codes such as MESA the relation $2E_{\mathrm{int}} + E_{\mathrm{grav}} \approx 0$ is a standard consistency check, and its negative-heat-capacity logic explains everything from pre-main-sequence contraction along the Hayashi track to the runaway helium flash in degenerate cores, where the thermostat fails because pressure is nearly temperature-independent. Its most dramatic modern face is the $\gamma = 4/3$ instability. In stars above roughly 100 solar masses, the core is hot enough that photon–photon collisions make electron–positron pairs; the energy sunk into rest mass softens the equation of state below $\gamma = 4/3$ and the threshold is crossed. The result is a pair-instability supernova — the core loses support, collapses, ignites explosively, and in extreme cases is unbound with no remnant. Candidate events such as the superluminous SN 2006gy, and the "black-hole mass gap" LIGO–Virgo–KAGRA sees in merging binaries, are read as fingerprints of this instability. JWST's hunt for the first very massive Population III stars is partly a search for where $\gamma \to 4/3$ physics — the relation Clausius wrote in 1870 — shaped the first stellar deaths.

A star is a self-gravitating ball of gas settled into balance. Gravity pulls every parcel inward; pressure — sustained by the gas's internal heat — pushes outward. Hydrostatic equilibrium says these cancel locally at every radius. The virial theorem says something stronger and global: in equilibrium the two are locked into a fixed energetic ratio, with exactly half the gravitational binding present at any instant as usable internal heat. Squeeze the star and both energies grow in lockstep; let it expand and both shrink — but the 2-to-1 accounting never changes.

The most startling consequence hides in that minus sign. If the star radiates, $E_{\mathrm{tot}}$ drops; since $E_{\mathrm{tot}} = \tfrac{1}{2}E_{\mathrm{grav}}$ the gravitational energy falls too — the star contracts — but $E_{\mathrm{int}} = -E_{\mathrm{tot}}$, so the internal energy rises and the interior gets hotter. A star that loses energy warms up. This negative heat capacity is why a contracting protostar heats its core until hydrogen ignites, and why a burning star self-corrects: overproduce energy and the core expands and cools, throttling the reactions back down. The theorem is also scale-free — it knows nothing about the density profile, only integrated energies, which is why it is the astrophysicist's favorite envelope tool: feed it $M$ and $R$ and read off a central temperature good to order unity. The one number that does control everything is the adiabatic index $\gamma$. As support shifts from ordinary gas ($\gamma = 5/3$) toward radiation or relativistic degeneracy ($\gamma = 4/3$), the virial coefficient slides to zero, and at $\gamma = 4/3$ the star becomes marginally bound — the physics behind the collapse of the most massive stars and of near-Chandrasekhar cores.

The whole theorem is an integrated consequence of one local statement — hydrostatic equilibrium. We start there, promote it to a global energy relation, and read off the interior conditions it predicts. Every equation below is exact for a star with negligible surface pressure unless a scaling symbol ($\sim$) marks an estimate. The starting point is the balance of pressure gradient against gravity at each radius, taken here as given from the hydrostatic-equilibrium article.

where: $P$ — local pressure; $r$ — radial coordinate; $m(r)$ — enclosed mass within $r$; $\rho(r)$ — local mass density; $G$ — Newton's constant. Reading: at every layer the outward push of the pressure drop exactly supports the weight above. Multiplying equation (1) by the interior volume and integrating over the whole star turns this local balance into a global one; the profile disappears, leaving the pressure integral and the gravitational energy.

where: $V$ — total stellar volume; $E_{\mathrm{grav}}$ — total gravitational potential energy, negative for a bound star. Reading: three times the mean pressure times the volume equals the binding energy — the pressure form of the virial theorem, valid for any pressure source. The gravitational energy is the work to assemble the star shell by shell from infinite dispersal — always negative.

where: $M$ — total stellar mass; $\mathrm{d}m$ — mass of a shell at $r$. Reading: deeper shells are bound more tightly; for a uniform sphere the integral is $-\tfrac{3}{5}GM^2/R$, larger for centrally concentrated stars. To turn the pressure integral into a statement about energy we need the equation of state. For a monatomic ideal gas the internal energy density is $\tfrac{3}{2}P$, so it becomes twice $E_{\mathrm{int}}$.

where: $E_{\mathrm{int}}$ — total internal (thermal) energy of the gas. Reading: the internal energy is exactly minus half the gravitational energy — the form usually meant by "the virial theorem for stars". Rearranging equation (4) exposes the loaded conclusion: the total energy is negative, half the (negative) gravitational energy.

where: $E_{\mathrm{tot}}$ — total (internal plus gravitational) energy. Reading: a bound star sits in a well half as deep as its binding; radiating makes $E_{\mathrm{tot}}$ more negative, forcing contraction and heat. The factor $\tfrac{1}{2}$ is special to $\gamma = 5/3$. For a general index it becomes a simple function of $\gamma$ whose sign decides stability.

where: $\gamma$ — adiabatic index of the dominant pressure source. Reading: $E_{\mathrm{tot}}$ is negative (bound, stable) only for $\gamma > 4/3$; at $\gamma = 4/3$ it vanishes and the star is marginally bound. Finally, the pressure form yields interior conditions directly. With $\langle P\rangle = -E_{\mathrm{grav}}/3V$ and $|E_{\mathrm{grav}}|\sim GM^2/R$, the pressure and — via the ideal-gas law — the temperature follow.

where: $\langle P\rangle$ — mean interior pressure; $R$ — stellar radius; $T_c$ — central temperature; $\mu$ — mean molecular weight (in units of $m_{\mathrm{H}}$); $m_{\mathrm{H}}$ — hydrogen mass; $k$ — Boltzmann constant. Reading: pressure climbs steeply with mass and shrinking radius, and temperature scales as $M/R$ — the reason massive, compact stars run hot.

The virial coefficient linking total energy to gravitational energy is set entirely by the adiabatic index $\gamma$ of the dominant pressure source.

1 · Central temperature of the Sun

Estimate the Sun's central temperature from the virial theorem, using $M_\odot = 2.0\times10^{30}\ \mathrm{kg}$, $R_\odot = 7.0\times10^8\ \mathrm{m}$, $\mu = 0.6$. Compare with the detailed model value.

Step 1. Equate internal energy to half the binding energy. With $N = M/(\mu m_{\mathrm{H}})$ particles each carrying $\tfrac{3}{2}kT_c$, the virial theorem sets $\tfrac{3}{2}\,\dfrac{M}{\mu m_{\mathrm{H}}}\,kT_c \sim \tfrac{1}{2}\,\dfrac{GM^2}{R}$.

Step 2. Solve for the temperature. Rearranging gives $T_c \sim \dfrac{1}{3}\dfrac{GM\mu m_{\mathrm{H}}}{kR}$, the prefactor order unity and profile-dependent.

Step 3. Insert the numbers ($G = 6.67\times10^{-11}$, $m_{\mathrm{H}} = 1.67\times10^{-27}\ \mathrm{kg}$, $k = 1.38\times10^{-23}\ \mathrm{J\,K^{-1}}$):

Step 4. Evaluate. The bracket without the $\tfrac{1}{3}$ is $1.4\times10^7\ \mathrm{K}$; with it, $T_c \sim 4.6\times10^6\ \mathrm{K}$ — so one estimates a few $\times10^6$ to $10^7\ \mathrm{K}$. The detailed model gives $T_c \approx 1.5\times10^7\ \mathrm{K}$; the virial estimate lands within a factor of a few, the shortfall reflecting the central concentration the mean ignores. Either way the theorem correctly predicts a core hot enough for hydrogen fusion from $M$ and $R$ alone.

2 · Kelvin–Helmholtz timescale of the Sun

If the Sun shone only by gravitational contraction, how long could it last? Use $L_\odot = 3.8\times10^{26}\ \mathrm{W}$ with the same $M_\odot$, $R_\odot$.

Step 1. Write the available energy. Only half the gravitational energy is radiated, so the reservoir is $|E_{\mathrm{tot}}| \sim \tfrac{1}{2}\,GM^2/R$, and dividing by luminosity gives $t_{\mathrm{KH}} = |E_{\mathrm{tot}}|/L \sim \dfrac{GM^2}{2RL}$.

Step 2. Insert the numbers.

Step 3. Convert to years. Dividing by $3.16\times10^7\ \mathrm{s\,yr^{-1}}$ gives $t_{\mathrm{KH}} \approx 1.6\times10^7\ \mathrm{yr}$ — only sixteen million years, hundreds of times shorter than the radiometric age of the Earth ($\sim4.5\times10^9$ yr). That mismatch was the nineteenth-century crisis that proved gravitational contraction cannot power the Sun; a nuclear source is required. The same $t_{\mathrm{KH}}$ still sets the contraction time of pre-main-sequence stars, and the same energetics fix the interior pressure scale at $\langle P\rangle \sim GM^2/R^4 \sim 10^{14}\ \mathrm{Pa}$.

The virial theorem was born in thermodynamics, not astronomy. The word "virial" — Latin vis, force — was coined by Rudolf Clausius, who also gave us entropy, and only later became a tool of astrophysics.

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., §2.4 and §10.2 · Kippenhahn, Weigert & Weiss, Stellar Structure and Evolution, 2nd ed., Ch. 3 · Prialnik, Theory of Stellar Structure and Evolution, 2nd ed., Ch. 2 · Binney & Tremaine, Galactic Dynamics, 2nd ed., §4.8 (tensor virial theorem) · Hansen, Kawaler & Trimble, Stellar Interiors, 2nd ed., §1.3. The key results to keep in mind: the integration by parts yields $3\int P\,\mathrm{d}V = -E_{\mathrm{grav}}$ with a vanishing surface term; the closure $u = P/(\gamma-1) = \tfrac{3}{2}P$ at $\gamma = 5/3$; the general coefficient $(3\gamma-4)/(3(\gamma-1))$ gives $\tfrac{1}{2}$ at $\gamma = 5/3$ and $0$ at $\gamma = 4/3$. The solar estimates use $M_\odot = 2.0\times10^{30}\ \mathrm{kg}$, $R_\odot = 7.0\times10^8\ \mathrm{m}$, $L_\odot = 3.8\times10^{26}\ \mathrm{W}$, $\mu = 0.6$, and yield $t_{\mathrm{KH}} \approx 1.6\times10^7\ \mathrm{yr}$, $T_c \sim$ few $\times10^6$–$10^7\ \mathrm{K}$, $\langle P\rangle \sim 10^{14}\ \mathrm{Pa}$. All constants are CODATA/IAU values.

Prerequisites. Hydrostatic Equilibrium in Stars · Mass Continuity Equation · Equation of State for Stellar Matter. Leads to. Equations of Stellar Structure · Radiation Pressure in Stellar Interiors · The Energy Generation Equation.

Ask the tutor: Ask me to show how the same virial reasoning that heats a contracting protostar also weighs a cluster of galaxies and reveals its dark matter.