Improving the precision of pulsar timing through polarization statistics

At the highest levels of pulsar timing precision achieved to date, experiments are limited by noise intrinsic to the pulsar. This stochastic wideband impulse modulated self-noise (SWIMS) limits pulsar timing precision by randomly biasing the measured times of arrival and thus increasing the root mean square (rms) timing residual. We discuss an improved methodology of removing this bias in the measured times of arrival by including information about polarized radiation. Observations of J0437-4715 made over a one-week interval at the Parkes Observatory are used to demonstrate a nearly 40 per cent improvement in the rms timing residual with this extended analysis. In this way, based on the observations over a 64 MHz bandwidth centred at 1341 MHz with integrations over 16.78 s we achieve a 476 ns rms timing residual. In the absence of systematic error, these results lead to a predicted rms timing residual of 30 ns in one hour integrations; however the data are currently limited by variable Faraday rotation in the Earths ionosphere. The improvement demonstrated in this work provides an opportunity to increase the sensitivity in various pulsar timing experiments, for example pulsar timing arrays that pursue the detection of the stochastic background of gravitational waves. The fractional improvement is highly dependent on the properties of the pulse profile and the stochastic wideband impulse modulated self-noise of the pulsar in question.

High signal-to-noise ratio observations and the ultimate limits of precision pulsar timing

We demonstrate that the sensitivity of high-precision pulsar timing experiments will be ultimately limited by the broadband intensity modulation that is intrinsic to the pulsars stochastic radio signal. That is, as the peak flux of the pulsar approaches that of the system equivalent flux density, neither greater antenna gain nor increased instrumental bandwidth will improve timing precision. These conclusions proceed from an analysis of the covariance matrix used to characterise residual pulse profile fluctuations following the template matching procedure for arrival time estimation. We perform such an analysis on 25 hours of high-precision timing observations of the closest and brightest millisecond pulsar, PSR J0437-4715. In these data, the standard deviation of the post-fit arrival time residuals is approximately four times greater than that predicted by considering the system equivalent flux density, mean pulsar flux and the effective width of the pulsed emission. We develop a technique based on principal component analysis to mitigate the effects of shape variations on arrival time estimation and demonstrate its validity using a number of illustrative simulations. When applied to our observations, the method reduces arrival time residual noise by approximately 20%. We conclude that, owing primarily to the intrinsic variability of the radio emission from PSR J0437-4715 at 20 cm, timing precision in this observing band better than 30 - 40 ns in one hour is highly unlikely, regardless of future improvements in antenna gain or instrumental bandwidth. We describe the intrinsic variability of the pulsar signal as stochastic wideband impulse modulated self-noise (SWIMS) and argue that SWIMS will likely limit the timing precision of every millisecond pulsar currently observed by Pulsar Timing Array projects as larger and more sensitive antennae are built in the coming decades.