Millisecond and Recycled Pulsars

A millisecond pulsar looks, at first sight, like a paradox dressed as a clock. It is a neutron star spinning hundreds of times a second — the fastest known, PSR J1748−2446ad, turns 716 times per second — and rapid rotation ought to signal youth, since pulsars are born fast and slow down. Yet millisecond pulsars are the oldest neutron stars we can see, with ages of $10^{9}$ – $10^{10}$ yr and fields a thousand times weaker than an ordinary pulsar's. The resolution is recycling: an old, radio-dead neutron star in a binary is spun back up by accreting matter and angular momentum from a companion that has swollen to fill its Roche lobe. The evolutionary chain ties together objects once studied in isolation. A pulsar is born with $P\sim 0.1$ s and $B\sim10^{12}$ G, spins down in $\sim10^{7}$ yr, and crosses the "death line" where its radio beam switches off. If a companion later overflows onto it, the system lights up as a low-mass X-ray binary; angular momentum from the inner disc edge torques the star up to milliseconds while its field decays or is buried. When accretion stops, a radio millisecond pulsar switches back on — now with a tiny field, a burnt-out white-dwarf companion, and a spin so steady it rivals an atomic clock. That stability is the payoff: pulsar timing turns these stars into instruments that test general relativity in the strong field, weigh neutron stars to constrain the dense-matter equation of state, and — as pulsar timing arrays — listen for the nanohertz gravitational-wave background.

Key takeaways

Millisecond pulsars (MSPs) are neutron stars with $P\lesssim30$ ms and weak fields $B\sim10^{8}$ – $10^{9}$ G, unlike young pulsars ( $P\sim1$ s, $B\sim10^{12}$ G).
They are old stars "recycled" by accretion in a binary: transferred matter carries angular momentum that spins the star up to near-millisecond periods while burying its field.
The signature is short $P$ with a tiny $\dot P\sim10^{-20}$ — together implying low $B$ — which is why the $P$ – $\dot P$ diagram cleanly separates the recycled population.
Recycling is confirmed by the observed link low-mass X-ray binary $\to$ accreting millisecond X-ray pulsar $\to$ radio MSP with a white-dwarf companion.
Because their spin-down is imperceptible, MSPs are macroscopic clocks measuring orbital dynamics, masses, and nanohertz gravitational waves.
Key idea: an old, radio-dead neutron star is reborn as a fast, stable clock — accretion in a binary transfers enough angular momentum to spin it to millisecond periods while burying its magnetic field, so a low $B$ and short $P$ together signal recycling rather than youth.

Where this lives today

Millisecond pulsars are working instruments. Their steadiest application is timing arrays: NANOGrav and its partners, now the International Pulsar Timing Array, monitor tens of the best MSPs for years and, in 2023, reported the Hellings–Downs correlation signalling a stochastic nanohertz gravitational-wave background — plausibly the hum of supermassive black-hole binaries. Individual systems weigh neutron stars through the Shapiro delay: PSR J0740+6620 at $\approx2.08\,M_\odot$ , combined with NICER's radius measurement, constrains the dense-matter equation of state and the TOV maximum mass. Relativistic binaries such as the double pulsar PSR J0737−3039 confirm general relativity to better than a part in a thousand. Globular clusters forge MSPs in their dense cores — Terzan 5 alone hosts more than forty. And the "spiders," black-widow and redback pulsars whose winds ablate their companions toward planetary masses, show recycling caught destroying the donor that made the pulsar.

The physical picture

Think of accretion as slowly pushing a merry-go-round. Each parcel of gas spiralling in from the companion arrives at the inner disc edge carrying orbital angular momentum $\sim\sqrt{GMr}$ per unit mass and hands it to the star; over $\sim10^{8}$ yr even a modest $\Delta M\sim0.1\,M_\odot$ delivers enough spin to reach milliseconds. The magnetic field is the governor: a strong field flings the flow out to where the disc orbits slowly and barely spins the star up, whereas a weak field lets the disc reach close in, where matter orbits fast, permitting the shortest periods. That single fact — weak field, deep disc, fast spin — is why recycled pulsars have both short $P$ and low $B$ . The first figure walks the life cycle stage by stage.

