Under the assumption that a dynamical scalar field is responsible for the current acceleration of the Universe, we explore the possibility of probing its physics in black hole merger processes with gravitational wave interferometers. Remaining agnostic about the microscopic physics, we use an effective field theory approach to describe the scalar dynamics. We investigate the case in which some of the higher derivative operators, that are highly suppressed on cosmological scales, instead become important on typical distances for black holes. If a coupling to the Gauss-Bonnet operator is one of them, a non-trivial background profile for the scalar field can be sourced in the surroundings of the black hole, resulting in a potentially large amount of "hair". In turn, this can induce sizeable modifications to the spacetime geometry or a mixing between the scalar and the gravitational perturbations. Both effects will ultimately translate into a modification of the quasi-normal mode spectrum in a way that is also sensitive to other operators besides the one sourcing the scalar background. The presence of deviations from the predictions of general relativity in the observed spectrum can therefore serve as a window onto dark energy physics.

With recent constraints on the propagation speed of gravitational waves, the class of scalar-tensor theories has significantly been reduced. We consider one of the surviving models still relevant for cosmology and investigate its radiative stability. The model contains operators with explicit breaking of the Galileon symmetry and we study whether they harm the re-organization of the effective field theory. Within the regime of validity we establish a non-renormalization theorem and show explicitly that the quantum corrections, to one-loop, do not detune the classical Lagrangian generating suppressed counterterms. This is striking since the non-renormalization theorem is established in the presence of a genuine Galileon symmetry breaking term.

Gravitational wave (GW) constraints have recently been used to significantly restrict models of dark energy and modified gravity. New bounds arising from GW decay and GW-induced dark energy instabilities are particularly powerful in this context, complementing bounds from the observed speed of GWs. We discuss the associated linear cosmology for Horndeski gravity models surviving these combined bounds and compute the corresponding cosmological parameter constraints, using CMB, redshift space distortion, matter power spectrum and BAO measurements from the Planck, SDSS/BOSS and 6dF surveys. The surviving theories are strongly constrained, tightening previous bounds on cosmological deviations from $\Lambda{}$CDM by over an order of magnitude. We also comment on general cosmological stability constraints and the nature of screening for the surviving theories, pointing out that a raised strong coupling scale can ensure compatibility with gravitational wave constraints, while maintaining a functional Vainshtein screening mechanism on solar system scales. Finally, we discuss the quasi-static limit as well as (constraints on) related observables for near-future surveys.

Positivity bounds - the consequences of requiring a unitary, causal, local UV completion - place strong restrictions on theories of dark energy and/or modified gravity. We derive and investigate such bounds for Horndeski scalar-tensor theories and for the first time pair these bounds with a cosmological parameter estimation analysis, using CMB, redshift space distortion, matter power spectrum and BAO measurements from the Planck, SDSS/BOSS and 6dF surveys. Using positivity bounds as theoretical priors, we show that their inclusion in the parameter estimation significantly improves the constraints on dark energy/modified gravity parameters. Considering as an example a specific class of models, which are particularly well-suited to illustrate the constraining power of positivity bounds, we find that these bounds eliminate over 60% of the previously allowed parameter space. We also discuss how combining positivity requirements with additional theoretical priors has the potential to further tighten these constraints: for instance also requiring a subluminal speed of gravitational waves eliminates all but 1% of the previously allowed parameter space.

Weyl (scale) invariant theories of scalars and gravity can generate all mass scales spontaneously. In this paper we study a particularly simple version -- scale invariant $R^2$ gravity -- and show that, during an inflationary period, it leads to fluctuations which, for a particular parameter choice, are almost indistinguishable from normal $R^2$ inflation. Current observations place tight constraints on the primary coupling constant of this theory and predict a tensor to scalar ratio, $0.0033 > r > 0.0026$, which is testable with the next generation of cosmic microwave background experiments.

We present new cosmological parameter constraints for general Horndeski scalar-tensor theories, using CMB, redshift space distortion, matter power spectrum and BAO measurements from the Planck, SDSS/BOSS and 6dF surveys. We focus on theories with cosmological gravitational waves propagating at the speed of light, $c_{\rm GW} = c$, implementing and discussing several previously unaccounted for aspects in the constraint derivation for such theories, that qualitatively affect the resulting constraints. In order to ensure our conclusions are robust, we compare results for three different parametrisations of the free functions in Horndeski scalar-tensor theories, identifying several parametrisation-independent features of the constraints. We also consider models, where $c_{\rm GW} \neq c$ in cosmological settings (still allowed after GW170817 for frequency-dependent $c_{\rm GW}$) and show how this affects cosmological parameter constraints.

