We present timing models for 20 millisecond pulsars in the Parkes Pulsar Timing Array. The precision of the parameter measurements in these models has been improved over earlier results by using ...longer data sets and modelling the non-stationary noise. We describe a new noise modelling procedure and demonstrate its effectiveness using simulated data. Our methodology includes the addition of annual dispersion measure (DM) variations to the timing models of some pulsars. We present the first significant parallax measurements for PSRs J1024−0719, J1045−4509, J1600−3053, J1603−7202, and J1730−2304, as well as the first significant measurements of some post-Keplerian orbital parameters in six binary pulsars, caused by kinematic effects. Improved Shapiro delay measurements have resulted in much improved pulsar mass measurements, particularly for PSRs J0437−4715 and J1909−3744 with M
p = 1.44 ± 0.07 and 1.47 ± 0.03 M⊙, respectively. The improved orbital period-derivative measurement for PSR J0437−4715 results in a derived distance measurement at the 0.16 per cent level of precision, D = 156.79 ± 0.25 pc, one of the most fractionally precise distance measurements of any star to date.
ABSTRACT We present a new model for the distribution of free electrons in the Galaxy, the Magellanic Clouds, and the intergalactic medium (IGM) that can be used to estimate distances to real or ...simulated pulsars and fast radio bursts (FRBs) based on their dispersion measure (DM). The Galactic model has an extended thick disk representing the so-called warm interstellar medium, a thin disk representing the Galactic molecular ring, spiral arms based on a recent fit to Galactic H ii regions, a Galactic Center disk, and seven local features including the Gum Nebula, Galactic Loop I, and the Local Bubble. An offset of the Sun from the Galactic plane and a warp of the outer Galactic disk are included in the model. Parameters of the Galactic model are determined by fitting to 189 pulsars with independently determined distances and DMs. Simple models are used for the Magellanic Clouds and the IGM. Galactic model distances are within the uncertainty range for 86 of the 189 independently determined distances and within 20% of the nearest limit for a further 38 pulsars. We estimate that 95% of predicted Galactic pulsar distances will have a relative error of less than a factor of 0.9. The predictions of YMW16 are compared to those of the TC93 and NE2001 models showing that YMW16 performs significantly better on all measures. Timescales for pulse broadening due to interstellar scattering are estimated for (real or simulated) Galactic and Magellanic Cloud pulsars and FRBs.
Signals from radio pulsars show a wavelength-dependent delay due to dispersion in the interstellar plasma. At a typical observing wavelength, this delay can vary by tens of microseconds on 5-yr ...time-scales, far in excess of signals of interest to pulsar timing arrays, such as that induced by a gravitational wave background. Measurement of these delay variations is not only crucial for the detection of such signals, but also provides an unparalleled measurement of the turbulent interstellar plasma at astronomical unit (au) scales.
In this paper we demonstrate that without consideration of wavelength-independent red noise, 'simple' algorithms to correct for interstellar dispersion can attenuate signals of interest to pulsar timing arrays. We present a robust method for this correction, which we validate through simulations, and apply it to observations from the Parkes Pulsar Timing Array. Correction for dispersion variations comes at a cost of increased band-limited white noise. We discuss scheduling to minimize this additional noise, and factors, such as scintillation, that can exacerbate the problem.
Comparison with scintillation measurements confirms previous results that the spectral exponent of electron density variations in the interstellar medium often appears steeper than expected. We also find a discrete change in dispersion measure of PSR J1603−7202 of ∼2 × 10−3 cm−3 pc for about 250 d. We speculate that this has a similar origin to the 'extreme scattering events' seen in other sources. In addition, we find that four pulsars show a wavelength-dependent annual variation, indicating a persistent gradient of electron density on an au spatial scale, which has not been reported previously.
ABSTRACT
Pulsar timing array projects measure the pulse arrival times of millisecond pulsars for the primary purpose of detecting nanohertz-frequency gravitational waves. The measurements include ...contributions from a number of astrophysical and instrumental processes, which can either be deterministic or stochastic. It is necessary to develop robust statistical and physical models for these noise processes because incorrect models diminish sensitivity and may cause a spurious gravitational wave detection. Here we characterize noise processes for the 26 pulsars in the second data release of the Parkes Pulsar Timing Array using Bayesian inference. In addition to well-studied noise sources found previously in pulsar timing array data sets such as achromatic timing noise and dispersion measure variations, we identify new noise sources including time-correlated chromatic noise that we attribute to variations in pulse scattering. We also identify ‘exponential dip’ events in four pulsars, which we attribute to magnetospheric effects as evidenced by pulse profile shape changes observed for three of the pulsars. This includes an event in PSR J1713+0747, which had previously been attributed to interstellar propagation. We present noise models to be used in searches for gravitational waves. We outline a robust methodology to evaluate the performance of noise models and identify unknown signals in the data. The detection of variations in pulse profiles highlights the need to develop efficient profile domain timing methods.
We present results of an all-sky search in the Parkes Pulsar Timing Array (PPTA) Data Release 1 data set for continuous gravitational waves (GWs) in the frequency range from 5 × 10−9 to 2 × 10−7 Hz. ...Such signals could be produced by individual supermassive binary black hole systems in the early stage of coalescence. We phase up the pulsar timing array data set to form, for each position on the sky, two data streams that correspond to the two GW polarizations and then carry out an optimal search for GW signals on these data streams. Since no statistically significant GWs were detected, we place upper limits on the intrinsic GW strain amplitude h
0 for a range of GW frequencies. For example, at 10−8 Hz our analysis has excluded with 95 per cent confidence the presence of signals with h
0 ≥ 1.7 × 10−14. Our new limits are about a factor of 4 more stringent than those of Yardley et al. based on an earlier PPTA data set and a factor of 2 better than those reported in the recent Arzoumanian et al. paper. We also present PPTA directional sensitivity curves and find that for the most sensitive region on the sky, the current data set is sensitive to GWs from circular supermassive binary black holes with chirp masses of 109 M⊙ out to a luminosity distance of about 100 Mpc. Finally, we set an upper limit of 4 × 10−3 Mpc−3 Gyr−1 at 95 per cent confidence on the coalescence rate of nearby (z ≲ 0.1) supermassive binary black holes in circular orbits with chirp masses of 1010 M⊙.
ABSTRACT
The main goal of pulsar timing array experiments is to detect correlated signals such as nanohertz-frequency gravitational waves. Pulsar timing data collected in dense monitoring campaigns ...can also be used to study the stars themselves, their binary companions, and the intervening ionized interstellar medium. Timing observations are extraordinarily sensitive to changes in path-length between the pulsar and the Earth, enabling precise measurements of the pulsar positions, distances and velocities, and the shapes of their orbits. Here we present a timing analysis of 25 pulsars observed as part of the Parkes Pulsar Timing Array (PPTA) project over time spans of up to 24 yr. The data are from the second data release of the PPTA, which we have extended by including legacy data. We make the first detection of Shapiro delay in four Southern pulsars (PSRs J1017−7156, J1125−6014, J1545−4550, and J1732−5049), and of parallax in six pulsars. The prominent Shapiro delay of PSR J1125−6014 implies a neutron star mass of Mp = 1.5 ± 0.2 M⊙ (68 per cent credibility interval). Measurements of both Shapiro delay and relativistic periastron advance in PSR J1600−3053 yield a large but uncertain pulsar mass of $M_p = 2.06^{+0.44}_{-0.41}$ M⊙ (68 per cent credibility interval). We measure the distance to PSR J1909−3744 to a precision of 10 lyr, indicating that for gravitational wave periods over a decade, the pulsar provides a coherent baseline for pulsar timing array experiments.
Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems would ...modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrained the characteristic amplitude of this background, Ac,yr, to be <1.0 × 10–15 with 95% confidence. This limit excludes predicted ranges for Ac,yr from current models with 91 to 99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments and that higher-cadence and shorter-wavelength observations would be more sensitive to gravitational waves.
High-sensitivity radio-frequency observations of millisecond pulsars usually show stochastic, broad-band, pulse-shape variations intrinsic to the pulsar emission process. These variations induce ...jitter noise in pulsar timing observations; understanding the properties of this noise is of particular importance for the effort to detect gravitational waves with pulsar timing arrays. We assess the short-term profile and timing stability of 22 millisecond pulsars that are part of the Parkes Pulsar Timing Array sample by examining intraobservation arrival time variability and single-pulse phenomenology. In 7 of the 22 pulsars, in the band centred at approximately 1400 MHz, we find that the brightest observations are limited by intrinsic jitter. We find consistent results, either detections or upper limits, for jitter noise in other frequency bands. PSR J1909−3744 shows the lowest levels of jitter noise, which we estimate to contribute ∼10 ns root mean square error to the arrival times for hour-duration observations. Larger levels of jitter noise are found in pulsars with wider pulses and distributions of pulse intensities. The jitter noise in PSR J0437−4715 decorrelates over a bandwidth of ∼2 GHz. We show that the uncertainties associated with timing pulsar models can be improved by including physically motivated jitter uncertainties. Pulse-shape variations will limit the timing precision at future, more sensitive, telescopes; it is imperative to account for this noise when designing instrumentation and timing campaigns for these facilities.
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, ...image, determined at the 1.5% level of precision (image 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 image allows a limit on the variation of Newton's gravitational constant, image yr super(-1). We also constrain any anomalous acceleration along the line of sight to the pulsar to image s super(-1) at 95% confidence, and derive a pulsar mass, image, one of the highest estimates so far obtained.
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.