Figure

The recycling pathway: a young pulsar (P ≈ 0.5 s, B ≈ 10¹² G) spins down and dies, then accretes ΔM ≈ 0.1 M⊙ from its companion at Ṁ ≈ 10⁻⁹ M⊙/yr over ≈ 10⁸ yr, and switches back on as a millisecond pulsar (P ≈ 3 ms, B ≈ 10⁸ G)

The decisive evidence lives on one plot. Because $B\propto\sqrt{P\dot P}$ and $\tau\sim P/\dot P$ , a scatter of period against period-derivative sorts every pulsar by field, age, and evolutionary state at a glance. Young pulsars pile up near $P\approx1$ s, $\dot P\approx10^{-15}$ ; recycled millisecond pulsars form a separate island at short $P$ and tiny $\dot P$ , bounded below-left by the "spin-up line." The second figure is the diagram a pulsar astronomer actually reaches for, with the field relation as diagonal rulings.

Figure

The P–Ṗ diagram, the master diagnostic: normal pulsars cluster near P ≈ 1 s, Ṗ ≈ 10⁻¹⁵ (B ≈ 10¹² G) while recycled millisecond pulsars sit at P ≈ 1.5–30 ms, Ṗ ≈ 10⁻²⁰ (B ≈ 10⁸–10⁹ G); the diagonal rulings are lines of constant surface field B = 3.2×10¹⁹ √(P Ṗ) G

The core equations

The whole subject rests on dipole spin-down plus the torque balance of accretion. The first equation is the workhorse: it converts the measured $P$ and $\dot P$ into an inferred surface field by equating dipole luminosity to the loss of rotational energy.

B_s = 3.2\times10^{19}\,\sqrt{P\,\dot P}\ \ \mathrm{G}\qquad (P\ \text{in s}) \qquad (1)

where: $B_s$ — surface polar magnetic field (Gauss); $P$ — spin period (s); $\dot P$ — period derivative (dimensionless). Reading it: a typical MSP with $P=3$ ms, $\dot P=10^{-20}$ gives $B_s\approx2\times10^{8}$ G — far below a young pulsar's field. The same spin-down defines a characteristic age. It is exact only for pure dipole braking ( $n=3$ ) from an initial spin much faster than today's, so for MSPs it is only an upper bound.

\tau_c \equiv \frac{P}{2\dot P} \qquad (2)

where: $\tau_c$ — characteristic (spin-down) age. Reading it: $P=3$ ms, $\dot P=10^{-20}$ gives $\tau_c\approx5\times10^{9}$ yr — comparable to the Galactic disc's age, flagging it as ancient. During accretion the spin runs up instead, stalling at an equilibrium period where the magnetospheric (Alfvén) radius equals the corotation radius — sharply dependent on field and $\dot M$ .

P_{\rm eq} \approx 2\pi\left(\frac{r_m^{3}}{GM}\right)^{1/2}\propto B^{6/7}\,\dot M^{-3/7}\,M^{-5/7}\,R^{18/7} \qquad (3)

where: $r_m$ — magnetospheric radius where magnetic stress balances the accreting ram pressure; $\dot M$ — accretion rate; $M,R$ — neutron-star mass and radius. Reading it: weak field and high $\dot M$ push $P_{\rm eq}$ toward a millisecond; strong field and low $\dot M$ leave it mildly recycled. The engine that changes the spin is the accretion torque — the rate at which inflowing matter delivers angular momentum at the inner edge.

\dot J = \dot M\,\sqrt{GM\,r_m} \qquad (4)

