Neutron Star Structure and Interior

A neutron star is what a massive star's iron core becomes when it loses the fight against gravity: the Chandrasekhar limit is breached, electrons are crushed into protons, and collapse halts only near the density inside an atomic nucleus. The result packs about $1.4\,M_\odot$ into a sphere roughly $12\,\mathrm{km}$ across — a mean density near $3.8\times10^{14}\,\mathrm{g\,cm^{-3}}$ , a surface gravity $\sim10^{11}$ times Earth's, and a compactness $GM/Rc^{2}\approx0.17$ so large that Newtonian gravity is simply wrong. Structure must be found by integrating the general-relativistic Tolman–Oppenheimer–Volkoff (TOV) equation, and that integration needs one input the terrestrial laboratory cannot supply: the equation of state $P(\rho)$ of cold matter above nuclear density. That single unknown is why the interior is the most interesting stratified object in physics. Reading inward, density climbs across fourteen decades and the composition changes at every landmark. A thin atmosphere and ocean give way to a solid outer crust — a Coulomb lattice of neutron-rich nuclei in a relativistic electron sea. At the neutron-drip point, $\rho_{\rm drip}\approx4\times10^{11}\,\mathrm{g\,cm^{-3}}$ , neutrons spill out of nuclei into a free superfluid; nuclei finally dissolve near nuclear saturation $\rho_0\approx2.8\times10^{14}\,\mathrm{g\,cm^{-3}}$ into a uniform fluid of neutrons, protons, electrons and muons, where neutrons pair into a $^{3}P_2$ superfluid and protons into a $^{1}S_0$ superconductor. The inner core, at a few $\rho_0$ , is the frontier: hyperons, meson condensates, or deconfined quark matter may live there, each rewriting the EOS in ways that measurements of the maximum mass, radius, and tidal deformability are only now beginning to decide. Neutron star structure is where nuclear physics, condensed matter, and relativistic gravity meet in one object.

Key takeaways

A neutron star is stratified by density: atmosphere and ocean, a solid crust of neutron-rich nuclei (free neutrons appear below the drip point), a fluid core of superfluid neutrons and superconducting protons, and an inner core of uncertain composition.
Density spans fourteen decades, from $\sim10^{6}\,\mathrm{g\,cm^{-3}}$ at the surface to $\gtrsim10^{15}\,\mathrm{g\,cm^{-3}}$ (several times $\rho_0\approx2.8\times10^{14}\,\mathrm{g\,cm^{-3}}$ ) at the centre.
Because $GM/Rc^{2}\approx0.17$ is order unity, structure obeys the relativistic TOV equation, not Newtonian hydrostatics; the supranuclear EOS is the central unknown fixing the mass–radius relation and maximum mass.
An ideal neutron gas caps at only $\approx0.7\,M_\odot$ ; the observed $\approx2\,M_\odot$ pulsars prove the repulsive nuclear force stiffens the EOS.
Superfluidity is a bulk property with observable consequences: quantized vortices ( $\kappa=h/2m_n$ ) pin to the crust and drive glitches, while pairing gaps regulate cooling.
Key idea: given a microphysical equation of state $P(\rho)$ , the Tolman–Oppenheimer–Volkoff equation integrates outward from a chosen central density to a surface where $P=0$ , uniquely fixing mass and radius; the layered composition of a neutron star is the coordinate-space image of that EOS across fourteen decades in density.

The core insight

Given a microphysical equation of state $P(\rho)$ , the Tolman–Oppenheimer–Volkoff equation integrates outward from a chosen central density to a surface where $P=0$ , uniquely determining mass and radius; the layered composition of a neutron star is the coordinate-space image of that EOS across fourteen decades in density.

Where this lives today

