Gravitational Collapse

Hydrostatic equilibrium, and how it ends

Every star is a balance between two forces. Gravity pulls every gram inward with acceleration GM(r)/r²; pressure pushes outward with a gradient dP/dr. Set them equal and you get the equation of hydrostatic equilibrium

dP/dr = − G M(r) ρ(r) / r²

A star sits where this equation is satisfied. The pressure is supplied by some microphysics — thermal pressure of an ideal gas, radiation pressure, electron degeneracy pressure, neutron degeneracy pressure, repulsion from the nuclear strong force. As long as some mechanism can deliver dP/dr at the right value, the configuration is stable.

Lose the mechanism and the balance breaks. Acceleration becomes GM(r)/r² inward, and the configuration falls. There is no slow leak; once the pressure support is gone the timescale is the free-fall time, which for ordinary stellar densities is minutes to hours. Gravitational collapse is therefore a runaway, not a creep.

Free-fall time and order-of-magnitude estimates

For an unconfined self-gravitating sphere of mean density ρ̄, the time for an outer shell to fall to the centre under pure Newtonian gravity is

τ_ff = (3π / (32 G ρ̄))^(1/2)
     ≈ √(R³ / (G M))      (within numerical factors of order unity).

Plug in a typical pre-collapse iron core: M = 10 M_☉, R = 10⁶ km = 10⁹ m. Mean density ρ̄ ≈ 4.7 × 10⁹ kg/m³.

τ_ff ≈ √(R³ / (G M))
    = √(10²⁷ / (6.674 × 10⁻¹¹ × 1.989 × 10³¹))
    = √(10²⁷ / 1.33 × 10²¹)
    = √(7.55 × 10⁵)
    ≈ 870 s
    ≈ 0.24 hr ≈ 14 min          (with the more accurate factor)
    ≈ 0.5 hr                     (the often-quoted half-hour order of magnitude).

The "half-hour" figure cited at the top is a round-number rule of thumb; the precise answer depends on the chosen radius (whether one means the iron-core radius or the original red-giant envelope) and on the (non-trivial) inward pull of the rest of the star. The headline lesson: stellar-mass collapse, once it starts, takes minutes to hours — not human-engineering timescales.

The Chandrasekhar limit — when collapse stops at a white dwarf

In 1931 Subrahmanyan Chandrasekhar, then a 19-year-old student on a steamship to England, computed the maximum mass of an electron-degenerate stellar remnant. Below this mass, the Pauli exclusion principle's degeneracy pressure of cold electrons in an ionised plasma can balance gravity. Above it, even relativistic degenerate electrons cannot. The result:

M_Ch = ω₃⁰ √(3π/2) · (ℏc/G)^(3/2) / (μ_e m_H)²  ≈ 1.456 (2/μ_e)² M_☉

For μ_e = 2 (typical for carbon-oxygen white dwarfs) the limit is ~ 1.4 M_☉. White dwarfs heavier than this collapse; lighter ones simply cool over Gyr without further collapse.

Type Ia supernovae are the cinema of this limit. A carbon-oxygen white dwarf in a close binary accretes mass from a companion. When its total mass approaches M_Ch, runaway carbon fusion ignites at the centre, deflagrates outward, and the entire star is unbound in a thermonuclear explosion — the bright Type Ia event used as a cosmological standard candle. The collapse-vs-explosion distinction here is critical: most Type Ia events do not produce a remnant, they vapourise the white dwarf.

The TOV limit — when collapse stops at a neutron star

Above M_Ch, electrons capture onto protons (p + e⁻ → n + ν_e), the gas becomes neutron-rich, and at densities of order 10¹⁷ kg/m³ the system stabilises on the degeneracy pressure of neutrons combined with the repulsive component of the nuclear strong force. The maximum mass of this configuration is the Tolman-Oppenheimer-Volkoff (TOV) limit. The exact value depends on the equation of state of nuclear matter at supranuclear density — a notoriously difficult quantum-chromodynamics problem.

Modern observations narrow the range:

• Pulsar mass measurements. Shapiro-delay-determined pulsar masses have produced solid measurements of 1.9–2.2 M_☉ neutron stars (PSR J1614-2230, PSR J0740+6620), proving the TOV limit lies above 2 M_☉.
• GW170817 (2017). The first observed binary neutron-star merger and its kilonova afterglow constrained the equation of state through the tidal deformability of the inspiralling neutron stars. Combined with the absence of a long-lived remnant (the post-merger object collapsed to a black hole within seconds), this anchors M_TOV ≈ 2.16 M_☉ in slow rotation; up to ~ 2.5 M_☉ with maximal rotation.
• NICER X-ray timing. The Neutron-star Interior Composition Explorer measures pulsar radii to ~ 10% accuracy, providing additional equation-of-state leverage.