where: $\dot J$ — spin-up torque; $\sqrt{GM\,r_m}$ — specific angular momentum at $r_m$ . Reading it: integrated over the accretion phase, $\dot J$ builds the final $I\Omega$ . That integration gives the recycling mass — the accreted mass needed to reach a given spin. Its smallness is the crux: a few percent to a few tenths of a solar mass is enough.

\Delta M \sim \frac{I\,\Omega}{\sqrt{GM\,r_m}}=\frac{2\pi I}{P\,\sqrt{GM\,r_m}} \qquad (5)

where: $\Delta M$ — accreted mass to reach spin period $P$ (order of magnitude). Reading it: for $P=3$ ms, $\Delta M\sim0.1\,M_\odot$ , delivered in $\sim10^{8}$ yr at LMXB rates. There is a hard floor on how fast any neutron star can spin: the centrifugal break-up period, where equatorial matter would be flung off. It depends on the equation of state via $R$ .

P_{\rm break} = 2\pi\left(\frac{R^{3}}{GM}\right)^{1/2}\approx0.5\text{–}1\ \mathrm{ms} \qquad (6)

where: $P_{\rm break}$ — minimum period before equatorial mass shedding. Reading it: for $M=1.4\,M_\odot$ , $R=10$ km this is $0.46$ ms; the fastest known pulsar at $1.396$ ms sits safely above, and how close pulsars approach it constrains the equation of state. Once accretion ends, the star radiates its rotational energy through its now-weak magnetosphere, powering the pulsar wind and high-energy emission.

\dot E_{\rm rot} = 4\pi^{2}\,I\,\frac{\dot P}{P^{3}} \qquad (7)

where: $\dot E_{\rm rot}$ — rate of rotational kinetic-energy loss. Reading it: an energetic MSP like B1937+21 radiates $\sim10^{36}$ erg s⁻¹ despite its tiny $\dot P$ , since $P^{-3}$ is enormous at ms spins. Finally, the braking index tests whether the spin-down is pure magnetic-dipole radiation or whether wind torques and field evolution intrude.

n \equiv \frac{\Omega\,\ddot\Omega}{\dot\Omega^{2}},\qquad n=3\ \text{(pure dipole)} \qquad (8)

where: $n$ — braking index. Reading it: a measured $n\neq3$ signals the dipole model is incomplete, which is why $\tau_c$ and $B_s$ are estimates, not precise measurements.

Derivation AThe inferred surface field from dipole spin-down

The goal is to build equation (1) from first principles, so that $B_s=3.2\times10^{19}\sqrt{P\dot P}$ is not a magic constant but a consequence of equating two energies: the rotational energy a pulsar loses as it slows, and the power a rotating magnetic dipole must radiate. This is the most-used relation in pulsar astronomy — every field and age here descends from it.

Step 1. Write the rotational energy $E=\tfrac12 I\Omega^{2}$ with $\Omega=2\pi/P$ , so the spin-down power is $\dot E=I\Omega\dot\Omega$ (negative, since the star slows). Why: rotation is the only large energy store an isolated old neutron star has.

Step 2. A magnetic dipole $m=B_sR^{3}$ inclined to the spin axis radiates $\dot E_{\rm dip}=\tfrac{2}{3c^{3}}B_s^{2}R^{6}\Omega^{4}\sin^{2}\alpha$ ; take the orthogonal rotator $\sin\alpha\to1$ . Why: an accelerating magnetic moment radiates, draining the reservoir.

Step 3. Equate them, $I\Omega|\dot\Omega|=\tfrac{2}{3c^{3}}B_s^{2}R^{6}\Omega^{4}$ , and solve: $B_s^{2}=\tfrac{3c^{3}I}{8\pi^{2}R^{6}}\,P\dot P$ , using $\Omega=2\pi/P$ , $|\dot\Omega|=2\pi\dot P/P^{2}$ . Why: the surviving observable is exactly $P\dot P$ , so constant-field lines are diagonals on the $P$ – $\dot P$ diagram.

