The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. ...We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the \(2015\) European Pulsar Timing Array data Desvignes et al. (in prep.). Our analysis of the GWB in the \(\sim 2 - 90\) nHz band shows consistency with isotropy, with the strain amplitude in \(l>0\) spherical harmonic multipoles \(\lesssim 40\%\) of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations.
Millisecond pulsars (MSPs) have been studied in detail since their discovery in 1982. The integrated pulse profiles of MSPs appear to be stable, which enables precision monitoring of the pulse times ...of arrival (TOAs). However, for individual pulses the shape and arrival phase can vary dramatically, which is known as pulse jitter. In this paper, we investigate the stability of integrated pulse profiles for 5 MSPs, and estimate the amount of jitter for PSR J0437-4715. We do not detect intrinsic profile shape variation based on integration times from ~10 to ~100 s with the provided instrumental sensitivity. For PSR J0437-4715 we calculate the jitter parameter to be f_J=0.067+-0.002, and demonstrate that the result is not significantly affected by instrumental TOA uncertainties. Jitter noise is also found to be independent of observing frequency and bandwidth around 1.4 GHz on frequency scales of <100 MHz, which supports the idea that pulses within narrow frequency scale are equally jittered. In addition, we point out that pulse jitter would limit TOA calculation for the timing observations with future telescopes like the Square Kilometre Array and the Five hundred metre Aperture Spherical Telescope. A quantitative understanding of pulse profile stability and the contribution of jitter would enable improved TOA calculations, which are essential for the ongoing endeavours in pulsar timing, such as the detection of the stochastic gravitational wave background.
Pulsar timing observations have revealed companions to neutron stars that include other neutron stars, white dwarfs, main-sequence stars, and planets. We demonstrate that the correlated and ...apparently stochastic residual times of arrival from the millisecond pulsar B1937+21 are consistent with the signature of an asteroid belt having a total mass less than approximately 0.05 Earth masses. Unlike the solar system's asteroid belt, the best fit pulsar asteroid belt extends over a wide range of radii, consistent with the absence of any shepherding companions. We suggest that any pulsar that has undergone accretion-driven spin-up and subsequently evaporated its companion may harbor orbiting asteroid mass objects. The resulting timing variations may fundamentally limit the timing precision of some of the other millisecond pulsars. Observational tests of the asteroid belt model include identifying periodicities from individual asteroids, which are difficult; testing for statistical stationarity that become possible when observations are conducted over a longer observing span; and searching for reflected radio emission.
We describe how to implement the spectral kurtosis method of interference removal (zapping) on a digitized signal of averaged power values. Spectral kurtosis is a hypothesis test, analogous to the ...t-test, with a null hypothesis that the amplitudes from which power is formed belong to a `good' distribution -- typically Gaussian with zero mean -- where power values are zapped if the hypothesis is rejected at a specified confidence level. We derive signal-to-noise ratios (SNRs) as a function of amount of zapping for folded radio pulsar observations consisting of a sum of signals from multiple telescopes in independent radio-frequency interference (RFI) environments, comparing four methods to compensate for lost data with coherent (tied-array) and incoherent summation. For coherently summed amplitudes, scaling amplitudes from non-zapped telescopes achieves a higher SNR than replacing zapped amplitudes with artificial noise. For incoherently summed power values, the highest SNR is given by scaling power from non-zapped telescopes to maintain a constant mean. We use spectral kurtosis to clean a tied-array radio pulsar observation by the Large European Array for Pulsars (LEAP): the signal from one telescope is zapped with time and frequency resolutions of 6.25 \(\mu\)s and 0.16 MHz, removing interference along with 0.27 per cent of `good' data, giving an uncertainty of 0.25 \(\mu\)s in pulse time of arrival (TOA) for PSR J1022+1001. We use a single-telescope observation to demonstrate recovery of the pulse profile shape, with 0.6 per cent of data zapped and a reduction from 1.22 to 0.70 \(\mu\)s in TOA uncertainty.
We have conducted radio timing observations of the eclipsing millisecond binary pulsar J2051-0827 with the European Pulsar Timing Array network of telescopes and the Parkes radio telescope, spanning ...over 13 years. The increased data span allows significant measurements of the orbital eccentricity, e = (6.2 {\pm} 1.3) {\times} 10^{-5} and composite proper motion, {\mu}_t = 7.3 {\pm} 0.4 mas/yr. Our timing observations have revealed secular variations of the projected semi-major axis of the pulsar orbit which are much more extreme than those previously published; and of the orbital period of the system. Investigation of the physical mechanisms producing such variations confirm that the variations of the semi-major axis are most probably caused by classical spin-orbit coupling in the binary system, while the variations in orbital period are most likely caused by tidal dissipation leading to changes in the gravitational quadrupole moment of the companion.
Direct detection of low-frequency gravitational waves (\(10^{-9} - 10^{-8}\) Hz) is the main goal of pulsar timing array (PTA) projects. One of the main targets for the PTAs is to measure the ...stochastic background of gravitational waves (GWB) whose characteristic strain is expected to approximately follow a power-law of the form \(h_c(f)=A (f/\hbox{yr}^{-1})^{\alpha}\), where \(f\) is the gravitational-wave frequency. In this paper we use the current data from the European PTA to determine an upper limit on the GWB amplitude \(A\) as a function of the unknown spectral slope \(\alpha\) with a Bayesian algorithm, by modelling the GWB as a random Gaussian process. For the case \(\alpha=-2/3\), which is expected if the GWB is produced by supermassive black-hole binaries, we obtain a 95% confidence upper limit on \(A\) of \(6\times 10^{-15}\), which is 1.8 times lower than the 95% confidence GWB limit obtained by the Parkes PTA in 2006. Our approach to the data analysis incorporates the multi-telescope nature of the European PTA and thus can serve as a useful template for future intercontinental PTA collaborations.
PSR J1518+4904 is one of only 9 known double neutron star systems. These systems are highly valuable for measuring the masses of neutron stars, measuring the effects of gravity, and testing ...gravitational theories. We determine an improved timing solution for a mildly relativistic double neutron star system, combining data from multiple telescopes. We set better constraints on relativistic parameters and the separate masses of the system, and discuss the evolution of PSR J1518+4904 in the context of other double neutron star systems. PSR J1518+4904 has been regularly observed for more than 10 years by the European Pulsar Timing Array (EPTA) network using the Westerbork, Jodrell Bank, Effelsberg and Nancay radio telescopes. The data were analysed using the updated timing software Tempo2. We have improved the timing solution for this double neutron star system. The periastron advance has been refined and a significant detection of proper motion is presented. It is not likely that more post-Keplerian parameters, with which the individual neutron star masses and the inclination angle of the system can be determined separately, can be measured in the near future. Using a combination of the high-quality data sets present in the EPTA collaboration, extended with the original GBT data, we have constrained the masses in the system to m_p<1.17 msun and m_c>1.55 msun (95.4% confidence), and the inclination angle of the orbit to be less than 47 degrees (99%). From this we derive that the pulsar in this system possibly has one of the lowest neutron star masses measured to date. From evolutionary considerations it seems likely that the companion star, despite its high mass, was formed in an electron-capture supernova.
We present results from the high precision timing analysis of the pulsar-white dwarf (WD) binary PSR J1012+5307 using 15 years of multi-telescope data. Observations were performed regularly by the ...European Pulsar Timing Array (EPTA) network, consisting of Effelsberg, Jodrell Bank, Westerbork and Nançay. All the timing parameters have been improved from the previously published values, most by an order of magnitude. In addition, a parallax measurement of \(\pi = 1.2(3)\) mas is obtained for the first time for PSR J1012+5307, being consistent with the optical estimation from the WD companion. Combining improved 3D velocity information and models for the Galactic potential the complete evolutionary Galactic path of the system is obtained. A new intrinsic eccentricity upper limit of \(e<8.4\times 10^{-7}\) is acquired, one of the smallest calculated for a binary system and a measurement of the variation of the projected semi-major axis also constrains the system's orbital orientation for the first time. It is shown that PSR J1012+5307 is an ideal laboratory for testing alternative theories of gravity. The measurement of the change of the orbital period of the system of \(\dot{P}_{b} = 5(1)\times 10^{-14}\) is used to set an upper limit on the dipole gravitational wave emission that is valid for a wide class of alternative theories of gravity. Moreover, it is shown that in combination with other binary pulsars PSR J1012+5307 is an ideal system to provide self-consistent, generic limits, based only on millisecond pulsar data, for the dipole radiation and the variation of the gravitational constant \(\dot{G}\).
High-precision pulsar timing relies on a solar-system ephemeris in order to convert times of arrival (TOAs) of pulses measured at an observatory to the solar system barycenter. Any error in the ...conversion to the barycentric TOAs leads to a systematic variation in the observed timing residuals; specifically, an incorrect planetary mass leads to a predominantly sinusoidal variation having a period and phase associated with the planet's orbital motion about the Sun. By using an array of pulsars (PSRs J0437-4715, J1744-1134, J1857+0943, J1909-3744), the masses of the planetary systems from Mercury to Saturn have been determined. These masses are consistent with the best-known masses determined by spacecraft observations, with the mass of the Jovian system, 9.547921(2)E-4 Msun, being significantly more accurate than the mass determined from the Pioneer and Voyager spacecraft, and consistent with but less accurate than the value from the Galileo spacecraft. While spacecraft are likely to produce the most accurate measurements for individual solar system bodies, the pulsar technique is sensitive to planetary system masses and has the potential to provide the most accurate values of these masses for some planets.