Crab Pulsar

A guest star, a thousand years ago

On July 4, 1054, court astronomers of the Song dynasty recorded a new star — a tiānguān kè xīng, "guest star at Tianguan" — bright enough to be visible in broad daylight for 23 days and seen at night for nearly two years. The Japanese diarist Fujiwara no Sadaie noted it as well; the Chaco Canyon petroglyphs in modern New Mexico may also record it. Then it faded from view, and for the next nine centuries nobody knew quite what they had seen.

What they had seen was a Type II core-collapse supernova in the constellation Taurus, about 2.0 kpc (6,500 light-years) away. The expanding ejecta is now catalogued as Messier 1, the Crab Nebula — the very first object Charles Messier added to his catalogue of "things that are not comets". The Crab Nebula's filaments are still flying outward at roughly 1,500 km/s, and at the dead centre of that expansion is the compact object that survived the explosion: a neutron star spinning so fast it would tear a normal star apart.

Discovery, 1968

The Crab Nebula had been studied for decades as a strong radio and X-ray source, but the existence of its central engine was only suspected. In late 1968 David Staelin and Edward Reifenstein, using the NRAO 300-foot transit telescope at Green Bank, found that the Crab radio emission contained sporadic, sharp pulses. By early 1969 the pulses had been timed to a precise periodicity: 33.3 milliseconds, the shortest pulsar period then known, and obviously associated with the nebula's central point source NP 0532 — soon designated PSR B0531+21 in the standard pulsar catalogue.

Within months, Cocke, Disney and Taylor at the Steward Observatory on Kitt Peak caught the same 33-ms pulse train in optical light — the first non-radio detection of any pulsar. The optical pulses were faint but unmistakable, identifying the south-western of two stars at the nebular centre as the pulsar itself. X-ray pulsations followed within the year, from rocket flights and the early X-ray satellites; gamma-ray pulsations were detected by SAS-2 and COS-B in the 1970s, by EGRET in the 1990s, by Fermi-LAT in the 2000s. In 2008 the MAGIC Cherenkov telescope on La Palma extended the pulsed emission above 100 GeV, where the standard outer-gap and slot-gap models of pulsar magnetospheres had predicted the spectrum should already have cut off. The Crab Pulsar's spectrum, in short, has broken essentially every prediction in turn — and there is still no completely consistent picture of how its magnetosphere makes the pulses we see.

Spin-down: the simplest possible engine

A rotation-powered pulsar is, at heart, a flywheel that radiates. Its observable parameters reduce to the spin period P and how fast P is lengthening, the spin-down rate Ṗ. For the Crab Pulsar:

P  = 33.3 ms        (and lengthening)
Ṗ  = 4.21 × 10⁻¹³ s/s

That Ṗ is the most important number. It means the period grows by about 36 nanoseconds per day — almost too small to imagine, until you compare it with the period itself and realise that an extraordinary torque is at work. The rotational kinetic energy stored in the neutron star is

E_rot = (1/2) I Ω²,    Ω = 2π/P

where I ≈ 10⁴⁵ g cm² is the standard estimate of the neutron-star moment of inertia. The rate at which that energy is being lost — the spin-down luminosity Ė — is

Ė = - dE_rot/dt = - I Ω Ω̇ = 4π² I Ṗ / P³

Plug in numbers and you get

Ė ≈ 4.6 × 10³⁸ erg/s
   ≈ 1.2 × 10⁵ L_☉

120,000 solar luminosities, almost all of it leaving the neutron star not as radiation from its surface but as a relativistic wind of electrons and positrons. This is the Crab Pulsar's engine output, and it matches — within a factor of order unity — the integrated radio-to-gamma-ray luminosity of the entire Crab Nebula. There is no other power source. Without the pulsar, the nebula would have faded into invisibility centuries ago.

Characteristic age vs true age

For a pulsar that has braked by a fixed mechanism (typically taken to be magnetic-dipole radiation in vacuum), the spin period evolves as Ω̇ = -K Ω^n, where n is the braking index. Integrating gives an estimate of the time elapsed since the pulsar was spinning much faster than now — the characteristic age

τ_c = P / [(n - 1) Ṗ]
    = P / (2 Ṗ)      for n = 3 (pure magnetic dipole)

For the Crab Pulsar with P = 33.3 ms and Ṗ = 4.21 × 10⁻¹³:

τ_c = 0.0333 / (2 × 4.21e-13) s
    ≈ 3.96 × 10¹⁰ s
    ≈ 1255 yr

So τ_c ≈ 1240 yr. But the true age, measured by counting back from 1054, is only ~970 yr. The discrepancy is real and informative. Three things conspire:

• The Crab Pulsar's measured braking index is n ≈ 2.51, not 3 — so the dipole-vacuum formula is wrong in detail.
• The pulsar's initial period P₀ at birth was almost certainly not negligible compared to today's; with P₀ ~ 19 ms, the integrated age comes out closer to 1000 yr.
• Pulsar wind torques add extra braking beyond pure dipole radiation, modifying n.

