skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Forecasting Pulsar Timing Array Sensitivity to Anisotropy in the Stochastic Gravitational Wave Background
Abstract Statistical anisotropy in the nanohertz-frequency gravitational wave background (GWB) is expected to be detected by pulsar timing arrays (PTAs) in the near future. By developing a frequentist statistical framework that intrinsically restricts the GWB power to be positive, we establish scaling relations for multipole-dependent anisotropy decision thresholds that are a function of the noise properties, timing baselines, and cadences of the pulsars in a PTA. We verify that (i) a larger number of pulsars, and (ii) factors that lead to lower uncertainty on spatial cross-correlation measurements between pulsars, lead to a higher overall GWB signal-to-noise ratio, and lower anisotropy decision thresholds with which to reject the null hypothesis of isotropy. Using conservative simulations of realistic NANOGrav data sets, we predict that an anisotropic GWB with angular power C l =1 > 0.3 C l =0 may be sufficient to produce tension with isotropy at the p = 3 × 10 −3 (∼3 σ ) level in near-future NANOGrav data with a 20 yr baseline. We present ready-to-use scaling relationships that can map these thresholds to any number of pulsars, configuration of pulsar noise properties, or sky coverage. We discuss how PTAs can improve the detection prospects for anisotropy, as well as how our methods can be adapted for more versatile searches.  more » « less
Award ID(s):
2146016 2020265 2007993
PAR ID:
10383488
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
940
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
173
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; ; ; ; et al (, The Astrophysical Journal)
    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
  2. ; ; (, The Astrophysical Journal)
    Abstract Pulsar timing arrays (PTAs) are sensitive to low-frequency gravitational waves (GWs), which induce correlated changes in millisecond pulsars’ timing residuals. PTA collaborations around the world have recently announced evidence of a nanohertz gravitational wave background (GWB), which may be produced by a population of supermassive black hole binaries (SMBHBs). The GWB is often modeled as following a power-law power spectral density (PSD); however, a GWB produced by a cosmological population of SMBHBs is expected to have a more complex power spectrum due to the discrete nature of the sources. In this paper, we investigate using at-process PSD to model the GWB, which allows us to fit for both the underlying power-law amplitude and spectral index as well as deviations from that power law, which may be produced by individual nearby binaries. We create simulated data sets based on the properties of the NANOGrav 15 yr data set, and we demonstrate that thet-process PSD can accurately recover the PSD when deviations from a power law are present. With longer timed data sets and more pulsars, we expect the sensitivity of our PTAs to improve, which will allow us to precisely measure the PSD of the GWB and study the sources producing it. 
    more » « less
  3. ; ; ; ; ; ; ; ; ; ; et al (, The Astrophysical Journal Letters)
    Abstract The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational-wave background (GWB) in its 15 yr data set. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy. By modeling the angular power distribution as a sum over spherical harmonics (where the coefficients are not bound to always generate positive power everywhere), we find that the Bayesian 95% upper limit on the level of dipole anisotropy is (Cl=1/Cl=0) < 27%. This is similar to the upper limit derived under the constraint of positive power everywhere, indicating that the dipole may be close to the data-informed regime. By contrast, the constraints on anisotropy at higher spherical-harmonic multipoles are strongly prior dominated. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15 yr data set and show that this data set has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data. 
    more » « less
  4. ; ; (, Classical and Quantum Gravity)
    Abstract As pulsar timing arrays (PTAs) transition into the detection era of the stochastic gravitational wave background (GWB), it is important for PTA collaborations to review and possibly revise their observing campaigns. The detection of a “single source” would be a boon for gravitational astrophysics, as such a source would emit gravitational waves for millions of years in the PTA frequency band. Here we present generic methods for studying the effects of various observational strategies, taking advantage of detector sensitivity curves, i.e., noise-averaged, frequency-domain detection statistics. The statistical basis for these methods is presented along with myriad examples of how to tune a detector towards single, deterministic signals or a stochastic background. We demonstrate that trading observations of the worst pulsars for high cadence campaigns on the best pulsars increases sensitivity to single sources at high frequencies while hedging losses in GWB and single source sensitivity at low frequencies. We also find that sky-targeted observing campaigns yield minimal sensitivity improvements compared with other PTA tuning options. Lastly, we show the importance of the uncorrelated half of the GWB, i.e. the pulsar-term, as an increasingly prominent sources of noise and show the impact of this emerging noise source on various PTA configurations. 
    more » « less
  5. ; ; ; ; ; ; ; ; ; ; et al (, The Astrophysical Journal)
    Abstract Evidence has emerged for a stochastic signal correlated among 67 pulsars within the 15 yr pulsar-timing data set compiled by the NANOGrav collaboration. Similar signals have been found in data from the European, Indian, Parkes, and Chinese pulsar timing arrays. This signal has been interpreted as indicative of the presence of a nanohertz stochastic gravitational-wave background (GWB). To explore the internal consistency of this result, we investigate how the recovered signal strength changes as we remove the pulsars one by one from the data set. We calculate the signal strength using the (noise-marginalized) optimal statistic, a frequentist metric designed to measure the correlated excess power in the residuals of the arrival times of the radio pulses. We identify several features emerging from this analysis that were initially unexpected. The significance of these features, however, can only be assessed by comparing the real data to synthetic data sets. After conducting identical analyses on simulated data sets, we do not find anything inconsistent with the presence of a stochastic GWB in the NANOGrav 15 yr data. The methodologies developed here can offer additional tools for application to future, more sensitive data sets. While this analysis provides an internal consistency check of the NANOGrav results, it does not eliminate the necessity for additional investigations that could identify potential systematics or uncover unmodeled physical phenomena in the data. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email