List of Figures

View Image Figure 1:
Venn diagram showing the numbers and locations of the various types of radio pulsars known as of January 2005. The large and small Magellanic clouds are denoted by LMC and SMC.
Watch/download Movie Figure 2: (Movie)
GIF movie showing the rotating neutron star (or “lighthouse”) model for pulsar emission. Animation designed by Michael Kramer.
View Image Figure 3:
The P -P diagram showing the current sample of radio pulsars. Binary pulsars are highlighted by open circles. Theoretical models [64] do not predict radio emission outside the dark blue region. Figure provided by Michael Kramer.
View Image Figure 4:
A variety of integrated pulse profiles taken from the available literature. References: Panels a, b, d, f [107], Panel c [23], Panels e, g, i [165], Panel h [29]. Each profile represents 360 degrees of rotational phase. These profiles are freely available from an on-line database [327].
View Image Figure 5:
Phenomenological models for pulse shape morphology produced by different line-of-sight cuts of the beam. Figure designed by Michael Kramer.
View Image Figure 6:
The sky distribution of pulsars in Galactic coordinates. The plane of the Galaxy is the central horizontal line. The Galactic centre is the midpoint of this line. Millisecond pulsars are indicated in red. Binary pulsars are highlighted by the open circles. The more isotropic distribution of the millisecond and binary pulsars reflects the differences in detecting them out to large distances cf. the normal population (see Section 3).
View Image Figure 7:
Cartoon showing various evolutionary scenarios involving binary pulsars.
View Image Figure 8:
Eccentricity versus orbital period for a sample of 21 low-mass binary pulsars which are not in globular clusters, with the triangles denoting three recently discovered systems [294]. The solid line shows the median of the predicted relationship between orbital period and eccentricity [249]. Dashed lines show 95% the confidence limit about this relationship. The dotted line shows 2 Pb oc e. Figure provided by Ingrid Stairs [294] using an adaptation of the orbital period-eccentricity relationship tabulated by Fernando Camilo.
View Image Figure 9:
Companion mass versus orbital period for binary pulsars showing the whole sample where, in the absence of mass determinations, statistical arguments based on a random distribution of orbital inclination angles (see Section 4.4) have been used to constrain the masses as shown (Panel a), and only those with well determined companion masses (Panel b). The dashed lines show the uncertainties in the predicted relation [307]. This relationship indicates that as these systems finished a period of stable mass transfer due to Roche-lobe overflow, the size and hence period of the orbit was determined by the mass of the evolved secondary star. Figure provided by Marten van Kerkwijk [330].
Watch/download Movie Figure 10: (Movie)
GIF movie showing a simulation following the motion of 100 pulsars in a model gravitational potential of our Galaxy for 200 Myr. The view is edge-on, i.e. the horizontal axis represents the Galactic plane (30 kpc across) while the vertical axis represents ±10 kpc from the plane. This snapshot shows the initial configuration of young neutron stars.
View Image Figure 11:
Left panel: The current sample of all known radio pulsars projected onto the Galactic plane. The Galactic centre is at the origin and the Sun is at (0,8.5) kpc. Note the spiral-arm structure seen in the distribution which is now required by the electron density model [73, 74]. Right panel: Cumulative number of pulsars as a function of projected distance from the Sun. The solid line shows the observed sample while the dotted line shows a model population free from selection effects.
View Image Figure 12:
Left panel: Pulse scattering caused by irregularities in the interstellar medium. The different path lengths and travel times of the scattered rays result in a “scattering tail” in the observed pulse profile which lowers its signal-to-noise ratio. Right panel: A simulation showing the percentage of Galactic pulsars that are likely to be undetectable due to scattering as a function of observing frequency. Low-frequency (<~ 1 GHz) surveys clearly miss a large percentage of the population due to this effect.
View Image Figure 13:
Left panel: A 22.5- min Arecibo observation of the binary pulsar B1913+16. The assumption that the pulsar has a constant period during this time is clearly inappropriate given the quadratic drifting in phase of the pulse during the observation (linear grey scale plot). Right panel: The same observation after applying an acceleration search. This shows the effective recovery of the pulse shape and a significant improvement in the signal-to-noise ratio.
View Image Figure 14:
Dynamic power spectra showing two recent pulsar discoveries in the globular cluster M62 showing fluctuation frequency as a function of time. Figure provided by Adam Chandler.
View Image Figure 15:
Small-number bias of the scale factor estimates derived from a synthetic population of sources where the true number of sources is known. Left panel: An edge-on view of a model Galactic source population. Right panel: The thick line shows NG, the true number of objects in the model Galaxy, plotted against Nobserved, the number detected by a flux-limited survey. The thin solid line shows Nest, the median sum of the scale factors, as a function of Nobs from a large number of Monte Carlo trials. Dashed lines show 25 and 75% percentiles of the N est distribution.
View Image Figure 16:
Beaming fraction plotted against pulse period for four different beaming models: Narayan & Vivekanand 1983 (NV83) [223], Lyne & Manchester 1988 (LM88) [202], Biggs 1990 (JDB90) [34] and Tauris & Manchester 1998 (TM88) [306].
View Image Figure 17:
Left panel: The corrected luminosity distribution (solid histogram with error bars) for normal pulsars. The corrected distribution before the beaming model has been applied is shown by the dot-dashed line. Right panel: The corresponding distribution for millisecond pulsars. In both cases, the observed distribution is shown by the dashed line and the thick solid line is a power law with a slope of -1. The difference between the observed and corrected distributions highlights the severe under-sampling of low-luminosity pulsars.
View Image Figure 18:
The relativistic binary merging plane. Top: Orbital eccentricity versus period for eccentric binary systems involving neutron stars. Bottom: Orbital period distribution for the massive white dwarf-pulsar binaries.
View Image Figure 19:
The current best empirical estimates of the coalescence rates of relativistic binaries involving neutron stars. The individual contributions from each known binary system are shown as dashed lines, while the solid line shows the total probability density function on a logarithmic and (inset) linear scale. The left panel shows DNS binaries [147, 148], while the right panel shows the equivalent results for NS-WD binaries [157].
View Image Figure 20:
Schematic showing the main stages involved in pulsar timing observations.
View Image Figure 21:
A 120 ms window centred on a coherently-dedispersed giant pulse from the Crab pulsar showing high-intensity nanosecond bursts. Figure provided by Tim Hankins [111].
View Image Figure 22:
Timing model residuals for PSR B1133+16. Panel a: Residuals obtained from the best-fitting model which includes period, period derivative, position and proper motion. Panel b: Residuals obtained when the period derivative term is set to zero. Panel c: Residuals showing the effect of a 1-arcmin position error. Panel d: Residuals obtained neglecting the proper motion. The lines in Panels b-d show the expected behaviour in the timing residuals for each effect. Data provided by Andrew Lyne.
View Image Figure 23:
Examples of timing residuals for a number of normal pulsars. Note the varying scale on the ordinate axis, the pulsars being ranked in increasing order of timing “activity”. Data taken from the Jodrell Bank timing program [281, 124]. Figure provided by Andrew Lyne.
View Image Figure 24:
The fractional stability of three millisecond pulsars compared to an atomic clock. Both PSRs B1855+09 and B1937+21 are comparable, or just slightly worse than, the atomic clock behaviour over timescales of a few years [214]. More recent timing of the millisecond pulsar J0437-4715 [126] indicates that it is more stable than the atomic clock.
View Image Figure 25:
Panel a: Keplerian orbital fit to the 669 -day binary pulsar J0407+1607 [190]. Panel b: Orbital fit in the period-acceleration plane for the globular cluster pulsar 47 Tuc S [102].
View Image Figure 26:
Orbital decay in the binary pulsar B1913+16 system demonstrated as an increasing orbital phase shift for periastron passages with time. The GR prediction due entirely to the emission of gravitational radiation is shown by the parabola. Figure provided by Joel Weisberg [336].
View Image Figure 27:
‘Mass-mass’ diagram showing the observational constraints on the masses of the neutron stars in the double pulsar system J0737-3039. Inset is an enlarged view of the small square encompassing the intersection of the tightest constraints. Figure provided by Michael Kramer [161].
View Image Figure 28:
Distribution of neutron star masses as inferred from timing observations of binary pulsars [292]. The vertical dotted line shows the canonical neutron star mass of 1.4Mo .. Figure provided by Ingrid Stairs.
View Image Figure 29:
Top panel: Observed timing residuals for PSR B1855+09. Bottom panel: Simulated timing residuals induced from a putative black hole binary in 3C66B. Figure provided by Rick Jenet [135].
View Image Figure 30:
Summary of the gravitational wave spectrum showing the location in phase space of the pulsar timing array (PTA) and its extension with the Square Kilometre Array (SKA). Figure updated by Michael Kramer [162] from an original design by Richard Battye.