The Tolman-Oppenheimer-Volkoff Limit

The last wall before a black hole

A neutron star is already absurd. Take 1.4 times the mass of the Sun, squeeze it into a sphere the size of a city — about 24 kilometres across — and you reach densities of several times 10¹⁷ kilograms per cubic metre, denser than an atomic nucleus. Nothing holds it up except quantum mechanics: the Pauli exclusion principle forbids two neutrons from occupying the same quantum state, so when you try to pack them tighter they push back with neutron degeneracy pressure. The repulsive part of the strong nuclear force adds even more stiffness. Together they brace the star against a gravitational field tens of billions of times stronger than Earth's.

But this defiance has a ceiling. Pile on more mass — from a companion star, a merger, or fallback after a supernova — and you eventually reach a point where degeneracy pressure and the strong force, even working together, cannot generate enough outward push to balance gravity. There is no heavier, denser remnant left to become. The core simply falls in on itself. That ceiling is the Tolman-Oppenheimer-Volkoff limit, and crossing it is the most decisive event in the life of compact matter: it is the boundary between a star and a black hole.

The TOV limit is the neutron-star sibling of the Chandrasekhar limit. Where Chandrasekhar's 1.4 M☉ marks where a white dwarf's electrons can no longer hold it up, the TOV limit — somewhere around 2.2 to 2.3 solar masses — marks where a neutron star's neutrons can no longer hold it up. Both are maximum masses; both end with gravity winning.

The physics: the TOV equation

A neutron star is a relativistic object — its surface gravity bends light, and the spacetime curvature inside is severe — so Newtonian hydrostatic equilibrium is not good enough. The correct description comes from general relativity. In 1939, Richard Tolman wrote down the static, spherically symmetric solution to Einstein's field equations, and Robert Oppenheimer and George Volkoff applied it to a ball of degenerate neutrons. The governing relation, the Tolman-Oppenheimer-Volkoff equation of hydrostatic equilibrium, is:

dP/dr = − G [ ρ + P/c² ] [ m(r) + 4π r³ P/c² ]
        ────────────────────────────────────────
                 r² [ 1 − 2 G m(r) / (r c²) ]

with   dm/dr = 4π r² ρ

Compare this to the Newtonian version, dP/dr = − G m ρ / r². Every extra factor in the TOV equation makes gravity stronger, not weaker, at high density:

• ρ + P/c² — in relativity, pressure itself has mass-energy and therefore gravitates. The harder the star pushes back, the more it weighs.
• m + 4π r³ P/c² — pressure acts as an additional source of gravity, on top of the enclosed mass.
• 1 / (1 − 2Gm/rc²) — the denominator blows up as the radius approaches the Schwarzschild radius r_s = 2GM/c², amplifying the gravitational pull.

This is the cruel twist that makes a maximum mass inevitable. For a white dwarf or low-mass star, adding pressure helps fight gravity. For a neutron star near the limit, adding pressure makes the star heavier and more compact, which strengthens gravity faster than the pressure can resist. Beyond a critical central density there is no stable solution: every configuration is unstable to collapse.

To actually integrate the TOV equation you need one more ingredient — the equation of state, the relation P(ρ) between pressure and density for cold matter above nuclear density. That single function is where all the uncertainty lives.

Why the equation of state sets the number

Above the nuclear saturation density ρ₀ ≈ 2.7 × 10¹⁷ kg/m³ (the density inside ordinary nuclei), matter is in a regime no laboratory can reproduce. We do not know with certainty what it is made of: it may be neutrons and protons with a sprinkling of electrons and muons, or it may transition to hyperons, a kaon or pion condensate, or even deconfined quark matter. Each possibility gives a different stiffness.

