Neutron Star and Black Hole Formation

The last hours of a massive star are a slow accounting that ends in sudden bankruptcy. Advanced burning has stacked the star into an onion of shells with an iron core, and iron is the end of the line: ${}^{56}\mathrm{Fe}$ sits at the peak of the binding-energy-per-nucleon curve, so no fusion can extract further energy from it. The inert core is now a white dwarf in all but name, held by relativistic electron degeneracy pressure and creeping toward the Chandrasekhar mass. Two thieves — photodisintegration, which shatters iron back into helium and free nucleons, and electron capture, which swallows the electrons providing pressure — quietly lower the threshold until the core is over it. Hydrostatic equilibrium fails, and roughly a solar mass falls almost freely inward in less than a second. At nuclear saturation density the strong force turns repulsive and the equation of state stiffens; the inner core overshoots, rebounds, and launches a shock. That prompt shock stalls within milliseconds, and the modern picture is that a fraction of the enormous neutrino luminosity from the newborn proto-neutron star is reabsorbed behind it and — helped by convection and hydrodynamic instabilities — revives it to unbind the envelope. About $3\times10^{53}\,\mathrm{erg}$ of gravitational binding energy is released, and $\sim99\%$ leaves as neutrinos. What is left behind is a compact remnant, and which one is decided by a single competition against the Tolman–Oppenheimer–Volkoff maximum mass. Below it, nucleon degeneracy and nuclear forces win and the object cools into a neutron star: $\sim1.4\,M_\odot$ packed into $\sim12\,\mathrm{km}$ , denser than an atomic nucleus. Above it — by an intrinsically heavy core or by fallback accretion — nothing can resist gravity, and the remnant collapses through its own event horizon to a stellar-mass black hole. This is where stellar evolution, nuclear physics, and general relativity meet, and the birthplace of the pulsars and gravitational-wave sources we observe.

Key takeaways

When a massive star ( $M_{\rm ZAMS}\gtrsim8\,M_\odot$ ) builds an iron core, further fusion is endothermic; the core loses pressure support and collapses on a millisecond timescale.
Photodisintegration of iron and electron capture accelerate the implosion until the inner core hits nuclear density $\rho_{\rm nuc}\approx2.7\times10^{14}\,\mathrm{g\,cm^{-3}}$ , where the stiff nuclear equation of state halts it and launches a bounce shock.
The gravitational binding energy released, $\sim3\times10^{53}\,\mathrm{erg}$ , escapes almost entirely as neutrinos; only $\sim10^{51}\,\mathrm{erg}$ becomes the visible explosion.
The remnant is a neutron star if the proto-neutron star stays below the TOV maximum ( $M_{\rm TOV}\approx2.0\text{–}2.3\,M_\odot$ ); above it, or with heavy fallback, collapse continues to a stellar-mass black hole.
Neutron stars ( $R\sim11\text{–}13\,\mathrm{km}$ ) are held by neutron degeneracy plus repulsive nuclear forces; a black hole is pure spacetime with horizon $r_s=2GM/c^2\approx2.95\,\mathrm{km}\,(M/M_\odot)$ .
Key idea: Core collapse is a gravitational implosion halted — or not — by degeneracy and nuclear forces; below the TOV mass the remnant is a neutron star, above it a black hole, and the released binding energy leaves overwhelmingly as neutrinos.

The core insight

Core collapse is a gravitational implosion halted (or not) by degeneracy and nuclear forces: below the TOV mass the remnant is a neutron star supported by nucleon degeneracy at nuclear density; above it, gravity wins and forms a black hole. The released binding energy escapes overwhelmingly as neutrinos.

Where this lives today

