Essential Radio Astronomy

Chapter 6 Pulsars

6.1 Pulsar Properties

Pulsars are magnetized neutron stars that appear to emit periodic short pulses of radio radiation with periods between 1.4 ms and 8.5 s. The radical proposal that neutron stars even exist was made with trepidation by Baade & Zwicky in 1934 [5]: “With all reserve we advance the view that a supernova represents the transition of an ordinary star into a new form of star, the neutron star, which would be the endpoint of stellar evolution. Such a star may possess a very small radius and an extremely high density. As neutrons can be packed much more closely than ordinary nuclei and electrons, the ‘gravitational packing’ energy in a cold neutron star may become very large, and, under certain circumstances, may far exceed the ordinary nuclear packing fractions.”

The name pulsar blends “pulse” and “star,” but pulsars are not pulsating stars. Like lighthouses, they continuously emit rotating beams of radiation and appear to flash each time the beam sweeps across the observer’s line of sight. The pulse periods are quite stable because they equal the rotation periods of massive neutron stars. Even though their radio emission mechanism is not well understood, pulsars have become uniquely valuable astrophysical tools:

  1. 1.

    Neutron stars are physics laboratories sampling extreme conditions—deep gravitational potentials GM/(rc2)1, densities ρ1014gcm-3 exceeding those in atomic nuclei, and magnetic field strengths as high as B1014 or even 1015 gauss—not reproducible on Earth.

  2. 2.

    Pulse periods can be timed with fractional errors as small as 10-16. Accurate pulsar timing permits exquisitely sensitive measurements of quantities such as the power of gravitational radiation emitted by a pulsar in a binary system, neutron-star masses, general relativistic effects in strong gravitational fields, orbital perturbations from binary companions as light as planets, accurate pulsar positions and proper motions, and potentially the distortions of interstellar space produced by long-wavelength gravitational radiation from the mergers of supermassive black holes throughout the universe.

Lorimer and Kramer [70] and Lyne and Graham-Smith [71] have written excellent reference books about pulsars and their astrophysical applications.

6.1.1 Discovery

Pulsars were discovered [50] serendipitously in 1967 on chart recordings (Figure 6.1) obtained during a low-frequency (ν=81 MHz) search for compact extragalactic radio sources that scintillate in the interplanetary plasma, just as stars twinkle in the Earth’s atmosphere.

Figure 6.1: The pulsar CP1919 (for Cambridge Pulsar at right ascension α=19h 19m) first appeared as “scruff” barely distinguishable from interference on a scintillation survey chart record (top). On a higher-speed chart recording of CP1919 (bottom), dips in the upper trace spaced by P1.3s showed that the “scruff” was actually a series of periodic pulses. Image credit: Jocelyn Bell Burnell and Antony Hewish.

This history of this important discovery is a warning against overprocessing data before looking at them, ignoring unexpected signals, and failing to explore the observational “parameter space” covered by the data—here the relevant parameter being time. The strong pulsar PSR B0329+54 was clearly visible in 1954 survey data from the Jodrell Bank 250-foot radio telescope [71], but nobody paid any attention to it at the time. The X-ray pulsar in the Crab Nebula (Figure 8.10) was present in data taken several months before the discovery of radio pulsars, but only after the radio pulsar in the Crab Nebula was announced were the X-ray pulses extracted [39]. As radio instrumentation and data-processing software become more sophisticated, more data are “cleaned up” automatically before they reach the astronomer. Matched filtering that brings out the expected signal usually suppresses the unexpected. Thus, clipping circuits or software remove the strong impulses that are usually caused by terrestrial interference, and integrators smooth out fluctuations shorter than the integration time. Most pulses seen by radio astronomers are just artificial interference from radar, electric cattle fences, etc., and short pulses from sources at interstellar distances imply unexpectedly high brightness temperatures Tb>1025K, which is much higher than the 1012 K upper limit for incoherent electron-synchrotron radiation imposed by synchrotron self-Compton cooling (Section 5.5.3).

Cambridge University graduate student Jocelyn Bell recognized that pulsars are astronomical sources where others had failed because she noticed that some “scruffy” pulses in her chart-recorder data (Figure 6.1) didn’t look like other forms of interference and they reappeared exactly once per sidereal day, indicating an origin outside the Solar System. Fortunately, she and her supervisor, Antony Hewish, “decided initially not to computerize the output because until we were familiar with the behavior of our telescope and receivers we thought it better to inspect the data visually, and because a human can recognize signals of different character whereas it is difficult to program a computer to do so.” Read for the full discovery story in her own words.

6.1.2 Neutron Star Masses and Densities

The sources of the pulses were originally unknown, and even intelligent transmissions by “LGM” (Little Green Men) were seriously suggested as explanations for pulsars. Their short periods imply very compact sources such as white dwarf stars, black holes, and neutron stars; their stable periods rule out black holes. Astronomers were familiar with slowly varying or pulsating emission from stars, but the natural period of a radially pulsating star depends on its mean density ρ and is typically days, not seconds. There is a comparable lower limit to the rotation period P of a gravitationally bound star, set by the requirement that the centrifugal acceleration at its equator not exceed the gravitational acceleration. If a nearly spherical star of mass M and radius R rotates with angular velocity Ω2π/P,

Ω2R<GMR2 (6.1)


P2>(4πR33)3πGM. (6.2)

In terms of the mean density

ρM(4πR33)-1, (6.3)
P>(3πGρ)1/2 (6.4)


ρ>3πGP2. (6.5)

Equation 6.5 gives a conservative lower limit to the mean density because a rapidly spinning star is oblate, which increases the centrifugal acceleration and decreases the gravitational acceleration at its equator.

