An official website of the United States government
Abstract The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational-wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within 1σ. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we “extended” each PTA’s data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings–Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA’s Data Release 3, which will involve not just adding in additional pulsars but also including data from all three PTAs where any given pulsar is timed by more than a single PTA.
more »
« less
Rapid construction of joint pulsar timing array data sets: the Lite method
ABSTRACT The International Pulsar Timing Array (IPTA)’s second data release (IPTA DR2) combines decades of observations of 65 millisecond pulsars from 7 radio telescopes. IPTA data sets should be the most sensitive data sets to nanohertz gravitational waves (GWs), but take years to assemble, often excluding valuable recent data. To address this, we introduce the IPTA ‘Lite’ analysis, where a Figure of Merit is used to select an optimal PTA data set to analyse for each pulsar, enabling immediate access to new data and preliminary results prior to full combination. We test the capabilities of the Lite analysis using IPTA DR2, finding that ‘DR2 Lite’ can be used to detect the common red noise process with an amplitude of $$A = 4.8^{+1.8}_{-1.8} \times 10^{-15}$$ at $$\gamma = 13/3$$. This amplitude is slightly large in comparison to the combined analysis, and likely biased high as DR2 Lite is more sensitive to systematic errors from individual pulsars than the full data set. Furthermore, although there is no strong evidence for Hellings-Downs correlations in IPTA DR2, we still find the full data set is better at resolving Hellings-Downs correlations than DR2 Lite. Alongside the Lite analysis, we also find that analysing a subset of pulsars from IPTA DR2, available at a hypothetical ‘early’ stage of combination (EDR2), yields equally competitive results as the full data set. Looking ahead, the Lite method will enable rapid synthesis of the latest PTA data, offering preliminary GW constraints before the superior full data set combinations are available.
more »
« less
- PAR ID:
- 10635460
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 542
- Issue:
- 4
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 3028-3048
- Size(s):
- p. 3028-3048
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Pulsar distances are notoriously difficult to measure, and play an important role in many fundamental physics experiments, such as pulsar timing arrays. Here, we perform a cross-match between International PTA pulsars (IPTA) and Gaia's Data Release 2 (DR2) and Data Release 3 (DR3). We then combine the IPTA pulsar’s parallax with its binary companion’s parallax, found in Gaia, to improve the distance measurement to the binary. We find seven cross-matched IPTA pulsars in Gaia DR2, and when using Gaia DR3 we find six IPTA pulsar cross-matches but with seven Gaia objects. Moving from Gaia DR2 to Gaia DR3, we find that the Gaia parallaxes for the successfully cross-matched pulsars improved by 53%, and pulsar distances improved by 29%. Finally, we find that binary companions with a <3.0σdetection are unreliable associations, setting a high bar for successful cross-matches.more » « less
-
ABSTRACT The search for gravitational waves using Pulsar Timing Arrays (PTAs) is a computationally expensive complex analysis that involves source-specific noise studies. As more pulsars are added to the arrays, this stage of PTA analysis will become increasingly challenging. Therefore, optimizing the number of included pulsars is crucial to reduce the computational burden of data analysis. Here, we present a suite of methods to rank pulsars for use within the scope of PTA analysis. First, we use the maximization of the signal-to-noise ratio as a proxy to select pulsars. With this method, we target the detection of stochastic and continuous gravitational wave signals. Next, we present a ranking that minimizes the coupling between spatial correlation signatures, namely monopolar, dipolar, and Hellings & Downs correlations. Finally, we also explore how to combine these two methods. We test these approaches against mock data using frequentist and Bayesian hypothesis testing. For equal-noise pulsars, we find that an optimal selection leads to an increase in the log-Bayes factor two times steeper than a random selection for the hypothesis test of a gravitational wave background versus a common uncorrelated red noise process. For the same test but for a realistic European PTA (EPTA) data set, a subset of 25 pulsars selected out of 40 can provide a log-likelihood ratio that is 89 % of the total, implying that an optimally selected subset of pulsars can yield results comparable to those obtained from the whole array. We expect these selection methods to play a crucial role in future PTA data combinations.more » « less
-
ABSTRACT 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.more » « less
-
Abstract Recently, many different pulsar timing array (PTA) collaborations have reported strong evidence for a common stochastic process in their data sets. The reported amplitudes are in tension with previously computed upper limits. In this paper, we investigate how using a subset of a set of pulsars biases Bayesian upper limit recovery. We generate 500 simulated PTA data sets, based on the NANOGrav 11 yr data set with an injected stochastic gravitational-wave background (GWB). We then compute the upper limits by sampling the individual pulsar likelihoods, and combine them through a factorized version of the PTA likelihood to obtain upper limits on the GWB amplitude, using different numbers of pulsars. We find that it is possible to recover an upper limit (95% credible interval) below the injected value, and that it is significantly more likely for this to occur when using a subset of pulsars to compute the upper limit. When picking pulsars to induce the maximum possible bias, we find that the 95% Bayesian upper limit recovered is below the injected value in 10.6% of the realizations (53 of 500). Further, we find that if we choose a subset of pulsars in order to obtain a lower upper limit than when using the full set of pulsars, the distribution of the upper limits obtained from these 500 realizations is shifted to lower-amplitude values.more » « less