Compact-remnant formation is now an observational science on three live frontiers. Gravitational-wave astronomy is the most spectacular: since GW150914 the LIGO–Virgo–KAGRA network has cataloged the masses of black holes ( $\sim5\text{–}100\,M_\odot$ ) and neutron stars from core collapse, and GW170817 — a binary neutron-star merger seen also in light — measured the tidal deformability of nuclear matter, constraining the equation of state that sets the TOV limit. Radio pulsar timing is complementary: Shapiro-delay masses like PSR J0740+6620 at $\sim2.08\,M_\odot$ push the observed maximum toward the ceiling, while NASA's NICER mission pins neutron-star radii on the mass–radius curve. The open questions stay sharp: the mass gap between the heaviest neutron stars ( $\sim2\,M_\odot$ ) and lightest black holes ( $\sim5\,M_\odot$ ), and the fraction of massive stars that collapse quietly — the "unnova" channel, hinted at by disappearing red supergiants.

The physical picture

Picture the iron core as a weight on a spring — the spring being relativistic electron degeneracy pressure. For an ultrarelativistic gas the pressure grows as $P\propto\rho^{4/3}$ , exactly the power gravity demands, so the spring can never gain on gravity once the load exceeds a critical value: the Chandrasekhar mass. Now let two processes soften the spring — photodisintegration absorbs thermal energy by tearing iron apart, and electron capture removes the degenerate electrons outright, lowering $Y_e$ and with it the supportable mass ( $M_{\rm Ch}\propto Y_e^2$ ). The core slips past its limit and falls; the first figure traces that fall.

Figure

Anatomy of core collapse in three acts: an inert iron core near the Chandrasekhar mass loses pressure support and free-falls at v ≈ 0.2–0.3 c, the inner core rebounds at nuclear density ρ_nuc = 2.7×10¹⁴ g cm⁻³ launching a shock, and the newborn proto-neutron star radiates E_bind ≈ 3×10⁵³ erg — about 99% of it — as neutrinos

The energy scale is set by gravity alone. Shrinking $1.4\,M_\odot$ from an Earth-sized core to a $\sim12\,\mathrm{km}$ ball releases $\sim GM^2/R\sim3\times10^{53}\,\mathrm{erg}$ — about a tenth of the remnant's rest-mass energy — and since even neutrinos are briefly trapped, this budget can only leak out by neutrino diffusion over $\sim10\,\mathrm{s}$ . Whether the proto-neutron star survives or crosses its horizon is decided by one number — its mass against the TOV limit — which the second figure plots on a mass–radius plane with the black-hole boundary.

Figure

The compact-object mass–radius plane: neutron stars occupy a narrow band near R ≈ 11–13 km rising to the TOV maximum M ≈ 2.2 M⊙, while the Schwarzschild line R = 2GM/c² ≈ 2.95 km per M⊙ bounds the black-hole region that no static star can enter

The core equations

We take the equations in order: a critical mass, the two agents that lower it, the energy released, the relativistic limit on the remnant, and the timescales that make it all happen in a heartbeat. The iron core is a white dwarf in disguise, held by relativistically degenerate electrons — an $n=3$ polytrope whose mass is fixed by fundamental constants and the electron fraction alone.

M_{\rm Ch} = \frac{\omega_3^0\sqrt{3\pi}}{2}\left(\frac{\hbar c}{G}\right)^{3/2}\frac{1}{(\mu_e m_H)^2}\approx\frac{5.83}{\mu_e^2}\,M_\odot \qquad (1)

where: $\mu_e$ — mean molecular weight per electron ( $\approx2$ for iron-group matter); $m_H$ — hydrogen atom (nucleon) mass; $\omega_3^0\approx2.018$ — Lane–Emden constant for $n=3$ ; $\hbar,c,G$ — reduced Planck constant, speed of light, gravitational constant. Reading it: the mass depends on $\mu_e$ but not on radius, so it is a hard ceiling — for $\mu_e=2$ , $M_{\rm Ch}\approx1.44\,M_\odot$ . At $T\gtrsim5\times10^9\,\mathrm{K}$ the blackbody photon bath holds enough gamma rays above threshold to undo the nuclear burning, tearing iron back into alpha particles and neutrons — absorbing energy and robbing the core of pressure.