The first pulsar discovered has a period P=1.3 s, so its mean density is at least


This limit is just consistent with the known densities of white dwarf stars. But soon the faster (P=0.033 s) pulsar in the Crab Nebula was discovered, and its period implies a density much higher than any stable white-dwarf star supported by electron-degeneracy pressure. Also, the Crab Nebula (Figure 8.10) is the remnant of a supernova recorded by Chinese astronomers as a “guest star” in 1054 AD, so the discovery of this pulsar confirmed the suggestion by Baade and Zwicky [5] that neutron stars are the compact remnants of supernovae. The fastest known pulsar has P=1.4×10-3 s implying ρ>1014 g cm-3, the density of atomic nuclei.

A star whose mass is greater than the Chandrasekhar mass

MCh(cG)3/21mp21.4M (6.6)

cannot be supported by electron degeneracy pressure and will collapse to become a neutron star. Equation 6.2 implies the maximum radius

R<(GMP24π2)1/3 (6.7)

of a P=1.4×10-3 s pulsar with mass M1.4M is

R <[6.67×10-8dynecm2g-21.42.0×1033g(1.4×10-3s)24π2]1/3
2×106cm=20km. (6.8)

These numbers motivate the definition of the canonical neutron star as being a uniform density sphere with mass M1.4M, radius R10 km, and moment of inertia I=2MR2/51045gcm2. “Canonical” is perhaps too strong a term for what is really no more than a convention for the typical properties of neutron stars; individual neutron stars can have different masses and radii. The extreme density and pressure turn most of the star into a neutron superfluid that is a superconductor at temperatures up to T109 K. The masses of individual pulsars have been measured with varying degrees of accuracy, and many are close to the canonical 1.4M. The highest accurately measured pulsar masses are M2.0M, and they rule out all currently proposed hyperon or boson condensate equations of state and “free” quarks [34]. Any neutron star of significantly higher mass (M3M in standard models) must collapse and become a black hole.

6.1.3 Magnetic Fields

The Sun and many other stars are known to possess approximately dipolar magnetic fields. Stellar interiors are fully ionized and hence good electrical conductors. Charged particles are constrained to move along magnetic field lines, and magnetic field lines are tied to the charged particles. When a star collapses from a radius 106 km to 10 km, its cross-sectional area a is divided by 1010, its magnetic flux ΦBn^𝑑a (where n^ is the unit vector normal to each infinitesimal surface area da) is conserved, and the magnetic field strength is multiplied by 1010. An initial magnetic field strength B100 G becomes B1012 G after collapse, so young neutron stars should have very strong dipolar fields. The best models of the core-collapse process show that a dynamo effect can generate even stronger magnetic fields. Such dynamos may be able to produce the 10141015 G fields observed in magnetars, which are neutron stars having such strong magnetic fields that their radiation is powered by magnetic field decay. Conservation of angular momentum during collapse increases the rotation rate by about the same factor, 1010, yielding initial rotation periods P0 in the millisecond range. Thus young neutron stars should contain rapidly rotating magnetic dipoles (Figure 6.2).

Figure 6.2: The traditional magnetic dipole model of a pulsar. Electrons and positrons from particle cascades are accelerated in one or more of the “gap” regions of the magnetosphere. They stream along the open magnetic field lines and emit coherent radio emission, and in the highest-energy pulsars potentially X-rays and γ-rays as well. The details of these dynamic, relativistic plasma processes are not well understood. Figure credit: Lorimer and Kramer [70].

6.1.4 Magnetic Dipole Radiation

If a rotating magnetic dipole is inclined by some inclination angle α>0 from the rotation axis, it emits electromagnetic radiation at the rotation frequency. Rewriting the Larmor formula (Equation 2.143) in terms of power radiated by a rotating electric dipole gives

Prad=2q2v˙23c3=23(qr¨sinα)2c3=23(p¨)2c3, (6.9)

where p=qr is the electric dipole moment and p=psinα is its component perpendicular to the rotation axis. J. J. Thomson’s derivation of the Larmor formula in terms of electric field lines (Section 2.7) is also valid for magnetic field lines. The Gaussian CGS units for magnetic and electric field are the same, so the power of magnetic dipole radiation is

Prad=23(m¨)2c3, (6.10)

where m=msinα is the perpendicular component of the magnetic dipole moment. For a uniformly magnetized sphere with radius R and surface magnetic field strength B, the magnitude of the magnetic dipole moment is [56]

m=BR3. (6.11)

If the inclined magnetic dipole rotates with angular velocity Ω=2π/P, then

Prad=23m2Ω4c3=23c3(BR3sinα)2(2πP)4, (6.12)

where P is the pulsar period. This electromagnetic radiation will appear at the very low radio frequency ν=P-1<1 kHz, so low that it cannot propagate through the surrounding ionized nebula or ISM. Magnetic dipole radiation extracts rotational kinetic energy from the neutron star and causes the pulsar period to increase with time. The absorbed radiation deposits energy in the surrounding nebula, the Crab Nebula (Figure 8.10) being a prime example.

6.1.5 Spin-Down Luminosity

The rotational kinetic energy E of a spinning object is related to its moment of inertia I by

E=IΩ22=2π2IP2. (6.13)

The moment of inertia of a small mass element about any rotation axis is its mass multiplied by the square of its radial distance r from the rotation axis. The moment of inertia of a sphere with radius R, mass M, and uniform density ρ=3M/(4πR3) spinning about its z-axis is obtained by summing over its mass elements:

I=-RR(0(R2-z2)1/2ρr2 2πr𝑑r)𝑑z, (6.14)

