Neutron Star

A neutron star is the ultra-compact degenerate stellar remnant formed by the gravitational core collapse of a massive star—with a zero-age main-sequence mass between 8 M_☉ and 25–30 M_☉—during a Type II, Type Ib, or Type Ic core-collapse supernova. Packing a mass between 1.17 M_☉ and ~2.16–2.30 M_☉ into a sphere of coordinate radius R ≈ 10.5–12.5 km (roughly the diameter of a terrestrial city), a neutron star is the densest directly observable object in the cosmos prior to event-horizon formation at the Schwarzschild radius. Its mean mass density, ρ̄ ≈ 3.7–5.9 × 10¹⁷ kg m⁻³ (3.7–5.9 × 10¹⁴ g cm⁻³), exceeds the saturation density of an atomic nucleus (ρ₀ = 2.7 × 10¹⁷ kg m⁻³, or 0.16 nucleons fm⁻³) by up to two- to threefold, supported against gravitational collapse by quantum-mechanical neutron degeneracy pressure and short-range repulsive strong nuclear (nucleon–nucleon) forces.

The existence of neutron stars was predicted theoretically in December 1933 (published January 1934) by Walter Baade and Fritz Zwicky—less than two years after James Chadwick's 1932 discovery of the neutron—who proposed that supernovae represent the transition of an ordinary star into a compact star consisting mainly of neutrons, liberating gravitational binding energy E_B ≈ 3*GM*²/(5*R*) ≈ 3 × 10⁴⁶ J (~10% of its rest-mass energy Mc², carried away 99% by thermal neutrinos). In 1939, J. Robert Oppenheimer and George M. Volkoff (using Richard C. Tolman's relativistic metric solutions) derived the general relativistic equation of hydrostatic equilibrium for a cold degenerate Fermi gas of neutrons.

Neutron stars remained purely theoretical until November 1967, when the graduate student Jocelyn Bell Burnell and her advisor Antony Hewish at the Mullard Radio Astronomy Observatory in Cambridge discovered regular 81.5 MHz radio pulses with a period of P = 1.3373011 seconds from PSR B1919+21 (1974 Nobel Prize in Physics). In 1968, Thomas Gold and Franco Pacini independently identified pulsars as rapidly rotating, highly magnetized neutron stars—confirmed months later by the discovery of the 33.1-millisecond pulsar PSR B0531+21 at the center of the Crab Nebula (the remnant of SN 1054).

---

Iron Core Collapse Beyond the Chandrasekhar Limit

When a massive star (M_ZAMS ≳ 8–10 M_☉) exhausts silicon burning in its core, it accumulates an electron-degenerate core of iron-group nuclei (⁵⁶Fe, ⁵⁶Ni) surrounded by concentric shells of Si, O, Ne, C, He, and H. Because ⁵⁶Fe and ⁵⁶Ni have among the highest nuclear binding energies per nucleon (~8.8 MeV nucleon⁻¹), endothermic fusion can no longer generate thermal pressure. Once the iron core exceeds its effective Chandrasekhar limit (M_Ch,Fe ≈ 1.26–1.40 M_☉), two catastrophic instabilities trigger gravitational collapse at ~0.23*c* within ~100 milliseconds:

1. Photodisintegration: At T ≳ 5 × 10⁹ K, thermal gamma-ray photons shatter iron nuclei endothermically into alpha particles and free nucleons (⁵⁶Fe + γ → 13 ⁴He + 4 n − 124.4 MeV), robbing the core of thermal energy.
2. Electron Capture (Neutronization): As core density exceeds 10¹² kg m⁻³, the electron Fermi energy (E_F,e > 10 MeV) surpasses the mass-energy threshold required to force protons inside nuclei to capture relativistic electrons via inverse beta decay:

p + e⁻ → n + ν_e

Removing degenerate electrons collapses the core until the inner core reaches nuclear saturation density (ρ₀ = 2.7 × 10¹⁷ kg m⁻³), where the strong nuclear force stiffens the equation of state, rebounding the infalling outer envelope in a neutrino-driven core-collapse supernova shock and leaving behind a proto-neutron star (T ~ 10¹¹ K) that cools rapidly via Urca neutrino emission (n → p + e⁻ + ν̄_e; p + e⁻ → n + ν_e).

