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.

The decay timescales of the quasinormal modes of a massive scalar field have an intriguing behavior: they either grow or decay with increasing angular harmonic numbers $\ell$, depending on whether the mass of the scalar field is small or large. We identify the properties of the effective potential of the scalar field that leads to this behavior and characterize it in detail. If the scalar field is non-minimally coupled, considered here, the scalar quasinormal modes will leak into the gravitational wave signal and will have decaying times that are comparable or smaller than those typical in General Relativity. Hence, these modes could be detectable in the future. Finally, we find that the anomalous behavior in the decay timescales of quasinormal modes is present in a much larger class of models beyond a simple massive scalar field.

We consider cosmological models where dark energy is described by a dynamical field equipped with the Chameleon screening mechanism, which serves to hide its effects in local dense regions and to conform to Solar System observations. In these models, there is no universal gravitational coupling and here we study the effective couplings that determine the force between massive objects, $G_N$, and the propagation of gravitational waves, $G_{gw}$. In particular, we revisit the Chameleon screening mechanism without neglecting the time dependence of the galactic environment where local regions are embedded in, and analyze the induced time evolution on $G_N$ and $G_{gw}$, which can be tested with Lunar Laser Ranging and direct gravitational waves observations. We explicitly show how and why these two couplings generically differ. We also find that due to the particular way the Chameleon screening mechanism works, their time evolutions are highly suppressed in the weak-field non-relativistic approximation.

In this paper we construct the first analytic examples of (3+1)-dimensional self-gravitating regular cosmic tube solutions which are superconducting, free of curvature singularities and with non trivial topological charge in the Einstein-SU(2) non-linear sigma-model. These gravitating topological solitons at large distance from the axis look like a (boosted) cosmic string with an angular defect given by the parameters of the theory, and near the axis the parameters of the solutions can be chosen so that the metric is singularity free and without angular defect. The curvature is concentrated on a tube around the axis. These solutions are similar to global strings but regular everywhere, and the non-linear sigma-model regularizes the global string in a similar way as a non-Abelian field regularizes the Dirac monopole. Also, these solutions can be promoted to self consistent solutions of the fully coupled Einstein Maxwell non-linear sigma-model in which the non-linear sigma model is minimally coupled both to the U(1) gauge field and to General Relativity. The analysis shows that these solutions behave as superconductors as they carry a persistent current even when the $U(1)$ field vanishes. Such persistent current cannot be continuously deformed to zero as it is tied to the topological charge of the solutions themselves. The peculiar features of the gravitational lensing of these gravitating solitons are shortly discussed.

Cosmological models with a dynamical dark energy field typically lead to a modified propagation of gravitational waves via an effectively time-varying gravitational coupling $G(t)$. The local variation of this coupling between the time of emission and detection can be probed with standard sirens. Here we discuss the role that Lunar Laser Ranging (LLR) and binary pulsar constraints play in the prospects of constraining $G(t)$ with standard sirens. In particular, we argue that LLR constrains the matter-matter gravitational coupling $G_N(t)$, whereas binary pulsars and standard sirens constrain the quadratic kinetic gravity self-interaction $G_{gw}(t)$. Generically, these two couplings could be different in alternative cosmological models, in which case LLR constraints are irrelevant for standard sirens. We use the Hulse-Taylor pulsar data and show that observations are highly insensitive to time variations of $G_{gw}(t)$ yet highly sensitive to $G_N(t)$. We thus conclude that future gravitational waves data will become the best probe to test $G_{gw}(t)$, and will hence provide novel constraints on dynamical dark energy models.

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 consider cosmological models with a dynamical dark energy field, and study the presence of three types of commonly found instabilities, namely ghost (when fields have negative kinetic energy), gradient (negative momentum squared) and tachyon (negative mass squared). In particular, we study the linear scalar perturbations of theories with two interacting scalar fields as a proxy for a dark energy and matter fields, and explicitly show how canonical transformations relate these three types of instabilities with each other. We generically show that low-energy ghosts are equivalent to tachyonic instabilities, and that high-energy ghosts are equivalent to gradient instabilities. Via examples we make evident the fact that whenever one of these fields exhibits an instability then the entire physical system becomes unstable, with an unbounded Hamiltonian. Finally, we discuss the role of interactions between the two fields, and show that whereas most of the time interactions will not determine whether an instability is present or not, they may affect the timescale of the instability. We also find exceptional cases in which the two fields are ghosts and hence the physical system is seemingly unstable, but the presence of interactions actually lead to stable solutions. These results are very important for assessing the viability of dark energy models that may exhibit ghost, gradient or tachyonic modes.

