Neutron Stars and Pulsars: Physics at the Edge of Matter

When a massive star exhausts its nuclear fuel and its core collapses, the collapse takes approximately 0.1 seconds. In that fraction of a second, a volume of matter roughly the size of Earth contracts to a sphere 20 kilometres in diameter. The density of the resulting object — a neutron star — reaches 5 × 10¹⁷ kg/m³, equivalent to squeezing the entire human population into a single cubic millimetre. Gravity at the surface is approximately 200 billion times Earth’s surface gravity. The escape velocity exceeds half the speed of light.

At these densities, individual atoms cannot exist. Electrons are forced into protons, producing neutrons and neutrinos in a process called neutronisation. The resulting stellar remnant is a quantum mechanical object governed by neutron degeneracy pressure — the Pauli exclusion principle applied at stellar scales — rather than thermal pressure. It is not a star in the conventional sense. It is an atomic nucleus the size of a city.

Understanding neutron stars requires physics at the boundary of three disciplines: nuclear physics at densities beyond anything achievable in terrestrial laboratories, general relativity in the strong-field regime, and condensed matter physics of a kind of matter that exists nowhere else in the universe. This is why they attract physicists as interesting objects in their own right — and why GW170817 produced 49 simultaneous journal papers.

Key parameters

Parameter	Value
Typical neutron star mass	1.4 M☉
Typical radius	10–12 km
Core density	~2–8 × nuclear density
Fastest pulsar (PSR J1748-2446ad)	716 Hz
Magnetar surface field	10¹⁴ – 10¹⁵ Gauss
GW170817 chirp mass	1.188 M☉

Formation: The Core Collapse Supernova

A star with initial mass above approximately 8–10 solar masses will end its life in a core collapse supernova. The endpoint depends on the initial mass:

• 8–20 M☉: Likely produces a neutron star
• 20–40 M☉: May produce a neutron star or black hole, depending on mass loss history and metallicity
• >40 M☉ or with low mass loss: Collapse directly to a black hole, possibly without a visible supernova

When the iron core — the final fusion product — reaches the Chandrasekhar mass (~1.4 M☉), electron degeneracy pressure can no longer support it. Core collapse begins. The inner core (approximately 0.5–0.6 M☉) collapses at roughly 0.25c, reaching nuclear density in 0.1 seconds. At nuclear density, the strong nuclear force provides repulsion, halting the collapse abruptly. The outer infalling material rebounds — the “bounce” — generating a shock wave that propagates outward through the star.

The shock initially stalls in the iron core as it loses energy to iron dissociation. The mechanism by which the shock is reinvigorated to produce the observed supernova explosion remains an active area of computational astrophysics research. Neutrino heating — the transfer of a small fraction (~1%) of the ~3 × 10⁵³ J of energy radiated as neutrinos in the first ~10 seconds — is the leading mechanism, but three-dimensional convection, standing accretion shock instabilities (SASI), and other multi-dimensional effects are critical.

The neutron star is born rotating rapidly — conservation of angular momentum from the progenitor’s rotation accelerates it dramatically, as a figure skater drawing in their arms — and intensely magnetised, with surface field strengths of 10⁸ to 10¹⁵ Tesla.

Pulsars: The Universe’s Most Precise Clocks

Jocelyn Bell Burnell, a graduate student at the University of Cambridge working with Antony Hewish, first detected pulsars in July 1967 using the Interplanetary Scintillation Array. The signal was a regular radio pulse repeating every 1.3373 seconds, precise to 10⁻⁷ — regular enough that it was briefly, half-seriously, labelled LGM-1 (Little Green Men 1) before the physical explanation became apparent.

Pulsars are rotating neutron stars emitting electromagnetic radiation from their magnetic poles — poles that are misaligned with the rotation axis, causing the beam to sweep across the sky like a lighthouse. When the beam passes Earth’s line of sight, we detect a pulse. The pulse period is the rotation period of the neutron star.

Newly formed pulsars rotate at hundreds of times per second (millisecond pulsars, recycled by mass accretion from a binary companion, can reach 716 Hz — PSR J1748-2446ad, the fastest known). They slow down as they lose rotational energy to electromagnetic radiation and pulsar wind. The slow-down rate (period derivative Ṗ) is measurable with extraordinary precision: for a pulsar with period P and derivative Ṗ, the derived “characteristic age” is τ = P/2Ṗ.

The Hulse-Taylor binary pulsar (PSR B1913+16), discovered by Russell Hulse and Joseph Taylor in 1974, provided the first indirect evidence for gravitational waves. The two neutron stars in this system orbit each other with a period of 7.75 hours. Taylor and colleagues measured the orbital decay over decades with pulsar timing — the period is decreasing precisely as general relativity predicts for energy loss to gravitational wave emission, at approximately 3.5 × 10⁻¹² seconds per orbit. This agreement between observation and prediction was accurate to one part in 10¹⁴, and earned Hulse and Taylor the 1993 Nobel Prize in Physics.