The Tolman–Oppenheimer–Volkoff (TOV) Equation and Maximum Mass

Because the gravitational compactness parameter of a neutron star is large (Ξ = GM/(Rc²) ≈ 0.15–0.30, so its coordinate radius R ≈ 11.5 km is only ~2.8 times its Schwarzschild radius r_s = 2*GM*/c² ≈ 4.13 km for a 1.4 M_☉ star), Newtonian gravity fails and hydrostatic structure is governed by the general relativistic Tolman–Oppenheimer–Volkoff (TOV) equation:

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

coupled with the enclosed gravitational mass equation dm(r)/dr = 4*πr*²ρ(r). Note that all three general-relativistic bracketed factors on the right-hand side exceed unity, meaning that pressure P(r) itself gravitates and spacetime curvature steepens the inward gravitational pull relative to Newtonian hydrostatic equilibrium. At the neutron star surface (r = R ≈ 11.5 km, M = 1.4 M_☉), the relativistic surface gravity and escape velocity reach:

g_surf = ( GM / R² ) ( 1 − 2*GM* / Rc² )^−1/2 ≈ 1.9 × 10¹² m s⁻² (~2 × 10¹¹ g_⊕)

v_esc = c √( 2*GM* / Rc² ) ≈ 0.58 c (~174,000 km s⁻¹)

When Oppenheimer and Volkoff (1939) solved the TOV equation assuming an ideal, non-interacting Fermi gas of neutrons, they obtained a maximum mass of only ~0.71 M_☉. In reality, at short inter-nucleon distances (r < 0.7 fm), repulsive vector-meson (ω, ρ) strong nuclear interactions dominate over degeneracy pressure alone, stiffening the nuclear equation of state P(ρ) and raising the Tolman–Oppenheimer–Volkoff (TOV) maximum mass limit to:

M_TOV ≈ 2.15–2.30 M_☉

Observationally, Shapiro time-delay and NASA NICER X-ray measurements have confirmed massive pulsars right up to this limit, including PSR J0740+6620 (M = 2.08 ± 0.07 M_☉, R = 12.39^+1.30_−0.98 km) and PSR J0348+0432 (M = 2.01 ± 0.04 M_☉), while multi-messenger constraints from the binary neutron star merger GW170817 bound the non-rotating ceiling to M_TOV ≲ 2.17–2.30 M_☉. Any compact remnant exceeding M_TOV collapses past its Schwarzschild radius into a stellar-mass black hole.

---

From surface to center, a neutron star is stratified into four distinct radial zones across 14 orders of magnitude in density:

