The International Pulsar Timing Array second data release: Search for an isotropic gravitational wave background

The subtle influence of gravitational waves on the apparent positioning of celestial bodies offers novel observational windows [1,2,3,4]. We calculate the expected astrometric signal induced by an isotropic Stochastic Gravitational Wave Background (SGWB) in the short distance limit. Our focus is on the resultant proper motion of Solar System objects, a signal on the same time scales addressed by Pulsar Timing Arrays (PTA). We derive the corresponding astrometric deflection patterns, finding that they manifest as distinctive dipole and quadrupole correlations or, in some cases, may not be present. Our analysis encompasses both Einsteinian and non-Einsteinian polarisations. We estimate the upper limits for the amplitude of SGWBs that could be obtained by tracking the proper motions of large numbers of solar system objects such as asteroids. We find that for SGWBs with negative spectral indices, such as that generated by Super Massive Black Hole Binaries (SMBHB), the constraints from these observations could rival those from PTAs. With the Gaia satellite and the Vera C. Rubin Observatory poised to track an extensive sample of asteroids — ranging from 𝒪(105) to 𝒪(106), we highlight the significant future potential for similar surveys to contribute to our understanding of the SGWB.

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 longterm sensitivity of the pulsar to a stochastic gravitational-wave (GW) background.

The highly stable spin of neutron stars can be exploited for a variety of (astro)physical investigations. In particular, arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such 'pulsar timing arrays' (PTAs) have been set up around the world over the past decades and collectively form the 'International' PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits

We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset 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$\%$ upper limit on the dimensionless strain amplitude $A$ of the background of $A<3.0\times 10^{-15}$ at a reference frequency of $1\mathrm{yr^{-1}}$ and a spectral index of $13/3$, corresponding to a background from inspiralling super-massive black hole binaries, constraining the GW energy density to $\Omega_\mathrm{gw}(f)h^2 < 1.1\times10^{-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 $\sim 5\times10^{-9}$~Hz. Finally we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95$\%$ upper limits on the string tension, $G\mu/c^2$, characterising 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\mu/c^2<1.3\times10^{-7}$, identical to that set by the {\it Planck} Collaboration, when combining {\it Planck} and high-$\ell$ 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 \times10^{-9}$, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.

We search for isotropic stochastic gravitational-wave background (SGWB) in the International Pulsar Timing Array second data release. By modeling the SGWB as a power-law, we find very strong Bayesian evidence for a common-spectrum process, and further this process has scalar transverse (ST) correlations allowed in general metric theory of gravity as the Bayes factor in favor of the ST-correlated process versus the spatially uncorrelated common-spectrum process is 30 ± 2. The median and the 90% equal-tail amplitudes of ST mode are ST=1.29−0.44+0.51×10−15 , or equivalently the energy density parameter per logarithm frequency is ΩGWST=2.31−1.30+2.19×10−9 , at frequency of 1 year−1. However, we do not find any statistically significant evidence for the tensor transverse (TT) mode and then place the 95% upper limits as TT<3.95×10−15 , or equivalently ΩGWTT<2.16×10−9 , at frequency of 1 year−1.

The observations from pulsar timing arrays (PTAs), led by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), have provided opportunities to constrain primordial gravitational waves at low frequencies. In this paper, we analyze the best-fit parameter values for different Gauss-Bonnet Inflationary Gravitational Wave (GB-IGW) models with the PTArcade program, and we compare the results with the observations of NANOGrav, European Pulsar Timing Array (EPTA), Parkes Pulsar Timing Array (PPTA), and International Pulsar Timing Array (IPTA). We find the potential parameter n derived from the GB-IGW model is not necessarily positive. Instead, PTA data suggest the possibility of new parameter ranges. Meanwhile, the other parameters' results, such as the spectral indices of the scalar and tensor perturbation also have differences with traditional Cosmic Microwave Background analyses from Planck 2018, showing the reason for the different phenomena of Gravitational Wave power spectra expected. Additionally, the GB-IGW models we calculated are closer to the central values of multiple PTA datasets compared to the standard inflation model, making the GB-IGW model more likely to be detected by future space-based gravitational wave observatories.

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released 11 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no evidence for a GWB, we place constraints on a population of inspiraling supermassive black hole (SMBH) binaries, a network of decaying cosmic strings, and a primordial GWB. For the first time, we find that the GWB constraints are sensitive to the solar system ephemeris (SSE) model used and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first pulsar-timing array (PTA) constraints that are robust against SSE errors. We thus place a 95% upper limit on the GW-strain amplitude of A GWB &lt; 1.45 × 10−15 at a frequency of f = 1 yr−1 for a fiducial f −2/3 power-law spectrum and with interpulsar correlations modeled. This is a factor of ∼2 improvement over the NANOGrav nine-year limit calculated using the same procedure. Previous PTA upper limits on the GWB (as well as their astrophysical and cosmological interpretations) will need revision in light of SSE systematic errors. We use our constraints to characterize the combined influence on the GWB of the stellar mass density in galactic cores, the eccentricity of SMBH binaries, and SMBH–galactic-bulge scaling relationships. We constrain the cosmic-string tension using recent simulations, yielding an SSE-marginalized 95% upper limit of Gμ &lt; 5.3 × 10−11—a factor of ∼2 better than the published NANOGrav nine-year constraints. Our SSE-marginalized 95% upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation universe) is ΩGWB(f) h 2 &lt; 3.4 × 10−10.