Neutron star structure has turned from theory into a measurement programme, and the target is always the equation of state at $(1\text{–}5)\,\rho_0$ . NASA's NICER X-ray telescope times the hot polar spots of pulsars, using relativistic light-bending of the pulse profile to solve simultaneously for mass and radius — PSR J0030+0451 near $1.4\,M_\odot$ and PSR J0740+6620 at $2.08\,M_\odot$ both land around $12\text{–}13\,\mathrm{km}$ , already excluding the softest and stiffest EOS. From the other side, the 2017 merger GW170817 measured the tidal deformability $\Lambda$ — how much each star's quadrupole is raised by its companion before merger — bounding the EOS from the stiff side and disfavouring very large radii. Radio timing supplies the high-mass anchor through the Shapiro delay, while pulsar glitches and the cooling curve of Cassiopeia A test the superfluid physics directly. The open frontier is the inner core: whether hyperons, a kaon condensate, or deconfined quark matter appears at a few $\rho_0$ is undecided, and the next generation of X-ray timing and louder merger signals is aimed squarely at it.

The physical picture

Imagine compressing ordinary matter until electrons are driven into nuclei, protons turn to neutrons, and the whole star becomes a single gravitationally bound nucleus — but one whose chemistry changes with depth, because pressure rises monotonically inward and each composition is favoured only over a specific pressure range. The radial cross-section is therefore a phase diagram of cold dense matter, traversed once from atmosphere to centre; the first figure lays it out, and each density boundary is a real number the EOS must reproduce.

Figure

Cross-section of a 1.4 M⊙ neutron star about 12 km in radius: density climbs inward across fourteen decades, from a thin atmosphere and ocean through the outer crust of neutron-rich nuclei, past the neutron-drip line at 4×10¹¹ g cm⁻³ into the inner crust, across nuclear saturation ρ₀ ≈ 2.8×10¹⁴ g cm⁻³ where nuclei dissolve, into the superfluid outer core and finally an inner core at 2 to 8 ρ₀ whose composition is unknown

Two different kinds of physics hold the star up: in the crust, the relativistic electron gas in a Coulomb lattice; in the core, neutron degeneracy pressure plus the stiff, short-range repulsion of the nuclear force, whose stiffness at $\rho\gg\rho_0$ is exactly what the maximum-mass measurement pins down. The outputs are not free — for a chosen EOS the TOV equation delivers one family of stars, one per central density, tracing a single curve in the mass–radius plane. The second figure shows what a representative EOS produces and where the observations, and the exclusions, fall.

Figure

The observable output of the TOV equation plus an equation of state is the mass–radius curve. A representative stiff EOS threads a typical 1.4 M⊙ star near 12.5 km, rises to a maximum around 2.2 M⊙ where dM/dρc = 0, and is capped from the upper left by general relativity as R approaches 2GM/c². An ideal non-interacting neutron gas would top out at only 0.7 M⊙, far below the measured 2.08 M⊙ of PSR J0740+6620

The core equations

The spine of the subject is the relativistic equation of hydrostatic equilibrium. For a static, spherical, perfect fluid, Einstein's equations reduce to the Tolman–Oppenheimer–Volkoff equation — the Newtonian balance $-Gm\rho/r^2$ decorated by three factors, each making gravity stronger.

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

where: $P$ — pressure at radius $r$ ; $\rho$ — total mass–energy density; $m=m(r)$ — enclosed gravitational mass; $G,\,c$ — gravitational constant and speed of light. Reading it: pressure itself gravitates ( $\rho\to\rho+P/c^2$ ), pressure adds to the active mass ( $m\to m+4\pi r^3P/c^2$ ), and curvature amplifies gravity through $(1-2Gm/rc^2)^{-1}$ . All three factors exceed unity, so relativity always makes a star harder to hold up. The TOV equation is closed by mass continuity, which converts the local density into the enclosed mass needed at each step of the outward integration.

\frac{\mathrm{d}m}{\mathrm{d}r}=4\pi r^{2}\rho \qquad (2)

