In this paper we study gravitational lensing of gravitational wave events. The main point of the present work is to introduce a semi-analytic approach so that each ingredient can be varied and tested individually. Our analytic model for the source population is motivated by a numerical study and we compare semi-analytical and numerical results. We determine the expected magnification for events seen at a given observed luminosity distance. We find that while the probability of significant magnification of the observed LIGO-Virgo events is very small, the probability distribution of the magnification has a significant tail to high $\mu$ such that e.g. the variance of the magnification is very large and even diverges in the geometric optics approximation. For the 10 binary black hole mergers observed by LIGO-Virgo in the O1+O2 observation campaigns, the probability that one of them has been magnified with magnification of 5 or bigger is $\mathcal{P}_{\text{obs}}(>5)\sim 0.01$ while the probability of magnification by 50 or bigger is $\mathcal{P}_{\text{obs}}(>50)\sim 0.005$.

Given the recent advances in gravitational-wave detection technologies, the detection and characterisation of gravitational-wave backgrounds (GWBs) with the Laser Interferometer Space Antenna (LISA) is a real possibility. To assess the abilities of the LISA satellite network to reconstruct anisotropies of different angular scales and in different directions on the sky, we develop a map-maker based on an optimal quadratic estimator. The resulting maps are maximum likelihood representations of the GWB intensity on the sky integrated over a broad range of frequencies. We test the algorithm by reconstructing known input maps with different input distributions and over different frequency ranges. We find that, in an optimal scenario of well understood noise and high frequency, high SNR signals, the maximum scales LISA may probe are $\ell_{\rm max} \lesssim 15$. The map-maker also allows to test the directional dependence of LISA noise, providing insight on the directional sky sensitivity we may expect.

We construct a model for the angular power spectrum of the instrumental noise in interferometer networks mapping gravitational wave backgrounds (GWBs) as a function of detector noise properties, network configuration and scan strategy. We use the model to calculate the noise power spectrum for current and future ground-based experiments, as well as for planned space missions. We present our results in a language similar to that used in cosmic microwave background and intensity mapping experiments, and connect the formalism with the sensitivity curves that are common lore in GWB analyses. Our formalism is implemented in a lightweight python module that we make publicly available at https://github.com/damonge/schNell.

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.

The spatial and temporal discreteness of gravitational wave sources leads to shot noise that may, in some regimes, swamp any attempts at measuring the anisotropy of the gravitational wave background. Cross-correlating a gravitational wave background map with a sufficiently dense galaxy survey can alleviate this issue, and potentially recover some of the underlying properties of the gravitational wave background. We quantify the shot noise level and we explicitly show that cross-correlating the gravitational wave background and a galaxy catalog improves the chances of a first detection of the background anisotropy with a gravitational wave observatory operating in the frequency range (10Hz,100Hz), given sufficient sensitivity.

It is standard practice to study the lensing of gravitational waves (GW) using the geometric optics regime. However, in many astrophysical configurations this regime breaks down as the wavelength becomes comparable to the Schwarzschild radius of the lens. We revisit the lensing of GW including corrections beyond geometric optics. We propose a perturbative method for calculating these corrections simply solving first order decoupled differential equations. We study the behavior of a single ray and find that the polarization plane defined in geometric optics is smeared due to diffraction effects, which leads to the rise of apparent vector and scalar polarization modes. We analyze how these modes depend on the observer choice, and we study the impact of diffraction on the pseudo-stress energy momentum tensor of the gravitational field.

We develop an approach to parametrize cosmological perturbations beyond linear order for general dark energy and modified gravity models characterized by a single scalar degree of freedom. We derive the full nonlinear action, focusing on Horndeski theories. In the quasi-static, non-relativistic limit, there are a total of six independent relevant operators, three of which start at nonlinear order. The new nonlinear couplings modify, beyond linear order, the generalized Poisson equation relating the Newtonian potential to the matter density contrast. We derive this equation up to cubic order in perturbations and, in a companion article, we apply it to compute the one-loop matter power spectrum. Within this approach, we also discuss the Vainshtein regime around spherical sources and the relation between the Vainshtein scale and the nonlinear scale for structure formation.