Observing and timing a group of millisecond pulsars with high rotational stability enables the direct detection of gravitational waves (GWs). The GW signals can be identified from the spatial correlations encoded in the times-of-arrival of widely spaced pulsar-pairs. The Chinese Pulsar Timing Array (CPTA) is a collaboration aiming at the direct GW detection with observations carried out using Chinese radio telescopes. This short article serves as a “table of contents” for a forthcoming series of papers related to the CPTA Data Release 1 (CPTA DR1) which uses observations from the Five-hundred-meter Aperture Spherical radio Telescope. Here, after summarizing the time span and accuracy of CPTA DR1, we report the key results of our statistical inference finding a correlated signal with amplitude logAc=−14.4−2.8+1.0 for spectral index in the range of α ∈ [ − 1.8, 1.5] assuming a GW background (GWB) induced quadrupolar correlation. The search for the Hellings–Downs (HD) correlation curve is also presented, where some evidence for the HD correlation has been found that a 4.6σ statistical significance is achieved using the discrete frequency method around the frequency of 14 nHz. We expect that the future International Pulsar Timing Array data analysis and the next CPTA data release will be more sensitive to the nHz GWB, which could verify the current results.

Pulsar timing arrays (PTAs) have revealed galaxy-size gravitational waves (GWs) in the form of a stochastic gravitational wave background (SGWB), correlating the radio pulses emitted by millisecond pulsars. This discovery naturally leads to the question of the origin and the nature of the SGWB; the latter is synonymous to testing how quadrupolar the inter-pulsar spatial correlation is. In this paper, we investigate the nature of the SGWB by considering correlations beyond the Hellings–Downs (HD) curve of Einstein’s general relativity. We scrutinize the HD and non-Einsteinian GW correlations with the North American Nanohertz Observatory for Gravitational Waves and the Chinese PTA data, and find that both data sets allow a graviton mass of mg ≲ 1.04 × 10−22 eV/c2 and subluminal traveling waves. We discuss gravitational physics scenarios beyond general relativity that could host non-Einsteinian GW correlations in the SGWB and highlight the importance of the cosmic variance inherited from stochastic variations across realizations in interpreting PTA observations.

Observations of low-frequency gravitational waves (GWs) will require the highest possible timing precision from an array of the most spin-stable pulsars. We can improve the sensitivity of a pulsar timing array (PTA) to different GW sources by observing pulsars with low timing noise over years to decades and distributed across the sky. We discuss observing strategies for a PTA focused on a stochastic GW background such as from unresolved supermassive black hole binaries as well as focused on single continuous-wave sources. First, we describe the method to calculate a PTA’s sensitivity to different GW-source classes. We then apply our method to the 45 pulsars presented in the North American Nanohertz Observatory for the GW 11 year data set. For expected amplitudes of the stochastic background, we find that all pulsars contribute significantly over the timescale of decades; the exception is for very pessimistic values of the stochastic-background amplitude. For individual single sources, we find that a number of pulsars contribute to the sensitivity of a given source, but that which pulsars contribute is different depending on the source, or versus an all-sky metric. Our results seem robust to the presence of red noise in pulsar arrival times. It is critical to obtain more robust pulsar-noise parameters as they heavily affect our results. Our results show that it is also imperative to locate and time as many high-precision pulsars as possible, as quickly as possible, to maximize the sensitivity of next-generation PTA detectors.