where r is the distance from the rotation axis and z is the height in cylindrical coordinates centered on the sphere. Then

I=πρ0R(R2-z2)2𝑑z=8πρR515=2MR25. (6.15)

The moment of inertia of a canonical neutron star is

I=2MR2521.42.0×1033g(106cm)251045gmcm2. (6.16)

The rotational kinetic energy of a canonical neutron star with the rotation period P=0.033 s of the Crab pulsar is


As magnetic dipole radiation extracts rotational energy, it slowly increases the period of a pulsar:

P˙dPdt>0. (6.17)

Note that the period derivative P˙ is a dimensionless (seconds per second) pure number. Combining the observed period P and period derivative P˙ yields an estimate of the rate E˙ at which the rotational energy is changing. The quantity

-E˙-dErotdt=-ddt(12IΩ2)=-IΩΩ˙ (6.18)

is called the spin-down luminosity. It is not a measured luminosity; it is the measured loss rate of rotational energy, which is presumed to equal the luminosity of magnetic dipole radiation. The spin-down luminosity is usually expressed in terms of the pulse period P:

Ω=2πP  so  Ω˙=2π(-P-2P˙), (6.19)

and Equation 6.18 becomes

-E˙=4π2IP˙P3. (6.20)

The Crab pulsar has P=0.033 s and P˙=10-12.4. If I=1045gcm2 (Equation 6.16), its spin-down luminosity is


If Prad-E˙105L, the luminosity of the low-frequency (ν=P-130Hz) magnetic dipole radiation from the Crab pulsar is comparable with the entire radio output of our Galaxy! It even exceeds the Eddington luminosity limit (Section 5.4.2) of the neutron star, which is possible because the energy source is not accretion. It greatly exceeds the average radio pulse luminosity of the Crab pulsar, 1030ergs-1. The long wavelength (λ=c/ν107 m) magnetic-dipole radiation is absorbed by and heats up the surrounding Crab Nebula (Figure 8.10), making it the “megawave oven” counterpart of a kitchen microwave oven. The observed bolometric luminosity of the Crab Nebula is comparable with the spin-down luminosity, supporting the model that rotational kinetic energy is converted to magnetic dipole radiation that is absorbed by the surrounding ionized nebula and later reradiated at radio through X-ray wavelengths.

6.1.6 Minimum Magnetic Field Strength

If -E˙Prad, Equations 6.12 and 6.20 can be combined to yield a lower limit to the magnetic field strength B>Bsinα at the neutron star surface, where α is the unknown inclination angle between the rotation and magnetic axes:

Prad=-E˙, (6.21)
23c3(BR3sinα)2(4π2P2)2=4π2IP˙P3, (6.22)
B2=3c3IPP˙24π2R6sin2α, (6.23)
B>(3c3I8π2R6)1/2(PP˙)1/2. (6.24)

Inserting the constants for the canonical pulsar in CGS units yields for the first factor

[3(3×1010cms-1)31045gcm28π2(106cm)6]1/23.2×1019, (6.25)

so the minimum magnetic field strength at the surface of a canonical pulsar is

(Bgauss)>3.2×1019(PP˙s)1/2. (6.26)

It is sometimes called the characteristic magnetic field of a pulsar.

The minimum magnetic field strength of the Crab pulsar (P=0.033 s, P˙=10-12.4) is


This is an amazingly strong magnetic field. Its energy density is

UB=B28π>6×1023ergcm-3. (6.27)

Just 1cm3 of this magnetic field contains over 6×1016J=6×1016Ws2GWyr of energy, the annual output of a large nuclear power station.

6.1.7 Characteristic Age

If the spin-down luminosity equals the magnetic dipole radiation luminosity and (Bsinα) doesn’t change significantly with time, a pulsar’s age τ can be estimated from PP˙ on the further assumption that the pulsar’s initial period P0 was much shorter than its current period. Solving Equation 6.23 for PP˙ shows that

PP˙=8π2R6(Bsinα)23c3I (6.28)

also doesn’t change with time. Rewriting the identity PP˙=PP˙ as PdP=PP˙dt and integrating over the pulsar’s age τ gives

P0PP𝑑P=0τ(PP˙)𝑑t=PP˙0τ𝑑t (6.29)

because PP˙ is constant. Integrating Equation 6.29 gives

P2-P022=PP˙τ. (6.30)

In the limit P02P2, the characteristic age of a pulsar defined by

τP2P˙ (6.31)

should be close to the actual age of the pulsar. The characteristic age depends only on the observables P and P˙; it does not depend on the unknown radius R, moment of inertia I, or perpendicular magnetic field Bsinα. However, a newly formed neutron star with small P0 may be quite oblate and will initially be slowed down by emitting quadrupole gravitational radiation. This can cause the characteristic age to be somewhat larger than the true age for young pulsars.

For example, the characteristic age of the Crab pulsar (P=0.033 s, P˙=10-12.4) is


It is slightly larger than the actual age of the Crab pulsar, which is known to be just under 1000yr because the Crab supernova was observed in 1054 AD.

