Tolman–Oppenheimer–Volkoff Limit

The Tolman–Oppenheimer–Volkoff limit (or TOV limit) is the fundamental upper bound on the gravitational mass of a cold, non-rotating neutron star in hydrostatic equilibrium, beyond which relativistic neutron degeneracy pressure and repulsive strong-force nucleon–nucleon interactions can no longer resist general-relativistic gravitational collapse into a stellar-mass black hole. Analogous to the Chandrasekhar limit (~1.44 M_☉) for electron-degenerate white dwarfs, the TOV limit was derived in 1939 by American physicists J. Robert Oppenheimer and George Volkoff, building on the relativistic static fluid sphere solutions formulated by Richard C. Tolman earlier that same year. While Oppenheimer and Volkoff's idealized calculation for a free, non-interacting Fermi gas of neutrons yielded a limit of only ~0.71 M_☉, modern nuclear equations of state (EoS) that incorporate repulsive short-range strong nuclear interactions and multi-messenger observations from pulsar timing and binary neutron star gravitational-wave mergers constrain the true non-rotating TOV limit to M_TOV ≈ 2.17–2.35 M_☉.

In Newtonian astrophysics, a spherically symmetric star of mass density ρ(r) and enclosed mass m(r) maintains hydrostatic equilibrium when the outward radial pressure gradient balances inward gravitational attraction:

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

Inside a neutron star, however, the gravitational potential parameter GM / (Rc²) approaches 0.2 to 0.3, and the internal pressure P becomes comparable to the rest-mass energy density ρc². Solving Einstein's field equations (G_μν = (8πG/c⁴)T_μν) for a static, spherically symmetric perfect fluid metric yields the Tolman–Oppenheimer–Volkoff equation:

dP / dr = − [G m(r) ρ(r) / r²] · [1 + P(r) / (ρ(r)c²)] · [1 + 4πr³P(r) / (m(r)c²)] · [1 − 2*Gm*(r) / (c²r)]⁻¹

coupled with the relativistic mass-continuity relation:

dm / dr = 4πr² ρ(r)

where ρ(r) = ε(r)/c² is the total mass-energy density (including rest mass plus internal kinetic and interaction energy). Compared to the Newtonian hydrostatic equation, every one of the three general-relativistic correction factors in brackets is greater than unity and acts to steepen the inward pressure gradient:

1. Inertial Mass of Pressure (1 + P / ρc²): In general relativity, enthalpy density (ρc² + P) rather than rest density alone responds to gravitational acceleration.
2. Active Gravitational Mass of Pressure (1 + 4πr³P / mc²): Because the stress-energy tensor T_μν includes spatial pressure components on its diagonal, internal pressure acts as an additional source of gravity ("pressure gravitates").
3. Schwarzschild Metric Curvature (1 − 2*Gm* / c²r)⁻¹: Spatial radial curvature intensifies the effective gravitational force as the radius r approaches the Schwarzschild radius R_S = 2*GM*/c².

Because increasing pressure eventually strengthens the inward gravitational pull faster than it can stiffen the core, there exists a maximum mass M_TOV at which dM/dρ_c = 0 (where ρ_c is the central density); any configuration on the branch dM/dρ_c < 0 is dynamically unstable to radial collapse.

To integrate the TOV differential equations from a chosen central density ρ_c outward to the stellar surface (P(R) = 0), one must specify a barotropic equation of state P(ρ) spanning densities from sub-nuclear crust matter up to 5–8 times the nuclear saturation density (ρ₀ ≈ 2.7 × 10¹⁴ g cm⁻³, or n₀ ≈ 0.16 nucleons fm⁻³):

Since 2010, three independent observational techniques have bracketed the empirical value of M_TOV from both below and above:

1. Lower Bound from Shapiro Delay in Heavy Millisecond Pulsars

When a radio millisecond pulsar orbits an inclined white dwarf companion nearly edge-on, radio pulses passing through the companion's gravitational potential well experience a general-relativistic Shapiro time delay that directly yields the pulsar's gravitational mass:

• PSR J1614−2230 (Demorest et al., 2010): Measured at M = 1.908 ± 0.016 M_☉, ruling out soft hyperon and kaon-condensate equations of state.
• PSR J0348+0432 (Antoniadis et al., 2013): Measured via white dwarf spectroscopy and orbital decay at M = 2.01 ± 0.04 M_☉.
• PSR J0740+6620 (Cromartie et al., 2020; Fonseca et al., 2021): Measured via Green Bank and CHIME Shapiro delay at M = 2.08 ± 0.07 M_☉ (with NICER X-ray pulse-profile modeling constraining its equatorial radius to R = 12.39^+1.30_−0.98 km), establishing a firm empirical lower bound of M_TOV > 2.01 M_☉.
• Black-Widow Pulsar PSR J0952−0607 (Romani et al., 2022): Optical radial-velocity spectroscopy of its irradiated companion yields M = 2.35 ± 0.17 M_☉ for a rapidly rotating (707 Hz) pulsar.

2. Upper Bound from Binary Neutron Star Merger GW170817

On August 17, 2017, the LIGO and Virgo gravitational-wave detectors observed GW170817, the inspiral and coalescence of two neutron stars with total binary mass M_tot ≈ 2.74 M_☉, followed 1.7 seconds later by short gamma-ray burst GRB 170817A and optical/infrared kilonova AT 2017gfo in galaxy NGC 4993. The electromagnetic kilonova and absence of an ultra-relativistic magnetar spindown wind indicated that the merger remnant temporarily formed a differentially rotating hypermassive neutron star (HMNS) that survived for ~100–1,000 milliseconds before collapsing into a black hole once neutrino cooling and magnetic braking removed its differential rotation. Combining the total remnant mass with rotational stability ratios constrains the non-rotating limit to M_TOV ≲ 2.17–2.30 M_☉.

Rapid rotation introduces centrifugal forces that partially counteract gravity, allowing a spinning neutron star to exceed the static limit M_TOV:

• Supramassive Neutron Stars (SMNS): A uniformly (rigidly) rotating neutron star spinning near its Keplerian mass-shedding frequency (f_K ≈ 1,000–1,500 Hz) can support up to M_max,rot ≈ 1.20 M_TOV (~2.5–2.8 M_☉). As magnetic dipole radiation gradually slows the star's spin over seconds to thousands of years, its centrifugal support decays until it crosses the stability boundary and collapses into a black hole (a proposed engine for Fast Radio Bursts and X-ray afterglow plateaus).
• Hypermassive Neutron Stars (HMNS): A differentially rotating merger remnant (in which the core spins faster than the outer envelope) can transiently support up to ~1.35–1.50 M_TOV (~3.0–3.4 M_☉), but survives for only tens to hundreds of milliseconds before magnetorotational viscosity redistributes angular momentum and triggers gravitational collapse.

1. Oppenheimer, J. R., & Volkoff, G. M. (1939). On massive neutron cores. Physical Review, 55(4), 374–381.
2. Tolman, R. C. (1939). Static solutions of Einstein's field equations for spheres of fluid. Physical Review, 55(4), 364–373.
3. Margalit, B., & Metzger, B. D. (2017). Constraining the maximum mass of neutron stars from multi-messenger observations of GW170817. The Astrophysical Journal Letters, 850(2), L19.
4. Fonseca, E., Cromartie, H. T., Pennucci, T. T., et al. (2021). Refined mass and geometric measurements of the high-mass PSR J0740+6620. The Astrophysical Journal Letters, 915(1), L12.