The noise properties of 42 millisecond pulsars from the European Pulsar Timing Array and their impact on gravitational wave searches

The sensitivity of Pulsar Timing Arrays to gravitational waves depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispersion and scattering effects, analysis errors and instrumental instabilities. We present the results from a noise analysis for 42 millisecond pulsars (MSPs) observed with the European Pulsar Timing Array. For characterising the low-frequency, stochastic and achromatic noise component, or "timing noise", we employ two methods, based on Bayesian and frequentist statistics. For 25 MSPs, we achieve statistically significant measurements of their timing noise parameters and find that the two methods give consistent results. For the remaining 17 MSPs, we place upper limits on the timing noise amplitude at the 95% confidence level. We additionally place an upper limit on the contribution to the pulsar noise budget from errors in the reference terrestrial time standards (below 1%), and we find evidence for a noise component which is present only in the data of one of the four used telescopes. Finally, we estimate that the timing noise of individual pulsars reduces the sensitivity of this data set to an isotropic, stochastic GW background by a factor of >9.1 and by a factor of >2.3 for continuous GWs from resolvable, inspiralling supermassive black-hole binaries with circular orbits.Comment: Accepted for publication by the Monthly Notices of the Royal Astronomical Societ

A millisecond pulsar in an extremely wide binary system

