No Arabic abstract
We use population inference to explore the impact that uncertainties in the distribution of binary black holes (BBH) have on the astrophysical gravitational-wave background (AGWB). Our results show that the AGWB monopole is sensitive to the nature of the BBH population (particularly the local merger rate), while the anisotropic $C_ell$ spectrum is only modified to within a few percent, at a level which is insignificant compared to other sources of uncertainty (such as cosmic variance). This is very promising news for future observational studies of the AGWB, as it shows that (i) the monopole can be used as a new probe of the population of compact objects throughout cosmic history, complementary to direct observations by LIGO and Virgo and (ii) we are able to make surprisingly robust predictions for the $C_ell$ spectrum, even with only very approximate knowledge of the black hole population. As a result, the AGWB anisotropies have enormous potential as a new probe of the large-scale structure of the Universe, and of late-Universe cosmology in general.
In the literature different approaches have been proposed to compute the anisotropies of the astrophysical gravitational wave background. The different expressions derived, although starting from our work Cusin, Pitrou, Uzan, Phys.Rev.D96, 103019 (2017) [1], seem to differ. This article compares the various theoretical expressions proposed so far and provides a separate derivation based on a Boltzmann approach. We show that all the theoretical formula in the literature are equivalent and boil down to the one of Ref. [1] when a proper matching of terms and integration by parts are performed. The difference between the various predictions presented for anisotropies in a cosmological context can only lie in the astrophysical modeling of sources, and neither in the theory nor in the cosmological description of the large scale structures. Finally we comment on the gauge invariance of expressions.
We calculate the noise induced in the anisotropies of the astrophysical gravitational-wave background by finite sampling of both the galaxy distribution and the compact binary coalescence event rate. This shot noise leads to a scale-invariant bias term in the angular power spectrum $C_ell$, for which we derive a simple analytical expression. We find that this bias dominates over the true cosmological power spectrum in any reasonable observing scenario, and that only with very long observing times and removal of a large number of foreground sources can the true power spectrum be recovered.
We investigate the discrepancy pointed out by Jenkins et al. in Ref. [1] between the predictions of anisotropies of the astrophysical gravitational wave (GW) background, derived using different methods in Cusin et al. [2] and in Jenkins et al. [3]. We show that this discrepancy is not due to our treatment of galaxy clustering, contrary to the claim made in Ref. [1] and we show that our modeling of clustering gives results in very good agreement with observations. Furthermore we show that the power law spectrum used in Refs. [1] and [3] to describe galaxy clustering is incorrect on large scales and leads to a different scaling of the multipoles $C_ell$. Moreover, we also explain that the analytic derivation of the gravitational wave background correlation function in Refs. [1] and [3] is mathematically ill-defined and predicts an amplitude of the angular power spectrum which depends on the (arbitrary) choice of a non-physical cut-off.
The formation of merging binary black holes can occur through multiple astrophysical channels such as, e.g., isolated binary evolution and dynamical formation or, alternatively, have a primordial origin. Increasingly large gravitational-wave catalogs of binary black-hole mergers have allowed for the first model selection studies between different theoretical predictions to constrain some of their model uncertainties and branching ratios. In this work, we show how one could add an additional and independent constraint to model selection by using the stochastic gravitational-wave background. In contrast to model selection analyses that have discriminating power only up to the gravitational-wave detector horizons (currently at redshifts $zlesssim 1$ for LIGO-Virgo), the stochastic gravitational-wave background accounts for the redshift integration of all gravitational-wave signals in the Universe. As a working example, we consider the branching ratio results from a model selection study that includes potential contribution from astrophysical and primordial channels. We renormalize the relative contribution of each channel to the detected event rate to compute the total stochastic gravitational-wave background energy density. The predicted amplitude lies below the current observational upper limits of GWTC-2 by LIGO-Virgo, indicating that the results of the model selection analysis are not ruled out by current background limits. Furthermore, given the set of population models and inferred branching ratios, we find that, even though the predicted background will not be detectable by current generation gravitational-wave detectors, it will be accessible by third-generation detectors such as the Einstein Telescope and space-based detectors such as LISA.
The black hole merging rates inferred after the gravitational-wave detection by Advanced LIGO/VIRGO and the relatively high mass of the progenitors are consistent with models of dark matter made of massive primordial black holes (PBH). PBH binaries emit gravitational waves in a broad range of frequencies that will be probed by future space interferometers (LISA) and pulsar timing arrays (PTA). The amplitude of the stochastic gravitational-wave background expected for PBH dark matter is calculated taking into account various effects such as initial eccentricity of binaries, PBH velocities, mass distribution and clustering. It allows a detection by the LISA space interferometer, and possibly by the PTA of the SKA radio-telescope. Interestingly, one can distinguish this background from the one of non-primordial massive binaries through a specific frequency dependence, resulting from the maximal impact parameter of binaries formed by PBH capture, depending on the PBH velocity distribution and their clustering properties. Moreover, we find that the gravitational wave spectrum is boosted by the width of PBH mass distribution, compared with that of the monochromatic spectrum. The current PTA constraints already rule out broad-mass PBH models covering more than three decades of masses, but evading the microlensing and CMB constraints due to clustering.