R. T. Edwards

tempo2, a new pulsar-timing package – I. An overview

Contemporary pulsar-timing experiments have reached a sensitivity level where systematic errors introduced by existing analysis procedures are limiting the achievable science. We have developed TEMPO2, a new pulsar-timing package that contains propagation and other relevant effects implemented at the 1-ns level of precision (a factor of ∼100 more precise than previously obtainable). In contrast with earlier timing packages, TEMPO2 is compliant with the general relativistic framework of the IAU 1991 and 2000 resolutions and hence uses the International Celestial Reference System, Barycentric Coordinate Time and up-to-date precession, nutation and polar motion models. TEMPO2 provides a generic and extensible set of tools to aid in the analysis and visualization of pulsar-timing data. We provide an overview of the timing model, its accuracy and differences relative to earlier work. We also present a new scheme for predictive use of the timing model that removes existing processing artefacts by properly modelling the frequency dependence of pulse phase.

tempo2, a new pulsar timing package – II. The timing model and precision estimates

Tempo2 is a new software package for the analysis of pulsar pulse times of arrival. In this paper we describe in detail the timing model used by tempo2, and discuss limitations on the attainable precision. In addition to the intrinsic slow-down behaviour of the pulsar, tempo2 accounts for the effects of a binary orbital motion, the secular motion of the pulsar or binary system, interstellar, Solar system and ionospheric dispersion, observatory motion (including Earth rotation, precession, nutation, polar motion and orbital motion), tropospheric propagation delay, and gravitational time dilation due to binary companions and Solar system bodies. We believe the timing model is accurate in its description of predictable systematic timing effects to better than one nanosecond, except in the case of relativistic binary systems where further theoretical development is needed. The largest remaining sources of potential error are measurement error, interstellar scattering, Solar system ephemeris errors, atomic clock instability and gravitational waves.

Upper Bounds on the Low-Frequency Stochastic Gravitational Wave Background from Pulsar Timing Observations: Current Limits and Future Prospects

Using a statistically rigorous analysis method, we place limits on the existence of an isotropic stochastic gravitational wave background using pulsar timing observations. We consider backgrounds whose characteristic strain spectra may be described as a power-law dependence with frequency. Such backgrounds include an astrophysical background produced by coalescing supermassive black-hole binary systems and cosmological backgrounds due to relic gravitational waves and cosmic strings. Using the best available data, we obtain an upper limit on the energy density per unit logarithmic frequency interval of Ω h2 ≤ 1.9 × 10-8 for an astrophysical background that is 5 times more stringent than the earlier limit of 1.1 × 10-7 found by Kaspi and colleagues. We also provide limits on a background due to relic gravitational waves and cosmic strings of Ω h2 ≤ 2.0 × 10-8 and Ω h2 ≤ 1.9 × 10-8, respectively. All of the quoted upper limits correspond to a 0.1% false alarm rate together with a 95% detection rate. We discuss the physical implications of these results and highlight the future possibilities of the Parkes Pulsar Timing Array project. We find that our current results can (1) constrain the merger rate of supermassive binary black hole systems at high redshift, (2) rule out some relationships between the black hole mass and the galactic halo mass, (3) constrain the rate of expansion in the inflationary era, and (4) provide an upper bound on the dimensionless tension of a cosmic string background.

Precision Timing of PSR J0437?4715: An Accurate Pulsar Distance, a High Pulsar Mass, and a Limit on the Variation of Newton's Gravitational Constant

Analysis of 10 years of high-precision timing data on the millisecond pulsar PSR J0437–4715 has resulted in a model-independent kinematic distance based on an apparent orbital period derivative, P_b, determined at the 1.5% level of precision (D_k = 157.0 ± 2.4 pc), making it one of the most accurate stellar distance estimates published to date. The discrepancy between this measurement and a previously published parallax distance estimate is attributed to errors in the DE200 solar system ephemerides. The precise measurement of P_b allows a limit on the variation of Newtons gravitational constant, |G/G| ≤ 23 × 10^−12 yr^−1. We also constrain any anomalous acceleration along the line of sight to the pulsar to |a⊙/c| ≤ 1.5 × 10^−18 s^−1 at 95% confidence, and derive a pulsar mass, m_(psr) = 1.76 ± 0.20 M⊙, one of the highest estimates so far obtained.

The subpulse modulation properties of pulsars at 92 cm and the frequency dependence of subpulse modulation

Context: A large sample of pulsars has been observed to study their subpulse modulation at an observing wavelength (when achievable) of both 21 and 92 cm using the Westerbork Synthesis Radio Telescope. In this paper we present the 92-cm data and a comparison is made with the already published 21-cm results. Aims: The main goals are to determine what fraction of the pulsars have drifting subpulses, whether those pulsars share some physical properties and to find out if subpulse modulation properties are frequency dependent. Methods: We analysed 191 pulsars at 92 cm searching for subpulse modulation using fluctuation spectra. The sample of pulsars is as unbiased as possible towards any particular pulsar characteristics. Results: For 15 pulsars drifting subpulses are discovered for the first time and 26 of the new drifters found in the 21-cm data are confirmed. We discovered nulling for 8 sources and 8 pulsars are found to intermittently emit single pulses that have pulse energies similar to giant pulses. Another pulsar was shown to exhibit a subpulse phase step. It is estimated that at least half of the total population of pulsars have drifting subpulses when observations with a high enough signal-to-noise ratio would be available. It could well be that the drifting subpulse mechanism is an intrinsic property of the emission mechanism itself, although for some pulsars it is difficult or impossible to detect. Drifting subpulses are in general found at both frequencies, although the chance of detecting drifting subpulses is possibly slightly higher at 92 cm. It appears that the youngest pulsars have the most disordered subpulses and the subpulses become more and more organized into drifting subpulses as the pulsar ages. The modulation indices measured at the two frequencies are clearly correlated, although at 92 cm they are on average possibly higher. At 92 cm the modulation index appears to be correlated with the characteristic age of the pulsar and the complexity parameters as predicted by three different emission models. The correlations with the modulation indices are argued to be consistent with the picture in which the radio emission can be divided in a drifting subpulse signal plus a quasi-steady signal which becomes, on average, stronger at high observing frequencies. The measured values of P3 at the two frequencies are highly correlated, but there is no evidence for a correlation with other pulsar parameters. Appendix A is only available in electronic form at http://www.aanda.org. Table [see full text] is also available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/469/607