Figure 6.3: PP˙ Diagram. The PP˙ diagram is useful for following the lives of pulsars, playing a role similar to the Hertzsprung–Russell diagram for ordinary stars. It encodes a tremendous amount of information about the pulsar population and its properties in terms of the observables P (the spin period) and P˙ (the time derivative of the spin period). Each pulsar’s characteristic age τ (Equation 6.31), minimum magnetic field strength B (Equation 6.26), and spin-down luminosity -E˙ (Equation 6.20) is determined by its location on the PP˙ diagram, as indicated by the contour lines for τ, B, and -E˙. Young pulsars in the upper middle of the diagram are often associated with supernova remnants (SNRs) and emit high-energy radiation. They slow down and evolve over millions of years to become one of the large population of P1 s “slow” pulsars, until their rate of rotation is too slow to power the radio emission mechanism, and they disappear from view. Rotating Radio Transients (RRATs) are pulsars that sporadically emit single pulses instead of continuous pulse trains. The millisecond pulsars in the lower left are mostly in binary systems and were “recycled” via accretion from their companion stars.

6.1.8 Braking Index

If magnetic dipole radiation is solely responsible for pulsar spin down, then Prad=-E˙ and Equations 6.12 and 6.18 together imply

Ω˙Ω3. (6.32)

The pulsar braking index n is defined by

Ω˙Ωn=CΩn, (6.33)

where C is the constant of proportionality. The braking index is completely determined by the observables P, P˙, and P¨, so comparing its observed value with the expected n=3 provides a useful check on pulsar spin-down models. The time derivative of Equation 6.33 is

Ω¨=CnΩn-1Ω˙=n(CΩnΩ)Ω˙=n(Ω˙Ω)Ω˙=nΩ˙2Ω, (6.34)


n=ΩΩ¨/Ω˙2. (6.35)

Converting from angular velocities to periods with

n=(2πP)2π(P¨P2+2P˙2P3)(-2πP-2P˙)-2, (6.36)
n=2-PP¨P˙2, (6.37)

the braking index in terms of the observables P, P˙, and P¨.

Braking indices in the range 1.4n<3 have been observed and used to investigate alternative spin-down mechanisms and to make more sophisticated estimates of pulsar ages and initial spin periods [60].

6.1.9 The Lives of Pulsars

Pulsars are born in supernovae and appear in the upper left corner of the pulsar PP˙ diagram (Figure 6.3). If B is conserved and they age as described above, they gradually move to the right and down, along lines of constant B and crossing lines of constant characteristic age. Pulsars with characteristic ages <105 yr are often found in or near recognizable supernova remnants (SNRs). Older pulsars are not, either because their SNRs have faded to invisibility or because the supernova explosions expelled the pulsars with enough speed that they have since escaped from their parent SNRs. The bulk of the pulsar population is older than 105 yr but much younger than the Galaxy (1010 yr). The observed distribution of pulsars in the PP˙ diagram indicates that something changes as pulsars age. One controversial possibility is that the magnetic fields of old pulsars must decay on timescales 107 yr, causing old pulsars to move almost straight down in the PP˙ diagram until either their magnetic field is too weak or their spin rate is too slow to produce radio emission via the normal, and as yet still highly uncertain, emission mechanism. Rotating Radio Transients (RRATs) are pulsars that emit so sporadically that they are more easily detected in searches for single pulses rather than for periodic pulse trains. RRATs with measurable periods usually have P>1s (Figure 6.3). Some RRATs may be old pulsars on the verge of becoming radio quiet.

Almost all short period (P<0.1 s) pulsars having fairly weak magnetic field strengths (B<1011 G) are in binary systems, as evidenced by periodic orbital variations in their observed pulse periods (Figure 6.4). These recycled pulsars have been spun up by accreting mass and angular momentum from their stellar companions, to the point that they emit radio pulses despite their relatively low magnetic field strengths B108 G. (Apparently, currents in the accreting plasma “bury” the magnetic field of the neutron star itself.) The magnetic fields of neutron stars funnel the ionized accreting material onto the magnetic polar caps, which become so hot that they, as well as the hot accretion disk, emit X-rays. We observe these systems as the Low Mass X-ray Binaries (LMXBs). As the neutron stars rotate, their inclined magnetic polar caps can appear and disappear from view, causing periodic fluctuations in X-ray flux, and making some of these systems detectable as X-ray pulsars.

Figure 6.4: Examples of Doppler variations observed in binary systems containing pulsars. Top: The Doppler variations of the globular cluster MSP J1748-2446N in Terzan 5. This pulsar is in a nearly circular orbit (eccentricity e=4.6×10-5) with a companion of minimum mass 0.47 M. The difference between the semimajor and semiminor axes for this orbit is only 51±4 cm! The thick lines show the periods as measured during GBT observations. Bottom: Similar Doppler variations from the highly eccentric binary MSP J0514-4002A in the globular cluster NGC 1851. This pulsar has one of the most eccentric orbits known (e=0.888) and a massive white-dwarf or neutron-star companion.

Millisecond pulsars (MSPs) with low mass (M0.11M) white-dwarf companions typically have orbits with small eccentricities owing to strong tidal dissipation during the accretion phase which spins up the pulsars. The eccentricity e of an elliptical orbit is defined as the ratio of the separation of the foci to the length of the major axis. It ranges between e=0 for a circular orbit and e=1 for a parabolic orbit. Pulsars with extremely eccentric orbits usually have neutron-star companions, indicating that these companions exploded as asymmetric supernovae and nearly disrupted the binary system. Stellar interactions in globular clusters cause a much higher fraction of recycled pulsars per unit mass than in the Galactic disk. These interactions can result in very strange systems such as pulsar–main-sequence-star binaries and MSPs in highly eccentric orbits. In both cases, the original low-mass companion star that recycled the pulsar was ejected in an interaction and replaced by another star.

About 15% of millisecond pulsars are isolated. They were probably recycled via the standard scenario in binary systems, but the energetic MSPs eventually ablated their companions away. Recently, many new “black widow” and “redback” systems have been detected which seem to be doing exactly this. They exhibit eclipses of their radio MSP emission, likely owing to free–free absorption by the ionized gas in the systems being blown off the companion stars.