The key property is stiffness: how rapidly pressure rises with density. A stiff equation of state generates a lot of pressure per unit density, supports more mass, and gives a larger, more massive maximum star. A soft equation of state gives less pressure, supports less mass, and lowers the TOV limit. Exotic phases (hyperons, quarks) generally soften the equation of state by opening new low-energy states for the particles to fall into, which is why the existence of 2-solar-mass pulsars is such a strong constraint — they rule out the softest models.

The result is that theory does not predict a single TOV mass but a range, roughly 2.0 to 2.9 M☉, with the favoured value near 2.2-2.3 M☉ once observations are folded in. There is, however, a firm upper bound from causality: the sound speed in the star cannot exceed the speed of light. The stiffest equation of state allowed by causality caps the non-rotating TOV mass at about 3 M☉ (specifically, no more than roughly 3.0-3.2 M☉ for plausible matching to known low-density physics). No neutron star, however exotic, can be heavier than that and still be a neutron star.

The key numbers

Quantity	Value	Notes
Original OV limit (1939)	0.71 M☉	Ideal free-neutron gas; no strong force
Modern TOV limit	≈ 2.2 – 2.3 M☉	With realistic stiff equation of state
Causality upper bound	≈ 3.0 M☉	Sound speed ≤ c
Heaviest measured pulsar	2.08 ± 0.07 M☉	PSR J0740+6620 (Shapiro delay)
Typical neutron star mass	1.1 – 1.5 M☉	Peaks near 1.4 M☉
Neutron star radius	11 – 13 km	NICER + GW constraints
Central density at TOV mass	~ 5 – 8 ρ₀	ρ₀ = 2.7 × 10¹⁷ kg/m³
Compactness GM/Rc²	~ 0.3	vs 0.5 at the horizon; ~3×10⁻⁶ for the Sun
Collapse timescale	< 1 ms	Free-fall across ~10 km

The compactness figure is worth dwelling on. The ratio GM/Rc² measures how close an object is to being a black hole; it is exactly 0.5 at the Schwarzschild radius. The Sun sits at about 2 × 10^-6; a white dwarf at about 3 × 10^-4; a neutron star near the TOV limit reaches about 0.3. A maximum-mass neutron star is, in a real sense, two-thirds of the way to being a black hole already. That is why such a small extra push — a few tenths of a solar mass — finishes the job.

How the limit is measured and pinned down

Because the TOV limit depends on physics we cannot compute from first principles, it has to be cornered observationally. Three independent lines of evidence converge on the same window.

1. Massive pulsars set a lower bound. The TOV mass must be at least as large as the heaviest neutron star we can weigh, because that star exists and has not collapsed. The most precise heavy-pulsar masses come from the Shapiro delay — the extra light-travel time as the pulsar's radio beam grazes the gravitational well of a white-dwarf companion. PSR J1614-2230 (1.97 M☉, 2010) was the first to definitively exceed 2 M☉ and kill off the softest equations of state. PSR J0740+6620, at 2.08 ± 0.07 M☉, is the current record, refined by NASA's NICER X-ray timing of its surface hot spots.

2. NICER measures radii. NASA's Neutron star Interior Composition Explorer, mounted on the International Space Station since 2017, watches the X-ray pulse profile of rotating hot spots warp as the star's own gravity bends the emitted light. Modelling that distortion yields both mass and radius simultaneously. For PSR J0740+6620 it found a radius near 12.4 km; a tight mass–radius pair directly constrains the equation of state and hence the TOV mass.

3. Gravitational waves set an upper bound. GW170817 (August 2017), the first binary neutron-star merger detected by LIGO and Virgo, had a total mass of about 2.7 M☉. The remnant did not promptly collapse — it briefly survived as a differentially rotating hypermassive star, evidenced by the kilonova AT2017gfo and its blue early emission — but it also did not survive as a stable star. Folding the merger dynamics, the tidal deformability measured from the inspiral waveform, and the lack of a long-lived remnant together brackets the non-rotating TOV mass to about 2.2-2.3 M☉. This is the single most powerful constraint to date.