Radiative stability places strong constraints on general dark energy and modified gravity theories. We consider Horndeski scalar-tensor theories with luminally propagating gravitational waves (as extensively discussed in the wake of GW170817) and show that generically there is a tension between obtaining observable deviations from General Relativity (GR) in cosmology and the requirement of radiative stability. Using this as a constraint, we discuss the subsets of theories that are capable of yielding observable, radiatively stable departures from GR. A key consequence are significantly tightened cosmological parameter constraints on dark energy and modified gravity parameters, which we explicitly compute using data from the Planck, SDSS/BOSS and 6dF surveys.

A scale-invariant universe can have a period of accelerated expansion at early times: inflation. We use a frame-invariant approach to calculate inflationary observables in a scale invariant theory of gravity involving two scalar fields - the spectral indices, the tensor to scalar ratio, the level of isocurvature modes and non-Gaussianity. We show that scale symmetry leads to an exact cancellation of isocurvature modes and that, in the scale-symmetry broken phase, this theory is well described by a single scalar field theory. We find the predictions of this theory strongly compatible with current observations.

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 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 analyse the late time cosmology and the gravitational properties of doubly coupled bigravity in the constrained vielbein formalism (equivalent to the metric formalism) when the mass of the massive graviton is of the order of the present Hubble rate. We focus on one of the two branches of background cosmology where the ratio between the scale factors of the two metrics is algebraically determined. The Universe evolves from a matter dominated epoch to a dark energy dominated era where the equation of state of dark energy can always be made close to -1 now by appropriately tuning the graviton mass. We also analyse the perturbative spectrum of the theory in the quasi static approximation well below the strong coupling scale where no instability is present and we show that there are five scalar degrees of freedom, two vectors and two gravitons. In a cosmological FRW background for both metrics, four of the five scalars are Newtonian potentials which lead to a modification of gravity on large scales. In this scalar sector, gravity is modified with effects on both the growth of structure and the lensing potential. In particular, we find that the $\Sigma$ parameter governing the Poisson equation of the weak lensing potential can differ from one in the recent past of the Universe. Overall, the nature of the modification of gravity at low energy, which reveals itself in the growth of structure and the lensing potential, is intrinsically dependent on the couplings to matter and the potential term of the vielbeins. We also find that the time variation of Newton's constant in the Jordan frame can easily satisfy the bound from solar system tests of gravity. Finally we show that the two gravitons present in the spectrum have a non-trivial mass matrix whose origin follows from the potential term of bigravity. This mixing leads to gravitational birefringence.

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.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

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.

In this paper we consider how the strong-coupling scale, or perturbative cutoff, in a multi-gravity theory depends upon the presence and structure of interactions between the different fields. This can elegantly be rephrased in terms of the size and structure of the `theory graph' which depicts the interactions in a given theory. We show that the question can be answered in terms of the properties of various graph-theoretical matrices, affording an efficient way to estimate and place bounds on the strong-coupling scale of a given theory. In light of this we also consider the problem of relating a given theory graph to a discretised higher dimensional theory, a la dimensional deconstruction.

Ghost-free massive gravity models generically have a strong coupling scale of $\Lambda_3 =(M_{\rm Pl} m^2)^{1/3}$. However, for one of these models - `minimal massive gravity' - it is not clear what this scale is in the subset of solutions with vanishing vector modes, since there are then no interactions at the scale $\Lambda_3$. We show that there are always scalar-tensor interactions at a scale arbitrarily close to $\Lambda_3$ around the Minkowski vacuum solution. This explicitly confirms and completes previous research establishing that $\Lambda_3$ effectively is the maximal strong coupling scale for all ghost-free (dRGT) massive gravity models (on Minkowski). In the process, we also revisit and clarify the construction of generic Lorentz-invariant massive gravity models, explicitly constructing an improved non-redundant expansion for these models.

We investigate matter couplings in massive bigravity. We find a new family of such consistent couplings, including and extending known consistent matter couplings, and we investigate their decoupling limits, ADM decompositions, Higuchi bounds and further aspects. We show that differences to previous known consistent couplings only arise beyond the $\Lambda_3$ decoupling limit and discuss the uniqueness of consistent matter couplings and how this is related to the so-called symmetric vielbein condition. Since we work in a vielbein formulation, these results easily generalise to multi-gravity.

