---
title: Molecular Clouds and Gravitational Collapse
draft: false
module: The Interstellar Medium
moduleNumber: 6
lessonNumber: 2
order: 602
summary: >
  Stars form in cold, dense molecular clouds when self-gravity overcomes thermal
  and magnetic support. The virial theorem fixes the Jeans mass and length at which
  a clump becomes unstable, the free-fall time sets how fast it collapses, and a
  fragmentation cascade — cut off at a minimum mass by the onset of opacity — turns
  one cloud into a whole cluster, imprinting the stellar initial mass function.
topics: [The Interstellar Medium]
sources:
  - book: Carroll & Ostlie
    ref: "Ch. 12 — The Interstellar Medium and Star Formation; §12.2 The Formation of Protostars"
  - book: Maoz
    ref: "Ch. 5 — Star Formation and the Interstellar Medium"
---

Stars form only in the coldest, densest phase of the interstellar medium, the
**giant molecular clouds** — self-gravitating complexes of $10^5$ to $10^6\,M_\odot$
of molecular hydrogen, tens of parsecs across, with mean densities of
$10^2$–$10^3\ \text{cm}^{-3}$ and temperatures near $10\ \text{K}$. At these low
temperatures thermal pressure is feeble, and the question of whether a region
collapses is a contest between self-gravity and the several forms of support:
thermal pressure, turbulence, and magnetic fields. This lesson derives the
threshold for collapse from the virial theorem, the timescale on which an unstable
clump falls together, the fragmentation cascade that divides a cloud into many stars,
and the magnetic physics that regulates the whole process.

## Giant molecular clouds and their support

A molecular cloud is bound by its own gravity, but it does not collapse freely. The
balance of forces is captured by the **virial theorem**, which for a cloud in
equilibrium relates the internal kinetic (thermal plus turbulent) energy $K$, the
gravitational potential energy $U$, and, in general, the magnetic and surface-pressure
terms. In the simplest thermal-only form,

$$
2K + U = 0,
\qquad
K = \frac{3}{2}\,\frac{M}{\mu m_{\mathrm H}}\,k T,
\qquad
U = -\frac{3}{5}\,\frac{G M^2}{R},
$$

for a uniform sphere of mass $M$, radius $R$, and mean molecular weight $\mu$
(about $2.3$ for molecular gas with helium). Equilibrium is the boundary case: if the
gravitational term dominates, $|U| > 2K$, the cloud is unbound against its own gravity
and collapses; if thermal energy dominates, it disperses or is confined only by
external pressure. The collapse condition $|U| > 2K$ is the physical content of the
Jeans criterion.[^co-collapse]

$$
% caption: A clump sits in an effective potential; if thermal support exceeds the
% well it is held, but past the Jeans threshold gravity wins and it collapses.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
% potential well
\draw[very thick] (0.4,4.4) .. controls (2.4,4.2) and (2.9,0.8) .. (4.6,0.8)
  .. controls (6.3,0.8) and (6.8,4.2) .. (8.8,4.4);
% supported marble high on the wall
\fill[black] (1.55,3.55) circle (5pt);
\draw[->, black] (1.55,3.9) -- (1.55,4.6);
\node[black!70, anchor=south] at (1.55,4.6) {supported};
% collapsing marble at the bottom
\fill[acc] (4.6,0.8) circle (5pt);
\draw[->, acc] (4.6,1.25) -- (4.6,0.35);
\node[acc, anchor=north] at (4.6,0.3) {collapse};
% axis label
\node[black!70, anchor=west] at (8.9,4.4) {radius};
\node[black!70, anchor=east, rotate=90] at (0.15,2.6) {energy};
\end{tikzpicture}
$$

Real clouds also carry supersonic turbulence and threading magnetic fields, both of
which add support and postpone collapse. Including them, the full virial balance reads
$2K + 2K_{\text{turb}} + U + M_{\text{mag}} = 0$, and observed giant molecular clouds
sit close to virial equilibrium, neither collapsing wholesale nor flying apart. Only
localized dense cores, where the balance tips, actually form stars.

## The Jeans mass and Jeans length