1. Atmosphere and Ocean (Thickness ~0.1–100 m; ρ < 10⁹ kg m⁻³): A millimeter-to-centimeter scale atmosphere of ionized hydrogen, helium, or carbon plasma overlying a liquid coulomb ocean of fully ionized ⁵⁶Fe nuclei.
2. Outer Solid Crust (Thickness ~0.3–0.5 km; 10⁹ ≤ ρ ≤ ρ_drip = 4.3 × 10¹⁴ kg m⁻³): Ions freeze into a solid body-centered cubic (bcc) Coulomb crystal lattice immersed in a relativistic degenerate Fermi sea of electrons. As density and electron Fermi energy E_F,e increase with depth, inverse beta decay shifts equilibrium from symmetric ⁵⁶Fe (Z = 26, N = 30) through increasingly neutron-rich magic-number nuclei (⁶²Ni, ⁸⁰Zr, ¹¹⁸Kr) that are prevented from beta-decaying because all electron momentum states below p_F,e are Pauli-blocked.
3. Inner Crust and "Nuclear Pasta" (Thickness ~0.5–1.0 km; ρ_drip < ρ < 0.5*ρ*₀ ≈ 1.4 × 10¹⁷ kg m⁻³): At the neutron drip density (ρ_drip = 4.3 × 10¹⁴ kg m⁻³), the chemical potential of neutrons inside nuclei reaches zero and unbound neutrons begin "dripping" out of nuclei to form an interstitial ¹S₀ superfluid neutron gas permeating the solid nuclear lattice. Near the base of the inner crust (0.2*ρ*₀ to 0.5*ρ*₀), frustration between short-range nuclear attraction and long-range Coulomb repulsion deforms spherical nuclei ("gnocchi") sequentially into cylindrical rods ("spaghetti"), planar slabs ("lasagna"), cylindrical holes ("tunnels"), and spherical bubbles ("Swiss cheese")—collectively termed nuclear pasta (Ravenhall, Pethick, & Wilson 1983)—with a shear modulus ~10¹⁰ times stronger than steel.
4. Outer Core (Radius ~9–10 km; 0.5*ρ*₀ ≤ ρ ≤ 2*ρ*₀): Nuclei dissolve completely into a homogeneous, electrically neutral beta-equilibrium (n–p–e–μ) quantum Fermi liquid satisfying chemical potential equilibrium μ_n = μ_p + μ_e = μ_p + μ_μ. It consists of ~90–95% neutrons (paired in a ³P₂ anisotropic superfluid state) and ~5–10% protons (paired in a ¹S₀ Type II superconductor state), balanced electrically by relativistic degenerate electrons and muons (μ⁻, which appear once μ_e exceeds the muon rest mass m_μc² = 105.66 MeV near ρ ≈ ρ₀).
5. Inner Core (Radius ~0–3 km; ρ > 2–3*ρ*₀): In the most massive neutron stars, central densities reach 5–8*ρ*₀, where the nucleon chemical potential may exceed strange baryon masses, potentially generating hyperons (Λ⁰, Σ⁻), Bose–Einstein condensates of pions (π⁻) or kaons (K⁻), or a phase transition to deconfined up–down–strange (u–d–s) color-superconducting quark matter.

---

Conservation of Angular Momentum and Magnetic Flux

As a stellar core of radius R_core ≈ 10⁴ km collapses to R_NS ≈ 10 km, conservation of rotational angular momentum (L = Iω ∝ MR²/P) amplifies spin frequency by (R_core/R_NS)² ≈ 10⁶, while conservation of magnetic flux (Φ_B = πR²B) amplifies the surface dipole magnetic field by the same factor of 10⁶:

• Rotation-Powered Pulsars (B ~ 10⁷–10⁹ T): Misalignment between the rotational axis and magnetic dipole axis causes the rotating magnetosphere to radiate electromagnetic dipole power Ė_rot = −*Iωω̇* = (2*π*⁴B²R⁶ sin²α)/(3*μ*₀c³P⁴), sweeping narrow relativistic beams of radio, optical, X-ray, and gamma-ray emission across Earth like a lighthouse. Old neutron stars in low-mass X-ray binaries accrete angular momentum from a companion star, spinning up into millisecond pulsars (P = 1.4–10 ms, B ~ 10⁴–10⁵ T); the fastest known spinner is PSR J1748−2446ad in the globular cluster Terzan 5 at 716 Hz (P = 1.396 ms, equatorial surface velocity ~0.24*c*), bounded below the centrifugal mass-shedding breakup period (P_min ≈ 0.5–0.8 ms). Sudden spin-up glitches (Δ*ν*/ν ~ 10⁻⁹ to 10⁻⁵) occur when quantized Feynman–Onsager superfluid neutron vortices unpin en masse from the inner crust lattice, transferring angular momentum to the solid crust.
• Magnetars (B ~ 10¹⁰–10¹¹ T, or 10¹⁴–10¹⁵ G): Proposed by Robert C. Duncan and Christopher Thompson (1992) to explain Soft Gamma Repeaters (SGRs) and Anomalous X-ray Pulsars (AXPs), magnetars possess magnetic fields amplified by convective proto-neutron star dynamos above the quantum electrodynamic (QED) Schwinger critical field:

B_QED = m_e²c² / ( eℏ ) = 4.414 × 10⁹ T (4.414 × 10¹³ G) at which the electron cyclotron energy ℏω_c equals the electron rest mass m_ec² = 511 keV, inducing vacuum birefringence and magnetic crust-cracking starquakes (such as the December 27, 2004 giant flare of SGR 1806−20, which released ~2 × 10³⁹ J in gamma rays in 0.2 s).

Gravitational Waves and r-Process Kilonova Nucleosynthesis