where: $m$ — enclosed gravitational mass; $\rho$ — local mass–energy density. Integrating (1) and (2) from a central density $\rho_c$ until $P=0$ returns one $(M,R)$ point; sweeping $\rho_c$ traces the whole mass–radius curve. The system is closed by a third relation — the equation of state $P(\rho)$ . A neutron star is cold on the nuclear scale (interior $\sim10^{8}\,\mathrm{K}$ against a Fermi temperature $\sim10^{12}\,\mathrm{K}$ ), so thermal pressure is negligible and the matter is barotropic, its composition fixed at each density by chemical equilibrium.

P=P(\rho)\quad\text{(cold, barotropic, }T\to0\text{)} \qquad (3)

where: $P(\rho)$ — the cold-catalyzed EOS. This one function, unknown above $\rho_0$ , is what the whole observational programme is trying to measure; every other quantity is a functional of it. At each density the proton fraction is set by beta equilibrium under $n\rightleftharpoons p+e^-+\bar\nu_e$ , together with charge neutrality. These conditions remove the compositional freedom and make the EOS a single-valued curve.

\mu_n=\mu_p+\mu_e,\qquad n_p=n_e+n_\mu \qquad (4)

where: $\mu_i$ — chemical potential of species $i$ ; $n_i$ — number density of species $i$ ( $n,p,e,\mu$ ). Muons switch on once the electron Fermi energy exceeds $m_\mu c^2$ ; the proton fraction stays small, a few to ten percent, across the outer core. The leading term in the core pressure is neutron degeneracy. Filling momentum states up to the Fermi momentum $p_F$ gives a stiff $\Gamma=5/3$ law non-relativistically, but once $p_F\gtrsim m_n c$ the gas turns ultra-relativistic and softens to $\Gamma=4/3$ — the marginal index for stability and the root of the maximum-mass problem.

P_{\rm deg}=\frac{(3\pi^{2})^{2/3}}{5}\frac{\hbar^{2}}{m_n}\,n^{5/3}\ (\text{non-rel})\ \xrightarrow{\ p_F\gtrsim m_nc\ }\ \frac{(3\pi^{2})^{1/3}}{4}\hbar c\,n^{4/3}\ (\text{ultra-rel}) \qquad (5)

where: $n$ — neutron number density; $m_n$ — neutron mass; $\hbar$ — reduced Planck constant. The $5/3$ law would hold up a Newtonian star; when it softens to $4/3$ , the GR terms in (1) tip the star over, capping the mass of a non-interacting gas. Two composition landmarks anchor the layered picture and recur throughout: neutron drip, where free neutrons first appear, and nuclear saturation, where nuclei dissolve into uniform matter.

\rho_{\rm drip}\approx4\times10^{11}\,\mathrm{g\,cm^{-3}},\qquad \rho_0\approx2.8\times10^{14}\,\mathrm{g\,cm^{-3}}\;(n_0\approx0.16\,\mathrm{fm^{-3}}) \qquad (6)

where: $\rho_{\rm drip}$ — density at which the neutron chemical potential reaches the continuum; $\rho_0$ — nuclear saturation, the density inside a heavy nucleus. The crust lives below $\rho_0$ , the core above, and the inner core reaches $(2\text{–}8)\,\rho_0$ . Finally, the core neutrons are a superfluid, and a superfluid cannot rotate rigidly — it threads the star with quantized vortices, each carrying a fixed circulation set by the Cooper-pair mass $2m_n$ .

\kappa=\frac{h}{2m_n}=\frac{\pi\hbar}{m_n}\approx1.98\times10^{-3}\,\mathrm{cm^{2}\,s^{-1}},\qquad n_v=\frac{2\Omega}{\kappa} \qquad (7)

where: $\kappa$ — circulation quantum per vortex line; $\Omega$ — stellar angular velocity; $n_v$ — areal vortex density. Matching the mean vorticity to solid-body rotation fixes $n_v$ ; pinning of these vortices to the crustal lattice is the mechanical model behind pulsar glitches.

Derivation AOrder-of-magnitude structure and why radius falls as mass rises

