Using a statistically rigorous analysis method, we place limits on the existence of an isotropic stochastic gravitational wave background using pulsar timing observations. We consider backgrounds whose characteristic strain spectra may be described as a power-law dependence with frequency. Such backgrounds include an astrophysical background produced by coalescing supermassive black-hole binary systems and cosmological backgrounds due to relic gravitational waves and cosmic strings. Using the best available data, we obtain an upper limit on the energy density per unit logarithmic frequency interval of SMBH g 1/(8 yr) ½ h 2 1:9 ; 10 À8 for an astrophysical background that is 5 times more stringent than the earlier limit of 1:1 ; 10 À7 found by Kaspi and colleagues. We also provide limits on a background due to relic gravitational waves and cosmic strings of relic g 1/(8 yr) ½ h 2 2:0 ; 10 À8 and cs g 1/(8 yr)½ h 2 1:9 ; 10 À8 , respectively. All of the quoted upper limits correspond to a 0.1% false alarm rate together with a 95% detection rate. We discuss the physical implications of these results and highlight the future possibilities of the Parkes Pulsar Timing Array project. We find that our current results can (1) constrain the merger rate of supermassive binary black hole systems at high redshift, (2) rule out some relationships between the black hole mass and the galactic halo mass, (3) constrain the rate of expansion in the inflationary era, and (4) provide an upper bound on the dimensionless tension of a cosmic string background.

General rightsIt is not permitted to download or to forward/distribute the text or part of it without the consent of the author(s) and/or copyright holder(s), other than for strictly personal, individual use, unless the work is under an open content license (like Creative Commons). Disclaimer/Complaints regulationsIf you believe that digital publication of certain material infringes any of your rights or (privacy) interests, please let the Library know, stating your reasons. In case of a legitimate complaint, the Library will make the material inaccessible and/or remove it from the website. Please Ask the Library: http://uba.uva.nl/en/contact, or a letter to: Library of the University of Amsterdam, Secretariat, Singel 425, 1012 WP Amsterdam, The Netherlands. You will be contacted as soon as possible. ABSTRACT Recent work has exploited pulsar survey data to identify temporally isolated, millisecond-duration radio bursts with large dispersion measures (DMs). These bursts have been interpreted as arising from a population of extragalactic sources, in which case they would provide unprecedented opportunities for probing the intergalactic medium; they may also be linked to new source classes. Until now, however, all so-called fast radio bursts (FRBs) have been detected with the Parkes radio telescope and its 13-beam receiver, casting some concern about the astrophysical nature of these signals. Here we present FRB 121102, the first FRB discovery from a geographic location other than Parkes. FRB 121102 was found in the Galactic anti-center region in the 1.4 GHz Pulsar Arecibo L-band Feed Array (ALFA) survey with the Arecibo Observatory with a DM = 557.4 ± 2.0 pc cm −3 , pulse width of 3.0 ± 0.5 ms, and no evidence of interstellar scattering. The observed delay of the signal arrival time with frequency agrees precisely with the expectation of dispersion through an ionized medium. Despite its low Galactic latitude (b = −0.• 2), the burst has three times the maximum Galactic DM expected along this particular line of sight, suggesting an extragalactic origin. A peculiar aspect of the signal is an inverted spectrum; we interpret this as a consequence of being detected in a sidelobe of the ALFA receiver. FRB 121102's brightness, duration, and the inferred event rate are all consistent with the properties of the previously detected Parkes bursts.

We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, and six high-timing-precision pulsars were observed weekly. All were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and (for binary pulsars) orbital parameters; time-variable dispersion delays; and parameters that 1 quantify pulse-profile evolution with frequency. The timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of significant orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. A companion paper uses these data to constrain the strength of the gravitational-wave background.

Here we present an analysis of high-precision pulsar timing data taken as part of the North American Nanohertz Observatory for Gravitational waves (NANOGrav) project. We have observed 17 pulsars for a span of roughly five years using the Green Bank and Arecibo radio telescopes. We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Sub-microsecond timing residuals are obtained in nearly all cases, and the best root-mean-square timing residuals in this set are ∼30-50 ns. We present methods for analyzing post-fit timing residuals for the presence of a gravitational wave signal with a specified spectral shape. These properly and optimally take into account the timing fluctuation power removed by the model fit, and can be applied to either data from a single pulsar, or to a set of pulsars to detect a correlated signal. We apply these methods to our dataset to set an upper limit on the strength of the nHz-frequency stochastic supermassive black hole gravitational wave background of h c (1 yr −1 ) < 7 × 10 −15 (95%). This result is dominated by the timing of the two best pulsars in the set, PSRs J1713+0747 and J1909−3744.

The direct detection of gravitational waves is a major goal of current astrophysics. We provide details of a new method for detecting a stochastic background of gravitational waves using pulsar timing data. Our results show that regular timing observations of 40 pulsars each with a timing accuracy of 100 ns will be able to make a direct detection of the predicted stochastic background from coalescing black holes within five years. With an improved pre-whitening algorithm, or if the background is at the upper end of the predicted range, a significant detection should be possible with only 20 pulsars.

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. Well-tested Bayesian techniques are used to set upper limits on the dimensionless strain amplitude (at a frequency of 1 yr −1 ) for a GWB from supermassive black hole binaries of <´-A 1.5 10 gw 15 . We also parameterize the GWB spectrum with a broken power-law model by placing priors on the strain amplitude derived from simulations of Sesana and McWilliams et al. Using Bayesian model selection we find that the data favor a broken power law to a pure power law with odds ratios of 2.2 and 22 to one for the Sesana and McWilliams prior models, respectively. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. Returning to a power-law model, we place stringent limits on the energy density of relic GWs, W <´-f h 4.2 10 gw 2 1 0 ( ) . Our limit on the cosmic string GWB, W <´-f h 2.2 10 gw 2 1 0 ( ) , translates to a conservative limit on the cosmic string tension with m <´-G 3.3 10 8 , a factor of four better than the joint Planck and high-lcosmic microwave background data from other experiments.

The analysis of high-precision timing observations of an array of ∼20 millisecond pulsars (a so-called 'timing array') may ultimately result in the detection of a stochastic gravitationalwave background. The feasibility of such a detection and the required duration of this type of experiment are determined by the achievable rms of the timing residuals and the timing stability of the pulsars involved. We present results of the first long-term, high-precision timing campaign on a large sample of millisecond pulsars used in gravitational-wave detection projects. We show that the timing residuals of most pulsars in our sample do not contain significant low-frequency noise that could limit the use of these pulsars for decade-long gravitational-wave detection efforts. For our most precisely timed pulsars, intrinsic instabilities of the pulsars or the observing system are shown to contribute to timing irregularities on a 5-year time-scale below the 100 ns level. Based on those results, realistic sensitivity curves for planned and ongoing timing array efforts are determined. We conclude that prospects for detection of a gravitational-wave background through pulsar timing array efforts within 5 years to a decade are good.

A 'pulsar timing array' (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of 'global' phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 ms pulsars is being observed at three radio-frequency bands, 50 cm (ß700 MHz), 20 cm (ß1400 MHz), and 10 cm (ß3100 MHz), with observations at intervals of two to three weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters, and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For 10 of the 20 pulsars, rms timing residuals are less than 1 μs for the best band after fitting for pulse frequency and its first time derivative. Significant 'red' timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array and a PTA based on the Square Kilometre Array. We also present an 'extended PPTA' data set that combines PPTA data with earlier Parkes timing data for these pulsars.