No Arabic abstract
The Vela pulsar is among a number of pulsars which show detectable optical pulsations. We performed optical observations of this pulsar in January and December 2009 with the Iqueye instrument mounted at the ESO 3.5 m New Technology Telescope. Our aim was to perform phase fitting of the Iqueye data, and to measure the optical pulse profile of the Vela pulsar at high time resolution, its absolute phase and rotational period. We calculated for the first time an independent optical timing solution and obtained the most detailed optical pulse profile available to date. Iqueye detected a distinct narrow component on the top of one of the two main optical peaks, which was not resolved in previous observations, and a third statistically significant optical peak not aligned with the radio one. The quality of the Iqueye data allowed us to determine the relative time of arrival of the radio-optical-gamma-ray peaks with an accuracy of a fraction of a millisecond. We compare the shape of the Iqueye pulse profile with that observed in other energy bands and discuss its complex multi-wavelength structure.
We extend profile domain pulsar timing to incorporate wide-band effects such as frequency-dependent profile evolution and broadband shape variation in the pulse profile. We also incorporate models for temporal variations in both pulse width and in the separation in phase of the main pulse and interpulse. We perform the analysis with both nested sampling and Hamiltonian Monte Carlo methods. In the latter case we introduce a new parameterisation of the posterior that is extremely efficient in the low signal-to-noise regime and can be readily applied to a wide range of scientific problems. We apply this methodology to a series of simulations, and to between seven and nine yr of observations for PSRs J1713$+$0747, J1744$-$1134, and J1909$-$3744 with frequency coverage that spans 700-3600MHz. We use a smooth model for profile evolution across the full frequency range, and compare smooth and piecewise models for the temporal variations in DM. We find the profile domain framework consistently results in improved timing precision compared to the standard analysis paradigm by as much as 40% for timing parameters. Incorporating smoothness in the DM variations into the model further improves timing precision by as much as 30%. For PSR J1713+0747 we also detect pulse shape variation uncorrelated between epochs, which we attribute to variation intrinsic to the pulsar at a level consistent with previously published analyses. Not accounting for this shape variation biases the measured arrival times at the level of $sim$30ns, the same order of magnitude as the expected shift due to gravitational-waves in the pulsar timing band.
We present a robust approach to incorporating models for the time-variable broadening of the pulse profile due to scattering in the ionized interstellar medium into profile-domain pulsar timing analysis. We use this approach to simultaneously estimate temporal variations in both the dispersion measure (DM) and scattering, together with a model for the pulse profile that includes smooth evolution as a function of frequency, and the pulsars timing model. We show that fixing the scattering timescales when forming time-of-arrival estimates, as has been suggested in the context of traditional pulsar timing analysis, can significantly underestimate the uncertainties in both DM, and the arrival time of the pulse, leading to bias in the timing parameters. We apply our method using a new, publicly available, GPU accelerated code, both to simulations, and observations of the millisecond pulsar PSR J1643$-$1224. This pulsar is known to exhibit significant scattering variability compared to typical millisecond pulsars, and we find including low-frequency ($< 1$ GHz) data without a model for these scattering variations leads to significant periodic structure in the DM, and also biases the astrometric parameters at the $4sigma$ level, for example, changing proper motion in right ascension by $0.50 pm 0.12$. If low frequency observations are to be included when significant scattering variations are present, we conclude it is necessary to not just model those variations, but also to sample the parameters that describe the variations simultaneously with all other parameters in the model, a task for which profile domain pulsar timing is ideally suited.
A pulsar timing array (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of global phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 millisecond pulsars is being observed at three radio-frequency bands, 50cm (~700 MHz), 20cm (~1400 MHz) and 10cm (~3100 MHz), with observations at intervals of 2 - 3 weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For ten of the 20 pulsars, rms timing residuals are less than 1 microsec for the best band after fitting for pulse frequency and its first time derivative. Significant red timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array (IPTA) and a PTA based on the Square Kilometre Array. We also present an extended PPTA data set that combines PPTA data with earlier Parkes timing data for these pulsars.
The Crab nebula pulsar was observed in 2009 January and December with a novel very fast optical photon counter, Iqueye, mounted at the ESO 3.5 m New Technology Telescope. Thanks to the exquisite quality of the Iqueye data, we computed accurate phase coherent timing solutions for the two observing runs and over the entire year 2009. Our statistical uncertainty on the determination of the phase of the main pulse and the rotational period of the pulsar for short (a few days) time intervals are $approx 1 , mu$s and ~0.5 ps, respectively. Comparison with the Jodrell Bank radio ephemerides shows that the optical pulse leads the radio one by ~240 $mu$s in January and ~160 $mu$s in December, in agreement with a number of other measurements performed after 1996. A third-order polynomial fit adequately describes the spin-down for the 2009 January plus December optical observations. The phase noise is consistent with being Gaussian distributed with a dispersion $sigma$ of $approx 15 , mu$s in most observations, in agreement with theoretical expectations for photon noise-induced phase variability.
Pulsars show two classes of rotational irregularities that can be used to understand neutron-star interiors and magnetospheres: glitches and timing noise. Here we present an analysis of the Vela pulsar spanning nearly 21 yr of observation and including 8 glitches. We identify the relative pulse number of all of the observations between glitches, with the only pulse-number ambiguities existing over glitch events. We use the phase coherence of the timing solution to simultaneously model the timing noise and glitches in a Bayesian framework, allowing us to select preferred models for both. We find the glitches can be described using only permanent and transient changes in spin frequency, i.e., no step changes in frequency derivative. For all of the glitches, we only need two exponentially decaying changes in spin frequency to model the transient components. In contrast to previous studies, we find that the dominant transient components decay on a common $approx$ 1300 d time scale, and that a larger fraction ( $gtrsim 25%$) of glitch amplitudes are associated with these transient components. We also detect shorter-duration transient components of $approx$ 25 d, as previously observed, but are limited in sensitivity to events with shorter durations by the cadence of our observations. The timing noise is well described by a steep power-law process that is independent of the glitches and subdominant to the glitch recovery. The braking index is constrained to be $<$ 8 with 95% confidence. This methodology can be used to robustly measure the properties of glitches and timing noise in other pulsars.