6.1.10 Emission Mechanisms

The radio pulses originate in the pulsar magnetosphere. Because the neutron star is a spinning magnetic dipole, it acts as a unipolar generator. The total Lorentz force acting on a charged particle is

F=q(E+v×Bc). (6.38)

Charges in the magnetic equatorial region redistribute themselves by moving along closed magnetic field lines until they build up an electrostatic field large enough to cancel the magnetic force and give |F|=0. The induced voltage is about 1016V in MKS units. However, the corotating field lines emerging from the polar caps cross the light cylinder (Figure 6.2), an imaginary cylinder centered on the pulsar and aligned with the rotation axis at whose radius the corotating speed equals the speed of light, so these field lines cannot close. Electrons in the polar cap are magnetically accelerated to very high energies along the open but curved field lines, where the acceleration resulting from the curvature causes them to emit curvature radiation that is strongly polarized in the plane of curvature. As the radio beam sweeps across the line of sight, the plane of polarization is observed to rotate by up to 180 degrees, a purely geometrical effect. Normal “slow” pulsars usually have quite large amounts of typically linear polarization. There are several relativistic effects in the rapidly rotating magnetospheres of millisecond pulsars that complicate this emission-geometry-based polarization, yet millisecond pulsars are typically highly polarized as well. High-energy photons produced by curvature radiation interact with the magnetic field and lower-energy photons to produce electron–positron pairs that radiate additional high-energy photons. The final results of this cascade process are bunches of charged particles that emit at radio wavelengths. Pulsars “die” in the lower right corner of the PP˙ diagram (Figure 6.3) when they have sufficiently low B and high P that the curvature radiation near the polar surface is no longer capable of generating particle cascades.

The extremely high brightness temperatures of pulsars are explained by coherent radiation. The electrons do not radiate as independent charges e; instead bunches of N electrons in volumes whose dimensions are less than a wavelength emit in phase as charges Ne. Larmor’s formula indicates that the power radiated by a charge q is proportional to q2, so the radiation intensity can be N times brighter than incoherent radiation from the same total number N of electrons. Because the coherent volume is smaller at shorter wavelengths, most pulsars have extremely steep radio spectra. Typical pulsar spectral indices (defined using the positive convention, see Equation 4.55) are α-1.7 (i.e., Sν-1.7), although some can be much steeper (α<-3), and a handful are almost flat (α>-0.5).

6.2 Pulsars and the Interstellar Medium

With their sharp and short-duration pulse profiles allowing continuous “on–off” differential measurements, small sizes, and high brightness temperatures, pulsars are unique probes of the interstellar medium (ISM). This section closely follows the discussion in Lorimer and Kramer [70], and Equations 6.39 through 6.42 are derived in Appendix D.

The electrons in the ISM make up a cold plasma whose refractive index is

μ=[1-(νpν)2]1/2, (6.39)

where ν is the frequency of the radio waves, νp is the plasma frequency

νp=(e2neπme)1/28.97kHz(necm-3)1/2, (6.40)

and ne is the electron number density. For a typical ISM value ne0.03cm-3, νp1.5 kHz. If ν<νp then μ is imaginary and radio waves cannot propagate through the plasma.

For propagating radio waves, μ<1 and the group velocity

vg=μc (6.41)

of pulses is less than the vacuum speed of light. For most radio observations νpν, so

vgc(1-νp22ν2). (6.42)

A broadband pulse moves through a plasma more slowly at lower frequencies than at higher frequencies. If the distance to the source is d, the dispersion delay t at frequency ν is

t =0ddlvg-dc (6.43)
=1c0d(1+νp22ν2)𝑑l-dc=(e22πmec)ν-20dne𝑑l. (6.44)

In astronomically convenient units the dispersion delay is

(tsec)4.149×103(DMpccm-3)(νMHz)-2, (6.45)
DM0dne𝑑l (6.46)

in units of pc cm-3 is called the dispersion measure and represents the integrated column density of electrons between the observer and the pulsar.

Figure 6.5: Pulsar dispersion. The gray scale shows the uncorrected dispersive delay (ν-2) as a function of frequency from an observation of the recently discovered pulsar J1800+5034. If the dispersion measure (DM) were higher or the pulse period shorter, then the pulses would wrap in pulse phase (potentially many times) because the data are folded (i.e., averaged) modulo the pulse period. The “de-dispersed” band-integrated pulse profile is shown at the top.

Because pulsar observations almost always cover a wide bandwidth, uncorrected differential delays across the band cause dispersive smearing of the pulsed signal when integrated across the band. If uncorrected, the delays would smear out the integrated pulse in time and make most pulsars undetectable (see Figure 6.5). For pulsar searches, the DM is unknown a priori and is a search parameter much like the pulsar spin frequency. This extra search dimension is one of the primary reasons that pulsar searches are computationally intensive. For high-precision timing observations of pulsars with known DMs, if Nyquist-sampled voltage data (often called “baseband” data) are available, the dispersion may be completely removed from the data by a technique known as coherent dedispersion. The dispersed data are de-convolved using the complex conjugate of the complex transfer function of the ISM that caused the dispersion in the first place.

Dispersion measures can be used to provide distance estimates to pulsars. Crude distances can be calculated for pulsars near the Galactic plane by assuming that ne0.03cm-3. However, several sophisticated models of the Galactic electron-density distribution now exist (e.g., NE2001; Cordes and Lazio [33]) that provide much better (Δd/d30% or less) distance estimates.

