We study how to probe bispectra of stochastic gravitational waves with pulsar timing arrays. The bispectrum is a key to probe the origin of stochastic gravitational waves. In particular, the shape of the bispectrum carries valuable information of inflation models. We show that an appropriate filter function for three point correlations enables us to extract a specific configuration of momentum triangles in bispectra. We also calculate the overlap reduction functions and discuss strategy for detecting the bispectrum with multiple pulsars.
The detection of a stochastic background of low-frequency gravitational waves by pulsar-timing and astrometric surveys will enable tests of gravitational theories beyond general relativity. These theories generally permit gravitational waves with non-Einsteinian polarization modes, which may propagate slower than the speed of light. We use the total-angular-momentum wave formalism to derive the angular correlation patterns of observables relevant for pulsar timing arrays and astrometry that arise from a background of subluminal gravitational waves with scalar, vector, or tensor polarizations. We find that the pulsar timing observables for the scalar longitudinal mode, which diverge with source distance in the luminal limit, are finite in the subluminal case. Furthermore, we apply our results to $f(R)$ gravity, which contains a massive scalar degree of freedom in addition to the standard transverse-traceless modes. The scalar mode in this $f(R)$ theory is a linear combination of the scalar-longitudinal and scalar-transverse modes, exciting only the monopole and dipole for pulsar timing arrays and only the dipole for astrometric surveys.
We introduce the concept of stationary graviton non-Gaussianity (nG), an observable that can be probed in terms of 3-point correlation functions of a stochastic gravitational wave (GW) background. When evaluated in momentum space, stationary nG corresponds to folded bispectra of graviton nG. We determine 3-point overlap functions for testing stationary nG with pulsar timing array GW experiments, and we obtain the corresponding optimal signal-to-noise ratio. For the first time, we consider 3-point overlap functions including scalar graviton polarizations (which can be motivated in theories of modified gravity); moreover, we also calculate 3-point overlap functions for correlating pulsar timing array with ground based GW detectors. The value of the optimal signal-to-noise ratio depends on the number and position of monitored pulsars. We build geometrical quantities characterizing how such ratio depends on the pulsar system under consideration, and we evaluate these geometrical parameters using data from the IPTA collaboration. We quantitatively show how monitoring a large number of pulsars can increase the signal-to-noise ratio associated with measurements of stationary graviton nG.
Gravitational wave burst is a catch-all category for signals whose durations are shorter than the observation period. We apply a method new to gravitational wave data analysis --- Bayesian non-parameterics --- to the problem of gravitational wave detection, with an emphasis on pulsar timing array observations. In Bayesian non-parametrics, constraints are set on the function space that may be reasonably thought to characterize the range of gravitational-wave signals. This differs from the approaches currently employed or proposed, which focus on introducing parametric signal models or looking for excess power as evidence of the presence of a gravitational wave signal. Our Bayesian nonparametrics analysis method addresses two issues: (1) investigate if a gravitational wave burst is present in the data; (2) infer the sky location of the source and the duration of the burst. Compared with the popular method proposed by Finn & Lommen, our method improves in two aspects: (1) we can estimate the burst duration by adding the prior that the gravitational wave signals are smooth, while Finn & Lommen ignored this important point; (2) we perform a full Bayesian analysis by marginalizing over all possible parameters and provide robust inference on the presence of gravitational waves, while Finn & Lommen chose to optimize over parameters, which would increase false alarm risk and also underestimate the parameter uncertainties.
Cosmic strings are potential gravitational wave (GW) sources that can be probed by pulsar timing arrays (PTAs). In this work we develop a detection algorithm for a GW burst from a cusp on a cosmic string, and apply it to Parkes PTA data. We find four events with a false alarm probability less than 1%. However further investigation shows that all of these are likely to be spurious. As there are no convincing detections we place upper limits on the GW amplitude for different event durations. From these bounds we place limits on the cosmic string tension of G mu ~ 10^{-5}, and highlight that this bound is independent from those obtained using other techniques. We discuss the physical implications of our results and the prospect of probing cosmic strings in the era of Square Kilometre Array (SKA).
Recent years have seen a burgeoning interest in using pulsar timing arrays (PTAs) as gravitational-wave (GW) detectors. To date, that interest has focused mainly on three particularly promising source types: supermassive--black-hole binaries, cosmic strings, and the stochastic background from early-Universe phase transitions. In this paper, by contrast, our aim is to investigate the PTA potential for discovering unanticipated sources. We derive significant constraints on the available discovery space based solely on energetic and statistical considerations: we show that a PTA detection of GWs at frequencies above ~3.e-5 Hz would either be an extraordinary coincidence or violate cherished beliefs; we show that for PTAs GW memory can be more detectable than direct GWs, and that, as we consider events at ever higher redshift, the memory effect increasingly dominates an events total signal-to-noise ratio. The paper includes also a simple analysis of the effects of pulsar red noise in PTA searches, and a demonstration that the effects of periodic GWs in the 10^-8 -- 10^-4.5 Hz band would not be degenerate with small errors in standard pulsar parameters (except in a few narrow bands).