We consider gravitational wave (GW) sources with an associated electromagnetic (EM) counterpart, and analyze the time delay between both signals in the presence of lensing. If GWs have wavelengths comparable to the Schwarzschild radius of astrophysical lenses, they must be treated with wave optics, whereas EM waves are typically well within the approximation of geometric optics. With concrete examples, we confirm that the GW signal never arrives before its EM counterpart, if both are emitted at the same time. However, during the inspiral of a binary, peaks of the GW waveform can arrive before their EM counterpart. We stress this is only an apparent superluminality since the GW waveform is both distorted and further delayed with respect to light. In any case, measuring the multi-messenger time delay and correctly interpreting it has important implications for unveiling the distribution of lenses, testing the nature of gravity, and probing the cosmological expansion history.

Gravitational wave (GW) detections of binary black holes (BBHs) from the first two observing runs of advanced LIGO and Virgo have shown evidence for a maximum component black hole mass of $\lesssim50 M_\odot$. This is consistent with expectations of a gap in the mass spectrum of black holes due to the existence of pair-instability supernovae (PISN). This upper mass gap is expected to exist in the range $\sim50M_\odot$-$120M_\odot$, beyond which a population of intermediate mass black holes (IMBHs) may arise. We argue that ground-based GW detectors are expected to be sensitive to BBHs with masses above the mass gap. In the absence of detections, two years at upgraded sensitivity (A+) would constrain the local merger rate of these "far side of the PISN gap" black holes to be lower than $0.01\,\text{yr}^{-1}\text{Gpc}^{-3}$. Alternatively, with a few tens of events it may be possible to constrain the minimum IMBH mass to the percent level, if the edge of the distribution is sufficiently sharp. These far side BBHs would also be observed by future third generation detectors such as the Einstein Telescope (ET) as well as space-based GW interferometers such as LISA, and may dominate the fraction of multi-band events. We show that by comparing the number of detections observed from ground and space it is possible to constrain the merger rate history. Moreover, we find that the location of the high end of the PISN mass gap leaves an imprint on the spectral shape of the stochastic background of unresolved binaries, and that this may be accessible with A+ sensitivity. Finally, we show that by exploiting the upper edge of the gap, these high-mass BBHs can be used as GW standard sirens to constrain the cosmic expansion at redshifts of $\sim0.4$, $0.8$, and $1.5$ with LISA, LIGO, and ET, respectively. These far-side black holes would be the most massive black holes LIGO could detect.

The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology.

Gravitational wave (GW) oscillations occur whenever there are additional tensor modes interacting with the perturbations of the metric coupled to matter. These extra modes can arise from new spin-2 fields (as in e.g. bigravity theories) or from non-trivial realisations of the cosmological principle induced by background vector fields with internal symmetries (e.g. Yang-Mills, gaugids or multi-Proca). We develop a general cosmological framework to study such novel features due to oscillations. The evolution of the two tensor modes is described by a linear system of coupled second order differential equations exhibiting friction, velocity, chirality and mass mixing. We follow appropriate schemes to obtain approximate solutions for the evolution of both modes and show the corresponding phenomenology for different mixings. Observational signatures include modulations of the wave-form, oscillations of the GW luminosity distance, anomalous GW speed and chirality. We discuss the prospects of observing these effects with present and future GW observatories such as LIGO/VIRGO and LISA.

The curvature perturbations produced during an early era of inflation are known to have quasi-Gaussian distribution functions close to their maximum, where they are well constrained by measurements of the cosmic microwave background anisotropies and of the large-scale structures. In contrast, the tails of these distributions are poorly known, although this part is the relevant one for rare, extreme objects such as primordial black holes. We show that these tails are highly non-Gaussian, and cannot be described with standard non-Gaussian expansions, that are designed to approximate the distributions close to their maximum only. Using the stochastic-$\delta N$ formalism, we develop a generic framework to compute the tails, which are found to have an exponential, rather than Gaussian, decay. These exponential tails are inevitable, and do not require any non-minimal feature as they simply result from the quantum diffusion of the inflaton field along its potential. We apply our formalism to a few relevant single-field, slow-roll inflationary potentials, where our analytical treatment is confirmed by comparison with numerical results. We discuss the implications for the expected abundance of primordial black holes in these models, and highlight that it can differ from standard results by several orders of magnitude. In particular, we find that potentials with an inflection point overproduce primordial black holes, unless slow roll is violated.

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.

