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.
We extend the formalisms developed in Gair et al. and Cornish and van Haasteren to create maps of gravitational-wave backgrounds using a network of ground-based laser interferometers. We show that in contrast to pulsar timing arrays, which are insensitive to half of the gravitational-wave sky (the curl modes), a network of ground-based interferometers is sensitive to both the gradient and curl components of the background. The spatial separation of a network of interferometers, or of a single interferometer at different times during its rotational and orbital motion around the Sun, allows for recovery of both components. We derive expressions for the response functions of a laser interferometer in the small-antenna limit, and use these expressions to calculate the overlap reduction function for a pair of interferometers. We also construct maximum-likelihood estimates of the + and x-polarization modes of the gravitational-wave sky in terms of the response matrix for a network of ground-based interferometers, evaluated at discrete times during Earths rotational and orbital motion around the Sun. We demonstrate the feasibility of this approach for some simple simulated backgrounds (a single point source and spatially-extended distributions having only grad or curl components), calculating maximum-likelihood sky maps and uncertainty maps based on the (pseudo)inverse of the response matrix. The distinction between this approach and standard methods for mapping gravitational-wave power is also discussed.
We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar System ephemeris errors, obtaining a robust 95$%$ upper limit on the dimensionless strain amplitude $A$ of the background of $A<3.0times 10^{-15}$ at a reference frequency of $1mathrm{yr^{-1}}$ and a spectral index of $13/3$, corresponding to a background from inspiralling super-massive black hole binaries, constraining the GW energy density to $Omega_mathrm{gw}(f)h^2 < 1.1times10^{-9}$ at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of $sim 5times10^{-9}$~Hz. Finally we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95$%$ upper limits on the string tension, $Gmu/c^2$, characterising a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu-Goto field theory cosmic string network, we set a limit $Gmu/c^2<1.3times10^{-7}$, identical to that set by the {it Planck} Collaboration, when combining {it Planck} and high-$ell$ Cosmic Microwave Background data from other experiments. For a stochastic relic background we set a limit of $Omega^mathrm{relic}_mathrm{gw}(f)h^2<1.2 times10^{-9}$, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.