Worked example: where does the OV limit come from?

You can recover the right order of magnitude for a degeneracy-supported maximum mass with a back-of-the-envelope argument — the same logic Chandrasekhar used. Balance the relativistic degeneracy energy of N fermions against the gravitational binding energy of a star of radius R.

For ultra-relativistic degenerate fermions, the Fermi energy per particle scales as E_F ~ ħc n^(1/3) ~ ħc N^(1/3) / R, so the total internal energy is

E_deg ~ N · ħc N^(1/3) / R = ħc N^(4/3) / R

The gravitational energy of a star of mass M = N m_n (with m_n the neutron mass) is

E_grav ~ − G M² / R = − G (N m_n)² / R

Both scale as 1/R, so the radius cancels: the star is in marginal balance only for one special particle number. Setting the magnitudes equal gives the maximum number of neutrons:

ħc N^(4/3) / R ~ G N² m_n² / R
→ N_max ~ ( ħc / G m_n² )^(3/2)
        ~ ( 1.7 × 10³⁸ )^(3/2)
        ~ 2.2 × 10⁵⁷ neutrons

Multiplying by the neutron mass m_n = 1.675 × 10⁻²⁷ kg:

M_max ~ N_max · m_n ~ 2.2 × 10⁵⁷ × 1.675 × 10⁻²⁷ kg
      ~ 3.7 × 10³⁰ kg
      ~ 1.8 M☉        (M☉ = 1.989 × 10³⁰ kg)

This dimensional estimate lands within a factor of two of the true answer — it captures the essential physics that the maximum mass is set entirely by the fundamental constants ħ, c, G, and the particle mass, with no free parameters. Order-unity prefactors (which the toy model drops) shift the precise number, and the full Oppenheimer-Volkoff integration of the relativistic equation, with an ideal-gas equation of state, gives 0.71 M☉. Adding the repulsive strong force — which the toy model omits — stiffens the matter and pushes the realistic value up to 2.2-2.3 M☉.

History: Tolman, Oppenheimer, Volkoff, and 1939

The story runs through a remarkable few years of 1930s physics. In 1930, on his ship voyage to England, the 19-year-old Subrahmanyan Chandrasekhar derived that a white dwarf supported by relativistic electron degeneracy has a maximum mass near 1.4 M☉. The natural question followed: what holds up a star heavier than that? The neutron had only just been discovered, by James Chadwick in 1932. In 1934, Walter Baade and Fritz Zwicky proposed that supernovae produce neutron stars — objects supported by degenerate neutrons.

In 1939, Richard Tolman, a physical chemist and relativist at Caltech, published the general-relativistic equations for a static fluid sphere. The same year, J. Robert Oppenheimer — later director of the Manhattan Project — and his student George Volkoff at Berkeley combined Tolman's equations with an ideal neutron-gas equation of state and integrated them numerically. Their landmark paper, "On Massive Neutron Cores," found a maximum mass of about 0.7 M☉. Oppenheimer, with Hartland Snyder, then published a companion paper the same year describing what happens when a star exceeds the limit: continued gravitational contraction — the first modern description of black-hole formation.

Their 0.7 M☉ figure was famously too low, because the ideal-gas model ignored the strong nuclear force. But the conceptual achievement stands: they proved that quantum degeneracy has a hard ceiling, that there is a maximum mass for cold matter, and that exceeding it leads to unstoppable collapse. The first actual neutron star, the Crab pulsar, would not be discovered until 1968 — nearly thirty years later.

TOV limit vs Chandrasekhar limit