Modifications of General Relativity leave their imprint both on the cosmic expansion history through a non-trivial dark energy equation of state, and on the evolution of cosmological perturbations in the scalar and in the tensor sectors. In particular, the modification in the tensor sector gives rise to a notion of gravitational-wave (GW) luminosity distance, different from the standard electromagnetic luminosity distance, that can be studied with standard sirens at GW detectors such as LISA or third-generation ground based experiments. We discuss the predictions for modified GW propagation from some of the best studied theories of modified gravity, such as Horndeski or the more general degenerate higher order scalar-tensor (DHOST) theories, non-local infrared modifications of gravity, bigravity theories and the corresponding phenomenon of GW oscillation, as well as theories with extra or varying dimensions. We show that modified GW propagation is a completely generic phenomenon in modified gravity. We then use a simple parametrization of the effect in terms of two parameters $(\Xi_0,n)$, that is shown to fit well the results from a large class of models, to study the prospects of observing modified GW propagation using supermassive black hole binaries as standard sirens with LISA. We construct mock source catalogs and perform detailed Markov Chain Monte Carlo studies of the likelihood obtained from LISA standard sirens alone, as well as by combining them with CMB, BAO and SNe data to reduce the degeneracies between cosmological parameters. We find that the combination of LISA with the other cosmological datasets allows one to measure the parameter $\Xi_0$ that characterizes modified GW propagation to the percent level accuracy, sufficient to test several modified gravity theories. [Abridged]

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

Gravitational waves (GWs) provide a new tool to probe the nature of dark energy (DE) and the fundamental properties of gravity. We review the different ways in which GWs can be used to test gravity and models for late-time cosmic acceleration. Lagrangian-based gravitational theories beyond general relativity (GR) are classified into those breaking fundamental assumptions, containing additional fields and massive graviton(s). In addition to Lagrangian based theories we present the effective theory of DE and the $\mu$-$\Sigma$ parametrization as general descriptions of cosmological gravity. Multi-messenger GW detections can be used to measure the cosmological expansion (standard sirens), providing an independent test of the DE equation of state and measuring the Hubble parameter. Several key tests of gravity involve the cosmological propagation of GWs, including anomalous GW speed, massive graviton excitations, Lorentz violating dispersion relation, modified GW luminosity distance and additional polarizations, which may also induce GW oscillations. We summarize present constraints and their impact on DE models, including those arising from the binary neutron star merger GW170817. Upgrades of LIGO-Virgo detectors to design sensitivity and the next generation facilities such as LISA or Einstein Telescope will significantly improve these constraints in the next two decades.

Primordial black-holes (PBH) can be produced in single-field models of inflation with a quasi-inflection point in the potential. In these models, a large production of PBHs requires a deviation from the slow-roll (SR) trajectory. In turn, this SR violation can produce an exponential growth of quantum fluctuations. We study the back-reaction of these quantum modes on the inflationary dynamics using stochastic inflation in the Hamilton-Jacobi formalism. We develop a methodology to solve quantum diffusion beyond SR in terms of the statistical moments of the probability distribution. We apply these techniques to a toy model potential with a quasi-inflection point. We find that there is an enhancement of the power spectrum due to the dominance of the stochastic noise in the phase beyond SR. Moreover, non-Gaussian corrections become as well relevant with a large positive kurtosis. Altogether, this produces a significant boost of PBH production. We discuss how our results extend to other single-field models with similar dynamics. We conclude that the abundance of PBHs in this class of models should be revisited including quantum diffusion.

Primordial Black Holes (PBH) arise naturally from high peaks in the curvature power spectrum of near-inflection-point single-field inflation, and could constitute today the dominant component of the dark matter in the universe. In this letter we explore the possibility that a broad spectrum of PBH is formed in models of Critical Higgs Inflation (CHI), where the near-inflection point is related to the critical value of the RGE running of both the Higgs self-coupling $\lambda(\mu)$ and its non-minimal coupling to gravity $\xi(\mu)$. We show that, for a wide range of model parameters, a half-domed-shaped peak in the matter spectrum arises at sufficiently small scales that it passes all the constraints from large scale structure observations. The predicted cosmic microwave background spectrum at large scales is in agreement with Planck 2015 data, and has a relatively large tensor-to-scalar ratio that may soon be detected by B-mode polarization experiments. Moreover, the wide peak in the power spectrum gives an approximately lognormal PBH distribution in the range of masses $0.01 - 100\,M_\odot$, which could explain the LIGO merger events, while passing all present PBH observational constraints. The stochastic background of gravitational waves coming from the unresolved black-hole-binary mergers could also be detected by LISA or PTA. Furthermore, the parameters of the CHI model are consistent, within $2\sigma$, with the measured Higgs parameters at the LHC and their running. Future measurements of the PBH mass spectrum could allow us to obtain complementary information about the Higgs couplings at energies well above the EW scale, and thus constrain new physics beyond the Standard Model.