{}^{56}_{26}\mathrm{Fe}+\gamma\rightarrow13\,{}^{4}_{2}\mathrm{He}+4\,n,\qquad Q=-124.4\,\mathrm{MeV} \qquad (2)

where: $\gamma$ — high-energy thermal photon in the $\sim10^{10}\,\mathrm{K}$ core; $Q$ — reaction energy, negative because energy is absorbed. Reading it: the reaction softens the gas when the core most needs pressure, and it runs away — contraction raises $T$ , making more threshold photons. The second thief removes the electrons themselves: degenerate electrons capture onto protons and neutron-rich nuclei, emitting neutrinos that carry off energy and lepton number.

e^-+p\rightarrow n+\nu_e,\qquad e^-+(Z,A)\rightarrow(Z-1,A)+\nu_e \qquad (3)

where: $e^-$ — degenerate electron providing pressure; $\nu_e$ — electron neutrino, escaping until trapping sets in; $(Z,A)$ — heavy nucleus becoming more neutron-rich. Reading it: since $M_{\rm Ch}\propto Y_e^2$ , dropping $Y_e$ from $\sim0.46$ toward $\sim0.3$ pushes the supportable mass below the core mass — the trigger for collapse. The gravitational energy freed as the core shrinks to nuclear size fuels everything that follows; a uniform-sphere estimate already gives the right number.

E_{\rm bind}\simeq\frac{3}{5}\frac{GM^2}{R}\approx3\times10^{53}\,\mathrm{erg}\left(\frac{M}{1.4\,M_\odot}\right)^2\left(\frac{10\,\mathrm{km}}{R}\right) \qquad (4)

where: $M$ — remnant mass; $R$ — remnant radius ( $\sim11\text{–}13\,\mathrm{km}$ ). Reading it: this is $\sim0.1\,Mc^2$ , of which $\sim99\%$ departs as neutrinos, $\sim10^{51}\,\mathrm{erg}$ becomes explosion kinetic energy, and only $\sim10^{49}\,\mathrm{erg}$ is light. Whether the remnant survives is set by general relativity, where pressure itself gravitates — so adding pressure to resist collapse also deepens the well.

\frac{\mathrm{d}P}{\mathrm{d}r}=-\frac{G\left(\rho+P/c^2\right)\left(m+4\pi r^3P/c^2\right)}{r^2\left(1-2Gm/rc^2\right)} \qquad (5)

where: $P(r)$ — pressure; $\rho(r)$ — mass–energy density; $m(r)$ — enclosed mass–energy; $c$ — speed of light. Reading it: all three relativistic factors steepen the gradient relative to the Newtonian $\mathrm{d}P/\mathrm{d}r=-Gm\rho/r^2$ , so any realistic nuclear equation of state yields a maximum mass $M_{\rm TOV}\approx2.0\text{–}2.3\,M_\odot$ above which no static solution exists. If the remnant exceeds that maximum, it collapses inside its own event horizon at the Schwarzschild radius.

r_s=\frac{2GM}{c^2}\approx2.95\,\mathrm{km}\left(\frac{M}{M_\odot}\right) \qquad (6)

where: $r_s$ — event-horizon radius; $M$ — black-hole mass. Reading it: a $2\,M_\odot$ neutron star of radius $\sim12\,\mathrm{km}$ has $r_s\approx5.9\,\mathrm{km}$ , so $r_s/R\sim0.5$ — it sits within a factor of a few of being a black hole, which is why GR is essential here. Once support is gone the core falls on the free-fall timescale set by its mean density.

t_{\rm ff}=\sqrt{\frac{3\pi}{32\,G\bar{\rho}}} \qquad (7)

where: $\bar\rho$ — mean density of the collapsing core; $t_{\rm ff}$ — characteristic collapse time. Reading it: at iron-core density $\sim10^9\,\mathrm{g\,cm^{-3}}$ this is $\sim0.07\,\mathrm{s}$ , plunging to sub-millisecond values near nuclear density — too fast for burning to respond. Finally, the trapped binding energy leaks out as a diffusive neutrino glow.

