European Pulsar Timing Array limits on an isotropic stochastic gravitational-wave background

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 Ω^(relic)_(gw)(f)h^2 <1.2×10^(−9), a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array

Fisher Formalism For Anisotropic Gravitational-Wave Background Searches With Pulsar Timing Arrays

Pulsar timing arrays (PTAs) are currently the only experiments directly sensitive to gravitational waves with decade-long periods. Within the next five to ten years, PTAs are expected to detect the stochastic gravitational-wave background (SGWB) collectively sourced by inspiraling supermassive black hole binaries. It is expected that this background is mostly isotropic, and current searches focus on the monopole part of the SGWB. Looking ahead, anisotropies in the SGWB may provide a trove of additional information on both known and unknown astrophysical and cosmological sources. In this paper, we build a simple yet realistic Fisher formalism for anisotropic SGWB searches with PTAs. Our formalism is able to accommodate realistic properties of PTAs and allows simple and accurate forecasts. We illustrate our approach with an idealized PTA consisting of identical, isotropically distributed pulsars. In a companion paper, we apply our formalism to current PTAs and show that it can be a powerful tool to guide and optimize real data analysis

Analysis of the first IPTA Mock Data Challenge by the EPTA timing data analysis working group

This is a summary of the methods we used to analyse the first IPTA Mock Data Challenge (MDC), and the obtained results. We have used a Bayesian analysis in the time domain, accelerated using the recently developed ABC-method which consists of a form of lossy linear data compression. The TOAs were first processed with Tempo2, where the design matrix was extracted for use in a subsequent Bayesian analysis. We used different noise models to analyse the datasets: no red noise, red noise the same for all pulsars, and individual red noise per pulsar. We sampled from the likelihood with four different samplers: "emcee", "t-walk", "Metropolis-Hastings", and "pyMultiNest". All but emcee agreed on the final result, with emcee failing due to artefacts of the high-dimensionality of the problem. An interesting issue we ran into was that the prior of all the 36 (red) noise amplitudes strongly affects the results. A flat prior in the noise amplitude biases the inferred GWB amplitude, whereas a flat prior in log-amplitude seems to work well. This issue is only apparent when using a noise model with individually modelled red noise for all pulsars. Our results for the blind challenges are in good agreement with the injected values. For the GWB amplitudes we found h_c = 1.03 +/- 0.11 [10^{-14}], h_c = 5.70 +/- 0.35 [10^{-14}], and h_c = 6.91 +/- 1.72 [10^{-15}], and for the GWB spectral index we found gamma = 4.28 +/- 0.20, gamma = 4.35 +/- 0.09, and gamma = 3.75 +/- 0.40. We note that for closed challenge 3 there was quite some covariance between the signal and the red noise: if we constrain the GWB spectral index to the usual choice of gamma = 13/3, we obtain the estimates: h_c = 10.0 +/- 0.64 [10^{-15}], h_c = 56.3 +/- 2.42 [10^{-15}], and h_c = 4.83 +/- 0.50 [10^{-15}], with one-sided 2 sigma upper-limits of: h_c <= 10.98 [10^{-15}], h_c <= 60.29 [10^{-15}], and h_c <= 5.65 [10^{-15}]

Insights Into Searches For Anisotropies In The Nanohertz Gravitational-Wave Background

Within the next several years pulsar timing arrays (PTAs) are positioned to detect the stochastic gravitational-wave background (GWB) likely produced by the collection of inspiralling supermassive black holes binaries, and potentially constrain some exotic physics. So far most of the pulsar timing data analysis has focused on the monopole of the GWB, assuming it is perfectly isotropic. The natural next step is to search for anisotropies in the GWB. In this paper, we use the recently developed PTA Fisher matrix to gain insights into optimal search strategies for GWB anisotropies. For concreteness, we apply our results to the European Pulsar Timing Array (EPTA) data, using realistic noise characteristics of its pulsars. We project the detectability of a GWB whose angular dependence is assumed to be a linear combination of predetermined maps, such as spherical harmonics or coarse pixels. We find that the GWB monopole is always statistically correlated with these maps, implying a loss of sensitivity to the monopole when searching simultaneously for anisotropies. We then derive the angular distributions of the GWB intensity to which a PTA is most sensitive, and illustrate how one may use these “principal maps” to approximately reconstruct the angular dependence of the GWB. Since the principal maps are neither perfectly anisotropic nor uncorrelated with the monopole, we also develop a frequentist criterion to specifically search for anisotropies in the GWB without any prior knowledge about their angular distribution. Lastly, we show how to recover existing EPTA results with our Fisher formalism, and clarify their meaning. The tools presented here will be valuable in guiding and optimizing the computationally demanding analyses of pulsar timing data

A data analysis library for gravitational wave detection

One of the main goals of Pulsar Timing Arrays (PTAs) is the direct detection of gravitational waves (GWs). A first detection will be a major leap for astronomy and substantial effort is currently going into timing as many pulsars as possible, with the highest possible accuracy. As part of the individual PTA projects, several groups are developing data analysis methods for the final stage of a gravitational-waves search pipeline: the analysis of the timing residuals. Here we report the progress of on-going work to develop, within a Bayesian framework, a comprehensive and user friendly analysis library to search for gravitational waves in PTA data

Observing the dynamics of super-massive black hole binaries with Pulsar Timing Arrays

Pulsar Timing Arrays are a prime tool to study unexplored astrophysical regimes with gravitational waves. Here we show that the detection of gravitational radiation from individually resolvable super-massive black hole binary systems can yield direct information about the masses and spins of the black holes, provided that the gravitational-wave induced timing fluctuations both at the pulsar and at the Earth are detected. This in turn provides a map of the non-linear dynamics of the gravitational field and a new avenue to tackle open problems in astrophysics connected to the formation and evolution of super-massive black holes. We discuss the potential, the challenges and the limitations of these observations.Comment: 5 pages, 1 figur

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 ∼60 per cent compared to a model that includes DM variations and spin-noise only

The NANOGrav Nine-year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background

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_(gw) < 1.5 x 10^(-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, Ω_(gw)(f)h^2 < 4.2 x 10^(-10). Our limit on the cosmic string GWB, Ω_(gw)(f)h^2 < 2.2 x 10^(-10), translates to a conservative limit on the cosmic string tension with Gµ < 3.3 x 10^(-8), a factor of four better than the joint Planck and high-l cosmic microwave background data from other experiments