Pulsar timing arrays (PTAs) — networks of millisecond pulsars monitored for years by radio telescopes — are sensitive to gravitational waves at nanohertz frequencies from supermassive black hole binaries. The NANOGrav collaboration announced in 2023 evidence for a gravitational wave background at these frequencies, consistent with the expected signal from the cosmic population of inspiralling supermassive black hole pairs. This represents a third gravitational wave detection method, complementary to LIGO’s 100 Hz band and LISA’s millihertz band.

The Equation of State: Dense Matter Physics

The central unsolved problem in neutron star physics is the equation of state (EOS) — the relationship between pressure and density inside the star. Above approximately 2–3 times nuclear saturation density (ρ₀ = 2.3 × 10¹⁷ kg/m³), the composition and behaviour of neutron star matter is genuinely unknown.

Proposed phases include:

• Hadronic matter: A densely packed fluid of neutrons and protons
• Hyperonic matter: Additional strange baryons (Λ, Σ, Ξ) appearing at high density as it becomes energetically favourable to create strange quarks
• Kaon condensate: A Bose-Einstein condensate of kaons (K⁻) at sufficiently high density
• Quark matter: Deconfined quarks forming colour-superconducting phases (CFL, 2SC) at extreme densities possibly present in the innermost core

The observational constraint is the neutron star mass-radius relationship. A stiff EOS (strong repulsion at high density) produces larger, higher-maximum-mass stars; a soft EOS produces smaller stars with lower maximum mass. Observed neutron star masses provide a lower bound on the maximum mass supported by the EOS.

The discovery of a 2 solar mass neutron star (PSR J1614-2230, Demorest et al. 2010; subsequently PSR J0952-0607 at 2.35 M☉, Romani et al. 2022) ruled out soft EOSs that predicted maximum masses below 2 M☉ — excluding many hyperonic and kaon condensate models. The equation of state must be stiff enough to support 2+ solar mass stars.

GW170817 provided complementary information via tidal deformability — how much each neutron star distorts in the gravitational field of its companion during inspiral. More compact, stiffer stars deform less; the phase evolution of the gravitational wave signal encodes this effect. The LIGO/Virgo constraint on tidal deformability Λ was consistent with neutron star radii of approximately 11–13 km, narrowing the EOS significantly and disfavouring both the very stiff and very soft extremes.

Magnetars: When Magnetic Fields Become Astrophysical Weapons

A magnetar is a neutron star with a surface magnetic field of 10¹³–10¹⁵ Tesla — approximately 1,000 times stronger than a typical pulsar. At these field strengths, the magnetic energy density exceeds the rest mass energy density of electrons; the quantum electrodynamic vacuum is no longer linear; X-ray photons split and merge. The physics requires non-perturbative QED.

On 27 December 2004, the magnetar SGR 1806-20 released more energy in 0.2 seconds than the Sun emits in 250,000 years. The gamma-ray burst briefly ionised Earth’s upper atmosphere from 50,000 light-years away — the most intense astrophysical event observed in historical records to have physically affected Earth. Had this magnetar been within 10 light-years, the effects would have been catastrophic.

Magnetars are responsible for the “soft gamma repeaters” (SGRs) and “anomalous X-ray pulsars” (AXPs) — energetically extreme transients powered by sudden releases of magnetically stored energy in the neutron star crust (“starquakes”) rather than by rotation or accretion.

From Stellar Core to Gold

The astrophysical r-process — rapid neutron capture nucleosynthesis — produces approximately half of all elements heavier than iron in the universe, including virtually all of the gold, platinum, and uranium. For decades, the primary candidate site was core collapse supernovae. GW170817 changed this.

The optical transient associated with GW170817 — the kilonova — showed spectroscopic evidence of r-process nucleosynthesis at luminosities consistent with the production of a substantial mass (0.03–0.05 M☉) of heavy elements. Strontium was identified directly; the broader light curve shape was consistent with a mixture of “light” r-process elements (first r-process peak, A ~ 90) and “heavy” elements (third r-process peak, A ~ 195, including gold and platinum).

The merger of two neutron stars — systems that take billions of years to spiral together after formation — is now confirmed as a major r-process site. Every kilogram of gold on Earth was synthesised in a neutron star collision billions of years before the solar system formed, ejected into the interstellar medium, and eventually incorporated into the molecular cloud that became the Sun and planets.

For how we detect the gravitational wave signals that revealed this process, see gravitational waves: how LIGO listens to spacetime.