Step 4. Insert $I=10^{45}$ g cm², $R=10$ km, $c=3\times10^{10}$ cm s⁻¹; the prefactor collapses to $B_s\approx3.2\times10^{19}\sqrt{P\dot P}$ G with $P$ in seconds. Why: the " $3.2\times10^{19}$ " just encodes the fiducial $I$ and $R$ , so $B_s$ is model-dependent at tens of percent. *Verified: $P=3$ ms, $\dot P=10^{-20}$ give $B_s\approx1.8\times10^{8}$ G, and the Crab's $P=33$ ms, $\dot P=4.2\times10^{-13}$ give $B_s\approx3.8\times10^{12}$ G — both matching standard tabulated values.*

Go deeper · Derivation B · 1 minAngular momentum budget for spinning up a dead pulsar

Derivation BAngular momentum budget for spinning up a dead pulsar

The goal is to show quantitatively that recycling is easy — that only a small fraction of a solar mass, supplied by an ordinary low-mass companion over a normal binary lifetime, carries enough angular momentum to spin a slow, dead neutron star up to milliseconds. We treat the star as a rigid flywheel gaining angular momentum from disc matter at its inner edge, compute the mass required, and check the time is short — the calculation that makes recycling compelling rather than merely plausible.

Step 1. The final angular momentum is $J=I\Omega=2\pi I/P$ ; for $I=10^{45}$ g cm² and $P=3$ ms, $J\approx2.1\times10^{48}$ g cm² s⁻¹. Why: this is the target the accretion torque must add up to.

Step 2. Each gram arrives at the inner edge with specific angular momentum $\ell=\sqrt{GM\,r_m}$ ; taking $r_m\approx R\approx10^{6}$ cm, $M=1.4\,M_\odot$ , $\ell\approx\sqrt{GMR}\approx1.4\times10^{16}$ cm² s⁻¹. Why: matter falling closer orbits faster and carries more spin per gram, so a deep disc (weak field) recycles far more efficiently.

Step 3. Divide: $\Delta M=J/\ell\approx1.5\times10^{32}$ g $\approx0.08\,M_\odot$ . Losses and imperfect coupling raise the practical figure to the standard $0.1$ – $0.2\,M_\odot$ . Why: a $\sim1\,M_\odot$ companion has mass to spare, so recycling is not resource-limited.

Step 4. Time check: at a typical LMXB rate $\dot M\sim10^{-9}\,M_\odot\,\mathrm{yr^{-1}}$ , delivering $0.08\,M_\odot$ takes $t\sim8\times10^{7}$ yr, well short of the donor's lifetime. Why: recycling must fit inside the binary's life, and it does with room to spare. *Verified: $\ell=\sqrt{GMR}=1.36\times10^{16}$ cm² s⁻¹, $\Delta M\approx0.08\,M_\odot$ , $t\approx8\times10^{7}$ yr — consistent with the $0.1$ – $0.2\,M_\odot$ range quoted by Bhattacharya and van den Heuvel (1991).*

Go deeper · Derivation C · 1 minThe spin-up line that bounds the recycled population

Derivation CThe spin-up line that bounds the recycled population

The goal is to derive the locus in the $P$ – $\dot P$ diagram below which no accreting pulsar can be spun — the "spin-up line" — and show that observed millisecond pulsars respect it, a striking confirmation of recycling. The idea: accretion drives the spin to equilibrium ( $r_m=r_c$ ); translating that equilibrium period at the maximum accretion rate into an inferred field gives a line the population must sit above and to the right of.

Step 1. The magnetospheric radius, from balancing dipole magnetic pressure against accretion ram pressure, is $r_m=\xi(\mu^{4}/2GM\dot M^{2})^{1/7}$ , with $\mu=B R^{3}$ and disc-coupling factor $\xi\approx0.5$ . Why: inside $r_m$ the field takes control and funnels matter along field lines.