We investigate a large class of infinitesimal, but fully nonlinear in the field, transformations of the Galileon and search for extended symmetries. The transformations involve powers of the coordinates $x$ and the field $\pi$ up to any finite order $N$. Up to quadratic order the structure of these symmetry transformations is the unique generalisation of both the infinitesimal version of the standard Galileon shift symmetry as well as a recently discovered infinitesimal extension of this symmetry. The only higher-order extensions of this symmetry we recover are (`Galileon dual' versions of) symmetries of the standard kinetic term.

The only known fully ghost-free and consistent Lorentz-invariant kinetic term for a graviton (or indeed for any spin-2 field) is the Einstein-Hilbert term. Here we propose and investigate a new family of candidate kinetic interactions and their extensions to derivative interactions involving several spin-2 fields. These new terms generically break diffeomorphism invariance(s) and as a result can lead to the propagation of 5 degrees of freedom for a single spin-2 field - analogous to ghost-free Massive Gravity. We discuss under what circumstances these new terms can be used to build healthy effective field theories and in the process establish the `Jordan' and `Einstein' frame pictures for Massive-, Bi- and Multi-Gravity.

In this paper we investigate the decoupling limit of a particular class of multi-gravity theories, i.e. of theories of interacting spin-2 fields. We explicitly compute the interactions of helicity-0 modes in this limit, showing that they take on the form of multi-Galileons and dual forms. In the process we extend the recently discovered Galileon dualities, deriving a set of new multi-Galileon dualities. These are also intrinsically connected to healthy, but higher-derivative, multi-scalar field theories akin to `beyond Horndeski' models.

In this paper we construct a family of ways in which matter can couple to one or more `metrics'/spin-2 fields in the vielbein formulation. We do so subject to requiring the weak equivalence principle and the absence of ghosts from pure spin-2 interactions generated by the matter action. Results are presented for Massive, Bi- and Multi-Gravity theories and we give explicit expressions for the effective matter metric in all of these cases.

We revisit and extend the `Effective field theory for massive gravitons' constructed by Arkani-Hamed, Georgi and Schwartz in the light of recent progress in constructing ghost-free theories with multiple interacting spin-2 fields. We show that there exist several dual ways of restoring gauge invariance in such multi-gravity theories, find a generalised Fierz-Pauli tuning condition relevant in this context and highlight subtleties in demixing tensor and scalar modes. The generic multi-gravity feature of scalar mixing and its consequences for higher order interactions are discussed. In particular we show how the decoupling limit is qualitatively changed in theories of interacting spin-2 fields. We relate this to dRGT (de Rham, Gabadadze, Tolley) massive gravity, Hassan-Rosen bigravity and the multi-gravity constructions by Hinterbichler and Rosen. As an additional application we show that EBI (Eddington-Born-Infeld) bigravity and higher order generalisations thereof possess ghost-like instabilities.

Relativistic scalar fields are ubiquitous in modified theories of gravity. An important tool in understanding their impact on structure formation, especially in the context of N-body simulations, is the quasi-static approximation in which the time evolution of perturbations in the scalar fields is discarded. We show that this approximation must be used with some care by studying linearly perturbed scalar field cosmologies and quantifying the errors that arise from taking the quasi-static limit. We focus on f(R) and chameleon models and link the accuracy of the quasi-static approximation to the fast/slow-roll behaviour of the background and its proximity to {\Lambda}CDM. Investigating a large range of scales, from super- to sub-horizon, we find that slow-rolling ({\Lambda}CDM-like) backgrounds generically result in good quasi-static behaviour, even on (super-)horizon scales. We also discuss how the approximation might affect studying the non-linear growth of structure in numerical N-body simulations.

We calculate the perturbed action, at second and third order, for a massive three-form field minimally coupled to gravity, and use it to explore the observational predictions of three-form inflation. One intriguing result is that the value of the spectral index is nearly independent of the three-form potential, being fixed solely by the number of e-folds of inflation, with n_s=0.97 for the canonical number of 60. Considering the bispectrum, we employ standard techniques to give explicit results for two models, one of which produces a large non-Gaussianity. Finally, we confirm our results by employing a duality relating the three-form theory to a non-canonical scalar field theory and explicitly re-computing results in this dual picture.

We present constraints on how far single field inflation may depart from the familiar slow-roll paradigm. Considering a fast-roll regime while requiring a (near)-scale-invariant power spectrum introduces large self-interactions for the field and consequently large and scale-dependent non-Gaussianities. Employing this signal, we use the requirement of weak-coupling together with WMAP constraints to derive bounds on generic $P(X,\phi)$ theories of single field inflation.

We consider generalized chameleon models where the conformal coupling between matter and gravitational geometries is not only a function of the chameleon field \phi, but also of its derivatives via higher order co-ordinate invariants. Specifically we consider the first such non-trivial conformal factor A(\phi,X), where X is the canonical kinetic term for \phi. The associated phenomenology is investigated and we show that such theories have a new generic mass-altering mechanism, potentially assisting the generation of a sufficiently large chameleon mass in dense environments. The most general effective potential is derived for such derivative chameleon setups and explicit examples are given. Interestingly this points us to the existence of a purely derivative chameleon protected by a shift symmetry for \phi. We also discuss potential ghost-like instabilities associated with mass-lifting mechanisms and find another, mass-lowering and instability-free, branch of solutions. This suggests that, barring fine-tuning, stable derivative models are in fact typically anti-chameleons that suppress the field's mass in dense environments. Furthermore we investigate modifications to the thin-shell regime and prove a no-go theorem for chameleon effects in non-conformal geometries of the disformal type.

We investigate non-Gaussianity in general single field models without assuming slow-roll conditions or the exact scale-invariance of the scalar power spectrum. The models considered include general single field inflation (e.g. DBI and canonical inflation) as well as bimetric models. We compute the full non-Gaussian amplitude, its size fnl, its shape, and the running with scale n_{NG}. In doing so we show that observational constraints allow significant violations of slow roll conditions and we derive explicit bounds on slow-roll parameters for fast-roll single field scenarios. A variety of new observational signatures is found for models respecting these bounds. We also explicitly construct concrete model implementations giving rise to this new phenomenology.

The minimal bimetric theory employing a disformal transformation between matter and gravity metrics is known to produce exactly scale-invariant fluctuations. It has a purely equilateral non-Gaussian signal, with an amplitude smaller than that of DBI inflation (with opposite sign) but larger than standard inflation. We consider non-minimal bimetric models, where the coupling $B$ appearing in the disformal transformation ${\hat g}_{\mn}= g_{\mn} -B\partial_\mu\phi\partial_\nu\phi$ can run with $\phi$. For power-law $B(\phi)$ these models predict tilted spectra. For each value of the spectral index, a distinctive distortion to the equilateral property can be found. The constraint between this distortion and the spectral index can be seen as a "consistency relation" for non-minimal bimetric models.

We consider general mixtures of isocurvature and adiabatic cosmological perturbations. With a minimal assumption set consisting of the linearized Einstein equations and a primordial perfect fluid we derive the second-order action and its curvature variables. We also allow for varying equation of state and speed of sound profiles. The derivation is therefore carried out at the same level of generality that has been achieved for adiabatic modes before. As a result we find a new conserved super-horizon quantity and relate it to the adiabatically conserved curvature perturbation. Finally we demonstrate how the formalism can be applied by considering a Chaplygin gas-like primordial matter model, finding two scale-invariant solutions for structure formation.

We revisit the question of whether fluctuations in hydrodynamical, adiabatical matter could explain the observed structures in our Universe. We consider matter with variable equation of state $w=p_0/\ep_0$ and a concomitant (under the adiabatic assumption) density dependent speed of sound, $c_s$. We find a limited range of possibilities for a set up when modes start inside the Hubble radius, then leaving it and freezing out. For expanding Universes, power-law $w(\ep_0)$ models are ruled out (except when $c_s^2\propto w \ll 1$, requiring post-stretching the seeded fluctuations); but sharper profiles in $c_s$ do solve the horizon problem. Among these, a phase transition in $c_s$ is notable for leading to scale-invariant fluctuations if the initial conditions are thermal. For contracting Universes all power-law $w(\ep_0)$ solve the horizon problem, but only one leads to scale-invariance: $w\propto \ep_0^2$ and $c_s\propto \ep_0$. This model bypasses a number of problems with single scalar field cyclic models (for which $w$ is large but constant).