Through numerical simulations of a minimally coupled massive Klein-Gordon scalar field, we show that it is possible to grow hair on a Schwarzschild black hole if one assumes an initial periodically time-varying but spatially homogeneous scalar background. By "hair", we mean a non-trivial profile in the scalar field. We find that this profile emerges on a timescale related to the mass of the black hole, with features related to the mass of the scalar particle. We undertake simulations with and without backreaction on the metric and see that the essential, qualitative features remain consistent. We also contrast the results from higher mass scalars to the case of a low mass with a large Compton wavelength. The results are particularly relevant for scalar-tensor theories of gravity and dark matter models consisting of a massive scalar, e.g. axions.

We discuss the difference between various gauge-invariant quantities typically used in single-field inflation, namely synchronous $\zeta_s$, comoving $\zeta_c$, and unitary $\zeta_u$ curvatures. We show that conservation of $\zeta_c$ outside the horizon is quite restrictive on models as it leads to conservation of $\zeta_s$ and $\zeta_u$, whereas the reverse does not hold. We illustrate the consequence of these differences with two inflationary models: ultra-slow-roll (USR) and braiding-ultra-slow-roll (BUSR). In USR, we show that out of the three curvatures, only $\zeta_s$ is conserved outside the horizon, and we connect this result to the concepts of separate universe and the usage of the $\delta N$ formalism. We find that even though $\zeta_s$ is conserved, there is still a mild violation of the separate universe approximation in the continuity equation. Nevertheless, the $\delta N$ formalism can still be applied to calculate the primordial power spectrum of some gauge-invariant quantities such as $\zeta_u$, although it breaks down for others such as the uniform-density curvature. In BUSR, we show that both $\zeta_u$ and $\zeta_s$ are conserved outside the horizon, but take different values. Additionally, since $\zeta_u\not=\zeta_c$ we find that the prediction for observable curvature fluctuations after inflation does not reflect $\zeta_c$ at horizon crossing during inflation and moreover involves not just $\zeta_u$ at that epoch but also the manner in which the braiding phase ends.

We consider the effect of a time-varying Planck mass on the propagation of gravitational waves (GWs). A running Planck mass arises naturally in several modified gravity theories, and here we focus on those that carry an additional dark energy field responsible for the late-time accelerated expansion of the universe, yet--like general relativity (GR)--propagate only two GW polarizations, both traveling at the speed of light. Because a time-varying Planck mass affects the amplitude of the GWs and therefore the inferred distance to the source, standard siren measurements of $H_0$ are degenerate with the parameter $c_M$ characterizing the time-varying Planck mass, where $c_M=0$ corresponds to GR with a constant Planck mass. The effect of non-zero $c_M$ will have a noticeable impact on GWs emitted by binary neutron stars (BNSs) at the sensitivities and distances observable by ground-based GW detectors such as advanced LIGO and A+, implying that standard siren measurements can provide joint constraints on $H_0$ and $c_M$. Assuming a $\Lambda$CDM evolution of the universe and taking Planck's measurement of $H_0$ as a prior, we find that GW170817 constrains $c_M = -9^{+21}_{-28}$ ($68.3\%$ credibility). We also discuss forecasts, finding that if we assume $H_0$ is known independently, then 100 BNS events detected by advanced LIGO can constrain $c_M$ to within $\pm0.9$. This is comparable to the current best constraints from cosmology. Similarly, for 100 LIGO A+ BNS detections, it is possible to constrain $c_M$ to $\pm0.5$. When analyzing joint $H_0$ and $c_M$ constraints we find that $\sim 400$ LIGO A+ events are needed to constrain $H_0$ to $1\%$ accuracy. Finally, we discuss the possibility of a nonzero value of $c_M$ biasing standard siren $H_0$ measurements from 100 LIGO A+ detections, and find that $c_M=+1.35$ could bias $H_0$ by 3-4$\sigma$ too low if we incorrectly assume $c_M=0$.