Step 2. The corotation radius is $r_c=(GMP^{2}/4\pi^{2})^{1/3}$ . Spin-up proceeds while $r_m<r_c$ ; at $r_m=r_c$ the star is at equilibrium; if $r_m>r_c$ the centrifugal barrier ejects matter (the propeller regime). Why: corotation is the gate — inside it the disc deposits angular momentum, outside it the magnetosphere flings matter away.

Step 3. Set $r_m=r_c$ and solve: $P_{\rm eq}\propto\mu^{6/7}\dot M^{-3/7}$ . Capping $\dot M$ at the Eddington rate $\dot M_{\rm Edd}\approx1.5\times10^{-8}\,M_\odot\,\mathrm{yr^{-1}}$ fixes the shortest period at a given field: $P_{\rm min}\simeq1.9\ \mathrm{ms}\,(B/10^{8}\,\mathrm{G})^{6/7}(\dot M/\dot M_{\rm Edd})^{-3/7}$ . Why: the maximum torque gives the lower-left envelope.

Step 4. Convert the field to $\dot P$ with equation (1); the result is a line every radio MSP lies above and to the right of. Why: a population respecting a line derived purely from accretion physics is a fingerprint of spin-up. *Verified: $P_{\rm min}(B=2\times10^{8}\,\mathrm{G})\approx3.4$ ms at Eddington, consistent with the standard spin-up line of Lorimer and Kramer.*

Limiting cases

The equilibrium spin period responds sharply to field and accretion rate, and the two extremes of the population — mildly recycled and fully recycled — fall straight out of equation (3).

Case	Result	Meaning
Strong field / low accretion	$P_{\rm eq}\propto B^{6/7}\dot M^{-3/7}$ large	Magnetosphere extends far out; only mild spin-up. Yields "slow recycled" pulsars, $P\sim20$ – $100$ ms, typically with heavier white-dwarf companions.
Weak field / high accretion	$P_{\rm eq}\sim$ a few ms	Disc penetrates near the surface where matter orbits fast; full recycling to $P\sim1$ – $3$ ms.
Break-up limit	$P_{\rm break}=2\pi\sqrt{R^{3}/GM}\sim0.5$ – $1$ ms	No star can spin faster than centrifugal mass-shedding; the fastest known, PSR J1748−2446ad at $1.396$ ms, sits safely above it.
Post-accretion coasting	$\dot P\propto B^{2}/P^{3}\to$ tiny	Once accretion halts, the low field gives negligible spin-down; the pulsar becomes an ultra-stable clock with $\dot P\sim10^{-20}$ .

Common trapA short spin period alone makes a millisecond pulsar

It does not — and this is the trap the whole $P$ – $\dot P$ diagram exists to defeat. The Crab pulsar spins at $\sim33$ ms, faster than plenty of objects loosely called "millisecond pulsars," yet it is a textbook young pulsar: born in SN 1054, it has a huge $\dot P\approx4\times10^{-13}$ and $B\approx4\times10^{12}$ G. What defines an MSP is short $P$ together with a minuscule $\dot P\sim10^{-20}$ , implying $B\sim10^{8}$ – $10^{9}$ G and hence recycling. A related trap: $\tau_c=P/2\dot P$ is a reasonable order-of-magnitude age for an MSP but is not the recycling time or true age — the real history is dominated by pre-accretion cooling and the accretion phase, which $\tau_c$ ignores.

Worked examples

1 · Field, age, and power of the first millisecond pulsar (fully worked)

PSR B1937+21, the first MSP ever found, has $P=1.5578$ ms and $\dot P=1.05\times10^{-19}$ . Take $I=10^{45}$ g cm². Find its inferred surface field, characteristic age, and spin-down luminosity, and say what each tells you.

Step 1. Field from equation (1): $P\dot P=1.5578\times10^{-3}\times1.05\times10^{-19}=1.64\times10^{-22}$ , so $B_s=3.2\times10^{19}\times\sqrt{1.64\times10^{-22}}\approx\mathbf{4.1\times10^{8}}$ G — four orders of magnitude below a young pulsar, the signature of a buried field.