The goal is to recover the whole scale of a neutron star — its kilometre radius and supranuclear central density — from nothing but a balance between gravity and non-relativistic neutron degeneracy pressure, and in the process to expose the counter-intuitive signature of every degenerate star: heavier means smaller. Treat the star as a uniform ball of $N=M/m_n$ neutrons and match the two competing pressures at the virial scale.

Step 1. Estimate the number density as $n\sim N/R^3=M/(m_n R^3)$ , then insert it into the Fermi-gas pressure (5): $P_{\rm deg}\sim(\hbar^2/m_n)[M/(m_n R^3)]^{5/3}$ . Why: degeneracy pressure is purely kinematic, depending on nothing but the density, which makes it the clean counter-force to gravity.

Step 2. Estimate the gravitational pressure as $P_{\rm grav}\sim GM^2/R^4$ , the virial pressure of a self-gravitating ball (energy $\sim GM^2/R$ over volume $\sim R^3$ ). Why: this is the pressure gravity can supply, and equilibrium demands it balance the degeneracy term.

Step 3. Set $P_{\rm deg}\sim P_{\rm grav}$ ; the $R^{-5}$ versus $R^{-4}$ powers give $R\sim\hbar^2/(G\,m_n^{8/3}\,M^{1/3})\propto M^{-1/3}$ . Why: the inverse-cube-root dependence is the hallmark of a degenerate star — adding mass compresses it — shared with the white dwarf, differing only in which fermion sets $m_n^{8/3}$ .

Step 4. Put in $M=1.4\,M_\odot$ : the bare expression gives $R\approx3\,\mathrm{km}$ , and the order-unity structural constant lifts it to $\sim10\text{–}12\,\mathrm{km}$ , with mean density a few $\rho_0$ . Why: recovering a kilometre radius at nuclear density from a one-line balance confirms the physical picture before any relativistic refinement.

Go deeper · Derivation B · 1 minThe maximum mass of an ideal neutron gas

Derivation BThe maximum mass of an ideal neutron gas

The goal is to show why a neutron star built from a non-interacting neutron gas cannot exceed about $0.7\,M_\odot$ , why this Oppenheimer–Volkoff bound sits absurdly below the $\approx2\,M_\odot$ stars we observe, and what that discrepancy says about the strong force. The argument combines the relativistic softening of the EOS with the destabilizing pressure terms of the TOV equation to produce a genuine maximum of mass versus central density.

Step 1. As central density rises the Fermi momentum grows until $p_F\gtrsim m_n c$ and the EOS softens from $P\propto\rho^{5/3}$ to the ultra-relativistic $P\propto\rho^{4/3}$ of equation (5). Why: at $\Gamma=4/3$ pressure and self-gravity scale the same way with compression, so the gas loses the stiffness margin it relies on to resist collapse.

Step 2. Feed this into the TOV terms $\rho+P/c^2$ and $m+4\pi r^3P/c^2$ : because pressure now gravitates strongly, raising $\rho_c$ boosts effective gravity faster than pressure support. Why: a $\Gamma=4/3$ gas is only marginally stable even in Newtonian gravity, and the relativistic pressure terms are strictly destabilizing, tipping the marginal case over.

Step 3. Integrate TOV with the free-neutron EOS and follow $M(\rho_c)$ : it rises, flattens, and turns over at $\rho_c\approx5\times10^{15}\,\mathrm{g\,cm^{-3}}$ where $\mathrm{d}M/\mathrm{d}\rho_c=0$ , giving $M_{\max}^{\rm OV}\approx0.71\,M_\odot$ at $R\approx9.6\,\mathrm{km}$ . Why: the turning point is a true maximum — stars beyond it are unstable and collapse — and it lies far below the measured $2.08\,M_\odot$ , falsifying the ideal-gas model.

Step 4. Restore the omitted physics: the nucleon–nucleon force is strongly repulsive at short range, stiffening the EOS above $\rho_0$ and pushing $M_{\max}\gtrsim2\,M_\odot$ while enlarging $R$ . Why: the gap between $0.7$ and $2\,M_\odot$ is direct evidence that dense matter is not an ideal gas — the maximum mass measures the stiffness of the strong force.