Setting $|U| = 2K$ at the boundary and solving for the mass gives the minimum mass a
region of temperature $T$ and density $\rho$ must exceed to collapse. Eliminating the
radius through $R = (3M/4\pi\rho)^{1/3}$ yields the **Jeans mass**,

$$
M_J = \left( \frac{5 k T}{G \mu m_{\mathrm H}} \right)^{3/2}
      \left( \frac{3}{4\pi\rho} \right)^{1/2},
$$

and the corresponding length scale, the **Jeans length**, is the size below which a
perturbation is stabilized by pressure and above which it collapses,

$$
\lambda_J = \left( \frac{15 k T}{4\pi G \mu m_{\mathrm H}\,\rho} \right)^{1/2}
= c_s \left( \frac{\pi}{G\rho} \right)^{1/2},
\qquad
c_s = \left( \frac{k T}{\mu m_{\mathrm H}} \right)^{1/2},
$$

with $c_s$ the isothermal sound speed. The Jeans length is, up to a factor of order
unity, the distance a sound wave travels in a free-fall time: a perturbation collapses
only if gravity acts faster than pressure can respond across it.

> **Worked example.** For a molecular clump at $T = 10\ \text{K}$ and
> $n(\mathrm{H}_2) = 10^3\ \text{cm}^{-3}$ ($\rho \approx 4.7\times10^{-18}\
> \text{kg}\,\text{m}^{-3}$), the isothermal sound speed is
> $c_s \approx 0.19\ \text{km}\,\text{s}^{-1}$, giving
> $$
> \lambda_J \approx 0.6\ \text{pc},
> \qquad
> M_J \approx 8\ M_\odot.
> $$
> A clump above about eight solar masses at this density and temperature is
> gravitationally unstable. Because $M_J \propto T^{3/2}\rho^{-1/2}$, the Jeans mass
> falls as the gas gets colder or denser, so cooling and compression both promote
> collapse and fragmentation into ever-smaller pieces.

$$
% caption: The Jeans mass falls with increasing density and rises with temperature,
% so colder, denser gas fragments into progressively lower-mass pieces.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (0,0) -- (10.4,0) node[right, black!70] {log density};
\draw[->, black] (0,0) -- (0,5.4) node[above, black!70] {log Jeans mass};
% M_J ~ rho^{-1/2}: falling straight lines; warmer sits above cooler
\draw[very thick] (0.6,4.8) -- (9.8,1.4);
\node[anchor=south west] at (7.3,2.2) {T = 20 K};
\draw[very thick, densely dashed] (0.6,3.9) -- (9.8,0.5);
\node[anchor=north east] at (6.6,1.6) {T = 10 K};
\node[black, anchor=west, font=\scriptsize] at (2.2,4.4) {slope = one half};
\end{tikzpicture}
$$

## The free-fall time

Once a clump exceeds the Jeans mass and pressure support becomes irrelevant, it
collapses under gravity alone. The **free-fall time** is the time for a
pressure-free, uniform sphere to collapse to a point. Each mass shell falls under the
mass interior to it; integrating the equation of motion for a shell starting from rest
gives

$$
t_{\text{ff}} = \left( \frac{3\pi}{32\, G \rho} \right)^{1/2},
$$

which depends only on the mean density, not on the size of the cloud — all shells
arrive at the centre simultaneously for a uniform sphere. For the clump above,
$t_{\text{ff}} \approx 10^6\ \text{yr}$, short compared to the million-to-ten-million
year lifetimes of the clouds themselves. That collapse does not happen everywhere at
this rate is direct evidence that additional support — turbulence and magnetic fields
— slows the process; the observed star-formation rate of the Galaxy is roughly a
hundred times below what unimpeded free-fall of all molecular gas would give.

The density dependence has a runaway character. As a region collapses its density
rises, so $t_{\text{ff}}$ shrinks and the densest sub-regions collapse fastest. A
uniform cloud is therefore unstable to fragmentation: small overdensities run ahead of
the global collapse and separate out as independent collapsing centres.

## The fragmentation cascade and the minimum mass

