No Arabic abstract
We describe an alternative approach to the analysis of gravitational-wave backgrounds, based on the formalism used to characterise the polarisation of the cosmic microwave background. In contrast to standard analyses, this approach makes no assumptions about the nature of the background and so has the potential to reveal much more about the physical processes that generated it. An arbitrary background can be decomposed into modes whose angular dependence on the sky is given by gradients and curls of spherical harmonics. We derive the pulsar timing overlap reduction functions for the individual modes, which are given by simple combinations of spherical harmonics evaluated at the pulsar locations. We show how these can be used to recover the components of an arbitrary background, giving explicit results for both isotropic and anisotropic uncorrelated backgrounds. We also find that the response of a pulsar timing array to curl modes is identically zero, so half of the gravitational-wave sky will never be observed using pulsar timing, no matter how many pulsars are included in the array. An isotropic, unpolarised and uncorrelated background can be accurately represented using only three modes, and so a search of this type will be only slightly more complicated than the standard cross-correlation search using the Hellings and Downs overlap reduction function. However, by measuring the components of individual modes of the background and checking for consistency with isotropy, this approach has the potential to reveal much more information. Each individual mode on its own describes a background that is correlated between different points on the sky. A measurement of the components that indicates the presence of correlations in the background on large angular scales would suggest startling new physics.
We extend our previous work on applying CMB techniques to the mapping of gravitational-wave backgrounds to backgrounds which have non-GR polarisations. Our analysis and results are presented in the context of pulsar-timing array observations, but the overarching methods are general, and can be easily applied to LIGO or eLISA observations using appropriately modified response functions. Analytic expressions for the pulsar-timing response to gravitational waves with non-GR polarisation are given for each mode of a spin-weighted spherical-harmonic decomposition of the background, which permit the signal to be mapped across the sky to any desired resolution. We also derive the pulsar-timing overlap reduction functions for the various non-GR polarisations, finding analytic forms for anisotropic backgrounds with scalar-transverse (breathing) and vector-longitudinal polarisations, and a semi-analytic form for scalar-longitudinal backgrounds. Our results indicate that pulsar-timing observations will be completely insensitive to scalar-transverse mode anisotropies in the polarisation amplitude beyond dipole, and anisotropies in the power beyond quadrupole. Analogously to our previous findings that pulsar-timing observations lack sensitivity to tensor-curl modes for a transverse-traceless tensor background, we also find insensitivity to vector-curl modes for a vector-longitudinal background.
Within the next several years, pulsar-timing array programs will likely usher in the next era of gravitational-wave astronomy through the detection of a stochastic background of nanohertz-frequency gravitational waves, originating from a cosmological population of inspiraling supermassive binary black holes. While the source positions will likely be isotropic to a good approximation, the gravitational-wave angular power distribution will be anisotropic, with the most massive and/or nearby binaries producing signals that may resound above the background. We study such a realistic angular power distribution, developing fast and accurate sky-mapping strategies to localize pixels and extended regions of excess power while simultaneously modeling the background signal from the less massive and more distant ensemble. We find that power anisotropy will be challenging to discriminate from isotropy for realistic gravitational-wave skies, requiring SNR $>10$ in order to favor anisotropy with $10:1$ posterior odds in our case study. Amongst our techniques, modeling the population signal with multiple point sources in addition to an isotropic background provides the most physically-motivated and easily interpreted maps, while spherical-harmonic modeling of the square-root power distribution, $P(hatOmega)^{1/2}$, performs best in discriminating from overall isotropy. Our techniques are modular and easily incorporated into existing pulsar-timing array analysis pipelines.
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).
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.
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has recently reported strong statistical evidence for a common-spectrum red-noise process for all pulsars, as seen in their 12.5-yr analysis for an isotropic stochastic gravitational-wave background. However, there is currently very little evidence for quadrupolar spatial correlations across the pulsars in the array, which is needed to make a confident claim of detection of a stochastic background. Here we give a frequentist analysis of a very simple signal+noise model showing that the current lack of evidence for spatial correlations is consistent with the magnitude of the correlation coefficients for pairs of Earth-pulsar baselines in the array, and the fact that pulsar timing arraysbare most-likely operating in the intermediate-signal regime. We derive analytic expressions that allow one to compare the expected values of the signal-to-noise ratios for both the common-spectrum and cross-correlation estimators.