Neutron Stars
Superfluid Vortex Creep and the Pulsar-Glitch Mechanism
On a single day in March 1969, the Vela pulsar suddenly spun faster — its rotation frequency jumped by about 20 microhertz, a fractional leap of roughly 2 parts in a million — even though pulsars are supposed to do nothing but slow down. That instantaneous spin-up, over in less than a minute, was the first recorded glitch, and explaining it forced astrophysicists to conclude that a neutron star is not a solid body at all but a frictionless neutron superfluid hiding beneath a rigid crust.
Superfluid vortex creep is the theory of how that hidden reservoir stores and releases angular momentum. A rotating superfluid cannot rotate as a rigid body; instead it threads itself with a dense forest of quantized vortex lines. When these vortices pin to the crustal lattice, the superfluid stops slowing with the crust, builds up a rotational lag, and then catastrophically dumps its stored spin — a glitch. Between glitches the vortices leak, or creep, thermally over their pinning barriers, and it is this creep that governs the slow post-glitch recovery.
- TypeNeutron-star interior dynamics / superfluid angular-momentum transfer
- RegimeInner crust, ρ ≈ 10¹¹–10¹⁴ g/cm³, T ≈ 10⁷–10⁸ K
- First observedVela glitch, March 1969 (Radhakrishnan & Manchester; Reichley & Downs)
- Typical glitch sizeΔν/ν ≈ 10⁻⁹ (small) to ~10⁻⁶ (Vela giant glitches)
- Key relationCreep rate ∝ exp(−E_p/kT); glitch when lag ω reaches critical ω_cr
- Observed in≳200 pulsars, ~600 glitches; Vela, Crab, PSR J0537−6910
Interactive visualization
Press play, or step through manually. The visualization is yours to drive — try it before reading on.
Watch the 60-second explainer
A condensed visual walkthrough — narrated, captioned, under a minute.
What It Is: A Superfluid That Cannot Spin Smoothly
Roughly a year after a neutron star is born, its interior cools below the critical temperature (T_c ~ 10⁹–10¹⁰ K) at which neutrons pair up via the attractive nuclear interaction and condense into a superfluid — a quantum fluid with zero viscosity. A superfluid cannot rotate like an ordinary fluid. Instead, by the Onsager–Feynman quantization condition, its circulation comes in fixed quanta: each quantized vortex line carries κ = h/2m_n ≈ 2 × 10⁻³ cm²/s of circulation (the factor of 2 because neutrons pair).
To mimic solid-body rotation, the superfluid threads itself with an array of these vortices whose surface density is n_v = 2Ω/κ ≈ 10⁴ vortices per cm² for a pulsar spinning at Ω ~ 70 rad/s (the Vela rate). The star's spin rate is set entirely by how many vortices thread it, and the fluid can only slow down by letting vortices migrate outward and annihilate at the boundary.
- Inner crust: a lattice of neutron-rich nuclei bathed in the free neutron superfluid — the glitch reservoir.
- Pinning: vortices lower their energy by sitting on lattice nuclei, so they get stuck.
The Mechanism: Pinning, Lag, Unpinning, and Creep
The physics was laid out by Anderson & Itoh (1975), building on Baym, Pethick & Pines (1969). The rigid crust, coupled to the star's magnetic dipole, is braked by electromagnetic radiation and steadily spins down. But the crustal superfluid's vortices are pinned to nuclei and cannot move outward, so the superfluid keeps its old, faster spin. A rotational lag ω = Ω_super − Ω_crust builds up.
Each pinned vortex feels a Magnus force per unit length F_M = ρ_s κ ω pushing it outward. When the lag reaches a critical value ω_cr, the Magnus force overwhelms the pinning force and a large population of vortices unpins at once, avalanching outward, transferring angular momentum to the crust — the glitch spins the star up.
Vortex creep (Alpar, Anderson, Pines & Shaham 1984) is the between-glitch behavior: even below ω_cr, thermal fluctuations let vortices hop the pinning barrier E_p at a rate ∝ exp(−E_p/kT). Because this is exponential in the lag, the response is intensely nonlinear — the creep self-regulates the lag and produces the slow post-glitch relaxation.
Key Quantities and a Worked Estimate
The angular-momentum bookkeeping is simple. If a reservoir of moment of inertia I_s, lagging by Δω, suddenly recouples to the crust (moment of inertia I_c), conservation gives the glitch size:
Δν/ν ≈ (I_s/I_c) · (Δω/Ω)
For a Vela giant glitch, Δν/ν ≈ 10⁻⁶. Observed glitch activity implies the coupled superfluid fraction must satisfy ΔI/I ≳ 1.6%, and once entrainment (the crustal lattice dragging superfluid neutrons, effective mass m*/m ~ 5) is included, the requirement rises to ≳ 7%.
- Pinning energy: E_p ~ 1–5 MeV per pinning site.
- Critical lag: ω_cr ~ 10⁻²–10⁻¹ rad/s (a few to tens of turns of accumulated lag).
- Vortex spacing: ~10⁻² cm for Vela.
- Recurrence: Vela glitches every ~2–3 years; the star ratchets its lag back to ω_cr in that time.
The large size demands more superfluid than the thin crust holds — the so-called glitch crisis — a live constraint on the nuclear equation of state.
How It's Observed: Pulsar Timing and Recovery Curves
Glitches are detected through pulsar timing: the arrival times of radio (or gamma-ray) pulses are tracked to sub-microsecond precision over years, and the rotation phase is fit to a smooth spin-down model. A glitch appears as an abrupt, permanent offset in that phase — a step increase in ν, usually accompanied by a step in the spin-down rate Δν̇/ν̇ ~ 10⁻².
The diagnostic signature of vortex creep is the post-glitch recovery: after the jump, part of the glitch relaxes exponentially back over timescales that, for Vela, span minutes to hundreds of days (classic components near ~6 and ~60 days). The fraction that recovers is the healing parameter Q = Δν_decay/Δν_total; small Q (~0.006 for the 2016 Vela glitch) means most of the spin-up is permanent, implying weak recoupling of the reservoir.
- Landmark data: Vela (PSR J0835−4510), the Crab, and PSR J0537−6910 (the most frequent glitcher, ~3 per year).
- The 2016 Vela glitch was caught mid-rise, resolving an overshoot and rapid ~minute-scale coupling.
Vortex Creep Versus Its Cousins
Vortex creep is one of two historic explanations for glitches; distinguishing it from rivals matters.
- Starquake model (Ruderman 1969; Baym & Pines): as the star spins down it becomes less oblate, the rigid crust cracks, and the moment of inertia drops, spinning the star up. Starquakes explain the tiny, rare Crab glitches but fail badly for Vela — the strain energy simply can't recharge fast enough for glitches every few years. Vortex unpinning is now favored for giant glitches.
- Two-component / vortex-avalanche models: the crust plus a tightly-coupled charged component is the 'normal' part; the pinned neutron superfluid is the reservoir. Modern hydrodynamic simulations treat glitches as scale-invariant self-organized-criticality avalanches, matching the observed power-law glitch-size distributions in most pulsars (Vela being a notable regular exception).
- Core superfluid: unlike the crustal creep picture, the neutron-star core hosts a proton superconductor and neutron superfluid coupled on much shorter timescales via electron scattering off magnetized vortices.
Significance, Famous Cases, and Open Questions
Glitches are the only direct probe of neutron-star interior dynamics. Because glitch sizes and recovery times depend on the superfluid moment of inertia, the pinning strength, and the internal temperature, they constrain the dense-matter equation of state, the crust's thickness, and even the star's mass.
- The Vela glitches (first seen 1969; giant events in 2000, 2016, 2024) remain the definitive testbed; the 2016 event's resolved rise time is among the cleanest data on superfluid coupling.
- The glitch crisis: if entrainment is as strong as microscopic calculations suggest, the crust alone cannot supply enough angular momentum, forcing the core superfluid to participate — still unresolved.
- Predicting glitches: the nonlinear creep model even predicts the interval to the next Vela glitch, and 2024 gravitational-wave searches (LIGO/Virgo) looked for a transient signal from a suddenly-deformed star.
Open questions include the true pinning energy, whether glitches emit detectable gravitational waves, and why Vela and PSR J0537−6910 glitch so quasi-periodically while most pulsars glitch randomly.
| Property | Between glitches (vortex creep) | During a glitch (vortex unpinning) |
|---|---|---|
| Vortex motion | Slow thermally-activated hopping over barriers | Sudden collective outward avalanche |
| Angular momentum flow | Gradual, keeps lag near steady state | Rapid dump from superfluid to crust |
| Timescale | Days to years (creep relaxation) | < 1 minute spin-up (often < seconds) |
| Driving factor | exp(−E_p/kT) Boltzmann activation | Lag ω exceeds critical ω_cr ~ 10⁻²–10⁻¹ rad/s |
| Observable | Slow ν̇ recovery, exponential decay | Step jump Δν/ν, spin-down step Δν̇/ν̇ ≈ 10⁻² |
| Benchmark: Vela | Recovery on ~min, days, ~60–100 d scales | Δν/ν ≈ 10⁻⁶, ~3-yr recurrence |
Frequently asked questions
What actually causes a pulsar glitch?
A glitch is a sudden transfer of angular momentum from an interior neutron superfluid to the star's solid crust. The superfluid's quantized vortices pin to crustal nuclei and stop slowing down with the crust, building a rotational lag. When the lag exceeds a critical value, the vortices unpin en masse and dump their stored spin into the crust in under a minute, spinning the star up by Δν/ν of order 10⁻⁹ to 10⁻⁶.
What is vortex creep, as opposed to vortex unpinning?
Vortex unpinning is the sudden avalanche that produces the glitch itself. Vortex creep is the slow, thermally-activated leakage of vortices over their pinning barriers that happens continuously between glitches. Because the creep rate depends exponentially on temperature and lag (∝ exp(−E_p/kT)), it is highly nonlinear and governs the post-glitch relaxation you see in timing data.
Why did the Vela glitch prove neutron stars contain a superfluid?
A glitch is a spin-up, but pulsars only spin down electromagnetically, so a rigid body cannot suddenly rotate faster on its own. The jump, plus a slow multi-day recovery, requires a decoupled interior reservoir that stored angular momentum and released it — exactly the behavior of a pinned superfluid weakly coupled to the crust. The two-timescale recovery matched Baym–Pethick–Pines superfluid predictions.
What is the 'glitch crisis'?
Giant Vela glitches require a superfluid reservoir with at least ~1.6% of the star's moment of inertia. When entrainment — the crustal lattice increasing the neutrons' effective mass by a factor of ~5 — is accounted for, the needed fraction rises to about 7%, more than the thin inner crust can hold. This forces the neutron-star core superfluid to participate, and remains an unsolved constraint on dense-matter physics.
How big and how frequent are glitches?
Fractional spin jumps range from Δν/ν ~ 10⁻¹¹ for tiny events up to ~3 × 10⁻⁶ for the largest Vela and Crab glitches, usually with a spin-down rate step Δν̇/ν̇ ~ 10⁻². Vela glitches every 2–3 years, the Crab far less often and more weakly, and PSR J0537−6910 is the champion at roughly three glitches per year, making it a favorite for gravitational-wave follow-up.
Do pulsar glitches emit gravitational waves?
Possibly, but none have been confirmed. A sudden crustal rearrangement or a transient stellar deformation during a glitch could radiate a short gravitational-wave burst or a longer quasi-monochromatic signal. LIGO and Virgo searched the 2016 and 2024 Vela glitches and set upper limits on the energy released, but no detection has been made — the strain is likely below current sensitivity.