The narrow window between Chandrasekhar (1.4 M_☉) and TOV (~ 2.16 M_☉) is where neutron stars exist as cold remnants. Outside it, the available physics cannot stop collapse.

Above the TOV limit — black-hole formation

If the post-collapse remnant exceeds the TOV mass, no quantum-mechanical mechanism provides enough pressure. The configuration collapses past its Schwarzschild radius r_s = 2GM/c² in roughly the free-fall time of the original core. The Penrose singularity theorem (1965, awarded the 2020 Nobel Prize) shows that, given mild positive-energy conditions, the formation of a closed trapped surface implies a curvature singularity inside the event horizon — a genuine prediction of general relativity, not a feature of the specific Schwarzschild solution.

In astrophysical practice, black holes form when a stellar progenitor with initial main-sequence mass > ~ 25 M_☉ collapses and either:

• Direct collapse / failed supernova. The shock fails to revive; the entire core falls into a black hole over seconds. Inferred from the disappearance of N6946-BH1 in 2009 — a 25 M_☉ red supergiant that vanished without a bright explosion.
• Fallback supernova. A weak supernova ejects the outer envelope, but a significant fraction of the inner ejecta falls back onto the proto-neutron star, pushing it past M_TOV and triggering collapse to a black hole on a timescale of tens of seconds.
• Pair-instability collapse. For very massive stars (helium core 65–135 M_☉), e+e− pair production drains photon pressure, the star pulsates and ultimately collapses; if the helium core exceeds 135 M_☉ a single black hole of mass 50–250 M_☉ forms directly.

Worked example: collapse of a 1.4 M_☉ iron core

Take a pre-collapse iron core: M = 1.4 M_☉, initial R_0 ≈ 3000 km, final R_NS ≈ 12 km. Compute infall time, kinetic energy released, and the resulting neutron-star compactness.

Free-fall time:
τ_ff ≈ √(R₀³ / (G M))
     = √((3 × 10⁶)³ / (6.67e−11 · 2.78e30))
     = √(2.7e19 / 1.85e20)
     = √(0.146)
     ≈ 0.38 s

Gravitational binding energy:
E_bind ≈ G M² / R_NS
       = 6.67e−11 · (2.78e30)² / 1.2e4
       = 6.67e−11 · 7.73e60 / 1.2e4
       = 4.3e46 J.

99% goes into neutrinos, ~1% into kinetic energy of ejecta:
E_ν       ≈ 4.3 × 10⁴⁶ J.   (in 10s of seconds, all flavours)
E_kinetic ≈ 4 × 10⁴⁴ J.     (a Type II supernova "Bethe", 10⁵¹ erg.)

Compactness:
GM / (R_NS c²) ≈ 0.17        (NS sits at 17% of Schwarzschild — relativistic).

The half-second figure is the dynamical timescale of the actual core collapse — the iron core falls from a few thousand km to roughly 100 km in less than a second, at which point neutrino pressure briefly halts the bounce and launches the supernova shock. SN 1987A — the closest naked-eye supernova since 1604, in the Large Magellanic Cloud — was the first confirmation: Kamiokande, IMB, and Baksan detected 24 neutrinos within 13 seconds, three hours before the optical brightening, consistent with this picture.

Endpoint vs progenitor mass

Initial M (main-seq)	Remnant M	Type	Endpoint	Notable example
0.08 – 8 M_☉	0.5 – 1.3 M_☉	Planetary nebula	White dwarf (CO/ONe)	Sirius B, 1.02 M_☉
8 – 20 M_☉	1.2 – 1.8 M_☉	Type II / Ib/c SN	Neutron star	Crab pulsar, 1.4 M_☉
20 – 25 M_☉	1.8 – 2.5 M_☉	Type II SN + fallback	Heavy NS or BH	PSR J0740+6620, 2.08 M_☉
25 – 65 M_☉	3 – 15 M_☉	Failed / fallback SN	Stellar black hole	Cyg X-1, ~21 M_☉
65 – 135 M_☉ (He core)	0 (pulsation)	Pulsational pair-instability	No remnant (mass gap)	BH "mass gap" 50–135 M_☉
> 135 M_☉ (He core)	50 – 260 M_☉	Direct collapse	Intermediate-mass BH	GW190521 (85+66→142 M_☉)
Z = 0 Pop III stars	10² – 10³ M_☉	Direct collapse	SMBH seed	JWST high-z AGN

Famous gravitational-collapse events