International audienceWe report on 22 yrs of radio timing observations of the millisecond pulsar J1024−0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869−0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, [M/H] = −1.0, T eff = 4050 ± 50 K) and that its position, proper motion and distance are consistent with those of PSR J1024−0719. We conclude that PSR J1024−0719 and 2MASS J10243869−0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives , which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide (P b > 200 yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of 384 ± 45 km s −1 with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024−0719 in light of its inclusion in pulsar timing arrays

Validation of heliospheric modeling algorithms through pulsar observations II: simulations with EUHFORIA

In space weather studies and forecasting we employ magnetohydrodynamic (MHD) simulations which can provide rather accurate reconstruction of the solar wind dynamics and its evolution. However, all MHD simulations are restricted by the input data and the modelled solar wind characteristics need to be validated with different types of observations. That is very difficult, in particular for the solar wind characteristics close to the Sun, since the majority of in-situ observations are taken in the vicinity of the Earth. This is why all alternative methods for estimation of solar wind plasma characteristics are very important. In this study we utilise low radio frequency observations of pulsars to probe the total electron content along the line of sight. For the first time, we compare density estimates from pulsars with predictions from the 3D MHD modelling code; the EUropean Heliospheric FORecasting Information Asset (EUHFORIA). We find a very good correlation for the solar wind density along a given line of sight obtained by EUHFORIA and pulsar observations. We also demonstrate that the pulsar observations can be very useful not only for the model validation but also for understanding its limitations.Comment: Published in Journal of Advances in Space Researc

Discovery and modelling of broad-scale plasma lensing in black-widow pulsar J2051−-0827

We report on an unusually bright observation of PSR J2051$-$0827 recorded during a regular monitoring campaign of black-widow pulsar systems with the Effelsberg 100-m telescope. Through fortunate coincidence, a particularly bright scintillation maximum is simultaneous with the eclipse by the companion, enabling precise measurements of variations in the flux density, dispersion measure (DM), and scattering strength throughout the eclipse. The flux density is highly variable throughout the eclipse, with a peak 1.7 times the average away from the eclipse, and yet does not significantly decrease on average. We recover the flux density variations from the measured DM variations using geometric optics, with a relative velocity as the only free parameter. We measure an effective velocity of (470 $\pm$ 10) km/s, consistent with the relative orbital motion of the companion, suggesting that the outflow velocity of the lensing material is low, or is directly along the line of sight. The 2 per cent uncertainty on the effective velocity is a formal error; systematics related to our current model are likely to dominate, and we detail several extensions to the model to be considered in a full treatment of lensing. This is a demonstration of the causal link between DM and lensing; the flux density variations can be predicted directly through the derivatives of DM. Going forward, this approach can be applied to investigate the dynamics of other eclipsing systems, and to investigate the physical nature of scintillation and lensing in the ionized interstellar medium.Comment: 12 pages, 8 figures, typos corrected, references update

High-precision timing of 42 millisecond pulsars with the European Pulsar Timing Array

International audienceWe report on timing, flux density, and polarimetric observations of the transient magnetar and 5.54 s radio pulsar XTE J1810-197 using the GBT, Nancay, and Parkes radio telescopes beginning in early 2006, until its sudden disappearance as a radio source in late 2008. Repeated observations through 2016 have not detected radio pulsations again. The torque on the neutron star, as inferred from its rotation frequency derivative f-dot, decreased in an unsteady manner by a factor of 3 in the first year of radio monitoring. In contrast, during its final year as a detectable radio source, the torque decreased steadily by only 9%. The period-averaged flux density, after decreasing by a factor of 20 during the first 10 months of radio monitoring, remained steady in the next 22 months, at an average of 0.7+/-0.3 mJy at 1.4 GHz, while still showing day-to-day fluctuations by factors of a few. There is evidence that during this last phase of radio activity the magnetar had a steep radio spectrum, in contrast to earlier behavior. There was no secular decrease that presaged its radio demise. During this time the pulse profile continued to display large variations, and polarimetry indicates that the magnetic geometry remained consistent with that of earlier times. We supplement these results with X-ray timing of the pulsar from its outburst in 2003 up to 2014. For the first 4 years, XTE J1810-197 experienced non-monotonic excursions in f-dot by at least a factor of 8. But since 2007, its f-dot has remained relatively stable near its minimum observed value. The only apparent event in the X-ray record that is possibly contemporaneous with the radio shut-down is a decrease of ~20% in the hot-spot flux in 2008-2009, to a stable, minimum value. However, the permanence of the high-amplitude, thermal X-ray pulse, even after the radio demise, implies continuing magnetar activity

From spin noise to systematics:stochastic processes in the first International Pulsar Timing Array data release

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 pulsar timing array (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 ionized 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 pulsar's 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 similar to 60 per cent compared to a model that includes DM variations and spin-noise only

The International Pulsar Timing Array: First data release

International audienceThe highly stable spin of neutron stars can be exploited for a variety of (astro)physical investigations. In particular, arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such 'pulsar timing arrays' (PTAs) have been set up around the world over the past decades and collectively form the 'International' PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available on-line) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limit

Validation of heliospheric modeling algorithms through pulsar observations I: Interplanetary scintillation-based tomography

Solar-wind 3-D reconstruction tomography based on interplanetary scintillation (IPS) studies provides fundamental information for space-weather forecasting models, and gives the possibility to determine heliospheric column densities. Here we compare the time series of Solar-wind column densities derived from long-term observations of pulsars, and the Solar-wind reconstruction provided by the UCSD IPS tomography. This work represents a completely independent comparison and validation of these techniques to provide this measurement, and it strengthens confidence in the use of both in space-weather analyses applications.Comment: Published in Journal of Advances in Space Researc

The International Pulsar Timing Array:First data release

The highly stable spin of neutron stars can be exploited for a variety of (astro)physical investigations. In particular, arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such 'pulsar timing arrays' (PTAs) have been set up around the world over the past decades and collectively form the 'International' PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of similar to 2 and provides a 2 sigma limit on the dimensionless amplitude of a stochastic gravitational-wave background of 1.7 x 10(-15) at a frequency of 1 yr(-1). This is 1.7 times less constraining than the limit placed by Shannon et al., due mostly to the more recent, high-quality data they used

The Scintillating Tail of Comet C/2020 F3 (Neowise)

Context. The occultation of a radio source by the plasma tail of a comet can be used to probe structure and dynamics in the tail. Such occultations are rare, and the occurrence of scintillation, due to small-scale density variations in the tail, remains somewhat controversial. Aims. A detailed observation taken with the Low-Frequency Array (LOFAR) of a serendipitous occultation of the compact radio source 3C196 by the plasma tail of comet C/2020 F3 (Neowise) is presented. 3C196 tracked almost perpendicularly behind the tail, providing a unique profile cut only a short distance downstream from the cometary nucleus itself. Methods. Interplanetary scintillation (IPS) is observed as the rapid variation of the intensity received of a compact radio source due to density variations in the solar wind. IPS in the signal received from 3C196 was observed for five hours, covering the full transit behind the plasma tail of comet C/2020 F3 (Neowise) on 16 July 2020, and allowing an assessment of the solar wind in which the comet and its tail are embedded. Results. The results reveal a sudden and strong enhancement in scintillation which is unequivocally attributable to the plasma tail. The strongest scintillation is associated with the tail boundaries, weaker scintillation is seen within the tail, and previously-unreported periodic variations in scintillation are noted, possibly associated with individual filaments of plasma. Furthermore, contributions from the solar wind and comet tail are separated to measure a sharp decrease in the velocity of material within the tail, suggesting a steep velocity shear resulting in strong turbulence along the tail boundaryComment: Accepted for publication in Astronomy and Astrophysics, 8 pages, 9 figure