L_\nu\approx\frac{E_{\rm bind}}{\tau_{\rm diff}}\sim\frac{3\times10^{53}\,\mathrm{erg}}{10\,\mathrm{s}}\sim3\times10^{52}\,\mathrm{erg\,s^{-1}} \qquad (8)

where: $L_\nu$ — total neutrino luminosity over all flavors; $\tau_{\rm diff}$ — neutrino diffusion (Kelvin–Helmholtz cooling) time $\sim10\,\mathrm{s}$ . Reading it: this briefly rivals the photon output of every star in the observable universe, and the $\sim20$ neutrinos from SN 1987A fixed the numerator.

Derivation AWhy the effective Chandrasekhar mass fails at collapse

The goal is to show that the iron core cannot avoid collapse once burning ceases, by tracking how its stability threshold moves relative to its actual mass. Electron capture drags the Chandrasekhar-type limit ( $M_{\rm Ch}\propto Y_e^2$ ) down while silicon shell burning pushes the mass up, so the two are guaranteed to cross.

Step 1. Start from the effective critical mass with a thermal correction, $M_{\rm Ch}^{\rm eff}=5.83\,Y_e^2\left[1+(s_e/\pi Y_e)^2\right]M_\odot$ , evaluated cold for fresh iron: $Y_e=26/56\approx0.464$ gives $M_{\rm Ch}\approx1.26\,M_\odot$ . Why: this fixes the baseline the core must beat; the entropy term $s_e$ only nudges it upward, so thermal support is a minor player.

Step 2. Turn on electron capture. As electrons are captured onto protons and heavy nuclei, $Y_e:0.46\rightarrow0.42\rightarrow\sim0.3$ , and because $M_{\rm Ch}^{\rm eff}\propto Y_e^2$ the supportable mass drops sharply toward $\sim1.0\text{–}1.1\,M_\odot$ . Why: deleptonization is the dominant destabilizer, attacking the pressure source directly.

Step 3. Grow the core. Silicon shell burning keeps dumping fresh iron on, raising the mass to $\sim1.3\text{–}1.5\,M_\odot$ — now above the shrinking $M_{\rm Ch}^{\rm eff}$ . Why: the mass rising while the threshold falls makes the crossing inevitable; collapse is scheduled, not chance.

Step 4. Seal it with photodisintegration. Contraction raises the temperature, endothermic iron photodisintegration switches on, and the last thermal pressure is drained. Why: once support is lost the free fall of equation (7) takes over and nothing slower can intervene in time. *Why the numbers work out: $M_{\rm Ch}=5.83/\mu_e^2\,M_\odot\approx1.44\,M_\odot$ for $\mu_e=2$ , and the cold-iron baseline $5.83\times0.464^2\approx1.26\,M_\odot$ confirms the $Y_e^2$ scaling that makes deleptonization so destructive.*

Go deeper · Derivation B · 1 minCore bounce and the stalled shock

Derivation BCore bounce and the stalled shock

The goal is to explain what stops a free-falling solar mass and why the shock it launches does not, by itself, blow up the star. The collapse splits into a subsonic homologous inner core and a supersonic outer core that lose causal contact; when the inner core hits nuclear saturation density the abruptly stiffening equation of state rebounds it, launching a shock that then stalls as it photodisintegrates iron and sheds neutrinos.

Step 1. Separate the flow. The inner $\sim0.6\,M_\odot$ collapses homologously ( $v\propto r$ ) and subsonically, while the outer core falls supersonically at $v\sim0.2\text{–}0.3\,c$ . Why: moving faster than its own sound speed, the outer core gets no pressure signal about the coming halt ahead of a shock — it does not know to stop.

Step 2. Stiffen the equation of state. At $\rho_{\rm nuc}\approx2.7\times10^{14}\,\mathrm{g\,cm^{-3}}$ the incompressibility jumps ( $P\propto\rho^\Gamma$ , $\Gamma>2$ ) as the repulsive core of the nuclear force engages. Why: this sudden hardening is the wall the inner core slams into, converting infall into a rebound.