• SN 1054 (Crab Nebula, 1054 AD). Recorded by Chinese astronomers; produced the Crab pulsar, a textbook young neutron star with P = 33 ms, B = 4 × 10⁸ T. Visible during daytime for 23 days.
• SN 1987A. In the LMC, 168 kly away. First neutrino detection from a stellar collapse: 24 events in three detectors, 13 s burst. The compact remnant was finally directly imaged by JWST in 2024 — a neutron star, not a black hole.
• N6946-BH1 (2009 disappearance). A 25 M_☉ red supergiant in NGC 6946 brightened briefly in 2009 and disappeared in subsequent imaging — a failed supernova, the cleanest observed example of direct black-hole formation.
• GW150914 (2015). First detected gravitational-wave event. A binary black-hole merger of 36 + 29 M_☉ stellar-mass holes, themselves the products of earlier gravitational collapses.
• GW170817 (2017). First binary neutron-star merger detected in GW + EM; the post-merger object collapsed to a black hole within seconds, anchoring the TOV limit.
• GW190521 (2019). A 85 + 66 → 142 M_☉ binary black-hole merger; primary mass falls in the pair-instability mass gap, hinting at hierarchical mergers or unusual collapse channels.

Variants and other contexts

• Jeans collapse. A gas cloud of mass M and radius R collapses if its thermal pressure cannot resist self-gravity. The Jeans mass M_J ≈ (k_B T / (G m_H))^(3/2) / √ρ separates collapsing from stable clouds. Star formation in molecular clouds proceeds by sequential Jeans-mass fragmentation.
• Cluster core collapse. A globular cluster reaches gravothermal collapse on a relaxation timescale (~ 10⁹ yr), when its central region develops negative heat capacity and contracts; halted by binary-binary scattering and stellar-mass-segregated formation.
• Galaxy formation. Cold dark matter haloes collapse from primordial density fluctuations; baryons cool inside and form stars. The Press-Schechter formalism (1974) computes the halo mass function from collapsed-fraction statistics.
• Primordial black holes. Sufficiently dense regions in the very early universe (radiation-dominated era) collapse directly to black holes if their density contrast exceeds δ_c ≈ 0.45. Constraints from microlensing, CMB, evaporation, and lensing of GRBs bound the PBH abundance.
• Cooling-flow clusters. The cores of relaxed galaxy clusters lose energy via X-ray emission; pure cooling would lead to runaway collapse and massive starburst at the cluster centre. Active-galactic-nucleus feedback (Perseus, M87) is the suppression mechanism.

Where gravitational collapse shows up

• Type II/Ib/Ic supernovae. Galactic SN rates ≈ 1–3 per century. Crab (1054), Cas A (~1680), SN 1987A, SN 2023ixf are recent benchmarks.
• Neutron star formation. ~10⁹ neutron stars in the Milky Way today; about 3000 pulsars catalogued; ~ 50 X-ray binaries with confirmed NS accretors.
• Stellar-mass black holes. ~10⁷ in the Milky Way (theoretical); ~ 30 confirmed via X-ray binaries (Cyg X-1, GRS 1915+105, V404 Cyg); > 200 confirmed via LIGO/Virgo gravitational-wave detections.
• Direct-collapse SMBH seeds. Required to explain z > 7 quasars (e.g., J1342+0928 at z = 7.54 with M = 8 × 10⁸ M_☉) which had only 700 Myr after the Big Bang to grow. Atomically cooled gas in pristine haloes may collapse directly to 10⁴–10⁵ M_☉ seeds.
• Star formation. Jeans-collapse fragmentation in molecular clouds produces 1–2 M_☉ stars in the Milky Way at ~ 1 M_☉/yr globally.

Common pitfalls

• Equating collapse with supernova. They are correlated but distinct. Core collapse to a neutron star drives the supernova explosion; direct collapse to a black hole can occur with no electromagnetic signature at all (failed supernova).
• Treating the Chandrasekhar limit as a "death zone". It is a stability limit for the cold electron-degenerate state. Hot or rotating white dwarfs can exceed M_Ch transiently; SN Ia progenitors do reach it.
• Ignoring rotation. A rotating collapsing star has angular momentum that can prevent direct collapse, force an accretion disk to form (BH binary inspirals, GRBs), or break the configuration into a binary. Spin parameter a/M up to ~ 0.7 is observationally common.
• Forgetting the neutrino burst. 99% of the gravitational binding energy released in a core-collapse SN is carried away by neutrinos, not photons. SN 1987A's 24 detected neutrinos hours before the optical brightening proved the mechanism in real time.
• Treating "1.4 M_☉" as the neutron-star mass. 1.4 is the canonical low-mass value; the observed distribution extends to ~ 2.1 M_☉ (PSR J0740+6620) and is determined by the equation of state, not by collapse physics alone.