Go deeper · Derivation C · 1 minThe onset of neutron drip

Derivation CThe onset of neutron drip

This derivation tracks how the outer crust's composition evolves under compression until, at a sharply defined density near $4\times10^{11}\,\mathrm{g\,cm^{-3}}$ , neutrons begin to leak from nuclei and form a free gas — the transition that creates the two-phase inner crust and the superfluid reservoir behind pulsar glitches. The controlling idea is that rising electron pressure forces electron capture, driving nuclei ever more neutron-rich until the last neutron level is unbound.

Step 1. As density climbs, the degenerate electrons' large Fermi energy makes electron capture $p+e^-\to n+\nu_e$ favourable, marching equilibrium nuclei from $^{56}$ Fe toward neutron-rich species such as $^{118}$ Kr. Why: minimizing total energy at fixed baryon number trades energetic electrons for bound neutrons, following a pressure-driven valley of stability, not the terrestrial one.

Step 2. Each added neutron is less bound than the last, so the neutron chemical potential $\mu_n$ rises toward the free-neutron rest energy. Why: neutron-rich nuclei sit near the edge of the binding surface, so the least-bound level climbs toward zero binding as more neutrons are packed in.

Step 3. Drip occurs when $\mu_n=m_n c^2$ : the least-bound neutrons reach the continuum and spill out. Crust models (Baym–Pethick–Sutherland) place this at $\rho_{\rm drip}\approx4.3\times10^{11}\,\mathrm{g\,cm^{-3}}$ . Why: once $\mu_n$ reaches the rest energy there is no barrier to a neutron leaving the nucleus, so free neutrons appear in equilibrium with the lattice.

Step 4. Above drip the matter is two-phase — a lattice bathed in a $^{1}S_0$ neutron superfluid — sculpted into nuclear pasta near the crust base by Coulomb–surface competition. Why: this dripped superfluid is the angular-momentum reservoir whose vortices pin to nuclear clusters, making neutron drip the microscopic origin of glitches.

Limiting cases

Case	Result	Meaning
Neutron drip	$\rho_{\rm drip}\approx4\times10^{11}\,\mathrm{g\,cm^{-3}}$	Neutron chemical potential reaches the continuum; free neutrons appear, marking the outer/inner crust boundary.
Nuclear saturation	$\rho_0\approx2.8\times10^{14}\,\mathrm{g\,cm^{-3}}\;(n_0\approx0.16\,\mathrm{fm^{-3}})$	Density inside heavy nuclei; near this value nuclei dissolve and the uniform core begins (around $0.5\rho_0$ ).
Core densities	$\rho_c\sim(2\text{–}8)\,\rho_0$	Central density of typical to massive stars; the regime where exotic phases may onset.
Neutrons at $\rho_0$	$p_Fc\approx331\,\mathrm{MeV},\ p_F/m_nc\approx0.35$	Mildly relativistic even at saturation, rising to $\approx0.60$ at $5\rho_0$ — relativity is unavoidable in the core.
Ideal neutron gas limit	$M_{\max}^{\rm OV}\approx0.7\,M_\odot$	Maximum mass of a non-interacting degenerate neutron gas — far below observed masses, proving interactions matter.
Newtonian limit	$GM/Rc^2\to0$	Recovers $\mathrm{d}P/\mathrm{d}r=-Gm\rho/r^2$ ; fine for white dwarfs, wrong here where $GM/Rc^2\approx0.17$ .

Common trapA neutron star is just a giant nucleus held up by neutron degeneracy pressure