In 1974, Russell A. Hulse and Joseph H. Taylor Jr. discovered the first binary pulsar, PSR B1913+16 (P_orb = 7.75 hours, M₁ = 1.4398 M_☉, M₂ = 1.3886 M_☉), whose periastron advances by 4.226598° yr⁻¹ (35,000× Mercury's relativistic perihelion precession) and whose orbital period decays at Ṗ_orb = −2.4025 × 10⁻¹² s s⁻¹—matching Albert Einstein's general relativistic quadrupole gravitational-wave emission prediction to within 0.16% (1993 Nobel Prize in Physics).

On August 17, 2017, the LIGO and Virgo laser interferometers detected GW170817, the gravitational-wave chirp from the inspiral and merger of two neutron stars (M_tot = 2.74 M_☉) in the galaxy NGC 4993 (d = 40 Mpc), followed 1.74 seconds later by the short gamma-ray burst GRB 170817A and the optical/infrared kilonova AT 2017gfo. Spectroscopic modeling of AT 2017gfo confirmed that tidal and shock decompression of ~0.05 M_☉ of ultra-neutron-rich neutron-star matter drives rapid neutron-capture (r-process) nucleosynthesis, synthesizing several Earth masses of heavy lanthanides and third-peak actinide/platinum-group elements—including cosmic gold (¹⁹⁷Au), platinum, europium, and uranium.

---

Why don't the free neutrons inside a neutron star beta-decay into protons and electrons?

In free space, an isolated neutron decays via beta decay (n → p + e⁻ + ν̄_e) with a half-life of 10.2 minutes (τ = 878.4 s), releasing a maximum electron energy of 0.782 MeV. Inside a neutron star core, however, degenerate electrons fill all quantum momentum states up to a Fermi energy E_F,e ≈ 100–200 MeV. By the Pauli exclusion principle, a decaying neutron cannot emit a 0.782 MeV electron because every available final electron state is already occupied (Pauli blocking), stabilizing neutrons in beta equilibrium (μ_n = μ_p + μ_e).

What is "nuclear pasta" inside a neutron star's inner crust?

At the base of the inner crust (ρ ≈ 0.5–1.4 × 10¹⁷ kg m⁻³), atomic nuclei are squeezed so closely together that they nearly touch. Competition between short-range attractive strong nuclear forces and long-range repulsive Coulomb forces causes spherical nuclei to fuse into extended geometric phases resembling Italian pasta—progressing with depth from spheres ("gnocchi") to cylindrical rods ("spaghetti"), flat parallel sheets ("lasagna"), and inverted bubble structures ("Swiss cheese") before dissolving into the uniform liquid core.

How did the 2017 neutron star merger GW170817 explain the origin of cosmic gold and platinum?

Core-collapse supernovae often lack the extreme neutron-to-proton ratios required to synthesize the heaviest elements beyond barium. When the two neutron stars in GW170817 collided, tidal tails and disk winds ejected ~0.05 M_☉ of decompressed neutron-rich matter. Free neutrons were captured by seed nuclei vastly faster than beta decay could occur (r-process), producing a radioactive kilonova (AT 2017gfo) whose infrared spectrum confirmed the synthesis of thousands of Earth masses of lanthanides, gold, platinum, and uranium.

---

1. Baade, W., & Zwicky, F. (1934). "Remarks on super-novae and cosmic rays." Physical Review, 46(1), 76–77. https://doi.org/10.1103/PhysRev.46.76.2
2. Oppenheimer, J. R., & Volkoff, G. M. (1939). "On massive neutron cores." Physical Review, 55(4), 374–381. https://doi.org/10.1103/PhysRev.55.374
3. Hewish, A., Bell, S. J., Pilkington, J. D. H., Scott, P. F., & Collins, R. A. (1968). "Observation of a rapidly pulsating radio source." Nature, 217(5130), 709–713. https://doi.org/10.1038/217709a0
4. Abbott, B. P., et al. (LIGO Scientific Collaboration and Virgo Collaboration) (2017). "GW170817: Observation of gravitational waves from a binary neutron star inspiral." Physical Review Letters, 119(16), 161101. https://doi.org/10.1103/PhysRevLett.119.161101