We develop an approach to compute observables beyond the linear regime of dark matter perturbations for general dark energy and modified gravity models. We do so by combining the Effective Field Theory of Dark Energy and Effective Field Theory of Large-Scale Structure approaches. In particular, we parametrize the linear and nonlinear effects of dark energy on dark matter clustering in terms of the Lagrangian terms introduced in a companion paper, focusing on Horndeski theories and assuming the quasi-static approximation. The Euler equation for dark matter is sourced, via the Newtonian potential, by new nonlinear vertices due to modified gravity and, as in the pure dark matter case, by the effects of short-scale physics in the form of the divergence of an effective stress tensor. The effective fluid introduces a counterterm in the solution to the matter continuity and Euler equations, which allows a controlled expansion of clustering statistics on mildly nonlinear scales. We use this setup to compute the one-loop dark-matter power spectrum.

The Universe is not perfectly homogeneous, the large scale structure forms overdense regions and voids. In this paper, we consider the possibility that we occupy a special position in our Universe, close to the center of a local underdense region that we model as an LTB void embedded in a homogeneous and isotropic Universe. The CMB sky measured by an off-center observer in this void is not statistically isotropic. In addition to the non-stochastic CMB anisotropies due to the geometry of the model, we also observe a lensing-like distortion of the CMB anisotropies. In this article, we propose a framework to forecast the precision with which we can measure the amplitude of the lensing-like deformation of the CMB temperature anisotropies. For illustrative purposes, we apply this method to a couple of large-scale void models differing over the matter density profile and we show that the CMB temperature data from the Planck satellite is potentially capable to detect the effect for the large voids chosen here. A companion paper will be dedicated to a systematic exploration of different realistic void models.

We show that the anisotropies of the astrophysical stochastic gravitational wave background in the mHz band have a strong dependence on the modelling of galactic and sub-galactic physics. We explore a wide range of self-consistent astrophysical models for stellar evolution and for the distribution of orbital parameters, all calibrated such that they predict the same number of resolved mergers to fit the number of detections during LIGO/Virgo O1+O2 observations runs. We show that different physical choices for the process of black hole collapse and cut-off in the black hole mass distribution give fractional differences in the angular power spectrum of anisotropies up to 50\% on all angular scales. We also point out that the astrophysical information which can be extracted from anisotropies is complementary to the isotropic background and individual mergers. These results underline the interest in the anisotropies of the stochastic gravitational wave background as a new and potentially rich field of research, at the cross-road between astrophysics and cosmology.

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model.

This article explores the properties (amplitude and shape) of the angular power spectrum of the anisotropies of the astrophysical gravitational wave background (AGWB) focusing on the signatures of the astrophysical models describing sub-galactic physics. It demonstrates that while some parameters have negligible impact others, and in particular the stellar evolution models, the metallicity and the merger time delay distribution can result in relative differences of order 40% in the angular power spectrum of anisotropies in both the LIGO/Virgo and LISA frequency bands. It is also shown that the monopole and the anisotropic components of the AGWB are complementary and sensitive to different astrophysical parameters. It follows that AGWB anisotropies are a new observable with the potential to provide new astrophysical information that can not be accessed otherwise.

This Astro-2020 White Paper deals with what we might learn from future gravitational wave observations about the early universe phase transitions and their energy scales, primordial black holes, Hubble parameter, dark matter and dark energy, modified theories of gravity and extra dimensions.

We study the cosmological evolution of some nonlocal gravity models, when the initial conditions are set during a phase of primordial inflation. We examine in particular three models, the so-called RT, RR and $\Delta_4$ models, previously introduced by our group. We find that the RR and $\Delta_4$ models have a stable evolution also during inflation. The RT model has an apparent instability, but we show that, because of the smallness of the scale associated to the nonlocal term compared to the inflationary scale, this instability is innocuous and also the RT model has a viable evolution even when its initial conditions are set during a phase of primordial inflation.

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 geometric optics approximation traditionally used to study the propagation of gravitational waves on a curved background, breaks down in the vicinity of compact and extended astrophysical objects, where wave-like effects like diffusion and generation of polarization occur. We provide a framework to study the generation of polarization of a stochastic background of gravitational waves propagating in an inhomogeneous universe. The framework is general and can be applied to both cosmological and astrophysical gravitational wave backgrounds in any frequency range. We derive an order of magnitude estimate of the amount of polarization generated for cosmological and astrophysical backgrounds, in the frequency range covered by present and planned gravitational wave experiments. For an astrophysical background in the PTA and LISA band, the amount of polarization generated is suppressed by a factor $10^{-4}$ $(10^{-5})$ with respect to anisotropies. For a cosmological background we get an additional $10^{-2}$ suppression. We speculate on using our approach to map the distribution of (unresolvable) structures in the Universe.