The Crab is therefore one of the very few pulsars where both a characteristic age and an independent true age are known. The 30% mismatch is not a problem — it is the constraint that lets us calibrate spin-down theory.

Magnetic field and the braking torque

Under the magnetic-dipole assumption, the surface dipole field strength B_s can be extracted from P and Ṗ via

B_s ≈ 3.2 × 10¹⁹ (P Ṗ)^(1/2)  G    (P in seconds)
    ≈ 3.78 × 10¹² G for the Crab

About four trillion gauss at the neutron-star surface — a trillion times stronger than the surface field of the Sun. The torque that brakes the rotation comes from electromagnetic radiation by this rotating dipole, augmented by the relativistic particle outflow it sustains; in reality the wind torque is comparable to the dipole torque, which is part of why the measured n falls below 3.

Energy budget: where the 4.6 × 10³⁸ goes

Almost none of the Crab Pulsar's spin-down power emerges as pulsed photons from the magnetosphere itself. The energy budget breaks down roughly as follows:

Channel	Fraction of Ė	Where you see it
Pulsed magnetospheric emission (radio→TeV)	~10⁻³ – 10⁻²	The 33-ms pulse train; main pulse and interpulse
Pulsar wind (cold relativistic e⁻/e⁺)	~0.99	Not directly visible; carries the energy outward
Synchrotron nebula (radio→hard X-ray)	~0.30	The full Crab Nebula glow
Inverse-Compton (GeV–TeV gamma)	~0.01	Imaging Cherenkov telescopes (HESS, MAGIC, VERITAS)

The bulk of the energy travels outward as a magnetised, ultrarelativistic wind (Lorentz factor γ ~ 10⁴ to 10⁶) that is invisible until it terminates in a shock about 0.1 pc from the pulsar — visible as the inner ring of the Chandra X-ray image of the Crab. Downstream of the termination shock, the freshly accelerated electrons and positrons radiate the broadband nebular spectrum.

Glitches and the interior superfluid

Every few years, the Crab Pulsar's spin period suddenly decreases — a glitch — by a few parts in 10⁸, followed by a slow relaxation back toward (but not all the way to) the pre-glitch trend. The Crab has been seen to glitch dozens of times since its discovery; the Vela pulsar is the other archetypal glitcher.

The standard interpretation invokes superfluid neutrons in the inner crust. The crust itself is a rigid crystalline lattice that brakes electromagnetically along with the surface. But the bulk of the neutrons inside the crust pair up into a superfluid, which (because vortex lines must pin to nuclei to transfer angular momentum) does not brake at the same rate. A reservoir of superfluid angular momentum slowly accumulates relative to the crust. When some critical pinning stress is exceeded, vortices unpin in an avalanche; angular momentum jumps from the superfluid to the crust; and we see a sudden spin-up. The slow recovery reflects the re-pinning of the vortex array.

Glitches are a unique probe of the interior physics of neutron stars — they constrain the mass fraction of the superfluid component, the strength of crust–core coupling, and the equation of state. The Crab is the most-monitored glitcher because its young age means relatively frequent activity and a very clean spin-down baseline against which to measure each event.

Pulsations from radio to TeV

The Crab Pulsar is detected as a pulsed source across an almost unbroken stretch of the electromagnetic spectrum. The pulse profile is bimodal: a main pulse (P1) and a secondary "interpulse" (P2) separated by about 145° in phase, both visible at almost every wavelength but with subtly shifting shape and phase. Spectral milestones include:

• Radio (100 MHz – 30 GHz). Three components — precursor, main pulse, interpulse — plus occasional giant pulses with peak flux densities orders of magnitude above the mean, lasting only nanoseconds. The brightest giant pulses are the most luminous coherent radio emission in the observable universe.
• Optical. The first non-radio pulsations ever found in a pulsar (1969); a faint star at V ≈ 16.5, modulated at 33 ms with high amplitude. The Crab is one of only a handful of optically pulsed pulsars.
• X-ray. Pulsed flux dominates the magnetospheric emission; both Chandra and NuSTAR have measured detailed pulse profiles and polarisation.
• GeV gamma. Fermi-LAT detects pulsations cleanly out past 30 GeV with sharp main and interpulse features.
• TeV gamma. Detection above 100 GeV by MAGIC (2008) and VERITAS (2011) ruled out simple outer-gap and slot-gap models, since both predicted a sharp super-exponential cutoff in this range. Pulsed emission has now been reported up to ~1.5 TeV.

No completely accepted theoretical framework yet explains all of this. The current best candidates involve a combination of curvature, synchrotron and inverse-Compton emission from particles accelerated in the current sheets just outside the light cylinder — but the geometry remains under active modelling.

The pulsar wind nebula

The pulsar wind expands relativistically outward through the slower supernova ejecta and is observed to terminate in a standing shock at roughly r_ts ≈ 0.1 pc from the pulsar — about 12 arcseconds at the Crab's distance. Inside r_ts is the cold, fast wind, dim and difficult to observe directly. Outside r_ts is the synchrotron-emitting nebula proper.