MEASURING THE MASS OF SOLAR SYSTEM PLANETS USING PULSAR TIMING

High-precision pulsar timing relies on a solar system ephemeris in order to convert times of arrival (TOAs) of pulses measured at an observatory to the solar system barycenter. Any error in the conversion to the barycentric TOAs leads to a systematic variation in the observed timing residuals; specifically, an incorrect planetary mass leads to a predominantly sinusoidal variation having a period and phase associated with the planets orbital motion about the Sun. By using an array of pulsars (PSRs J0437–4715, J1744–1134, J1857+0943, J1909–3744), the masses of the planetary systems from Mercury to Saturn have been determined. These masses are consistent with the best-known masses determined by spacecraft observations, with the mass of the Jovian system, 9.547921(2) ×10–4 M ☉, being significantly more accurate than the mass determined from the Pioneer and Voyager spacecraft, and consistent with but less accurate than the value from the Galileo spacecraft. While spacecraft are likely to produce the most accurate measurements for individual solar system bodies, the pulsar technique is sensitive to planetary system masses and has the potential to provide the most accurate values of these masses for some planets.

tempo2: a new pulsar timing package – III. Gravitational wave simulation

Analysis of pulsar timing data sets may provide the first direct detection of gravitational waves. This paper, the third in a series describing the mathematical framework implemented into the tempo2 pulsar timing package, reports on using tempo2 to simulate the timing residuals induced by gravitational waves. The tempo2 simulations can be used to provide upper bounds on the amplitude of an isotropic, stochastic, gravitational wave background in our Galaxy and to determine the sensitivity of a given pulsar timing experiment to individual, supermassive, binary black hole systems.

Discovery and timing of the first 8gr8 Cygnus survey pulsars

Context. Since 2004 we have been carrying out a pulsar survey of the Cygnus region with the Westerbork Synthesis Radio Telescope (WSRT) at a frequency of 328 MHz. The survey pioneered a novel interferometric observing mode, termed 8gr8 (eight-grate), whereby multiple simultaneous digital beams provide high sensitivity over a large field of view. Aims. Since the Cygnus region is known to contain OB associations, it is likely that pulsars are formed here. Simulations have shown that this survey could detect 70 pulsars, which would increase our understanding of the radio pulsar population in this region. We also aim to expand the known population of intermittent and rotating radio transient (RRAT)-like pulsars. Methods. In this paper we describe our methods of observation, processing and data analysis, and we present the first results. Our observing method exploits the way a regularly spaced, linear array of telescopes yields a corresponding regularly spaced series of so-called “grating” beams on the sky. By simultaneously forming a modest number (eight) of offset digital beams, we can utilize the entire field of view of each WSRT dish, but retain the coherently summed sensitivity of the entire array. For the processing we performed a large number of trial combinations of period and dispersion measure (DM) using a computer cluster. Results. In the first processing cycle of the WSRT 8gr8 Cygnus Survey, we have discovered three radio pulsars, with spin periods of 1.657, 1.099 and 0.445 s. These pulsars have been observed on a regular basis since their discovery, both in a special follow-up programme as well as in the regular timing programme. The timing solutions are presented in this paper. We also discuss this survey method in the context of the SKA and its pathfinders. Conclusions. We have found three new pulsars using the WSRT. Reprocessing and further analysis of the data will reveal dimmer pulsars, and RRAT-like or intermittent pulsars.

TEMPO2: a New Pulsar Timing Package

TEMPO2 is a new pulsar timing package that provides the precision necessary for modern millisecond timing projects and should supersede all existing pulsar timing packages such as TEMPO and PSRTIME. As TEMPO2 is the only program available that can analyse multiple pulsar datasets simultaneously it will become an integral part of pulsar timing array projects which aim to detect the signatures of gravitational radiation in pulsar timing residuals. In this paper we describe the basic functionality of TEMPO2.

Measuring the mass of solar system planets using pulsar timing

High‐precision pulsar timing relies on a solar system ephemeris in order to convert times of arrival (TOAs) of pulses measured at an observatory to the solar system barycenter. Any error in the conversion to the barycentric TOAs leads to a systematic variation in the observed timing residuals; specifically, an incorrect planetary mass leads to a predominantly sinusoidal variation having a period and phase associated with the planet’s orbital motion about the Sun. By using an array of pulsars (PSRs J0437−4715, J1744−1134, J1857+0943, J1909−3744), the masses of the planetary systems from Mercury to Saturn have been determined. These masses are consistent with the best‐known masses determined by spacecraft observations, with the mass of the Jovian system, 9.547921(2)×10−4M⊙, being significantly more accurate than the mass determined from the Pioneer and Voyager spacecraft, and consistent with but less accurate than the value from the Galileo spacecraft. While spacecraft are likely to produce the most accurat...