Step 3. Overshoot and rebound. The inner core compresses past equilibrium to $\sim2\text{–}4\,\rho_{\rm nuc}$ , then springs back, driving a pressure wave into the supersonic infall where it steepens into a shock. Why: energy stored in overcompressed nuclear matter gives the shock its initial push.

Step 4. Watch it stall. Plowing outward, the shock spends $\sim8\,\mathrm{MeV}$ per nucleon dissociating the iron it sweeps up and loses more to neutrinos, halting at $\sim100\text{–}200\,\mathrm{km}$ within milliseconds. Why: this is why the prompt bounce fails and delayed neutrino heating is needed to revive the shock.

Go deeper · Derivation C · 1 minThe neutron-star radius from degeneracy pressure

Derivation CThe neutron-star radius from degeneracy pressure

The goal is to see why a neutron star is only $\sim10\,\mathrm{km}$ across — a factor $\sim1000$ smaller than a white dwarf of similar mass — using nothing more than white-dwarf physics with electrons swapped for neutrons. To first approximation it is one gigantic nucleus held by non-relativistic neutron degeneracy pressure, and the same polytropic balance, with $m_n$ replacing $\mu_e m_H$ , shrinks it by the mass ratio to a high power.

Step 1. Write the pressure. A non-relativistic degenerate neutron gas has $P\sim\hbar^2n^{5/3}/m_n$ with number density $n$ and neutron mass $m_n$ — an $n=3/2$ polytrope. Why: neutrons resist compression by the same Pauli exclusion as electrons, so white-dwarf structure carries over unchanged in form.

Step 2. Balance against gravity. Equating the polytropic pressure gradient to self-gravity gives $R\propto M^{-1/3}$ , identical in scaling to a white dwarf but with the internal mass factor now $m_n$ rather than $\mu_e m_H$ . Why: heavier degenerate particles pack into a smaller radius, so swapping electrons for neutrons ( $\sim2000\times$ heavier) shrinks the star.

Step 3. Put in numbers. The scaling $R\sim\hbar^2/(G m_n^{8/3}M^{1/3})\approx10\,\mathrm{km}\,(M/1.4\,M_\odot)^{-1/3}$ ; realistic nuclear forces stiffen the equation of state and raise this to $R\approx12\,\mathrm{km}$ . Why: the $m_n^{-8/3}$ factor is the $\sim(m_n/m_e)$ shrink from $\sim10^4$ to $\sim10\,\mathrm{km}$ , and the same forces lift the maximum mass toward TOV.

Limiting cases

The remnant type is decided by the progenitor mass and the core structure it builds:

Case	Result	Meaning
Neutron star	$M_{\rm core}\lesssim M_{\rm TOV}\approx2\text{–}2.3\,M_\odot$	Proto-neutron star cools to a stable neutron star; typical gravitational mass $\sim1.2\text{–}1.6\,M_\odot$ ; progenitors $\sim8\text{–}20\,M_\odot$
Black hole via fallback	$M_{\rm remnant}>M_{\rm TOV}$ after accretion	A neutron star forms briefly, accretes fallback, and collapses; weak or failed explosion; progenitors $\sim20\text{–}40\,M_\odot$
Direct collapse	$M_{\rm ZAMS}\gtrsim40\,M_\odot$ (low $Z$ )	Core too massive; collapse straight to a black hole with little or no explosion (an unnova)
Lower boundary	$M_{\rm ZAMS}\approx8\,M_\odot$	Below this, stars end as white dwarfs (or electron-capture supernovae near the edge), not iron-core collapse

Common trapThe bounce shock does not power the supernova, and the Chandrasekhar mass is not the TOV mass

