The virial theorem provides a fundamental link between the internal pressure (or kinetic energy) and the gravitational binding energy of a self-gravitating system. For a star, galaxy, or molecular cloud in steady equilibrium, it can be derived directly from the condition of hydrostatic equilibrium.

Starting from hydrostatic equilibrium,

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

we multiply both sides by the spherical volume element \(4\pi r^3 dr\) and integrate from the center (\(r=0\)) to the surface (\(r=R\)):

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

—

Left-hand side (pressure term):

Integrate by parts: $$ \int_0^R 4\pi r^3 dP = \Big[4\pi r^3 P\Big]_0^R - \int_0^R 12\pi r^2 P\, dr. $$

At the surface \(P(R)\approx 0\), and at the center \(r=0\), the first term vanishes, leaving: $$ \int_0^R 4\pi r^3 dP \approx -3\!\int_0^R 4\pi r^2 P\, dr = -3\!\int P\, dV. $$

—

Right-hand side (gravitational term):

Using \(dM = 4\pi r^2 \rho\, dr\), $$ -\int_0^R 4\pi r^3 \rho \frac{G M}{r^2} dr = -\int_0^R \frac{G M(r)}{r}\, dM(r). $$

But this is precisely the gravitational potential energy of the system, $$ \Omega = -\int_0^R \frac{G M(r)}{r}\, dM(r). $$

—

Equating both sides gives the pressure–gravity balance:

\begin{equation} 3\!\int P\, dV + \Omega = 0. \end{equation}

For an ideal gas, where the internal (thermal) kinetic energy is related to pressure by \( K = \tfrac{3}{2}\!\int P\, dV \), this can be written as:

\begin{equation} 2K + \Omega = 0. \end{equation}

This is the virial theorem — it expresses the global equilibrium between the total kinetic energy and the gravitational potential energy of a bound system.

—

A simple example of the virial theorem in action is a satellite orbiting a planet of mass \(M\) at radius \(r\). Its kinetic and potential energies are:

$$ E_k = \tfrac{1}{2} m v^2 = \tfrac{1}{2}\frac{G M m}{r}, \qquad E_p = -\frac{G M m}{r}. $$

Hence,

\begin{equation} E_k = -\tfrac{1}{2} E_p, \end{equation}

which satisfies \(2K + \Omega = 0\).

—

Now imagine a star with no internal nuclear energy source. It behaves similarly: as it radiates energy from its surface, its total energy \(E = K + \Omega\) decreases. To restore equilibrium, the star contracts, making \(\Omega\) more negative and \(K\) larger, because energy conservation requires that the potential energy released during contraction appears as heat.

As the star shrinks:

• \(E_p\) becomes more negative (stronger gravity),
• \(E_k\) increases (particles move faster),
• and the temperature rises.

Thus, when the star loses energy, it actually gets hotter — a hallmark of negative specific heat. Removing energy makes it hotter because more gravitational binding energy is converted into kinetic energy. This leads to a quasi-stable equilibrium during stellar formation or collapse.

During stellar birth, the collapsing gas cloud spins faster and heats up until internal pressure and rotation counteract gravity. During stellar death, as nuclear fuel runs out, the core again contracts and heats, triggering fusion of heavier elements until no further fusion is possible.

When nuclear reactions in the core produce energy at exactly the rate it is lost from the surface, the star remains in true equilibrium — this is hydrostatic equilibrium, maintained by the virial balance.

—

The virial theorem can also be used to estimate the total mass of a galaxy cluster from the observed velocities of its member galaxies.

For a system of \(N\) galaxies, each of mass \(m\), the virial condition is:

$$ 2 \sum_i \frac{1}{2} m v_i^2 - \sum_{i\neq j} \frac{G m^2}{r_{ij}} = 0. $$

Rewriting and simplifying for identical masses: $$ m \sum_i v_i^2 - G m^2 \sum_{i\neq j} \frac{1}{r_{ij}} = 0. $$

Multiply the first term by \(N/N\) and the second term by \(N^2/N^2\):

$$ N m \left[\frac{1}{N}\sum_i v_i^2 \right] - G \frac{(N m)^2}{2} \left[\frac{1}{N(N-1)/2} \sum_{i\neq j} \frac{1}{r_{ij}}\right] = 0. $$

If \(N(N-1)\approx N^2\) and \(M = N m\) is the total mass, we obtain:

\begin{equation} M \langle v_i^2 \rangle_{av} - \frac{G M^2}{2} \langle r_{ij}^{-1} \rangle_{av} = 0. \end{equation}

Hence, the virial mass of the cluster is:

\begin{equation} M = \frac{2 \langle v_i^2 \rangle_{av}}{G \langle r_{ij}^{-1} \rangle_{av}}. \end{equation}

—

In practice, we observe only the line-of-sight velocity component \(v_{i,\text{los}}\) via Doppler shift. Assuming the velocity distribution is isotropic:

$$ \langle v_i^2 \rangle_{av} = 3 \langle v_{i,\text{los}}^2 \rangle_{av}. $$

Thus, using the virial theorem, one can infer the total (including dark) mass of a cluster from measured velocity dispersions and projected separations between galaxies.

The virial theorem states that in any self-gravitating system in equilibrium, the total kinetic energy equals half the magnitude of the (negative) potential energy. It connects local pressure balance to global energetics.

• In stars, it links thermal pressure to gravitational binding energy: \(3\int P\,dV + \Omega = 0\).
• In galaxy clusters, it relates velocity dispersion to total mass, including unseen dark matter.
• It explains why stars heat up as they lose energy — a consequence of negative specific heat.
• It underlies stellar stability, cluster dynamics, and the estimation of galaxy masses through velocity dispersion.

1. Derive the step \(3\int P\, dV + \Omega = 0\) from the hydrostatic equilibrium equation, showing the integration by parts.
2. Explain why a system obeying \(2K + \Omega = 0\) can have a negative heat capacity (losing energy makes it hotter).
3. Using the virial theorem, derive the expression for the virial mass \(M = \frac{2\langle v^2 \rangle}{G \langle r^{-1}\rangle}\) and describe how \(v_{\text{los}}\) measurements yield total cluster mass.