No Arabic abstract
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.
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 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.
We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsars sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $sim 60%$ compared to a model that includes DM variations and spin-noise only.
A new Bayesian method for the analysis of folded pulsar timing data is presented that allows for the simultaneous evaluation of evolution in the pulse profile in either frequency or time, along with the timing model and additional stochastic processes such as red spin noise, or dispersion measure variations. We model the pulse profiles using `shapelets - a complete ortho-normal set of basis functions that allow us to recreate any physical profile shape. Any evolution in the profiles can then be described as either an arbitrary number of independent profiles, or using some functional form. We perform simulations to compare this approach with established methods for pulsar timing analysis, and to demonstrate model selection between different evolutionary scenarios using the Bayesian evidence. %s The simplicity of our method allows for many possible extensions, such as including models for correlated noise in the pulse profile, or broadening of the pulse profiles due to scattering. As such, while it is a marked departure from standard pulsar timing analysis methods, it has clear applications for both new and current datasets, such as those from the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA).
Pulsar timing is a technique that uses the highly stable spin periods of neutron stars to investigate a wide range of topics in physics and astrophysics. Pulsar timing arrays (PTAs) use sets of extremely well-timed pulsars as a Galaxy-scale detector with arms extending between Earth and each pulsar in the array. These challenging experiments look for correlated deviations in the pulsars timing that are caused by low-frequency gravitational waves (GWs) traversing our Galaxy. PTAs are particularly sensitive to GWs at nanohertz frequencies, which makes them complementary to other space- and ground-based detectors. In this chapter, we will describe the methodology behind pulsar timing; provide an overview of the potential uses of PTAs; and summarise where current PTA-based detection efforts stand. Most predictions expect PTAs to successfully detect a cosmological background of GWs emitted by supermassive black-hole binaries and also potentially detect continuous-wave emission from binary supermassive black holes, within the next several years.