No Arabic abstract
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).
Pulsar timing experiments are currently searching for gravitational waves, and this dissertation focuses on the development and study of the pulsar timing residual models used for continuous wave searches. The first goal of this work is to re-present much of the fundamental physics and mathematics concepts behind the calculations and theory used in pulsar timing. While there exist many reference sources in the literature, I try to offer a fully self-contained explanation of the fundamentals of this research which I hope the reader will find helpful. The next goal broadly speaking has been to further develop the mathematics behind the currently used pulsar timing models for detecting gravitational waves with pulsar timing experiments. I classify four regimes of interest, governed by frequency evolution and wavefront curvature effects incorporated into the timing residual models. Of these four regimes the plane-wave models are well established in previous literature. I add a new regime which I label Fresnel, as I show it becomes important for significant Fresnel numbers describing the curvature of the gravitational wavefront. Then I give two in-depth studies. The first forecasts the ability of future pulsar timing experiments to probe and measure these Fresnel effects. The second further generalizes the models to a cosmologically expanding universe, and I show how the Hubble constant can be measured directly in the most generalized pulsar timing residual model. This offers future pulsar timing experiments the possibility of being able to procure a purely gravitational wave-based measurement of the Hubble constant. The final chapter shows the initial steps taken to extend this work in the future toward Doppler tracking experiments.
In our previous paper cite{PTA1} we derived a generic expression for the pulse redshift the main observable for the Pulsar Timing Array (PTA) experiment for detection of gravitational waves for all possible polarizations induced by modifications of general relativity (GR). In this work we provide a generic expression of the overlap reduction function for PTA without using the short wavelength approximation for tensorial polarization. We are convinced, that the short wavelength approximation is not applicable to the overlap reduction function for PTAs, since the removal of the exponential terms in the integrand would lead to poles for $x, y$ and $l$ polarizations and discontinuities for $+$ and $times$. In this work we provide a series expansion to calculate the integral exactly and investigate the behaviour of the series for short wavelength values via numerical evaluation of the analytical series. We find a disagreement for the limit of co-located pulsars with the Hellings & Downs curve.
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.
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.
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.