The recent detection of GRB 170817A and GW170817 constrains the speed of gravity waves $c_T$ to be that of light, which severely restricts the landscape of modified gravity theories that impact the cosmological evolution of the universe. In this work, we investigate the presence of black hole hair in the remaining viable cosmological theories of modified gravity that respect the constraint $c_T=1$. We focus mainly on scalar-tensor theories of gravity, analyzing static, asymptotically flat black holes in Horndeski, Beyond Horndeski, Einstein-Scalar-Gauss-Bonnet, and Chern-Simons theories. We find that in all of the cases considered here, theories that respect $c_T=1$ do not allow for hair, or have negligible hair. We further comment on vector-tensor theories including Einstein Yang-Mills, Einstein-Aether, and Generalised Proca theories, as well as bimetric theories.

We use the covariant formulation proposed in Tattersall et al (2017) to analyse the structure of linear perturbations about a spherically symmetric background in different families of gravity theories, and hence study how quasi-normal modes of perturbed black holes may be affected by modifications to General Relativity. We restrict ourselves to single-tensor, scalar-tensor and vector-tensor diffeomorphism-invariant gravity models in a Schwarzschild black hole background. We show explicitly the full covariant form of the quadratic actions in such cases, which allow us to then analyse odd parity (axial) and even parity (polar) perturbations simultaneously in a straightforward manner.

We analyse cosmological perturbations around a homogeneous and isotropic background for scalar-tensor, vector-tensor and bimetric theories of gravity. Building on previous results, we propose a unified view of the effective parameters of all these theories. Based on this structure, we explore the viable space of parameters for each family of models by imposing the absence of ghosts and gradient instabilities. We then focus on the quasistatic regime and confirm that all these theories can be approximated by the phenomenological two-parameter model described by an effective Newton's constant and the gravitational slip. Within the quasistatic regime we pinpoint signatures which can distinguish between the broad classes of models (scalar-tensor, vector-tensor or bimetric). Finally, we present the equations of motion for our unified approach in such a way that they can be implemented in Einstein-Boltzmann solvers.

The detection of an electromagnetic counterpart (GRB 170817A) to the gravitational wave signal (GW170817) from the merger of two neutron stars opens a completely new arena for testing theories of gravity. We show that this measurement allows us to place stringent constraints on general scalar-tensor and vector-tensor theories, while allowing us to place an independent bound on the graviton mass in bimetric theories of gravity. These constraints severely reduce the viable range of cosmological models that have been proposed as alternatives to general relativistic cosmology.

We present a covariant formulation for constructing general quadratic actions for cosmological perturbations, invariant under a given set of gauge symmetries for a given field content. This approach allows us to analyse scalar, vector and tensor perturbations at the same time in a straightforward manner. We apply the procedure to diffeomorphism invariant single-tensor, scalar-tensor and vector-tensor theories and show explicitly the full covariant form of the quadratic actions in such cases, in addition to the actions determining the evolution of vector and tensor perturbations. We also discuss the role of the symmetry of the background in identifying the set of cosmologically relevant free parameters describing these classes of theories, including calculating the relevant free parameters for an axisymmetric Bianchi-I vacuum universe.

We implement the method developed in [1] to construct the most general parametrised action for linear cosmological perturbations of bimetric theories of gravity. Specifically, we consider perturbations around a homogeneous and isotropic background, and identify the complete form of the action invariant under diffeomorphism transformations, as well as the number of free parameters characterising this cosmological class of theories. We discuss, in detail, the case without derivative interactions, and compare our results with those found in massive bigravity.