Unresolved and resolved sources of gravitational waves are at the origin of a stochastic gravitational wave background. While the computation of its mean density as a function of frequency in a homogeneous and isotropic universe is standard lore, the computation of its anisotropies requires to understand the coarse graining from local systems, to galactic scales and then to cosmology. An expression of the gravitational wave energy density valid in any general spacetime is derived. It is then specialized to a perturbed Friedmann-Lema\^itre spacetime in order to determine the angular power spectrum of this stochastic background as well as its correlation with other cosmological probes, such as the galaxy number counts and weak lensing. Our result for the angular power spectrum also provides an expression for the variance of the gravitational wave background.

We study the impact on the galaxy correlation function of the presence of a vector component in the tracers' peculiar velocities, in the case in which statistical isotropy is violated. We present a general framework - based on the bipolar spherical harmonics expansion - to study this effect in a model independent way, without any hypothesis on the origin or the properties of these vector modes. We construct six new observables, that can be directly measured in galaxy catalogs in addition to the standard monopole, quadrupole and hexadecapole, and we show that they completely describe any deviations from isotropy. We then perform a Fisher analysis in order to quantify the constraining power of future galaxy surveys. As an example, we show that the SKA2 would be able to detect anisotropic rotational velocities with amplitudes as low as 1% of that of the vorticity generated during shell-crossing in standard dark matter scenarios.

We present the first predictions for the angular power spectrum of the astrophysical gravitational wave background constituted of the radiation emitted by all resolved and unresolved astrophysical sources. Its shape and amplitude depend on both the astrophysical properties on galactic scales and on cosmological properties. We show that the angular power spectrum behaves as $C_{\ell}\propto 1/{\ell}$ on large scales and that relative fluctuations of the signal are of order 30% at 100 Hz. We also present the correlations of the astrophysical gravitational wave background with weak-lensing and galaxy distribution. These numerical results pave the way to the study of a new observable at the crossroad between general relativity, astrophysics and cosmology.

The gravitational wave radiation emitted by all, resolved and unresolved, astrophysical sources in the observable universe generates a stochastic background. This background has a directional dependence inherited from the inhomogeneities of the matter distribution. This article proposes a new and independent derivation of the angular dependence of its energy density by focusing on the total gravitational wave signal produced by an ensemble of incoherent sources. This approach clarifies the origin of the angular correlation and the relation between the gravitational wave signal that can be measured by interferometers and the energy density of the stochastic background.

We compute the generation of vorticity from velocity dispersion in the dark matter fluid. For dark matter at zero temperature Helmholtz's theorem dictates that no vorticity is generated and we therefore allow the dark matter fluid to have a non-vanishing velocity dispersion. This implies a modification to the usual hydrodynamical system (continuity and Euler equations): we have to consider the Boltzmann hierarchy up to the second moment. This means that the Euler equation is modified with a source term that describes the effect of non-zero velocity dispersion. We write an equation for the Eulerian vorticity in Lagrangian coordinates and show that it has a growing mode already at second order in perturbation theory. We compute the power spectrum of the vorticity and the rotational velocity at second order in perturbation theory.

The properties of the cosmic microwave background (CMB) temperature and polarisation anisotropies measured by a static, off-centered observer located in a local spherically symmetric void, are described. In particular in this paper we compute, together with the standard 2- point angular correlation functions, the off-diagonal correlators, which are no more vanishing by symmetry. While the energy shift induced by the off-centered position of the observer can be suppressed by a proper choice of the observer velocity, a lensing-like effect on the CMB emission point remains. This latter effect is genuinely geometrical (e.g. non-degenerate with a boost) and reflects in the structure of the off-diagonal correlators. At lowest order in this effect, the temperature and polarisation correlation matrices have non-vanishing diagonal elements, as usual, and all the off-diagonal terms are excited. This particular signature of a local void model allows one, in principle, to disentangle geometrical effects from local kinematical ones in CMB observations.