Three traps recur. First, the prompt bounce shock does not blow up the star — it stalls within milliseconds, drained by photodisintegration and neutrino losses, and must be revived by delayed neutrino heating. Second, do not conflate the Chandrasekhar mass ( $\sim1.4\,M_\odot$ , the iron-core limit from relativistic electron degeneracy) with the TOV mass ( $\sim2\,M_\odot$ , the neutron-star limit from neutron degeneracy plus nuclear forces in GR); they come from different pressure sources. Third, the visible explosion ( $\sim10^{51}\,\mathrm{erg}$ ) is only about $1\%$ of the total budget ( $\sim10^{53}\,\mathrm{erg}$ ) — the supernova you see is the faint afterthought of a neutrino event.

Worked examples

1 · The energy budget of core collapse (fully worked)

A proto-neutron star of $M=1.4\,M_\odot$ settles to $R=10\,\mathrm{km}$ . Estimate the gravitational binding energy released, compare it to the rest-mass energy, and partition it into neutrinos, explosion, and light.

Step 1. Use the uniform-sphere estimate, equation (4): $E_{\rm bind}=\tfrac{3}{5}GM^2/R$ . With $M=1.4\times1.989\times10^{33}=2.79\times10^{33}\,\mathrm{g}$ , $R=10^6\,\mathrm{cm}$ , and $G=6.674\times10^{-8}$ , $E_{\rm bind}=0.6\times6.674\times10^{-8}\times(2.79\times10^{33})^2/10^6=\mathbf{3.1\times10^{53}\,erg}$ .

Step 2. Compare to rest mass: $Mc^2=2.79\times10^{33}\times(2.998\times10^{10})^2=2.5\times10^{54}\,\mathrm{erg}$ , so $E_{\rm bind}/Mc^2\approx\mathbf{0.12}$ — close to the full-GR value of $\sim0.1\,Mc^2$ ; collapse liberates roughly a tenth of the rest energy.

Step 3. Partition it: about $99\%$ ( $\sim3\times10^{53}\,\mathrm{erg}$ ) leaves as neutrinos over $\sim10\,\mathrm{s}$ ; $\sim10^{51}\,\mathrm{erg}$ becomes explosion kinetic energy; only $\sim10^{49}\,\mathrm{erg}$ emerges as light. Neutrinos carry ~10² times the visible explosion energy — a nearby core collapse is first of all a neutrino source.

2 · How fast is the fall?

Compare the free-fall time at the iron-core density $\bar\rho=10^9\,\mathrm{g\,cm^{-3}}$ with that at nuclear density $\bar\rho=2.7\times10^{14}\,\mathrm{g\,cm^{-3}}$ .

Step 1. At iron-core density, equation (7): $t_{\rm ff}=\sqrt{3\pi/(32G\bar\rho)}=\sqrt{9.42/(32\times6.674\times10^{-8}\times10^{9})}=\mathbf{0.066\,s}$ — under a tenth of a second.

Step 2. At nuclear density the same formula gives $t_{\rm ff}=\sqrt{9.42/(32\times6.674\times10^{-8}\times2.7\times10^{14})}=\mathbf{1.3\times10^{-4}\,s}$ — about $0.1\,\mathrm{ms}$ .

Step 3. Since $t_{\rm ff}\propto\bar\rho^{-1/2}$ , raising the density by $2.7\times10^5$ shortens the timescale by $\sim520$ . The bounce-and-shock drama plays out in milliseconds — far faster than any nuclear burning could respond, so collapse is a purely hydrodynamic-plus-gravity event.

3 · How relativistic is a neutron star?

Compute the Schwarzschild radius of a $2\,M_\odot$ neutron star and its compactness $r_s/R$ for $R=12\,\mathrm{km}$ .

Step 1. From equation (6), $r_s=2.95\,\mathrm{km}\times(M/M_\odot)=\mathbf{5.9\,km}$ , so the compactness is $r_s/R=5.9/12=\mathbf{0.49}$ — the star sits at roughly half its own horizon radius. A neutron star is not Newtonian-with-corrections: at $r_s/R\sim0.5$ general relativity shapes its structure and maximum mass at the tens-of-percent level, which is why the TOV equation (5), not Newtonian hydrostatics, sets the mass limit.

