Time-domain Implementation of the Optimal Cross-Correlation Statistic for Stochastic Gravitational-Wave Background Searches in Pulsar Timing Data

Supermassive black hole binaries, cosmic strings, relic gravitational waves from inflation, and first order phase transitions in the early universe are expected to contribute to a stochastic background of gravitational waves in the 10^(-9) Hz-10^(-7) Hz frequency band. Pulsar timing arrays (PTAs) exploit the high precision timing of radio pulsars to detect signals at such frequencies. Here we present a time-domain implementation of the optimal cross-correlation statistic for stochastic background searches in PTA data. Due to the irregular sampling typical of PTA data as well as the use of a timing model to predict the times-of-arrival of radio pulses, time-domain methods are better suited for gravitational wave data analysis of such data. We present a derivation of the optimal cross-correlation statistic starting from the likelihood function, a method to produce simulated stochastic background signals, and a rigorous derivation of the scaling laws for the signal-to-noise ratio of the cross-correlation statistic in the two relevant PTA regimes: the weak signal limit where instrumental noise dominates over the gravitational wave signal at all frequencies, and a second regime where the gravitational wave signal dominates at the lowest frequencies.

The Loudest Event Statistic: General Formulation, Properties and Applications

The use of the loudest observed event to generate statistical statements about rate and strength has become standard in searches for gravitational waves from compact binaries and pulsars. The Bayesian formulation of the method is generalized in this paper to allow for uncertainties both in the background estimate and in the properties of the population being constrained. The method is also extended to allow rate interval construction. Finally, it is shown how to combine the results from multiple experiments and a comparison is drawn between the upper limit obtained in a single search and the upper limit obtained by combining the results of two experiments each of half the original duration. To illustrate this, we look at an example case, motivated by the search for gravitational waves from binary inspiral.