This picture is doubly wrong. First, an ideal degenerate neutron gas gives an Oppenheimer–Volkoff maximum mass of only $\approx0.7\,M_\odot$ — you cannot even build a $1.4\,M_\odot$ star that way, let alone the observed $2.08\,M_\odot$ PSR J0740+6620. The strong short-range repulsion of the nuclear force does most of the real supporting work above $\rho_0$ ; degeneracy pressure is the seed, not the whole story, and the maximum mass is essentially a measurement of how stiff that force is. Second, never reach for the Newtonian hydrostatic equation. With compactness $GM/Rc^2\approx0.17$ , the three general-relativistic factors in the TOV equation are order-unity corrections — pressure gravitates, pressure adds to the source mass, and the metric factor $(1-2Gm/rc^2)^{-1}\approx1.5$ — and omitting them overestimates the maximum mass badly. A neutron star is a relativistic object made of strongly interacting matter, not a cold Fermi gas in Newtonian gravity.

Worked examples

1 · Why general relativity is not optional (fully worked)

Take a canonical neutron star, $M=1.4\,M_\odot$ and $R=12\,\mathrm{km}$ . Compute its compactness, the TOV metric factor, the surface gravitational redshift, the surface gravity relative to Earth, and the mean density in units of $\rho_0$ . Decide whether Newtonian hydrostatics could ever be adequate.

Step 1. Compactness. With $M=2.78\times10^{30}\,\mathrm{kg}$ and $R=1.2\times10^{4}\,\mathrm{m}$ , $GM/Rc^2=(6.674\times10^{-11})(2.78\times10^{30})/[(1.2\times10^{4})(2.998\times10^{8})^2]=\mathbf{0.172}$ — a fifth of the way to a black hole ( $GM/Rc^2=0.5$ ), so relativistic corrections are order-unity.

Step 2. Metric factor. The TOV curvature term is $(1-2GM/Rc^2)^{-1}=(1-0.345)^{-1}=\mathbf{1.53}$ at the surface — gravity amplified $50\%$ by spacetime curvature alone, before the pressure terms.

Step 3. Surface redshift. A photon climbing out is redshifted by $1+z=(1-2GM/Rc^2)^{-1/2}=1.24$ , so $\mathbf{z\approx0.24}$ — directly measurable in X-ray spectra and a clean handle on $M/R$ .

Step 4. Surface gravity. $g=GM/R^2=1.29\times10^{12}\,\mathrm{m\,s^{-2}}$ , which is $\mathbf{1.3\times10^{11}}$ times Earth's $9.81\,\mathrm{m\,s^{-2}}$ .

Step 5. Mean density. $\bar\rho=3M/(4\pi R^3)=3.85\times10^{14}\,\mathrm{g\,cm^{-3}}\approx\mathbf{1.4\,\rho_0}$ — nuclear density averaged over the whole star. A compactness of 0.17, a metric factor of 1.5, and a redshift of 0.24 all say the same thing: Newtonian hydrostatics is disqualified and the TOV equation is mandatory.

2 · Superfluid vortices in the Vela pulsar (fully worked)

Vela (PSR B0833-45) spins at $\nu=11.19\,\mathrm{Hz}$ . Compute the circulation quantum, the areal vortex density in the core, and the total number of vortices threading a $10\,\mathrm{km}$ core.

Step 1. Circulation quantum. $\kappa=h/2m_n=(6.626\times10^{-34})/(2\times1.675\times10^{-27})=1.98\times10^{-7}\,\mathrm{m^2\,s^{-1}}=\mathbf{1.98\times10^{-3}\,cm^2\,s^{-1}}$ , set entirely by the Cooper-pair mass $2m_n$ .

Step 2. Vortex density. With $\Omega=2\pi\nu=70.3\,\mathrm{rad\,s^{-1}}$ , equation (7) gives $n_v=2\Omega/\kappa=140.6/(1.98\times10^{-7})=7.1\times10^{8}\,\mathrm{m^{-2}}=\mathbf{7.1\times10^{4}\,cm^{-2}}$ .

Step 3. Total vortices. Over a core cross-section $\pi R_{\rm core}^2=\pi(10^6\,\mathrm{cm})^2=3.1\times10^{12}\,\mathrm{cm^2}$ , the count is $N\approx n_v\pi R_{\rm core}^2=(7.1\times10^4)(3.1\times10^{12})\approx\mathbf{2\times10^{17}\ vortices}$ . A macroscopic pulsar is threaded by hundreds of quadrillions of quantized vortex lines; pinning and collective unpinning of this array is the accepted engine of the Vela glitches.