The gravitational waves (GWs) from supermassive binary black holes (BBHs) have long been sought by pulsar timing array (PTA) experiments, in the forms of both a stochastic GW background (GWB) and individual sources. Evidence for a GWB was reported recently by several PTAs with origins to be determined. Here we use a BBH population synthesis model to investigate the detection probability of individual BBHs by the Chinese PTA (CPTA) and the constraint on the GWB origin that may be obtained by PTA observations of both GWB and individual BBHs. If the detected GWB signal is entirely due to BBHs, a significantly positive redshift evolution (∝ (1 + z)2.07) of the mass scaling relation between supermassive black holes and their host galaxies is required. In this case, we find that the detection probability of individual BBHs is ∼85% or 64% if using a period of 3.4 yr of CPTA observation data, with an expectation of ∼1.9 or 1.0 BBHs detectable with a signal-to-noise ratio ≥3 or 5, and it is expected to increase to >95% if the observation period is extended to 5 yr or longer. Even if the contribution from BBHs to the GWB power signal is as small as ∼10%, a positive detection of individual BBHs can still be expected within an observation period of ∼10 yr. A nondetection of individual BBHs within several years from now jointly with the detected GWB signal can put a strong constraint on the upper limit of the BBH contribution to the GWB signal and help identify/falsify a cosmological origin.

We study the detectability of circular polarization in a stochastic gravitational wave background from various sources such as supermassive black hole binaries, cosmic strings, and inflation in the early universe with pulsar timing arrays. We calculate generalized overlap reduction functions for the circularly polarized stochastic gravitational wave background. We find that the circular polarization can not be detected for an isotropic background. However, there is a chance to observe the circular polarization for an anisotropic gravitational wave background. We also show how to separate polarized gravitational waves from unpolarized gravitational waves.

We perform a direct search for an isotropic stochastic gravitational-wave background (SGWB) produced by cosmic strings in the Parkes Pulsar Timing Array (PPTA) Data Release 2 (DR2). We find no evidence for such an SGWB, and therefore place a 95% confidence level upper limit on the cosmic string tension, G μ, as a function of the reconnection probability, p, which can be less than 1 in the string-theory-inspired models or pure Yang–Mills theory. The upper bound on the cosmic string tension is G μ ≲ 5.1 × 10−10 for p = 1, which is about five orders of magnitude tighter than the bound derived from the null search of individual gravitational-wave bursts from cosmic string cusps in the PPTA DR2, and comparable to previous bounds derived from the null search of the SGWB from cosmic strings.

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 inspiralling 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 both on 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.

A population of supermassive black hole (SMBH) binaries is expected to generate a stochastic gravitational wave background (SGWB) in the pulsar timing array (PTA) frequency range of 10−9 to $10^{-7}\, {\rm Hz}$. Detection of this signal is a current observational goal and so predictions of its characteristics are of significant interest. In this work, we use SMBH binary mergers from the MassiveBlackII simulation to estimate the characteristic strain of the stochastic background. We examine both a gravitational wave (GW) driven model of binary evolution and a model which also includes the effects of stellar scattering and a circumbinary gas disc. Results are consistent with PTA upper limits and similar to estimates in the literature. The characteristic strain at a reference frequency of $1\, {\rm yr}^{-1}$ is found to be $A_{\rm {yr}^{-1}}= 6.9 \times 10^{-16}$ and $A_{\rm {yr}^{-1}}= 6.4 \times 10^{-16}$ in the GW-driven and stellar scattering/gas disc cases, respectively. Using the latter approach, our models show that the SGWB is mildly suppressed compared to the purely GW-driven model as frequency decreases inside the PTA frequency band.