Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of the Parkes Pulsar Timing Array (PPTA), which spans 18 yr. Using current Bayesian inference techniques we recover and characterize a common-spectrum noise process. Represented as a strain spectrum h c = A ( f / 1 yr − 1 ) α , we measure A = 3.1 − 0.9 + 1.3 × 10 − 15 and α = −0.45 ± 0.20, respectively (median and 68% credible interval). For a spectral index of α = −2/3, corresponding to an isotropic background of GWs radiated by inspiraling supermassive black hole binaries, we recover an amplitude of A = 2.04 − 0.22 + 0.25 × 10 − 15 . However, we demonstrate that the apparent signal strength is time-dependent, as the first half of our data set can be used to place an upper limit on A that is in tension with the inferred common-spectrum amplitude using the complete data set. We search for spatial correlations in the observations by hierarchically analyzing individual pulsar pairs, which also allows for significance validation through randomizing pulsar positions on the sky. For a process with α = −2/3, we measure spatial correlations consistent with a GWB, with an estimated false-alarm probability of p ≲ 0.02 (approx. 2σ). The long timing baselines of the PPTA and the access to southern pulsars will continue to play an important role in the International Pulsar Timing Array.

Bayesian parameter estimation is fast becoming the language of gravitational-wave astronomy. It is the method by which gravitational-wave data is used to infer the sources’ astrophysical properties. We introduce a user-friendly Bayesian inference library for gravitational-wave astronomy, Bilby. This Python code provides expert-level parameter estimation infrastructure with straightforward syntax and tools that facilitate use by beginners. It allows users to perform accurate and reliable gravitational-wave parameter estimation on both real, freely available data from LIGO/Virgo and simulated data. We provide a suite of examples for the analysis of compact binary mergers and other types of signal models, including supernovae and the remnants of binary neutron star mergers. These examples illustrate how to change the signal model, implement new likelihood functions, and add new detectors. Bilby has additional functionality to do population studies using hierarchical Bayesian modeling. We provide an example in which we infer the shape of the black hole mass distribution from an ensemble of observations of binary black hole mergers.

Astronomy & Astrophysics (A&A) is an international journal which publishes papers on all aspects of astronomy and astrophysics

We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power-law strain spectrum of the form $h_c = A(f/1\, \mathrm{yr}^{-1})^{\alpha }$, we found strong evidence for a spectrally similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times 10^{-15}$ and spectral index α = −0.5 ± 0.5, where the uncertainties represent 95 per cent credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of α = −2/3, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times 10^{-15}$. None the less, no significant evidence of the Hellings–Downs correlations that would indicate a gravitational-wave origin was found. We also analysed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability.

Pulsar timing arrays offer a probe of the low-frequency gravitational wave spectrum (1–100 nHz), which is intimately connected to a number of markers that can uniquely trace the formation and evolution of the Universe. We present the dataset and the results of the timing analysis from the second data release of the European Pulsar Timing Array (EPTA). The dataset contains high-precision pulsar timing data from 25 millisecond pulsars collected with the five largest radio telescopes in Europe, as well as the Large European Array for Pulsars. The dataset forms the foundation for the search for gravitational waves by the EPTA, presented in associated papers. We describe the dataset and present the results of the frequentist and Bayesian pulsar timing analysis for individual millisecond pulsars that have been observed over the last ~25 yr. We discuss the improvements to the individual pulsar parameter estimates, as well as new measurements of the physical properties of these pulsars and their companions. This data release extends the dataset from EPTA Data Release 1 up to the beginning of 2021, with individual pulsar datasets with timespans ranging from 14 to 25 yr. These lead to improved constraints on annual parallaxes, secular variation of the orbital period, and Shapiro delay for a number of sources. Based on these results, we derived astrophysical parameters that include distances, transverse velocities, binary pulsar masses, and annual orbital parallaxes.

A nanohertz-frequency stochastic gravitational-wave background can potentially be detected through the precise timing of an array of millisecond pulsars. This background produces low-frequency noise in the pulse arrival times that would have a characteristic spectrum common to all pulsars and a well-defined spatial correlation. Recently the North American Nanohertz Observatory for Gravitational Waves collaboration (NANOGrav) found evidence for the common-spectrum component in their 12.5 yr data set. Here we report on a search for the background using the second data release of the Parkes Pulsar Timing Array. If we are forced to choose between the two NANOGrav models—one with a common-spectrum process and one without—we find strong support for the common-spectrum process. However, in this paper, we consider the possibility that the analysis suffers from model misspecification. In particular, we present simulated data sets that contain noise with distinctive spectra but show strong evidence for a common-spectrum process under the standard assumptions. The Parkes data show no significant evidence for, or against, the spatially correlated Hellings–Downs signature of the gravitational-wave background. Assuming we did observe the process underlying the spatially uncorrelated component of the background, we infer its amplitude to be A = 2.2 − 0.3 + 0.4 × 10 − 15 in units of gravitational-wave strain at a frequency of 1 yr⁻¹. Extensions and combinations of existing and new data sets will improve the prospects of identifying spatial correlations that are necessary to claim a detection of the gravitational-wave background.