Step 2. Age from equation (2): $\tau_c=P/2\dot P=1.5578\times10^{-3}/(2.1\times10^{-19})\approx\mathbf{2.4\times10^{8}}$ yr. Old, though B1937+21 is one of the more energetic MSPs; most have $\tau_c\sim10^{9}$ – $10^{10}$ yr.

Step 3. Spin-down power from equation (7): $\dot E_{\rm rot}=4\pi^{2}I\dot P/P^{3}=39.5\times10^{45}\times1.05\times10^{-19}/(3.78\times10^{-9})\approx\mathbf{1.1\times10^{36}}$ erg s⁻¹. Despite the tiny $\dot P$ , the $P^{-3}$ factor makes this enough to light B1937+21 as a gamma-ray source. B_s ≈ 4×10⁸ G, τ_c ≈ 2×10⁸ yr, Ėrot ≈ 10³⁶ erg s⁻¹.

2 · How much mass to spin a dead star to 3 ms

A neutron star ( $M=1.4\,M_\odot$ , $R=10$ km, $I=10^{45}$ g cm²) is spun to $P=3$ ms by accretion delivered near the surface. Estimate the accreted mass and the time at $\dot M=10^{-9}\,M_\odot\,\mathrm{yr^{-1}}$ .

Step 1. Target angular momentum: $J=2\pi I/P=\mathbf{2.1\times10^{48}}$ g cm² s⁻¹.

Step 2. Specific angular momentum at the surface: with $GM=1.86\times10^{26}$ cm³ s⁻², $R=10^{6}$ cm, $\ell=\sqrt{GMR}=\mathbf{1.4\times10^{16}}$ cm² s⁻¹.

Step 3. Required mass: $\Delta M=J/\ell=1.5\times10^{32}$ g $\approx\mathbf{0.08\,M_\odot}$ , rising to $\sim0.1$ – $0.2\,M_\odot$ with realistic losses.

Step 4. Time: $t=\Delta M/\dot M\approx\mathbf{8\times10^{7}}$ yr — comfortably within an LMXB lifetime. A few 0.01–0.1 M⊙ over ≈ 10⁸ yr is enough to reach 3 ms.

3 · The spin-up line versus the break-up floor

For a recycled field $B=2\times10^{8}$ G accreting at the Eddington rate, what is the shortest period accretion can produce, and how does it compare to centrifugal break-up?

Step 1. Spin-up line (derivation C): $P_{\rm min}=1.9\ \mathrm{ms}\,(B/10^{8})^{6/7}=1.9\times2^{6/7}\approx\mathbf{3.4}$ ms — even maximal accretion stalls the star near $3$ ms at this field.

Step 2. Break-up floor from equation (6): $P_{\rm break}=2\pi\sqrt{R^{3}/GM}\approx\mathbf{0.46}$ ms — the equilibrium period is seven times longer than the absolute floor.

Step 3. Interpretation: sub-millisecond spins need a weaker field or super-Eddington torque, which is why the record ( $1.396$ ms) sits between the field-limited value and the break-up floor. The field, not break-up, sets where most MSPs pile up.

Physics of the portrait

The millisecond pulsar was stumbled upon more than predicted, and the person tied to both its discovery and its greatest application is Don Backer — who found the first one and turned the whole class into a gravitational-wave detector.

Origins · 1 minDonald Backer and how this was discovered

Donald Backer

1943–2010 · discovery of millisecond pulsars and pulsar timing

