Placing limits on the stochastic gravitational-wave background using European Pulsar Timing Array data

Abstract: Direct detection of low-frequency gravitational waves ($10^{-9} - 10^{-8}$ 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 $h_c(f)=A (f/\hbox{yr}^{-1})^{\alpha}$, where $f$ is the gravitational-wave 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 $\alpha$ with a Bayesian algorithm, by modelling the GWB as a random Gaussian process. For the case $\alpha=-2/3$, which is expected if the GWB is produced by supermassive black-hole binaries, we obtain a 95% confidence upper limit on $A$ of $6\times 10^{-15}$, which is 1.8 times lower than the 95% confidence GWB limit obtained by the Parkes PTA in 2006. Our approach to the data analysis incorporates the multi-telescope nature of the European PTA and thus can serve as a useful template for future intercontinental PTA collaborations.

• The paper presents a Bayesian analysis that constrains the stochastic gravitational-wave background amplitude to 6×10⁻¹⁵ with 95% confidence.
• It employs multi-telescope EPTA datasets to enhance pulsar timing cross-correlation precision in the 10⁻⁹–10⁻⁸ Hz range.
• The improved limits on gravitational waves challenge models of supermassive black-hole binaries and alternative sources like cosmic strings.

Analysis of Gravitational-Wave Background Constraints from the European Pulsar Timing Array

The paper "Placing limits on the stochastic gravitational-wave background using European Pulsar Timing Array data" by R. van Haasteren et al. presents a significant empirical study aimed at setting stringent upper limits on the stochastic gravitational-wave background (GWB) using data from the European Pulsar Timing Array (EPTA). This work utilizes pulsar timing, specifically with millisecond pulsars, as a tool to investigate the GWB within the low-frequency range of $10^{-9}$-$10^{-8}$ Hz. This research is entrenched in the context of testing general relativity and expanding the observational capabilities in the domain of gravitational-wave astronomy.

Methodological Framework

The study employs a Bayesian statistical framework to set the upper limit on the GWB amplitude. The GWB is modeled as a random Gaussian process, and the characteristic strain is approximated by a power-law form, $h_c(f)=A (f/\text{yr}^{-1})^{\alpha}$, where $f$ is the gravitational wave frequency. This model is motivated by the theoretical consideration that the GWB is likely generated predominantly by an ensemble of supermassive black-hole binaries, which predicts a spectral slope $\alpha = -2/3$.

The analysis incorporates multi-telescope datasets from the EPTA, which include several European observatories, offering a comprehensive combination of datasets characteristic of a collaborative PTA effort. The inclusion of multiple telescopes enhances the precision in the timing analysis by allowing for a more accurate cross-correlation of timing residuals observed across different locations.

Results

Key findings from the paper include:

• A 95% confidence upper limit on the GWB amplitude $A$ of $6 \times 10^{-15}$, assuming a spectral index $\alpha = -2/3$. This represents an improvement by a factor of 1.8 over the previously established limit by the Parkes Pulsar Timing Array in 2006.
• The Bayesian approach leveraging cross-correlations between different pulsar timing sites provided a robust framework adaptable for future intercontinental collaborations.

Implications

The results from this study impose increasingly stringent constraints on the characteristics of the GWB. The implications of placing such tight upper limits are multifaceted:

• The limits closely approach the predictions from theoretical models involving supermassive black-hole binaries. Therefore, they are beginning to constrain some of these models more tightly, providing valuable input to astrophysical scenarios of galaxy evolution and merger histories.
• The constraints derived also bear significance for alternative models of GWB sources, such as cosmic strings, which are speculative in nature but possess crucial theoretical implications. A failure to detect a GWB within the newly established limits narrows the parameter space for cosmic-string models particularly concerning the string tension $G\mu$.

Future Prospects

The methodology and improvement in sensitivity to the GWB as delineated in this paper present a promising stepping stone for future investigations. The continued improvement in timing precision and the extension of observational baselines will enhance the effectiveness of PTAs. Moreover, the prospect of an International Pulsar Timing Array (IPTA)—by unifying efforts globally—hints at the promisingly imminent ability to probe deeper into the frequencies conducive to gravitational-wave detection.

In conclusion, this paper illustrates the potential of PTAs as a collection of the most precise natural clocks in probing foundational aspects of gravitational-wave physics and the cosmos at large. Future advancements could usher in a new epoch of direct gravitational-wave detection at nanohertz frequencies, substantially enriching our understanding of the universe's structure and dynamics.