3 Tests of GR - 2 PulsarsObservations, and Timing2.2 Pulsar observations

2.3 Pulsar timing 

Once dispersion has been removed, the resultant time series is typically folded modulo the expected pulse period, in order to build up the signal strength over several minutes and to obtain a stable time-averaged profile. The pulse period may not be very easily predicted from the discovery period, especially if the pulsar happens to be in a binary system. The goal of pulsar timing is to develop a model of the pulse phase as a function of time, so that all future pulse arrival times can be predicted with a good degree of accuracy.

The profile accumulated over several minutes is compared by cross-correlation with the ``standard profile'' for the pulsar at that observing frequency. A particularly efficient version of the cross-correlation algorithm compares the two profiles in the frequency domain [130]. Once the phase shift of the observed profile relative to the standard profile is known, that offset is added to the start time of the observation in order to yield a ``Time of Arrival'' (TOA) that is representative of that few-minute integration. In practice, observers frequently use a time- and phase-stamp near the middle of the integration in order to minimize systematic errors due to a poorly known pulse period. As a rule, pulse timing precision is best for bright pulsars with short spin periods, narrow profiles with steep edges, and little if any profile corruption due to interstellar scattering.

With a collection of TOAs in hand, it becomes possible to fit a model of the pulsar's timing behaviour, accounting for every rotation of the neutron star. Based on the magnetic dipole model [108, 53], the pulsar is expected to lose rotational energy and thus ``spin down''. The primary component of the timing model is therefore a Taylor expansion of the pulse phase tex2html_wrap_inline2677 with time t :

  equation99

where tex2html_wrap_inline2681 and tex2html_wrap_inline2683 are a reference phase and time, respectively, and the pulse frequency tex2html_wrap_inline2685 is the time derivative of the pulse phase. Note that the fitted parameters tex2html_wrap_inline2685 and tex2html_wrap_inline2689 and the magnetic dipole model can be used to derive an estimate of the surface magnetic field tex2html_wrap_inline2691 :

  equation104

where tex2html_wrap_inline2693 is the inclination angle between the pulsar spin axis and the magnetic dipole axis, R is the radius of the neutron star (about tex2html_wrap_inline2697), and the moment of inertia is tex2html_wrap_inline2699 . In turn, integration of the energy loss, along with the assumption that the pulsar was born with infinite spin frequency, yields a ``characteristic age'' tex2html_wrap_inline2701 for the pulsar:   equation119

2.3.1 Basic transformation 

Equation (1Popup Equation) refers to pulse frequencies and times in a reference frame that is inertial relative to the pulsar. TOAs derived in the rest frame of a telescope on the Earth must therefore be translated to such a reference frame before Equation (1Popup Equation) can be applied. The best approximation available for an inertial reference frame is that of the Solar System Barycentre (SSB). Even this is not perfect; many of the tests of GR described below require correcting for the small relative accelerations of the SSB and the centre-of-mass frames of binary pulsar systems. But certainly for the majority of pulsars it is adequate. The required transformation between a TOA at the telescope tex2html_wrap_inline2703 and the emission time t from the pulsar is

  equation129