Donald Charles Backer was a radio astronomer at Berkeley whose career bracketed the entire arc of this article. In 1982, following up a steep-spectrum, unusually compact source, Backer, Shrinivas Kulkarni, Carl Heiles, Michael Davis and Miller Goss discovered PSR B1937+21 spinning at 1.558 ms — 642 times a second, twenty times faster than the Crab and unlike anything then known. It seemed to demand either an impossibly stiff equation of state or a new evolutionary channel. The channel had been sketched already — Bisnovatyi-Kogan and Komberg proposed accreting old neutron stars in 1974, and in the same 1982 volume of Nature Alpar, Cheng, Ruderman and Shaham, and independently Radhakrishnan and Srinivasan, laid out the recycling scenario in which mass transfer in a low-mass X-ray binary spins an old, field-decayed star up to milliseconds. Backer then made the larger move: recognizing that an MSP is a clock of near-atomic stability, he and Roger Foster proposed in 1990 that an array of such clocks could detect the correlated timing signature of passing gravitational waves — the pulsar timing array. Dismissed by some as hopelessly ambitious, that idea matured into NANOGrav and its global partners, which in 2023 reported the first evidence for a nanohertz gravitational-wave background. Backer died in 2010, but the two things he is remembered for — the first millisecond pulsar and the timing array — are the beginning and the current frontier of this field.

History

1967

— Bell and Hewish detect the first radio pulsar

1974

— Bisnovatyi-Kogan and Komberg propose recycling of old neutron stars by accretion

1982

— Backer and colleagues discover the first millisecond pulsar PSR B1937+21 at 1.558 ms

1998

— Wijnands and van der Klis catch a neutron star being spun up in the accreting X-ray pulsar SAX J1808.4−3658

2023

— pulsar timing arrays report evidence for a nanohertz gravitational-wave background

Going further

Texts. Carroll and Ostlie, An Introduction to Modern Astrophysics, 2nd ed., ch. 16 · Lorimer and Kramer, Handbook of Pulsar Astronomy (2005), chs. 3, 8 · Shapiro and Teukolsky, Black Holes, White Dwarfs, and Neutron Stars, ch. 10 · Bhattacharya and van den Heuvel, Physics Reports 203, 1 (1991) · Alpar, Cheng, Ruderman and Shaham, Nature 300, 728 (1982) · Backer et al., Nature 300, 615 (1982), first MSP · Wijnands and van der Klis, Nature 394, 344 (1998), first accreting ms X-ray pulsar.

Provenance. The field relation $B_s=3.2\times10^{19}\sqrt{P\dot P}$ and age $\tau_c=P/2\dot P$ are recomputed from equating dipole radiation to $I\Omega\dot\Omega$ with $I=10^{45}$ g cm², $R=10$ km (Carroll and Ostlie ch. 16; Lorimer and Kramer ch. 3); the recycling mass and $P_{\rm eq}\propto B^{6/7}\dot M^{-3/7}$ follow from $r_m=r_c$ (Bhattacharya and van den Heuvel 1991; Alpar et al. 1982). A verifier should confirm: (1) $B_s$ from the dipole balance; (2) $\Delta M\approx0.08\,M_\odot$ from $\ell=\sqrt{GMR}=1.4\times10^{16}$ cm² s⁻¹ — the base migration used $\ell\approx4\times10^{15}$ , an arithmetic slip that inflated $\Delta M$ to $0.2\,M_\odot$ ; (3) $P_{\rm break}=0.46$ ms; (4) B1937+21 giving $B_s\approx4\times10^{8}$ G, $\tau_c\approx2\times10^{8}$ yr, $\dot E_{\rm rot}\approx10^{36}$ erg s⁻¹. Masses (J0740+6620 $\approx2.08\,M_\odot$ ) and the fastest period (J1748−2446ad, $1.396$ ms) are standard measured values.

Prerequisites. Neutron Stars · Pulsars · Pulsar Magnetospheres and Spin-Down · White Dwarfs. Leads to. Equation of State of Dense Matter · Tolman–Oppenheimer–Volkoff Equation · Pulsar Glitches and Neutron Star Superfluidity.

Ask the tutor: Ask me why a millisecond pulsar's low field and short period together prove it was recycled — or to trace one star from young pulsar to timing-array clock.