Physics of the portrait

The structure of a neutron star was computed before anyone had seen one. In the 1930s it was a purely theoretical object — predicted, given a governing equation, and handed a paradox (the maximum mass) within a decade, then left waiting three more for the observation that proved it real.

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

J. Robert Oppenheimer

1904–1967 · relativistic structure of dense stars

Oppenheimer is remembered publicly for the Manhattan Project, but to compact-star physics he is the author, with George Volkoff, of the 1939 paper that founded the subject. The stage was set quickly: Chadwick found the neutron in 1932, Landau speculated about stars at nuclear density, and in 1934 Baade and Zwicky proposed that a supernova turns an ordinary star into a neutron star. Richard Tolman supplied the general-relativistic machinery, and Oppenheimer and Volkoff integrated what is now the Tolman–Oppenheimer–Volkoff equation for an ideal degenerate neutron gas. The result was startling: rather than a stiffer white dwarf, the ideal gas gave a maximum mass of only about $0.7\,M_\odot$ — below the Chandrasekhar mass, so such stars could not even form from the cores they were meant to describe. The resolution, unavailable to them, is that dense matter is not an ideal gas: the repulsive core of the nuclear force stiffens the EOS and lifts the maximum past $2\,M_\odot$ , and that the ideal-gas bound is wrong by a factor of three is exactly what makes the maximum mass such a sharp probe of the strong interaction. The object Oppenheimer described sat unconfirmed until Jocelyn Bell's pulsar in 1967 revealed it spinning in the sky.

History

1934

— Baade and Zwicky propose neutron stars as the remnants of supernovae

1939

— Oppenheimer and Volkoff integrate the TOV equation and find a maximum mass near 0.7 M⊙

1967

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

1971

— Baym, Pethick and Sutherland compute the crust equation of state and the neutron-drip density

2017

— GW170817 measures the tidal deformability, constraining the dense-matter equation of state

Going further

Texts. Carroll and Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 16 (The Degenerate Remnants of Stars) · Shapiro and Teukolsky, Black Holes, White Dwarfs, and Neutron Stars, Ch. 8–11 · Glendenning, Compact Stars: Nuclear Physics, Particle Physics, and General Relativity, 2nd ed. · Haensel, Potekhin and Yakovlev, Neutron Stars 1: Equation of State and Structure · Lattimer and Prakash, "Neutron Star Structure and the Equation of State," ApJ 550, 426 (2001).

Provenance. A verifier should reproduce: the TOV equation and its three GR factors from the interior Schwarzschild solution (Shapiro and Teukolsky; Glendenning); the Newtonian scaling $R\propto M^{-1/3}$ giving a kilometre-scale radius at $1.4\,M_\odot$ (Carroll and Ostlie Ch. 16); the Oppenheimer–Volkoff maximum mass $\approx0.7\,M_\odot$ for an ideal neutron gas (Oppenheimer and Volkoff 1939); the landmarks $\rho_{\rm drip}\approx4\times10^{11}$ and $\rho_0\approx2.8\times10^{14}\,\mathrm{g\,cm^{-3}}$ (Baym, Pethick and Sutherland 1971); the circulation quantum $\kappa=h/2m_n\approx1.98\times10^{-3}\,\mathrm{cm^2\,s^{-1}}$ with Vela vortex count $\sim10^{17}$ ; and the worked compactness $GM/Rc^2\approx0.17$ and redshift $z\approx0.24$ . All numbers are standard tabulated or elementary-calculation quantities; none were invented.

Prerequisites. Neutron Stars · Tolman–Oppenheimer–Volkoff Equation · Equation of State of Dense Matter · Electron Degeneracy Pressure · Chandrasekhar Mass Limit. Leads to. Pulsar Glitches and Neutron Star Superfluidity · Magnetars · Millisecond and Recycled Pulsars · Pulsars.

