No Arabic abstract
Within the next several years pulsar timing arrays (PTAs) are positioned to detect the stochastic gravitational-wave background (GWB) likely produced by the collection of inspiralling super-massive black holes binaries, and potentially constrain some exotic physics. So far most of the pulsar timing data analysis has focused on the monopole of the GWB, assuming it is perfectly isotropic. The natural next step is to search for anisotropies in the GWB. In this paper, we use the recently developed PTA Fisher matrix to gain insights into optimal search strategies for GWB anisotropies. For concreteness, we apply our results to EPTA data, using realistic noise characteristics of its pulsars. We project the detectability of a GWB whose angular dependence is assumed to be a linear combination of predetermined maps, such as spherical harmonics or coarse pixels. We find that the GWB monopole is always statistically correlated with these maps, implying a loss of sensitivity to the monopole when searching simultaneously for anisotropies. We then derive the angular distributions of the GWB intensity to which a PTA is most sensitive, and illustrate how one may use these principal maps to approximately reconstruct the angular dependence of the GWB. Since the principal maps are neither perfectly anisotropic nor uncorrelated with the monopole, we also develop a frequentist criterion to specifically search for anisotropies in the GWB without any prior knowledge about their angular distribution. Lastly, we show how to recover existing EPTA results with our Fisher formalism, and clarify their meaning. The tools presented here will be valuable in guiding and optimizing the computationally demanding analyses of pulsar timing data.
We present a new veto procedure to distinguish between continuous gravitational wave (CW) signals and the detector artifacts that can mimic their behavior. The veto procedure exploits the fact that a long-lasting coherent disturbance is less likely than a real signal to exhibit a Doppler modulation of astrophysical origin. Therefore, in the presence of an outlier from a search, we perform a multi-step search around the frequency of the outlier with the Doppler modulation turned off (DM-off), and compare these results with the results from the original (DM-on) search. If the results from the DM-off search are more significant than those from the DM-on search, the outlier is most likely due to an artifact rather than a signal. We tune the veto procedure so that it has a very low false dismissal rate. With this veto, we are able to identify as coherent disturbances >99.9% of the 6349 candidates from the recent all-sky low-frequency Einstein@Home search on the data from the Advanced LIGO O1 observing run [1]. We present the details of each identified disturbance in the Appendix.
In hierarchical searches for continuous gravitational waves, clustering of candidates is an important postprocessing step because it reduces the number of noise candidates that are followed-up at successive stages [1][7][12]. Previous clustering procedures bundled together nearby candidates ascribing them to the same root cause (be it a signal or a disturbance), based on a predefined cluster volume. In this paper, we present a procedure that adapts the cluster volume to the data itself and checks for consistency of such volume with what is expected from a signal. This significantly improves the noise rejection capabilities at fixed detection threshold, and at fixed computing resources for the follow-up stages, this results in an overall more sensitive search. This new procedure was employed in the first Einstein@Home search on data from the first science run of the advanced LIGO detectors (O1) [11].
We present an analysis of high-precision pulsar timing data taken as part of the North American Nanohertz Observatory for Gravitational waves (NANOGrav) project. We have observed 17 pulsars for a span of roughly five years using the Green Bank and Arecibo radio telescopes. We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Sub-microsecond timing residuals are obtained in nearly all cases, and the best root-mean-square timing residuals in this set are ~30-50 ns. We present methods for analyzing post-fit timing residuals for the presence of a gravitational wave signal with a specified spectral shape. These optimally take into account the timing fluctuation power removed by the model fit, and can be applied to either data from a single pulsar, or to a set of pulsars to detect a correlated signal. We apply these methods to our dataset to set an upper limit on the strength of the nHz-frequency stochastic supermassive black hole gravitational wave background of h_c (1 yr^-1) < 7x10^-15 (95%). This result is dominated by the timing of the two best pulsars in the set, PSRs J1713+0747 and J1909-3744.
The field of gravitational-wave astronomy has been opened up by gravitational-wave observations made with interferometric detectors. This review surveys the current state-of-the-art in gravitational-wave detectors and data analysis methods currently used by the Laser Interferometer Gravitational-Wave Observatory in the United States and the Virgo Observatory in Italy. These analysis methods will also be used in the recently completed KAGRA Observatory in Japan. Data analysis algorithms are developed to target one of four classes of gravitational waves. Short duration, transient sources include compact binary coalescences, and burst sources originating from poorly modelled or unanticipated sources. Long duration sources include sources which emit continuous signals of consistent frequency, and many unresolved sources forming a stochastic background. A description of potential sources and the search for gravitational waves from each of these classes are detailed.
We investigate the isotropic and anisotropic components of the Stochastic Gravitational Wave Background (SGWB) originated from unresolved merging compact binaries in galaxies. We base our analysis on an empirical approach to galactic astrophysics that allows to follow the evolution of individual systems. We then characterize the energy density of the SGWB as a tracer of the total matter density, in order to compute the angular power spectrum of anisotropies with the Cosmic Linear Anisotropy Solving System (CLASS) public code in full generality. We obtain predictions for the isotropic energy density and for the angular power spectrum of the SGWB anisotropies, and study the prospect for their observations with advanced Laser Interferometer Gravitational-Wave and Virgo Observatories and with the Einstein Telescope. We identify the contributions coming from different type of sources (binary black holes, binary neutron stars and black hole-neutron star) and from different redshifts. We examine in detail the spectral shape of the energy density for all types of sources, comparing the results for the two detectors. We find that the power spectrum of the SGWB anisotropies behaves like a power law on large angular scales and drops at small scales: we explain this behaviour in terms of the redshift distribution of sources that contribute most to the signal, and of the sensitivities of the two detectors. Finally, we simulate a high resolution full sky map of the SGWB starting from the power spectra obtained with CLASS and including Poisson statistics and clustering properties.