While the gas stays cold and optically thin, it collapses **isothermally**: the
compressional heat is radiated away efficiently, so $T$ stays near $10\ \text{K}$ and
the density climbs. Along an isothermal collapse the Jeans mass drops as
$M_J \propto \rho^{-1/2}$, so a clump that was marginally unstable as a whole soon
contains many new Jeans masses. Each fragments in turn, and the process repeats:
a **hierarchical fragmentation cascade** subdivides the cloud into smaller and smaller
bound pieces, which is why molecular clouds form gravitationally bound star clusters
rather than single stars.[^maoz-sf]

$$
% caption: Isothermal collapse lowers the Jeans mass as density rises, so a cloud
% subdivides hierarchically into clumps, cores, and finally individual protostars.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
% level 1: cloud
\draw[very thick] (0.6,2.3) ellipse (0.95 and 1.3);
\node at (0.6,2.3) {cloud};
% level 2: clumps
\foreach \y in {1.2,2.3,3.4}
  \draw[thick] (3.4,\y) ellipse (0.55 and 0.42);
\node[anchor=south] at (3.4,4.0) {clumps};
% level 3: cores
\foreach \y in {0.7,1.5,2.3,3.1,3.9}
  \draw[thick] (6.0,\y) circle (0.28);
\node[anchor=south] at (6.0,4.4) {cores};
% level 4: protostars
\foreach \y in {0.6,1.2,1.8,2.4,3.0,3.6,4.2}
  \fill[black] (8.6,\y) circle (3pt);
\node[anchor=south] at (8.6,4.6) {protostars};
% connectors
\draw[->, black] (1.6,2.3) -- (2.75,2.3);
\draw[->, black] (4.0,2.3) -- (5.6,2.3);
\draw[->, black] (6.35,2.3) -- (8.3,2.3);
\end{tikzpicture}
$$

The cascade cannot continue indefinitely. As fragments grow dense they become opaque
to their own thermal radiation. When the collapse can no longer radiate away its
compressional energy, it turns **adiabatic**: the temperature rises with further
compression, $T \propto \rho^{\gamma-1}$, and now $M_J \propto \rho^{(3\gamma-4)/2}$
increases with density for $\gamma > 4/3$. Fragmentation halts at the **opacity limit**,
the density at which the gas first traps its own radiation. A short calculation of when
the radiated luminosity of a collapsing fragment falls below its compressional heating
rate gives a minimum fragment mass of order

$$
M_{\min} \sim 0.01\ M_\odot,
$$

a few Jupiter masses. This sets the low-mass end of star formation: objects below the
opacity-limited minimum cannot form by direct fragmentation, and the smallest stellar
fragments are near the boundary between the lowest-mass stars and brown dwarfs.

## The initial mass function

The distribution of masses with which stars are born, the **initial mass function**
(IMF), is the statistical outcome of the fragmentation cascade. Empirically it is a
declining power law in mass: many low-mass stars form for every high-mass one. Salpeter
first fit the number of stars per unit mass as

$$
\xi(m)\,\d m \propto m^{-2.35}\,\d m
\qquad (m \gtrsim 0.5\,M_\odot),
$$

a slope that holds well above half a solar mass. At lower masses the function flattens;
the Kroupa parametrization uses a broken power law, with slope $-2.3$ above
$0.5\,M_\odot$, $-1.3$ between $0.08$ and $0.5\,M_\odot$, and a further flattening
into the brown-dwarf regime below the hydrogen-burning limit at $0.08\,M_\odot$. The
turnover near a few tenths of a solar mass marks the characteristic mass of star
formation, plausibly tied to the Jeans mass at the density where cloud cores become
opaque.

$$
% caption: The initial mass function is a declining power law: a steep Salpeter
% slope above half a solar mass, flattening to a turnover near a few tenths.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
\draw[->, black] (0,0) -- (10.4,0) node[right, black!70] {log mass};
\draw[->, black] (0,0) -- (0,5.4) node[above, black!70] {log number};
% rising to a turnover then steep Salpeter decline
\draw[acc, very thick] (0.6,2.4) .. controls (1.8,3.6) and (2.6,4.2) .. (3.6,4.2)
  -- (9.6,0.7);
