We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar data set spanning 18 yr of observations from the 2015 European Pulsar Timing Array data ...release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar system ephemeris errors, obtaining a robust 95 per cent upper limit on the dimensionless strain amplitude A of the background of A < 3.0 × 10−15 at a reference frequency of 1 yr−1 and a spectral index of 13/3, corresponding to a background from inspiralling supermassive black hole binaries, constraining the GW energy density to Ωgw(f)h
2 < 1.1 × 10−9 at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of ∼5 × 10−9 Hz. Finally, we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95 per cent upper limits on the string tension, Gμ/c
2, characterizing a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu–Goto field theory cosmic string network, we set a limit Gμ/c
2 < 1.3 × 10−7, identical to that set by the Planck Collaboration, when combining Planck and high-ℓ cosmic microwave background data from other experiments. For a stochastic relic background, we set a limit of
$\Omega ^\mathrm{relic}_\mathrm{gw}(f)h^2<1.2 \times 10^{-9}$
, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.
We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines ...range from 7–18 yr. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with temponest
yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semimajor axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz–Kelker bias) than the M2 and M3 models by Schnitzeler. However, we measure an average uncertainty of 80 per cent (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binary MSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600−3053 and J1918−0642, implying pulsar and companion masses
$m_{\rm p}=1.22_{-0.35}^{+0.5}\ {\rm M}_{{\odot }}$
,
$m_{\rm c} = 0.21_{-0.04}^{+0.06}\ {\rm M}_{{\odot } }$
and
$m_{\rm p}=1.25_{-0.4}^{+0.6}\ {\rm M}_{{\odot }}$
,
$m_{\rm c} = 0.23_{-0.05}^{+0.07}\ {\rm M}_{{\odot } }$
, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012+5307 and J1909−3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600−3053 and J1909−3744.
Direct detection of low-frequency gravitational waves (GWs,
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
, where f is the GW 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 α with a Bayesian algorithm, by modelling the GWB as a random Gaussian process. For the case α=−2/3, which is expected if the GWB is produced by supermassive black hole binaries, we obtain a 95 per cent confidence upper limit on A of 6 × 10−15, which is 1.8 times lower than the 95 per cent confidence GWB limit obtained by the Parkes PTA in 2006. Our approach to the data analysis incorporates the multitelescope nature of the European PTA and thus can serve as a useful template for future intercontinental PTA collaborations.
We present results from the high-precision timing analysis of the pulsar-white dwarf (WD) binary PSR J1012+5307 using 15 years of multitelescope 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 π= 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 × 10−7 is acquired, one of the smallest calculated for a binary system and a measurement of the variation of the projected semimajor 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 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 .
We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest European Pulsar Timing Array ...(EPTA) data set, which consists of ultraprecise timing data on 41-ms pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95 per cent upper limit on the sky-averaged strain amplitude lies in the range ... This limit varies by a factor of five, depending on the assumed source position and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit - obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is ... These limits, the most stringent to date at f < 10...nHz, exclude the presence of sub-centiparsec binaries with chirp mass ... out to a distance of about 25 Mpc, and with ... out to a distance of about 1Gpc (...). We show that state-of-the-art SMBHB population models predict <1 per cent probability of detecting a CGW with the current EPTA data set, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years. (ProQuest: ... denotes formulae/symbols omitted.)
The sensitivity of Pulsar Timing Arrays to gravitational waves (GWs) 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 characterizing 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 per cent 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 per cent), 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.
LEAP: the Large European Array for Pulsars Bassa, C. G; Janssen, G. H; Karuppusamy, R ...
Monthly Notices of the Royal Astronomical Society,
02/2016, Letnik:
456, Številka:
2
Journal Article
Recenzirano
Odprti dostop
The Large European Array for Pulsars (LEAP) is an experiment that harvests the collective power of Europe's largest radio telescopes in order to increase the sensitivity of high-precision pulsar ...timing. As part of the ongoing effort of the European Pulsar Timing Array, LEAP aims to go beyond the sensitivity threshold needed to deliver the first direct detection of gravitational waves. The five telescopes presently included in LEAP are the Effelsberg Telescope, the Lovell Telescope at Jodrell Bank, the Nançay Radio Telescope, the Sardinia Radio Telescope and the Westerbork Synthesis Radio Telescope. Dual polarization, Nyquist-sampled time series of the incoming radio waves are recorded and processed offline to form the coherent sum, resulting in a tied-array telescope with an effective aperture equivalent to a 195-m diameter circular dish. All observations are performed using a bandwidth of 128 MHz centred at a frequency of 1396 MHz. In this paper, we present the design of the LEAP experiment, the instrumentation, the storage and transfer of data and the processing hardware and software. In particular, we present the software pipeline that was designed to process the Nyquist-sampled time series, measure the phase and time delays between each individual telescope and a reference telescope and apply these delays to form the tied-array coherent addition. The pipeline includes polarization calibration and interference mitigation. We also present the first results from LEAP and demonstrate the resulting increase in sensitivity, which leads to an improvement in the pulse arrival times.
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 <, ~0.05 M sub(+ in circle). 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, which becomes possible when observations are conducted over a longer observing span; and searching for reflected radio emission.
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; this is known as pulse jitter. In this paper, we investigate the stability of integrated pulse profiles for five 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-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 a narrow frequency scale are equally jittered. In addition, we point out that pulse jitter would limit TOA calculation for 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 ongoing endeavours in pulsar timing such as detection of the stochastic gravitational wave background.
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 suboptimal. ...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. Our analysis of the GWB in the ~2-90 nHz band shows consistency with isotropy, with the strain amplitude in l>0 spherical harmonic multipoles ≲40% of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations.