An official website of the United States government
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
Quality over quantity: Optimizing pulsar timing array analysis for stochastic and continuous gravitational wave signals
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
- PAR ID:
- 10381483
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 518
- Issue:
- 2
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 1802-1817
- Size(s):
- p. 1802-1817
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Pulsar timing arrays (PTAs) are ensembles of millisecond pulsars observed for years to decades. The primary goal of PTAs is to study gravitational-wave astronomy at nanohertz frequencies, with secondary goals of undertaking other fundamental tests of physics and astronomy. Recently, compelling evidence has emerged in established PTA experiments for the presence of a gravitational-wave background. To accelerate a confident detection of such a signal and then study gravitational-wave emitting sources, it is necessary to observe a larger number of millisecond pulsars to greater timing precision. The SKA telescopes, which will be a factor of three to four greater in sensitivity compared to any other southern hemisphere facility, is poised to make such an impact. In this chapter, we motivate an SKAO pulsar timing array (SKAO PTA) experiment. We discuss the classes of gravitational waves present in PTA observations and how an SKAO PTA can detect and study them. We then describe the sources that can produce these signals. We discuss the astrophysical noise sources that must be mitigated to undertake the most sensitive searches. We then describe a realistic PTA experiment implemented with the SKA and place it in context alongside other PTA experiments likely ongoing in the 2030s. We describe the techniques necessary to search for gravitational waves in the SKAO PTA and motivate how very long baseline interferometry can improve the sensitivity of an SKAO PTA. The SKAO PTA will provide a view of the Universe complementary to those of the other large facilities of the 2030s.more » « less
-
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
-
Abstract Pulsar timing arrays (PTAs) are galactic-scale gravitational wave (GW) detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency GW signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15 yr data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white-noise parameters and two red-noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of 7 × 10−15at 5 nHz. A power-law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav’s 15 yr GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals.more » « less