We discuss a scale-free model of bigravity, in which the mass parameter of the standard bigravity potential is promoted to a dynamical scalar field. This modification retains the ghost-free bigravity structure, in particular it remains free of the Boulware-Deser ghost. We investigate the theory's interaction structure, focusing on its consistent scaling limits and strong coupling scales. Furthermore we explore the model's quadratic action, both around generic background configurations and paying special attention to cosmological backgrounds and to the associated background evolution. Finally we consider the possibility of realizing a phase of late-time acceleration as well as a quasi-de Sitter inflationary stage at early times, when the promoted "mass scalar" becomes the inflaton.

We introduce a new formalism to study perturbations of Hassan-Rosen bigravity theory, around general backgrounds for the two dynamical metrics. In particular, we derive the general expression for the mass term of the perturbations and we explicitly compute it for cosmological settings. We study tensor perturbations in a specific branch of bigravity using this formalism. We show that the tensor sector is affected by a late-time instability, which sets in when the mass matrix is no longer positive definite.

We study two nonlinear extensions of the nonlocal $R\,\Box^{-2}R$ gravity theory. We extend this theory in two different ways suggested by conformal symmetry, either replacing $\Box^{-2}$ with $(-\Box + R/6)^{-2}$, which is the operator that enters the action for a conformally-coupled scalar field, or replacing $\Box^{-2}$ with the inverse of the Paneitz operator, which is a four-derivative operator that enters in the effective action induced by the conformal anomaly. We show that the former modification gives an interesting and viable cosmological model, with a dark energy equation of state today $w_{\rm DE}\simeq -1.01$, which very closely mimics $\Lambda$CDM and evolves asymptotically into a de Sitter solution. The model based on the Paneitz operator seems instead excluded by the comparison with observations. We also review some issues about the causality of nonlocal theories, and we point out that these nonlocal models can be modified so to nicely interpolate between Starobinski inflation in the primordial universe and accelerated expansion in the recent epoch.

Recent work has shown that modifications of General Relativity based on the addition to the action of a non-local term $R\,\Box^{-2}R$, or on the addition to the equations of motion of a term involving $(g_{\mu\nu}\Box^{-1} R)$, produce dynamical models of dark energy which are cosmologically viable both at the background level and at the level of cosmological perturbations. We explore a more general class of models based on the addition to the action of terms proportional to $R_{\mu\nu}\,\Box^{-2}R^{\mu\nu}$ and $C_{\mu\nu\rho\sigma}\, \Box^{-2}C^{\mu\nu\rho\sigma}$, where $C_{\mu\nu\rho\sigma}$ is the Weyl tensor. We find that the term $R_{\mu\nu}\,\Box^{-2}R^{\mu\nu}$ does not give a viable background evolution. The non-local Weyl-square term, in contrast, does not contribute to the background evolution but we find that, at the level of cosmological perturbations, it gives instabilities in the tensor sector. Thus, only non-local terms which depend just on the Ricci scalar $R$ appear to be cosmologically viable. We discuss how these results can provide a hint for the mechanism that might generate these effective non-local terms from a fundamental local theory.

In this paper we study the generation of primordial perturbations in a cosmological setting of bigravity during inflation. We consider a model of bigravity which can reproduce the $\Lambda$CDM background and large scale structure and a simple model of inflation with a single scalar field and a quadratic potential. Reheating is implemented with a toy-model in which the energy density of the inflaton is entirely dissipated into radiation. We present analytic and numerical results for the evolution of primordial perturbations in this cosmological setting. We find that the amplitude of tensor perturbations generated during inflation is sufficiently suppressed to avoid the effects of the tensor instability discovered in Refs.[1,2] which develops during the cosmological evolution in the physical sector. We argue that from a pure analysis of the tensor perturbations this bigravity model is compatible with present observations. However, we derive rather stringent limits on inflation from the vector and scalar sectors.

In this paper we study gravitational wave perturbations in a cosmological setting of bigravity which can reproduce the {\Lambda}CDM background and large scale structure. We show that in general gravitational wave perturbations are unstable and only for very fine tuned initial conditions such a cosmology is viable. We quantify this fine tuning. We argue that similar fine tuning is also required in the scalar sector in order to prevent the tensor instability to be induced by second order scalar perturbations. Finally, we show that due to this power law instability, models of bigravity can lead to a large tensor to scalar ratio even for low scale inflation.