Property	Chandrasekhar limit	TOV limit
Object	White dwarf	Neutron star
Supporting pressure	Electron degeneracy	Neutron degeneracy + strong force
Maximum mass	1.4 M☉ (sharp)	≈ 2.2–2.3 M☉ (uncertain)
Physics needed	Special relativity + QM	General relativity + nuclear EoS
How sharp is it	Precise (clean physics)	A range (unknown dense matter)
Density at limit	~ 10⁹ – 10¹⁰ kg/m³	~ 10¹⁸ kg/m³ (several × ρ₀)
What happens above it	Collapse to neutron star or Type Ia SN	Collapse to a black hole
Derived by / year	Chandrasekhar, 1930–31	Oppenheimer & Volkoff, 1939

The two limits are rungs on the same ladder. A degenerate object hits one ceiling, transforms into the next-denser state, and is then held up against a fiercer gravity by the next quantum-mechanical pressure — until, at the TOV limit, there is no next state. Quark stars have been proposed as one more conceivable rung, but if they exist their maximum mass is similar, and the practical end of the line remains a black hole.

Rotation, supramassive stars, and related limits

The plain TOV limit assumes a non-rotating star. Real neutron stars spin, and rapid rotation provides extra centrifugal support. This raises the maximum mass by roughly 18-20% for a star spinning at its mass-shedding (Keplerian) limit, giving a maximum maximum mass of perhaps 2.6-2.8 M☉.

• Stable neutron star. Mass below the non-rotating TOV limit. Lives forever (barring accretion).
• Supramassive neutron star (SMNS). Mass above the TOV limit but below the maximum rotating mass. Held up only by rotation; as magnetic braking spins it down, it eventually crosses the TOV threshold and collapses to a black hole — sometimes long after a merger, possibly powering some short gamma-ray bursts and their X-ray plateaus.
• Hypermassive neutron star (HMNS). Even heavier, supported by differential rotation (the core spinning faster than the envelope). Extremely short-lived — milliseconds to seconds — as differential rotation is erased by viscosity and magnetic fields. The likely fate of the GW170817 remnant.
• The "mass gap." Below about 5 M☉ but above the TOV limit lies a regime once thought to be empty of compact objects. The 2.6 M☉ secondary in the GW190814 merger sits provocatively in this gap — it is either the heaviest neutron star or the lightest black hole ever seen, and we cannot yet tell which.

Common misconceptions and subtleties

• "The TOV limit is exactly 3 solar masses." Three solar masses is the causality upper bound — the absolute most a neutron star could be under the stiffest physically allowed matter. The realistic, observationally favoured value is closer to 2.2-2.3 M☉. The figure "2 to 3 M☉" you see quoted brackets the uncertainty, not a precise prediction.
• "Neutron degeneracy pressure alone holds the star up." No — for a realistic neutron star, the repulsive strong nuclear force provides most of the support near the limit. Pure neutron degeneracy (the Oppenheimer-Volkoff ideal-gas model) caps out at only 0.7 M☉, well below every observed neutron star.
• "Exceeding the TOV limit is like exceeding Chandrasekhar — you get a supernova." Not quite. Exceeding Chandrasekhar in a white dwarf can trigger a thermonuclear runaway (Type Ia) or core collapse. Exceeding the TOV limit in a neutron star produces direct collapse to a black hole, generally without a luminous explosion — the matter simply disappears behind a horizon, often dubbed a "failed supernova" or quiet collapse.
• "The collapse takes a while." Once stability is lost, collapse proceeds on the free-fall timescale across ~10 km — under a millisecond. There is no slow squeeze; the transition from neutron star to black hole is essentially instantaneous on astronomical scales.
• "It applies to the original supernova too." A core-collapse supernova whose iron core exceeds the effective TOV mass during collapse can skip the neutron-star stage and go straight to a black hole. This direct-collapse channel is one reason some massive stars vanish without a bright supernova, as seen in surveys monitoring red supergiants that simply disappear.
• "PSR J0952-0607 at 2.35 M☉ breaks the limit." It does not — it pushes the lower bound on the TOV mass upward and tightens constraints on the equation of state. The limit must be at least as large as the heaviest star that exists; a heavier pulsar means a stiffer equation of state, not a violated limit.