Physics of the portrait

The dividing line between a neutron star and a black hole was drawn on paper a quarter-century before either object was seen — by a man better known for something else entirely.

Origins · 1 minJ. Robert Oppenheimer and how this was discovered

J. Robert Oppenheimer

1904–1967 · neutron-star maximum mass and gravitational collapse

By the late 1930s the neutron was barely five years old and the neutron star was a speculation of Baade and Zwicky's, but J. Robert Oppenheimer — then building the theoretical-physics school at Berkeley — turned it into hard general relativity. In 1939, with his student George Volkoff and building on Richard Tolman's work, he integrated the equation now bearing all three names and found that a neutron star cannot be arbitrarily heavy: above a maximum mass (modern nuclear equations of state place it at $\sim2.0\text{–}2.3\,M_\odot$ ) no static solution exists. The same year, with Hartland Snyder, "On Continued Gravitational Contraction" described matter falling through what we now call an event horizon — the first modern account of black-hole formation, written before the word existed. Then the war took him to Los Alamos and he never returned; the 1939 papers were fully appreciated only after his death, once pulsars (1967) proved neutron stars real. Oppenheimer thus fixed the two poles of the subject — the mass above which a neutron star cannot exist, and the fate of anything heavier — leaving the surrounding cast (Chandrasekhar's electron-degeneracy limit, Baade and Zwicky's supernova link, Landau's neutron-core idea) to fill in the rest.

History

1931

— Chandrasekhar derives the maximum mass of a white dwarf

1934

— Baade and Zwicky propose that supernovae mark the birth of neutron stars

1939

— Oppenheimer and Volkoff compute the neutron-star maximum mass; Oppenheimer and Snyder model continued collapse to a black hole

1967

— Bell and Hewish detect the first pulsar, confirming neutron stars exist

1987

— Neutrinos from SN 1987A confirm proto-neutron-star formation and its energy budget

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 15–16 (the fate of massive stars; neutron stars and black holes) · Shapiro & Teukolsky, Black Holes, White Dwarfs, and Neutron Stars, Ch. 6, 9–11 · Bethe, Rev. Mod. Phys. 62, 801 (1990), "Supernova mechanisms" · Janka, Annu. Rev. Nucl. Part. Sci. 62, 407 (2012), "Explosion Mechanisms of Core-Collapse Supernovae" · Oppenheimer & Volkoff, Phys. Rev. 55, 374 (1939). Key numerical checks: the Chandrasekhar mass $5.83/\mu_e^2\,M_\odot\approx1.44\,M_\odot$ for $\mu_e=2$ ; the iron photodisintegration $Q=-124.4\,\mathrm{MeV}$ follows from binding energies as $13\times28.3-492.3$ ; the binding energy $\tfrac{3}{5}GM^2/R\approx3.1\times10^{53}\,\mathrm{erg}\approx0.12\,Mc^2$ ; the Schwarzschild radius $2GM/c^2\approx2.95\,\mathrm{km}\,(M/M_\odot)$ ; nuclear saturation density $2.7\times10^{14}\,\mathrm{g\,cm^{-3}}$ ; free-fall times $0.066\,\mathrm{s}$ at iron-core density and $1.3\times10^{-4}\,\mathrm{s}$ at nuclear density; and the TOV maximum $\sim2.0\text{–}2.3\,M_\odot$ .

Prerequisites. Advanced Burning Stages in Massive Stars · White Dwarfs and Electron Degeneracy Pressure · Nuclear Reaction Rates and the Gamow Peak · Helium Burning and the Triple-Alpha Process. Leads to. Pulsars and Neutron Star Observations · Compact Binary Mergers and Gravitational Waves · Type II Supernova Light Curves.

Ask the tutor: Ask me to walk through why lowering the electron fraction $Y_e$ shrinks the Chandrasekhar mass, or to trace a $25\,M_\odot$ star step by step from iron core to final remnant.

Why does the core collapse instead of settling into a bigger, cooler white dwarf?

Because its support has a hard ceiling that its mass exceeds. Relativistic electron degeneracy pressure grows as $P\propto\rho^{4/3}$ , the same power as gravity's demand, so beyond the Chandrasekhar mass no compression can restore balance. Worse, electron capture and photodisintegration actively remove pressure as the core contracts, lowering the limit even as the mass climbs. A white dwarf settles because it stays below its limit; an iron core, pushed above it and drained by two runaway processes, has no stable configuration and simply falls.

If the bounce shock stalls, what actually makes the supernova explode?

The delayed neutrino-heating mechanism. The prompt shock stalls within milliseconds because it spends about $8\,\mathrm{MeV}$ per nucleon breaking up the iron it sweeps up and loses more to neutrinos. But the newborn proto-neutron star radiates $\sim3\times10^{52}\,\mathrm{erg\,s^{-1}}$ in neutrinos, and a small fraction is reabsorbed in the hot, low-density region behind the stalled shock. Aided by convection that lets matter linger in the heating zone, this re-energizes the shock over a few hundred milliseconds until it breaks out and unbinds the envelope.

How can we possibly test the claim that 99% of the energy leaves as neutrinos?

We already did, once. On 23 February 1987 the water-Cherenkov detectors Kamiokande-II and IMB, with the Baksan scintillator, recorded about two dozen neutrinos over $\sim13$ seconds — hours before the optical brightening of SN 1987A in the Large Magellanic Cloud. The inferred total energy ( $\sim3\times10^{53}\,\mathrm{erg}$ ), the duration matching the Kelvin–Helmholtz cooling time, and the mean energy of $\sim10\text{–}20\,\mathrm{MeV}$ all matched core-collapse predictions — confirming that a proto-neutron star forms and cools by neutrino emission, and remaining the direct measurement of the energy budget.

What decides whether the remnant is a neutron star or a black hole?

A competition between the remnant's mass and the TOV maximum of $\sim2.0\text{–}2.3\,M_\odot$ . Below it, neutron degeneracy plus repulsive nuclear forces win and the proto-neutron star cools into a stable neutron star. If it is intrinsically heavier, or if fallback rains back and pushes it over the limit, no pressure can resist gravity and it collapses through its event horizon to a black hole — often with a weak or failed explosion. Progenitor mass sets the odds: $\sim8\text{–}20\,M_\odot$ stars tend to leave neutron stars, heavier ones black holes.

Why isn't the Chandrasekhar mass the same as the neutron-star maximum mass?

Because they are held up by different particles and different physics. The Chandrasekhar mass ( $\sim1.4\,M_\odot$ ) is the ceiling for relativistic electron degeneracy pressure and governs white dwarfs and iron cores. The TOV mass ( $\sim2\,M_\odot$ ) is the ceiling for neutron degeneracy stiffened by repulsive nuclear forces, computed in general relativity where pressure itself gravitates — the numbers differ because neutrons are far heavier than electrons and the nuclear force adds stiffness. Exceeding the Chandrasekhar mass triggers collapse; exceeding the TOV mass makes it end in a black hole rather than a neutron star.

Is a neutron star really just a giant nucleus?

As a first approximation, yes — it is bulk nuclear matter held together by gravity rather than by a single nucleus's surface tension, at densities near or above nuclear saturation ( $2.7\times10^{14}\,\mathrm{g\,cm^{-3}}$ ). That is why balancing neutron degeneracy pressure against gravity, exactly as for a white dwarf but with neutrons replacing electrons, already gives the right $\sim10\,\mathrm{km}$ radius. But the analogy is imperfect: a real neutron star has a crust of nuclei, a fluid neutron interior with a few percent protons, and possibly exotic phases (hyperons, quark matter) in the core — the uncertainties that mass and radius measurements aim to resolve.