The ionized ISM affects pulsar signals in several other ways besides dispersion. Inhomogeneities in the turbulent ISM result in both diffractive and refractive scintillation, or time- and frequency-dependent flux-density variations of the pulsar signal much like the “twinkling” of starlight by Earth’s turbulent atmosphere. Diffractive scintillations occur over typical timescales of minutes to hours and radio bandwidths of kHz to hundreds of MHz, and they can cause more than order-of-magnitude flux-density fluctuations. Refractive scintillations tend to be less than a factor of 2 in amplitude and occur on timescales of weeks.

A phenomenon related to scintillation is pulse broadening caused by scattering of radio waves. ISM inhomogeneities cause multipath propagation of a pulsar signal, whereby some rays travel longer physical distances because they do not follow straight lines to the observer. Such rays are delayed in time relative to those traveling more direct paths, and so cause a strongly frequency-dependent (typically ν-4) exponential-like scattering tail of the pulse (see Figure 6.6). This scatter broadening can greatly decrease both the observed pulsed flux density from a pulsar and its timing precision.

Figure 6.6: Pulse broadening caused by scattering. These data show the strong frequency dependence (typically ν-3 to ν-4.4 depending on the nature of the scattering media) of pulse broadening due to multipath propagation as exhibited by the Galactic Center magnetar J1745-2900 [80]. At each frequency the vertical dotted lines span the pulse phase range containing most of the pulse flux.

Finally, pulsars have broadband continuum spectra, so if there are gas clouds along our line of sight, pulsars can be used to probe the ISM via absorption by spectral lines of Hi or molecules. Such absorption spectra can be used to estimate pulsar distances.

6.3 Pulsar Timing

Pulsars are intrinsically interesting and exotic objects, but much of the best science based on pulsar observations has come from their use as tools via pulsar timing. Pulsar timing is the regular monitoring of the rotation of the neutron star by tracking (nearly exactly) the arrival times of the radio pulses. The key point to remember is that pulsar timing unambiguously accounts for every single rotation of the neutron star over long periods (years to decades) of time. This unambiguous and very precise tracking of rotation phase allows pulsar astronomers to probe the interior physics of neutron stars, make extremely accurate astrometric measurements, uniquely test gravitational theories in the strong-field regime, and possibly within the next few years, directly detect gravitational waves (GWs) (GWs are propagating distortions of space–time) from supermassive black hole binaries.

For pulsar timing, astronomers “fold” (average) the data from many pulses modulo the instantaneous pulse period P. The instantaneous pulse frequency is f=1/P, and the instantaneous pulse phase ϕ is defined by dϕ/dt=f. Pulse phase is usually measured in turns of 2π radians, so 0<ϕ<1.

Averaging over many pulses yields an average pulse profile. Although the shapes of individual pulses vary considerably because pulsar emission is intrinsically a noise process, the shape of the average profile is quite stable.

For timing, the average pulse profile is correlated with a template or model profile so that its phase offset can be determined. When multiplied by the instantaneous pulse period, that phase yields a time offset that can be added to a high-precision reference point on the profile (for example, the left edge of the profile based on the recorded time of the first sample of the observation) to create the time of arrival (TOA).

The precision with which a TOA can be determined is approximately equal to the duration of a sharp pulse feature (e.g., the leading edge) divided by the signal-to-noise ratio of the average profile. It is usually expressed in terms of the width Wf of the pulse features in units of the pulse period P, and the signal-to-noise ratio (S/N) such that σTOAWfP/(S/N). Strong, fast pulsars with narrow pulse profiles provide the most accurate arrival times.

In the nearly inertial frame of the Solar System barycenter (center of mass), the rotation period of a pulsar is nearly constant, so the time-dependent phase ϕ(t) of a pulsar can be approximated by the Taylor expansion

ϕ(t)=ϕ0+f(t-t0)+12f˙(t-t0)2+, (6.47)

where ϕ0 and t0 are arbitrary reference phases and times for each pulsar. The critical constraint for pulsar timing is that the observed rotational phase difference between each of the TOAs must contain an integer number of rotations. Each TOA corresponds to a different time t, so the correct fitting parameters (e.g., f and f˙) must result in a phase change between each pair of TOAs i and j that is an integer number of turns, or Δϕij=n turns. Because all measurements are made with regard to the integrated pulse phase rather than the instantaneous pulse period, the precision with which astronomers can make long-term timing measurements can be quite extraordinary.

With what precision can timing determine the spin frequency f of a pulsar? Because f=dϕ/dt when ϕ is measured in turns, the precision is based on how precisely a change in phase Δϕ can be measured over some time interval ΔT. Typically, ΔT is the length of time (up to several tens of years for many pulsars now) over which a pulsar’s phase has been tracked through regular monitoring. Δϕ is determined principally by the individual TOA precisions, although for some types of measurements a statistical component is important as well because precision improves as the number of measurements N-1/2 if the random errors are larger than the systematic errors. For example, the original millisecond pulsar B1937+21 has pulse period P0.00156 s and the TOA precision is σTOA1μs, which corresponds to a phase error of Δϕ6×10-4 turns. This pulsar has been timed for ΔT>25 years, so

ΔfΔϕΔT=6×10-425yrs107.5syr-1=8×10-13Hz. (6.48)

That pulsar’s spin frequency is 642 Hz, so the absolute frequency error Δf implies that the spin rate is known to 15 significant figures!