Multi-messenger gravitational wave (GW) astronomy has commenced with the detection of the binary neutron star merger GW170817 and its associated electromagnetic counterparts. The almost coincident observation of both signals places an exquisite bound on the GW speed $|c_g/c-1|\leq5\cdot10^{-16}$. We use this result to probe the nature of dark energy (DE), showing that a large class of scalar-tensor theories and DE models are highly disfavored. As an example we consider the covariant Galileon, a cosmologically viable, well motivated gravity theory which predicts a variable GW speed at low redshift. Our results eliminate any late-universe application of these models, as well as their Horndeski and most of their beyond Horndeski generalizations. Three alternatives (and their combinations) emerge as the only possible scalar-tensor DE models: 1) restricting Horndeski's action to its simplest terms, 2) applying a conformal transformation which preserves the causal structure and 3) compensating the different terms that modify the GW speed (to be robust, the compensation has to be independent on the background on which GWs propagate). Our conclusions extend to any other gravity theory predicting varying $c_g$ such as Einstein-Aether, Ho\v{r}ava gravity, Generalized Proca, TeVeS and other MOND-like gravities.

The direct detection of gravitational waves (GWs) is an invaluable new tool to probe gravity and the nature of cosmic acceleration. A large class of scalar-tensor theories predict that GWs propagate with velocity different than the speed of light, a difference that can be $\mathcal{O}(1)$ for many models of dark energy. We determine the conditions behind the anomalous GW speed, namely that the scalar field spontaneously breaks Lorentz invariance and couples to the metric perturbations via the Weyl tensor. If these conditions are realized in nature, the delay between GW and electromagnetic (EM) signals from distant events will run beyond human timescales, making it impossible to measure the speed of GWs using neutron star mergers or other violent events. We present a robust strategy to exclude or confirm an anomalous speed of GWs using eclipsing binary systems, whose EM phase can be exquisitely determined. he white dwarf binary J0651+2844 is a known example of such system that can be used to probe deviations in the GW speed as small as $c_g/c-1\gtrsim 2\cdot 10^{-12}$ when LISA comes online. This test will either eliminate many contender models for cosmic acceleration or wreck a fundamental pillar of general relativity.

We study the role of field redefinitions in general scalar-tensor theories. In particular, we first focus on the class of field redefinitions linear in the spin-2 field and involving derivatives of the spin-0 mode, generically known as disformal transformations. We start by defining the action of a disformal transformation in the tangent space. Then, we take advantage of the great economy of means of the language of differential forms to compute the full transformation of Horndeski's theory under general disformal transformations. We obtain that Horndeski's action maps onto itself modulo a reduced set of non-Horndeski Lagrangians. These new Lagrangians are found to be invariant under disformal transformation that depend only in the first derivatives of the scalar. Moreover, these combinations of Lagrangians precisely appear when expressing in our basis the constraints of the recently proposed Extended Scalar-Tensor (EST) theories. These results allow us to classify the different orbits of scalar-tensor theories invariant under particular disformal transformations, namely the special disformal, kinetic disformal and disformal Horndeski orbits. In addition, we consider generalizations of this framework. We find that there are possible well-defined extended disformal transformations that have not been considered in the literature. However, they generically cannot link Horndeski theory with EST theories. Finally, we study further generalizations in which extra fields with different spin are included. These field redefinitions can be used to connect different gravity theories such as multi-scalar-tensor theories, generalized Proca theories and bi-gravity. We discuss how the formalism of differential forms could be useful for future developments in these lines.

We use a description based on differential forms to systematically explore the space of scalar-tensor theories of gravity. Within this formalism, we propose a basis for the scalar sector at the lowest order in derivatives of the field and in any number of dimensions. This minimal basis is used to construct a finite and closed set of Lagrangians describing general scalar-tensor theories invariant under Local Lorentz Transformations in a pseudo-Riemannian manifold, which contains ten physically distinct elements in four spacetime dimensions. Subsequently, we compute their corresponding equations of motion and find which combinations are at most second order in derivatives in four as well as arbitrary number of dimensions. By studying the possible exact forms (total derivatives) and algebraic relations between the basis components, we discover that there are only four Lagrangian combinations producing second order equations, which can be associated with Horndeski's theory. In this process, we identify a new second order Lagrangian, named kinetic Gauss-Bonnet, that was not previously considered in the literature. However, we show that its dynamics is already contained in Horndeski's theory. Finally, we provide a full classification of the relations between different second order theories. This allows us to clarify, for instance, the connection between different covariantizations of Galileons theory. In conclusion, our formulation affords great computational simplicity with a systematic structure. As a first step we focus on theories with second order equations of motion. However, this new formalism aims to facilitate advances towards unveiling the most general scalar-tensor theories.