Why are millisecond pulsars old, when fast rotation usually means youth?

Because their speed was acquired late, not at birth. Among isolated pulsars fast does mean young. But a millisecond pulsar got its spin from accretion long after it had slowed down and gone radio-quiet: a companion evolved, spilled matter onto the star, and the infalling gas spun it back up. So the millisecond spin marks the recycling episode, while the star — and its low, decayed field — is billions of years old. The clean tell is the period derivative: a young fast pulsar has a large $\dot P$ ; an MSP's is a hundred million times smaller.

What is the recycling scenario, in one paragraph?

A neutron star is born as an ordinary pulsar, spins down over $\sim10^{7}$ yr, and dies. Much later its companion swells and overflows its Roche lobe; the transferred gas forms a disc and the system shines as a low-mass X-ray binary. Matter reaching the inner edge spins the star up while accretion decays or buries its field. After $\sim10^{8}$ yr the companion is exhausted — usually leaving a helium white dwarf — accretion stops, and the fast, weak-field star reappears as a radio millisecond pulsar. The observed sequence confirms it: low-mass X-ray binaries, then accreting millisecond X-ray pulsars caught mid-spin-up, then radio MSPs with white-dwarf companions.

Why do most millisecond pulsars have white-dwarf companions?

Because the companion is the burnt-out donor. Spinning the star up over $\sim10^{8}$ yr needs a long, stable phase of mass transfer, which a low-mass star climbing the giant branch supplies gently as it overflows its Roche lobe. By the time enough mass has crossed, the donor has shed its envelope and is left as its degenerate helium (or carbon–oxygen) core — a white dwarf. So the very process that recycles the pulsar consumes the star that fed it. Fully recycled, few-millisecond pulsars typically carry low-mass helium white dwarfs in nearly circular orbits, a direct prediction of recycling.

What sets the absolute fastest a pulsar can spin?

Two limits, with observations between them. The hard floor is centrifugal break-up, $P_{\rm break}\approx0.5$ – $1$ ms, depending on the equation of state through $R$ . But most MSPs never get that fast, because accretion stalls earlier at the field-set equilibrium period $P_{\rm eq}\propto B^{6/7}\dot M^{-3/7}$ — a few milliseconds for realistic fields. The fastest known, PSR J1748−2446ad at $1.396$ ms, lies safely above break-up; the absence of any sub-millisecond pulsar suggests either fields halt the spin-up first or gravitational-wave emission from a slight deformation caps it. Finding one would sharply constrain the equation of state.

Why are millisecond pulsars such good clocks, and how does that find gravitational waves?

A weak-field star spinning hundreds of times a second has an almost imperceptible spin-down, $\dot P\sim10^{-20}$ , so with a timing model you can predict pulse arrival times years ahead to better than a microsecond — rivaling atomic standards. A passing gravitational wave shifts those arrival times by nanoseconds. No single pulsar can be told from clock noise, but a wave imprints a correlated pattern across many pulsars that depends only on their angular separation — the Hellings–Downs curve. Searching an array of the best MSPs for that correlation is what pulsar timing arrays do, and in 2023 they reported the first evidence for a nanohertz gravitational-wave background.

What are black-widow and redback pulsars?

They are millisecond pulsars caught eating their companions — the "spider" binaries. The recycled pulsar's relativistic wind and radiation ablate the companion's surface, driving off gas that eclipses the radio pulses for part of each orbit. Black widows have very-low-mass, semi-degenerate companions ( $\lesssim0.05\,M_\odot$ ) being whittled toward planetary masses; redbacks have heavier, non-degenerate ones ( $\sim0.1$ – $0.4\,M_\odot$ ). They matter because some redbacks swing between accretion-powered X-ray and rotation-powered radio states on human timescales — the "transitional" millisecond pulsars that directly link the low-mass X-ray binary and radio-MSP phases.