Many corrections have to be applied to the observed TOAs before ϕ(t) can be expressed as a Taylor series (Equation 6.47). The arrival of a pulse at an observatory on Earth at topocentric (topocentric means measured from a fixed point on the Earth’s surface) time tt can be corrected to the time t in the nearly inertial Solar System barycentric frame, which we assume to be nearly the same as the time in the frame comoving with the pulsar. Note that the measured pulse rates will differ from the actual pulse rates in the pulsar frame by the Doppler factor resulting from the unknown line-of-sight pulsar velocity as well as unknown relativistic corrections from within the pulsar system itself.

The principal terms in the timing equation are

t=tt-t0+Δclock-ΔDM+ΔR+ΔE+ΔS+ΔR+ΔE+ΔS. (6.49)

As before, t0 is a reference epoch, Δclock represents a clock correction that accounts for differences between the observatory clocks and terrestrial time standards, and ΔDM is the frequency-dependent dispersion delay caused by the ISM. The other Δ terms are delays from within the Solar System and, if the pulsar is in a binary, from within its orbit. The Roemer delay ΔR is the classical light travel time across the Earth’s orbit. Its magnitude is 500cosβ s, where β is the ecliptic latitude of the pulsar (the angle between the pulsar and the ecliptic plane containing the Earth’s orbit around the Sun), and ΔR is the corresponding delay across the orbit of a pulsar in a binary or multiple system. The Einstein delay ΔE accounts for the time dilation from the moving pulsar (and observatory) and the gravitational redshift caused by the Sun and planets or the pulsar and any companion stars. The Shapiro delay ΔS is the extra time required by the pulses to travel through the curved space–time containing the Sun, planets, and pulsar companions. Errors in any of these parameters, as well as other parameters such as f, f˙, and proper motion, give very specific systematic signatures in plots of timing residuals (see Figure 6.7), which are simply the differences between the observed TOAs and the predicted TOAs based on the current timing model parameters.

Figure 6.7: Pulsar timing residuals. The top panel shows a “good” timing solution for a fairly average millisecond pulsar with an rms timing precision of about 5 μs over 4 years. The four remaining panels show how the timing residuals are affected by various timing parameter errors. From top to bottom, an error of 1% in the spin-down rate of the pulsar (causing a quadratic drift in pulsar phase), a position error in either right ascension or declination of only 50 mas (an annual sinusoid reflecting the Earth’s motion), a pulsar proper motion of 10 mas/yr (an annual sinusoid growing linearly with time), or the presence of a planet with the mass and orbital period of Mars around the pulsar.

A good example showing how pulsar timing can be extremely useful is timing-based pulsar astrometry. Pulsar positions on the sky are determined by timing a pulsar over the course of a year as the Earth orbits the Sun and tracking the changing Roemer delay.

The Roemer delay ΔR of a pulsar at ecliptic coordinates λ (ecliptic longitude) and β (ecliptic latitude) is

ΔR500scos(β)cos(θ(t)+λ), (6.50)

where θ(t) is the orbital phase of the Earth with respect to the vernal equinox (the intersection of the ecliptic plane and the celestial equator). Equation 6.50 is only approximate because the Earth’s orbit is not quite circular.

If there is an error in our position estimate, the individual position error components Δλ and Δβ cause a differential time delay Δτ to be present in the timing residuals with respect to the correct Roemer delay:

Δτ500s[cos(β+Δβ)cos(θ(t)+λ+Δλ)-cosβcos(θ(t)+λ)]. (6.51)

If the position errors are small enough that sinxx, cosx1, and ΔβΔλ0, we can use trigonometric angle-sum identities and then simplify to get

Δτ-500s[Δλcosβsin(θ(t)+λ)+Δβsinβcos(θ(t)+λ)]. (6.52)

Applying the trig identity sin(θ(t)+ϕ)=cosϕsinθ(t)+sinϕcosθ(t) to the equation for Δτ, we see that

Acosϕ=-500sΔλcosβ, (6.53)
Asinϕ=-500sΔβsinβ, (6.54)

and therefore

Δλ=-Acosϕ500scosβ, (6.55)
Δβ=-Asinϕ500ssinβ, (6.56)

where A and ϕ are the amplitude (in time units because it corresponds to light-travel time delay) and phase of an error sinusoid that we would observe in the timing residuals from the pulsar (see the middle panel in Figure 6.7).

If the pulsar is located near the ecliptic plane (β0), then cosβ1 and there is maximum timing leverage, and therefore minimum error, in the determination of λ. However, sinβ0 and so the errors on β are huge. A similar problem occurs near the ecliptic pole for λ. VLBI positions give better astrometric accuracy for pulsars near the ecliptic plane or ecliptic pole.

In a timing fit for position, the amplitude A of the error sinusoid in the timing residuals as described above will be determined to an absolute precision ΔA approximately equal to the TOA uncertainty divided by the square-root of the number of observations made of the pulsar (so long as there is at least one year of timing data).

Example. How precisely can we measure the position of a millisecond pulsar at an ecliptic latitude of β45, with typical TOA errors of 2 μs, and where there have been 16 observations over a year? With 2 μs TOAs and 16 independent observations, we should be able to constrain the amplitude of the sinusoidal position error amplitude A to be within about ΔA2μs/16=5×10-7 s. Because both cosβ=sinβ0.71, our errors on both λ and β should be similar: ΔλΔβΔA500scosβ=5×10-7s500s0.71=1.4×10-9radians. There are 206,265 arcseconds (as) per radian, so those errors correspond to errors in both directions of only 290 μas! Even normal pulsars with slow spin periods typically provide astrometric precisions of 0.1 as or better.

For binary pulsars, the time delays across the neutron-star orbit allow for similar high-precision measurements of the pulsar orbital parameters. The binary pulsar Roemer delays comprise up to five Keplerian parameters describing elliptical orbits: the projected semimajor axis x(a1sini)/c, the longitude of periastron ω, the time of periastron passage T0, the orbital period Pb, and the orbital eccentricity e.