Chandra resolves striking structures at this scale: an inner ring of bright X-ray knots tracing the termination shock, two opposing equatorial jets perpendicular to the rotation axis, and an outer torus of dense particle outflow in the equatorial plane. The whole inner nebula varies on timescales of weeks, with knots brightening, fading and migrating — direct evidence that we are watching particle acceleration in action.

Downstream of the shock, electrons and positrons radiate synchrotron from radio (the lowest-energy electrons) through optical, X-ray, and into the MeV gamma-ray band. The same population of relativistic particles inverse-Compton scatters synchrotron photons (the synchrotron self-Compton process) and the cosmic microwave background up to TeV energies. The inner pulsar wind nebula is the brightest steady source of TeV gamma rays in the sky and serves as the standard candle of TeV astronomy.

Why the Crab matters

The Crab Pulsar and its nebula are the closest thing astrophysics has to a complete textbook example of relativistic compact-object physics. Every step of the chain — supernova collapse to neutron star, rapid initial spin, magnetic braking, particle acceleration in the magnetosphere, relativistic wind launching, termination shock, downstream synchrotron and inverse-Compton emission — is on display in a single object, at known distance, with a known true age, and bright enough to be observed in every band of the electromagnetic spectrum.

It is also a working laboratory. The Crab is the calibration source for X-ray detectors (the "Crab unit", 1 Crab = 2.4 × 10⁻⁸ erg cm⁻² s⁻¹ in 2–10 keV, is still in active use). Its glitches probe nuclear superfluidity. Its giant radio pulses constrain coherent emission theory. Its TeV detection forced rewrites of magnetosphere models. And its 33-ms clock — known to a precision better than one part in 10¹⁵ over decadal baselines — is one of the steadiest oscillators in the known universe.

Comparison with other young pulsars

The Crab is one member of a small family of young, energetic, rotation-powered pulsars associated with their own supernova remnants. A handful of others have comparable spin-down luminosities and similar pulsar wind nebulae:

Pulsar	P (ms)	Ė (erg/s)	τ_c (kyr)	Host SNR
Crab (B0531+21)	33.3	4.6 × 10³⁸	1.24	M1 (SN 1054)
PSR B0540-69	50.5	1.5 × 10³⁸	1.6	SNR 0540-69.3 (LMC)
Vela (B0833-45)	89.3	7 × 10³⁶	11	Vela SNR
PSR J0537-6910	16.1	4.8 × 10³⁸	5	N157B (LMC)
PSR B1509-58	150.6	1.8 × 10³⁷	1.6	MSH 15-52

Only J0537-6910 and B0540-69 rival the Crab in spin-down luminosity. The Vela pulsar is the next most monitored glitcher. None has been observed across as wide a spectrum, and only the Crab benefits from a securely dated historical supernova.

Open questions

• What is the pulsar wind magnetisation? Theory predicts σ (the ratio of Poynting to particle energy flux) at the termination shock should be small, of order unity or less, to power efficient particle acceleration. Realistic models of how the wind reaches this state — the so-called σ-problem — remain unresolved.
• Where in the magnetosphere are the TeV photons produced? Detection of pulsations above 1 TeV requires emission well outside the light cylinder. Current-sheet acceleration models can in principle reach these energies, but a self-consistent treatment of geometry, particle content, and the observed spectral shape is still missing.
• What sets the giant-pulse phenomenon? Nanosecond-duration radio bursts with brightness temperatures exceeding 10³⁷ K demand coherent emission, but no single mechanism has been shown to reproduce all the observed statistics.
• How will the Crab age? Continued spin-down will drive the pulsar across the "death line" in the P–Ṗ diagram on timescales of ~10⁷ yr. Long before that, the wind power will fade and the nebula will dissipate, leaving a quiet, isolated neutron star.

Common pitfalls

• Treating τ_c as the true age. Characteristic age is exact only under specific assumptions (n = 3, P₀ ≪ P). For the Crab the assumptions fail at the ~30% level.
• Confusing the pulsar luminosity with the nebula luminosity. The pulsed magnetospheric flux is a small fraction (≲ 1%) of Ė. The remainder is in the wind and radiates downstream in the nebula.
• Identifying the Crab Pulsar with a magnetar. Magnetars are a different class — much higher B (~10¹⁴–10¹⁵ G), much longer P (~5–10 s), and are powered by magnetic-field decay, not rotation. The Crab is the canonical rotation-powered pulsar.
• Reading Ṗ as a measurement uncertainty. Ṗ is the actual rate of period change — a physical observable, not a measurement error. The Crab's Ṗ is so large that it is detected with ~15-digit precision after only weeks of timing.
• Forgetting the wind torque. Pure magnetic-dipole braking gives n = 3; including the relativistic wind torque drops n. The measured Crab braking index of ~2.51 reflects the wind contribution and is therefore expected, not anomalous.