Here tex2html_wrap_inline2707 accounts for the dispersive delay in seconds of the observed pulse relative to infinite frequency; the parameter D is derived from the pulsar's dispersion measure by tex2html_wrap_inline2711, with DM in units of tex2html_wrap_inline2713 and the observing frequency f in MHz. The Roemer term tex2html_wrap_inline2717 takes out the travel time across the solar system based on the relative positions of the pulsar and the telescope, including, if needed, the proper motion and parallax of the pulsar. The Einstein delay tex2html_wrap_inline2719 accounts for the time dilation and gravitational redshift due to the Sun and other masses in the solar system, while the Shapiro delay tex2html_wrap_inline2721 expresses the excess delay to the pulsar signal as it travels through the gravitational well of the Sun - a maximum delay of about tex2html_wrap_inline2723 at the limb of the Sun; see [11] for a fuller discussion of these terms. The terms tex2html_wrap_inline2725, tex2html_wrap_inline2727, and tex2html_wrap_inline2729 in Equation (4Popup Equation) account for similar ``Roemer'', ``Einstein'', and ``Shapiro'' delays within the pulsar binary system, if needed, and will be discussed in Section  2.3.2 below. Most observers accomplish the model fitting, accounting for these delay terms, using the program TEMPO  [110]. The correction of TOAs to the reference frame of the SSB requires an accurate ephemeris for the solar system. The most commonly used ephemeris is the ``DE200'' standard from the Jet Propulsion Laboratory [127]. It is also clear that accurate time-keeping is of primary importance in pulsar modeling. General practice is to derive the time-stamp on each observation from the Observatory's local time standard - typically a Hydrogen maser - and to apply, retroactively, corrections to well-maintained time standards such as UTC(BIPM), Universal Coordinated Time as maintained by the Bureau International des Poids et Mesures in Paris.

2.3.2 Binary pulsars 

The terms tex2html_wrap_inline2725, tex2html_wrap_inline2727, and tex2html_wrap_inline2729 in Equation (4Popup Equation), describe the ``Roemer'', ``Einstein'', and ``Shapiro'' delays within a pulsar binary system. The majority of binary pulsar orbits are adequately described by five Keplerian parameters: the orbital period tex2html_wrap_inline2737, the projected semi-major axis x, the eccentricity e, and the longitude tex2html_wrap_inline2743 and epoch tex2html_wrap_inline2745 of periastron. The angle tex2html_wrap_inline2743 is measured from the line of nodes tex2html_wrap_inline2749 where the pulsar orbit intersects the plane of the sky. In many cases, one or more relativistic corrections to the Keplerian parameters must also be fit. Early relativistic timing models, developed in the first years after the discovery of PSR B1913+16, either did not provide a full description of the orbit (see, e.g., [22]), or else did not define the timing parameters, in a way that allowed deviations from GR to be easily identified (see, e.g., [49, 58]). The best modern timing model [33Jump To The Next Citation Point In The Article, 133Jump To The Next Citation Point In The Article, 43Jump To The Next Citation Point In The Article] incorporates a number of ``post-Keplerian'' timing parameters which are included in the description of the three delay terms, and which can be fit in a completely phenomenological manner. The delays are defined primarily in terms of the phase of the orbit, defined by the eccentric anomaly u and true anomaly tex2html_wrap_inline2753, as well as tex2html_wrap_inline2743, tex2html_wrap_inline2737, and their possible time derivatives. These are related by

eqnarray171

where tex2html_wrap_inline2759 is the reference value of tex2html_wrap_inline2743 at time tex2html_wrap_inline2745 . The delay terms then become:

eqnarray187

Here tex2html_wrap_inline2625 represents the combined time dilation and gravitational redshift due to the pulsar's orbit, and r and s are, respectively, the range and shape of the Shapiro delay. Together with the orbital period derivative tex2html_wrap_inline2771 and the advance of periastron tex2html_wrap_inline2773, they make up the post-Keplerian timing parameters that can be fit for various pulsar binaries. A fuller description of the timing model also includes the aberration parameters tex2html_wrap_inline2775 and tex2html_wrap_inline2777, but these parameters are not in general separately measurable. The interpretation of the measured post-Keplerian timing parameters will be discussed in the context of double-neutron-star tests of GR in Section  4 .



3 Tests of GR - 2 PulsarsObservations, and Timing2.2 Pulsar observations

image Testing General Relativity with Pulsar Timing
Ingrid H. Stairs
http://www.livingreviews.org/lrr-2003-5
© Max-Planck-Gesellschaft. ISSN 1433-8351
Problems/Comments to livrev@aei-potsdam.mpg.de