Toward robust detections of nanohertz gravitational waves

Abstract

The recent observation of a common red-noise process in pulsar timing arrays (PTAs) suggests that the detection of nanohertz gravitational waves might be around the corner. However, in order to confidently attribute this red process to gravitational waves, one must observe the Hellings-Downs curve -- the telltale angular correlation function associated with a gravitational-wave background. This effort is complicated by the complex modelling of pulsar noise. Without proper care, mis-specified noise models can lead to false-positive detections. Background estimation using bootstrap methods such as sky scrambles and phase shifts, which use the data to characterize the noise, are therefore important tools for assessing significance. We investigate the ability of current PTA experiments to estimate their background with "quasi-independent" scrambles -- characterized by a statistical "match" below the fiducial value: $|M|<0.1$. We show that sky scrambling is affected by "saturation" after $O(10)$ quasi-independent realizations; subsequent scrambles are no longer quasi-independent. We show phase scrambling saturates after $O(100)$ quasi-independent realizations. With so few independent scrambles, it is difficult to make reliable statements about the $\gtrsim 5\sigma$ tail of the null distribution of the detection statistic. We discuss various methods by which one may increase the number of independent scrambles. We also consider an alternative approach wherein one re-frames the background estimation problem so that the significance is calculated using statistically dependent scrambles. The resulting $p$-value is in principle well-defined but may be susceptible to failure if assumptions about the data are incorrect.

Paper accepted for publication on ApJ on the 7th of August 2023

References in corpus (20)

• The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background
• Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array
• Searching for the nano-Hertz stochastic gravitational wave background with the Chinese Pulsar Timing Array Data Release I
• The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background
• The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background
• On the evidence for a common-spectrum process in the search for the nanohertz gravitational-wave background with the Parkes Pulsar Timing Array
• Common-red-signal analysis with 24-yr high-precision timing of the European Pulsar Timing Array: Inferences in the stochastic gravitational-wave background search
• Detection of 107 glitches in 36 southern pulsars
• The MeerKAT Telescope as a Pulsar Facility: System verification and early science results from MeerTime
• The Parkes Pulsar Timing Array Project: Second data release
• Identifying and mitigating noise sources in precision pulsar timing data sets
• The MeerKAT Pulsar Timing Array: First Data Release
• Time-domain Implementation of the Optimal Cross-Correlation Statistic for Stochastic Gravitational-Wave Background Searches in Pulsar Timing Data
• The disturbance of a millisecond pulsar magnetosphere
• Gravitational wave background from rotating neutron stars
• Detecting a Stochastic Gravitational-Wave Background: The Overlap Reduction Function
• Evidence for profile changes in PSR J1713+0747 using the uGMRT
• Systematic errors in estimation of gravitational-wave candidate significance
• Gravitational wave detection without boot straps: a Bayesian approach
• Bayesian search for gravitational wave bursts in pulsar timing array data

Cited by in corpus (8)

• The NANOGrav 15-year Gravitational-Wave Background Methods
• Pulsar Timing Arrays require hierarchical models
• Systematic errors in searches for nanohertz gravitational waves
• Reading signatures of supermassive binary black holes in pulsar timing array observations
• The NANOGrav 15 yr data set: Posterior predictive checks for gravitational-wave detection with pulsar timing arrays
• On the calculation of p-values for quadratic statistics in Pulsar Timing Arrays
• Measuring the crust-superfluid coupling time-scale for 105 UTMOST pulsars with a Kalman filter
• Model-agnostic search of gravitational wave echoes in LVK data