% turnover marker
\fill[black] (3.6,4.2) circle (2pt);
\node[anchor=south] at (3.6,4.3) {turnover};
% Salpeter slope label
\node[anchor=south west] at (6.4,2.1) {Salpeter slope};
% brown-dwarf region
\draw[black, densely dotted] (2.0,0) -- (2.0,4.0);
\node[black, anchor=north, font=\scriptsize] at (1.1,0.9) {brown dwarfs};
\node[black, anchor=north, font=\scriptsize] at (7.4,4.6) {massive stars rare};
\end{tikzpicture}
$$

## Magnetic support and ambipolar diffusion

Molecular clouds are threaded by magnetic fields of $10$–$100\ \mu\text{G}$, and the
field resists compression across its lines. The magnetic support of a cloud is set by
its **mass-to-flux ratio** $M/\Phi$, where $\Phi = B\,\pi R^2$ is the magnetic flux
through it. There is a critical value

$$
\left( \frac{M}{\Phi} \right)_{\text{crit}} = \frac{c_\Phi}{\sqrt{G}},
\qquad c_\Phi \approx 0.13,
$$

below which the magnetic tension can support any compression and the cloud is
**subcritical** (stable against collapse), and above which it is **supercritical** and
collapses regardless of field strength. Because flux is conserved as gas is compressed
along field lines, the mass-to-flux ratio is nearly invariant during collapse, so a
subcritical cloud stays subcritical.

The way a subcritical cloud forms stars is **ambipolar diffusion**. The magnetic field
is tied only to the ionized component; the neutral gas, which is the overwhelming
majority, feels the field only indirectly through ion–neutral collisions. In the
weakly ionized interior of a molecular cloud (ionization fraction $\sim 10^{-7}$,
maintained by cosmic rays), the neutrals slowly drift inward through the ions and the
field, concentrating mass at the centre while leaving the flux behind. Over a
characteristic time of $\sim 10\ \text{Myr}$ the central region becomes supercritical
and collapses dynamically. Ambipolar diffusion is the slow, rate-limiting step that
makes the observed star-formation efficiency low even though the free-fall time is
short.

$$
% caption: In a subcritical cloud the neutral gas drifts inward through the
% ion-anchored magnetic lines, concentrating mass until the core turns supercritical.
\begin{tikzpicture}[scale=1.0, font=\footnotesize]
\definecolor{acc}{HTML}{4A6FA5}
% vertical magnetic lines
\foreach \x in {1.2,2.6,4.0,5.4,6.8}
  \draw[black, thick] (\x,0.3) -- (\x,4.9);
\node[black, anchor=south] at (4.0,4.9) {magnetic lines};
% central dense core
\fill[acc!18] (4.0,2.6) circle (0.9);
\draw[acc, very thick] (4.0,2.6) circle (0.9);
\node[acc] at (4.0,2.6) {dense core};
% neutral drift arrows inward
\draw[->, thick] (1.5,2.6) -- (3.0,2.6);
\draw[->, thick] (6.5,2.6) -- (5.0,2.6);
\draw[->, thick] (4.0,0.7) -- (4.0,1.6);
\draw[->, thick] (4.0,4.5) -- (4.0,3.6);
\node[anchor=west] at (7.1,2.6) {neutrals drift in};
\end{tikzpicture}
$$

The Jeans mass and free-fall time fix when and how fast a cloud collapses, the
fragmentation cascade and its opacity-limited cutoff fix the range of masses produced,
and ambipolar diffusion sets the slow clock in magnetically supported clouds. What
happens to a single collapsing core once it becomes optically thick — the birth of a
protostar and its climb toward the main sequence — is the subject of the [next
lesson](/astrophysics-cosmology/ism-and-star-formation/protostars-and-the-pre-main-sequence).

[^co-collapse]: Carroll & Ostlie, §12.2 — The Formation of Protostars: the Jeans criterion, the Jeans mass and length, the free-fall time, and fragmentation.
[^maoz-sf]: Maoz, §5.3 — gravitational collapse, the Jeans instability, fragmentation, the opacity-limited minimum mass, and the initial mass function.