Ask the tutor: Ask me to integrate the TOV equation for a toy polytropic EOS and watch the maximum mass appear as the mass–radius curve turns over.

What actually sets the 1.4 M⊙ and 12 km scale?

The mass scale comes from stellar evolution: neutron stars are born from iron cores at the Chandrasekhar mass, $\approx1.4\,M_\odot$ , so birth masses cluster there. The radius comes from the order-of-magnitude balance — gravity against neutron degeneracy pressure near nuclear density — which fixes $R$ to about $10\text{–}13\,\mathrm{km}$ almost independently of mass, because the EOS is stiff enough to hold the radius nearly constant. The precise radius at a given mass is set by the EOS at $(1\text{–}2)\,\rho_0$ , which is why measuring $R$ to a kilometre with NICER so tightly constrains the pressure of dense matter.

What is the difference between the crust and the core?

The dividing line is nuclear saturation, $\rho_0\approx2.8\times10^{14}\,\mathrm{g\,cm^{-3}}$ . Below it, in the crust, matter is inhomogeneous: distinct nuclei on a Coulomb lattice, a relativistic electron gas throughout, and (below neutron drip) a free neutron superfluid between the nuclei. Above it, in the core, the nuclei have dissolved into a uniform fluid of neutrons, protons, electrons and muons. The crust is solid, roughly a kilometre thick, and stores the pinned vorticity that glitches release; the fluid core holds most of the mass and moment of inertia, and its EOS sets the maximum mass.

How do we actually measure the equation of state?

Through the mass–radius relation, because for a given EOS the TOV equation produces exactly one $M(R)$ curve — so measuring points on that curve inverts to the EOS. Three methods contribute: radio timing weighs binary pulsars via the Shapiro delay, delivering precise masses and the crucial $2\,M_\odot$ lower bound; NICER's X-ray timing models hot-spot pulse profiles to extract mass and radius together; and gravitational-wave detectors measure the tidal deformability in a merger, which depends steeply on radius. Each carves the allowed EOS band from a different side, closing in on the pressure–density curve at a few times nuclear density.

What are pulsar glitches, and what do they tell us about the interior?

A glitch is a sudden, tiny spin-up of an otherwise steadily slowing pulsar, seen most dramatically in Vela. In the standard model the crustal neutron superfluid rotates via quantized vortices pinned to the lattice; as the crust spins down, the pinned superfluid lags and stores angular momentum, until the vortices unpin catastrophically and dump that momentum onto the crust. The glitch size and slow post-glitch relaxation measure the fraction of the star's moment of inertia held in the superfluid and the strength of vortex pinning — a direct dynamical probe of interior superfluidity that no static structure model could reveal.

Could the inner core be quark matter, and how would we ever know?

It is a live possibility. Above a few $\rho_0$ the nucleons are packed so densely that their quarks may deconfine, and hyperons or a meson condensate could also appear. We cannot make this matter in the laboratory — heavy-ion collisions reach the required density only at high temperature, not the cold, dense corner a neutron star occupies — so the evidence must be astrophysical: each candidate phase softens or stiffens the EOS differently, changing the maximum mass, the radius at fixed mass, and the tidal deformability. A sufficiently precise mass–radius measurement, or a merger showing a phase-transition signature, could reveal a kink in the EOS that only deconfinement would produce. For now the question is open.

Why does the maximum mass matter so much?

Because it is the sharpest single number the EOS predicts, and it is cleanly measurable. A stiffer EOS resists compression more and supports a larger maximum mass. So every confirmed high-mass pulsar sets a hard lower bound on the maximum mass, and therefore on the stiffness of dense matter, ruling out any EOS too soft to reach it. The $2.08\,M_\odot$ of PSR J0740+6620 already eliminates a swathe of soft models, including any that let hyperons excessively soften the core. The maximum mass also fixes the boundary between the heaviest neutron star and the lightest black hole — the "mass gap" probed by gravitational-wave catalogues.