We present a method for parametrizing linear cosmological perturbations of theories of gravity, around homogeneous and isotropic backgrounds. The method is sufficiently general and systematic that it can be applied to theories with any degrees of freedom (DoFs) and arbitrary gauge symmetries. In this paper, we focus on scalar-tensor and vector-tensor theories, invariant under linear coordinate transformations. In the case of scalar-tensor theories, we use our framework to recover the simple parametrizations of linearized Horndeski and "Beyond Horndeski" theories, and also find higher-derivative corrections. In the case of vector-tensor theories, we first construct the most general quadratic action for perturbations that leads to second-order equations of motion, which propagates two scalar DoFs. Then we specialize to the case in which the vector field is time-like (\`a la Einstein-Aether gravity), where the theory only propagates one scalar DoF. As a result, we identify the complete forms of the quadratic actions for perturbations, and the number of free parameters that need to be defined, to cosmologically characterize these two broad classes of theories.

In this paper we consider the issue of paradigm evaluation by applying Bayes' theorem along the following nested hierarchy of progressively more complex structures: i) parameter estimation (within a model), ii) model selection and comparison (within a paradigm), iii) paradigm evaluation. In such a hierarchy the Bayesian evidence works both as the posterior's normalization at a given level and as the likelihood function at the next level up. Whilst raising no objections to the standard application of the procedure at the two lowest levels, we argue that it should receive a considerable modification when evaluating paradigms, when testability and fitting data are equally important. By considering toy models we illustrate how models and paradigms that are difficult to falsify are always favoured by the Bayes factor. We argue that the evidence for a paradigm should not only be high for a given dataset, but exceptional with respect to what it would have been, had the data been different. With this motivation we propose a measure which we term predictivity, as well as a prior to be incorporated into the Bayesian framework, penalising unpredictivity as much as not fitting data. We apply this measure to inflation seen as a whole, and to a scenario where a specific inflationary model is hypothetically deemed as the only one viable as a result of information alien to cosmology (e.g. Solar System gravity experiments, or particle physics input). We conclude that cosmic inflation is currently hard to falsify, but that this could change were external/additional information to cosmology to select one of its many models. We also compare this state of affairs to bimetric varying speed of light cosmology.

We study a previously largely unexplored branch of homogeneous and isotropic background solutions in ghost-free massive bigravity with consistent double matter coupling. For a certain family of parameters we find `self-inflated' FLRW cosmologies, i.e. solutions with an accelerated early-time period during the radiation-dominated era. In addition, these solutions also display an accelerated late-time period closely mimicking GR with a cosmological constant. Interestingly, within this family, the particular case of $\beta_1=\beta_3=0$ gives bouncing cosmologies, where there is an infinite contracting past, a non-zero minimum value of the scale factor at the bounce, and an infinite expanding future.

We present a comprehensive analysis of classical scalar, vector and tensor cosmological perturbations in ghost-free massive bigravity. In particular, we find the full evolution equations and analytical solutions in a wide range of regimes. We show that there are viable cosmological backgrounds but, as has been found in the literature, these models generally have exponential instabilities in linear perturbation theory. However, it is possible to find stable scalar cosmological perturbations for a very particular choice of parameters. For this stable subclass of models we find that vector and tensor perturbations have growing solutions. We argue that special initial conditions are needed for tensor modes in order to have a viable model.

We propose and develop a general algorithm for finding the action for cosmological perturbations which rivals the conventional, gauge-invariant approach and can be applied to theories with more than one metric. We then apply it to a particular case of bigravity, focusing on the Eddington- inspired Born-Infeld theory, and show that we can obtain a nearly scale-invariant power spectrum for both scalar and tensor primordial quantum perturbations. Unfortunately, in the case of the minimal Eddington-inspired Born-Infeld theory, we find that the tensor-to-scalar ratio of perturbations is unacceptably large. We discuss the applicability of our general method and the possibility of resurrecting the specific theory we have looked at.

We explore the observational implications of models of intermediate inflation driven by modified dispersion relations, specifically those representing the phenomenon of dimensional reduction in the ultraviolet limit. These models are distinct from the standard ones because they do not require violations of the strong energy condition, and this is reflected in their structure formation properties. We find that they can naturally accommodate deviations from exact scale-invariance. They also make clear predictions for the running of the spectral index and tensor modes, rendering the models straightforwardly falsifiable. We discuss the observational prospects for these models and the implications these may have for quantum gravity scenarios.