Relativistic binaries, particularly those with compact and elliptical orbits, may allow the measurement of up to five post-Keplerian (PK) parameters: the rate of periastron advance ω˙ (elliptical orbits do not close in relativistic theories), the orbital period decay Pb˙ caused by the emission of gravitational radiation, the relativistic γ term describing time dilation and gravitational redshift, and the Shapiro delay terms r (range) and s (shape). An example of a full “timing solution,” listing high-precision spin, astrometric, binary, and two post-Keplerian parameters (r,s), is shown in Figure 6.8.

Figure 6.8: An example millisecond pulsar timing ephemeris for the two-solar-mass neutron star system J1614-2230. Arrival times for this pulsar can be measured at the microsecond level, allowing extremely precise measurements of astrometric, spin, and orbital parameters, as well as a high-precision measurement of the relativistic Shapiro delay, which provided the masses of the neutron star and its white-dwarf companion. The numbers in parentheses for each parameter are the errors in the last digits of those parameters. From Demorest et al. [34].

In any theory of gravity, the five PK parameters are functions only of the pulsar mass m1, the companion mass m2, and the standard five Keplerian orbital parameters. For general relativity, the formulas are

ω˙ =3(Pb2π)-5/3(TM)2/3(1-e2)-1, (6.57)
γ =e(Pb2π)1/3T2/3M-4/3m2(m1+2m2), (6.58)
P˙b =-192π5(Pb2π)-5/3(1+7324e2+3796e4)(1-e2)-7/2T5/3m1m2M-1/3, (6.59)
r =Tm2, (6.60)
s =x(Pb2π)-2/3T-1/3M2/3m2-1. (6.61)

In these equations, TGM/c3=4.925490947μs is the solar mass in time units (which is known much more precisely than either G or M individually), m1, m2, and Mm1+m2 are in solar masses, and ssini (where i is the orbital inclination). If any two of these PK parameters are measured, the masses of the pulsar and its companion can be determined. If more than two are measured, each additional PK parameter yields a different test of a gravitational theory.

For the famous case of the Hulse–Taylor binary pulsar B1913+16, high-precision measurements of ω˙ and γ were first made to determine the masses of the two neutron stars accurately. The Nobel-prize-winning measurement came with the eventual detection of P˙b, which implied that the orbit was decaying in accordance with general relativity’s predictions for the the emission of gravitational radiation.

More recently, the double-pulsar system J0737-3039 was discovered, which contains two pulsar clocks in a more compact orbit (2.4 hrs compared to 7.7 hrs for PSR B1913+16), allowing the measurement of all five PK parameters as well as the pulsar mass ratio R and relativistic spin precession, giving a total of five tests of general relativity. Kramer et al. [62] showed that general relativity is correct at the 0.05% level and measured the masses of the two neutron stars to within 1 part in 104, and the measurements continue to improve (see Figure 6.9).

The double neutron-star systems have indirectly confirmed the existence of gravitational waves by matching the predicted orbital decays as orbital energy is carried away by gravitational radiation. But over the last decade, one of the driving efforts in pulsar astronomy has become the direct detection of GWs using pulsar timing arrays (PTAs). A PTA is an array of MSPs spread over the sky rather than an array of telescopes, and the goal is to use that pulsar array to detect correlated signals in the timing residuals of dozens of MSPs caused by the distortions of interstellar space from nanohertz (i.e., periods of years) GWs passing through our Galaxy.

Because GW emission and propagation is a quadrupolar process in general relativity, in contrast to the dipolar emission of electromagnetic waves from an accelerating electron, GWs cause specific angular correlations in the timing residuals between pairs of pulsars on the sky. Pulsars close together on the sky will be similarly affected by a passing GW, whereas those much farther apart on the sky will be uncorrelated or even negatively correlated by the same GW. That angular pattern is known as the Hellings and Downs curve [49], and detecting it would be the key to confirming that correlated signals in timing residuals are caused by GWs and not by other effects such as clock or planetary ephemeris errors.

The most likely sources of detectable nanohertz GWs are supermassive black hole binaries (with total masses of 1081010 M) and years-long orbital periods after their parent galaxies have merged. A single massive nearby system, or an ensemble of more distant, less massive systems (making a “stochastic background” of GWs), will cause tens-of-nanosecond spatially and temporally correlated systematics in PTA timing residuals. This level of long-term timing precision for the best MSPs is now being achieved.

Three PTA experiments have been working on this endeavor: NANOGrav in North America and the Parkes and European PTAs in Australia and Europe, respectively. Together, they are collaborating in the International Pulsar Timing Array (IPTA), and huge progress toward a detection has been made recently. Current PTA limits, based on upper limits for appropriately correlated low-frequency timing residuals, are beginning to constrain models of galaxy mergers throughout the universe. And given continued improvements in pulsar timing capability and many new high-precision MSPs from recent surveys, a direct detection of GWs seems possible or even likely within the next five years.

Figure 6.9: A recent PSR J0737-3039 mass versus mass diagram showing lines corresponding to the post-Keplerian parameters measured for the system. In this case seven(!) parameters have been measured, including the mass ratio R, since both neutron stars are (or were) pulsar clocks. The dark region in the inset plot is the only portion of parameter space consistent with general relativity, effectively testing it at the 0.05% level [19]. The parameter ΩB is a measurement of the relativistic spin precession of the “B” pulsar in the system from detailed modeling of the eclipses of the “A” pulsar by the magnetosphere of “B” [17].