We revisit the constraints on scalar tensor theories of modified gravity following the purge of GW170817. We pay particular attention to dynamical loopholes where the anomalous speed of propagation of the gravitational wave can vanish on-shell, when we impose the dynamical field equations. By working in the effective field theory formalism we are able to improve on previous analyses, scanning a much wider class of theories, including Beyond Horndeski and DHOST. Furthermore, the formalism is well adapted to consider the effect of inhomogeneous perturbations, where we explicitly take into account the fact that the galactic overdensities are pressureless to leading order.

We aim to provide spectroscopic evidence regarding the nature of the putative star cluster in Eridanus 2 and to place constraints on the mass and abundance of massive astrophysical compact halo objects (MACHOs) as a constituent of dark matter. Methods. We present spectroscopic observations of the central square arcminute of Eridanus 2 from MUSE-Faint, a survey of ultra-faint dwarf galaxies with the Multi Unit Spectroscopic Explorer on the Very Large Telescope. We derive line-of-sight velocities for possible member stars of the putative cluster and for stars in the centre of Eridanus 2. We discuss the existence of the cluster and determine new constraints for MACHOs using the Fokker-Planck diffusion approximation. Results. Out of 182 extracted spectra, we identify 26 member stars of Eridanus 2, seven of which are possible cluster members. We find intrinsic mean line-of-sight velocities of $79.7^{+3.1}_{-3.8}\,\mathrm{km}\,\mathrm{s}^{-1}$ and $76.0^{+3.2}_{-3.7}\,\mathrm{km}\,\mathrm{s}^{-1}$ for the cluster and the bulk of Eridanus 2, respectively, and intrinsic velocity dispersions of ${<}7.6\,\mathrm{km}\,\mathrm{s}^{-1}$ (68-$\%$ upper limit) and $10.3^{+3.9}_{-3.2}\,\mathrm{km}\,\mathrm{s}^{-1}$, respectively. This indicates the cluster most likely exists as a distinct dynamical population hosted by Eridanus 2, without surplus of dark matter over the background distribution. Among the member stars in the bulk of Eridanus 2, we find possible carbon stars, alluding to the existence of an intermediate-age population. We derive constraints on the fraction of dark matter that can consist of MACHOs with a given mass between $1$-$10^5\,M_\mathrm{sun}$. For dark matter consisting purely of MACHOs, the mass of the MACHOs must be less than ${\sim}7.6\,M_\mathrm{sun}$ and ${\sim}44\,M_\mathrm{sun}$ at a $68$- and $95$-$\%$ confidence level, respectively. (Abridged)

We show that a pair of field theory monodromies in which the shift symmetry is broken by small, well motivated deformations, naturally incorporates a mechanism for cancelling off radiative corrections to the cosmological constant. The lighter monodromy sector plays the role of inflation as well as providing a rigid degree of freedom that acts as a dynamical counterterm for the cosmological constant. The heavier monodromy sector includes a rigid dilaton that forces a global constraint on the system and the cancellation of vacuum energy loops occurs at low energies via the sequestering mechanism. This suggests that monodromy constructions in string theory could be adapted to incorporate mechanisms to stabilise the cosmological constant in their low energy descriptions.

We study quantum aspects of the galileon duality, especially in the case of a particular interacting galileon theory that is said to be dual to a free theory through the action of a simultaneous field and coordinate transformation. This would appear to map a theory with multiple vacua to one with a unique vacuum state. However, by regulating the duality transformation using external sources, we are able to preserve the full vacuum structure in the dual frame. By explicitly calculating the one-particle irreducible effective action on a maximally symmetric background, we identify a semi-classical contribution to the Wightman functions that has not been taken into account in previous analyses due to the singular point in the duality map. This may affect its spectral properties at high energy scales. These observations cast doubt on the main evidence in support of a non-local UV structure for galileons.

Galaxy mergers are key events in galaxy evolution, often causing massive starbursts and fueling active galactic nuclei (AGN). In these highly dynamic systems, it is not yet precisely known how much starbursts and AGN respectively contribute to the total luminosity, at what interaction stages they occur, and how long they persist. Here we estimate the fraction of the bolometric infrared (IR) luminosity that can be attributed to AGN by measuring and modeling the full ultraviolet to far-infrared spectral energy distributions (SEDs) in up to 33 broad bands for 24 merging galaxies with the Code for Investigating Galaxy Emission. In addition to a sample of 12 confirmed AGN in late-stage mergers, found in the $Infrared$ $Array$ $Satellite$ Revised Bright Galaxy Sample or Faint Source Catalog, our sample includes a comparison sample of 12 galaxy mergers from the $Spitzer$ Interacting Galaxies Survey, mostly early-stage. We perform identical SED modeling of simulated mergers to validate our methods, and we supplement the SED data with mid-IR spectra of diagnostic lines obtained with $Spitzer$ InfraRed Spectrograph. The estimated AGN contributions to the IR luminosities vary from system to system from 0% up to 91% but are significantly greater in the later-stage, more luminous mergers, consistent with what is known about galaxy evolution and AGN triggering.

In the first part of this paper we critically examine the ultra-violet implications of theories that exhibit Vainshtein screening, taking into account both the standard Wilsonian perspective as well as more exotic possibilities. Aspects of this discussion draw on results from the second part of the paper in which we perform a general study of derivatively coupled scalar theories using non-perturbative exact renormalisation group techniques, which are of interest independently of their application to modified gravity. In this context, we demonstrate the suppression of quantum corrections within the Vainshtein radius and discuss the potential relation with the classicalisation conjecture. We question whether the latter can be considered a realistic candidate for UV completion of large-scale modifications of gravity on account of a dangerously low classicalisation/strong coupling scale.

Gravitational models of self-tuning are those in which vacuum energy has no observable effect on spacetime curvature, even though it is a priori unsuppressed below the cut-off. We complement Weinberg's no go theorem by studying field theoretic completions of self-adjustment allowing for broken translations as well as other generalisations, and identify new obstructions. Our analysis uses a very general Kallen-Lehmann spectral representation of the exchange amplitude for conserved sources of energy-momentum and exploits unitarity and Lorentz invariance to show that a transition from self-tuning of long wavelength sources to near General Relativity on shorter scales is generically not possible. We search for novel ways around our obstructions and highlight two interesting possibilities. The first is an example of a unitary field configuration on anti-de Sitter space with the desired transition from self-tuning to GR. A second example is motivated by vacuum energy sequestering.

Recently two of the authors proposed a mechanism of vacuum energy sequester as a means of protecting the observable cosmological constant from quantum radiative corrections. The original proposal was based on using global Lagrange multipliers, but later a local formulation was provided. Subsequently other interesting claims of a different non-local approach to the cosmological constant problem were made, based again on global Lagrange multipliers. We examine some of these proposals and find their mutual relationship. We explain that the proposals which do not treat the cosmological constant counterterm as a dynamical variable require fine tunings to have acceptable solutions. Furthermore, the counterterm often needs to be retuned at every order in the loop expansion to cancel the radiative corrections to the cosmological constant, just like in standard GR. These observations are an important reminder of just how the proposal of vacuum energy sequester avoids such problems.

We recently formulated a local mechanism of vacuum energy sequester. This mechanism automatically removes all matter loop contributions to vacuum energy from the stress energy tensor which sources the curvature. Here we adapt the local vacuum energy sequestering mechanism to also cancel all the vacuum energy loops involving virtual gravitons, in addition to the vacuum energy generated by matter fields alone.

Recently it was argued that the conformal coupling of the chameleon to matter fields created an issue for early universe cosmology. As standard model degrees of freedom become non-relativistic in the early universe the chameleon is attracted towards a "surfing" solution, so that it arrives at the potential minimum with too large a velocity. This leads to rapid variations in the chameleon's mass and excitation of high energy modes, casting doubts on the classical treatment at Big Bang Nucleosynthesis. Here we present the DBI chameleon, a consistent high energy modification of the chameleon theory that dynamically renders it weakly coupled to matter during the early universe thereby eliminating the adverse effects of the `kicks'. This is done without any fine tuning of the coupling between the chameleon and matter fields and retains its screening ability in the solar system. We demonstrate this explicitly with a combination of analytic and numerical results.

We consider phase transitions and their contributions to vacuum energy in the manifestly local theory of vacuum energy sequestering. We demonstrate that the absence of instabilities imposes constraints on the couplings of gravitating and non-gravitating sectors, which can be satisfied in a large class of models. We further show by explicit construction that the vacuum energy contributions to the effective cosmological constant in the descendant vacua are generically strongly suppressed by the ratios of spacetime volumes of parent and descendant geometries. This means that the cosmological constant in de Sitter descendant vacua remains insensitive to phase transitions which may have occurred in the course of its cosmic history.

We present a manifestly local, diffeomorphism invariant and locally Poincare invariant formulation of vacuum energy sequestering. In this theory, quantum vacuum energy generated by matter loops is cancelled by auxiliary fields. The auxiliary fields decouple from gravity almost completely. Their only residual effect is an a priori arbitrary, finite contribution to the curvature of the background geometry, which is radiatively stable. Its value is to be determined by a measurement, like the finite part of any radiatively stable UV-sensitive quantity in quantum field theory.

We discuss unimodular gravity at a classical level, and in terms of its extension into the UV through an appropriate path integral representation. Classically, unimodular gravity is simply a gauge fixed version of General Relativity (GR), and as such it yields identical dynamics and physical predictions. We clarify this and explain why there is no sense in which it can "bring a new perspective" to the cosmological constant problem. The quantum equivalence between unimodular gravity and GR is more of a subtle question, but we present an argument that suggests one can always maintain the equivalence up to arbitrarily high momenta. As a corollary to this, we argue that whenever inequivalence is seen at the quantum level, that just means we have defined two different quantum theories that happen to share a classical limit.

Light degrees of freedom that modify gravity on cosmological scales must be "screened" on solar system scales in order to be compatible with data. The Vainshtein mechanism achieves this through a breakdown of classical perturbation theory, as large interactions involving new degrees of freedom become important below the so-called Vainshtein radius. We begin to develop an extension of the Parameterized Post-Newtonian (PPN) formalism that is able to handle Vainshteinian corrections. We argue that theories with a unique Vainshtein scale must be expanded using two small parameters. In this Parameterized Post-Newtonian-Vainshteinian (PPNV) expansion, the primary expansion parameter that controls the PPN order is, as usual, the velocity $v$. The secondary expansion parameter, $\alpha$, controls the strength of the Vainshteinian correction and is a theory-specific combination of the Schwarzschild radius and the Vainshtein radius of the source that is independent of its mass. We present the general framework and apply it to Cubic Galileon theory both inside and outside the Vainshtein radius. The PPNV framework can be used to determine the compatibility of such theories with solar system and other strong-field data.

Recently we proposed a mechanism for sequestering the Standard Model vacuum energy that predicts that the universe will collapse. Here we present a simple mechanism for bringing about this collapse, employing a scalar field whose potential is linear and becomes negative, providing the negative energy density required to end the expansion. The slope of the potential is chosen to allow for the expansion to last until the current Hubble time, about $10^{10}$ years, to accommodate our universe. Crucially, this choice is technically natural due to a shift symmetry. Moreover, vacuum energy sequestering selects {\it radiatively stable} initial conditions for the collapse, which guarantee that immediately before the turnaround the universe is dominated by the linear potential which drives an epoch of accelerated expansion for at least an efold. Thus a single, technically natural choice for the slope ensures that the collapse is imminent and is preceded by the current stage of cosmic acceleration, giving a new answer to the 'Why Now?' problem.

These lectures on the cosmological constant problem were prepared for the X Mexican School on Gravitation and Mathematical Physics. The problem itself is explained in detail, emphasising the importance of radiative instability and the need to repeatedly fine tune as we change our effective description. Weinberg's no go theorem is worked through in detail. I review a number of proposals including Linde's universe multiplication, Coleman's wormholes, the fat graviton, and SLED, to name a few. Large distance modifications of gravity are also discussed, with causality considerations pointing towards a global modification as being the most sensible option. The global nature of the cosmological constant problem is also emphasized, and as a result, the sequestering scenario is reviewed in some detail, demonstrating the cancellation of the Standard Model vacuum energy through a global modification of General Relativity.

Recently we suggested a reformulation of General Relativity which completely sequesters from gravity {\it all} of the vacuum energy from a protected matter sector, assumed to contain the Standard Model. Here we elaborate further on the mechanism, presenting additional details of how it cancels all loop corrections and renders all contributions from phase transitions automatically small. We also consider cosmological consequences in more detail and show that the mechanism is consistent with a variety of inflationary models that make a universe big and old. We discuss in detail the underlying assumptions behind the dynamics of our proposal, and elaborate on the relationship of the physical interpretation of divergent operators in quantum field theory and the apparent `acausality' which our mechanism seems to entail, which we argue is completely harmless. It is merely a reflection of the fact that any UV sensitive quantity in quantum field theory cannot be calculated from first principles, but is an input whose numerical value must be measured. We also note that since the universe should be compact in spacetime, and so will collapse in the future, the current phase of acceleration with $w_{DE}\approx-1$ is just a transient. This could be tested by future cosmological observations.

We propose a very simple reformulation of General Relativity, which completely sequesters from gravity {\it all} of the vacuum energy from a matter sector, including all loop corrections and renders all contributions from phase transitions automatically small. The idea is to make the dimensional parameters in the matter sector functionals of the 4-volume element of the universe. For them to be nonzero, the universe should be finite in spacetime. If this matter is the Standard Model of particle physics, our mechanism prevents any of its vacuum energy, classical or quantum, from sourcing the curvature of the universe. The mechanism is consistent with the large hierarchy between the Planck scale, electroweak scale and curvature scale, and early universe cosmology, including inflation. Consequences of our proposal are that the vacuum curvature of an old and large universe is not zero, but very small, that $w_{DE} \simeq -1$ is a transient, and that the universe will collapse in the future.

We have recently proposed a novel self tuning mechanism to alleviate the famous cosmological constant problem, based on the general scalar tensor theory proposed by Horndeski. The self-tuning model ends up consisting of four geometric terms in the action, with each term containing a free potential function of the scalar field; the four together being labeled as the Fab-Four. In this paper we begin the important task of deriving the cosmology associated with the Fab-Four Lagrangian. Performing a phase plane analysis of the system we are able to obtain a number of fixed points for the system, with some remarkable new solutions emerging from the trade-off between the various potentials. As well as obtaining inflationary solutions we also find conventional radiation/matter-like solutions, but in regimes where the energy density is dominated by a cosmological constant, and where we do not have any explicit forms of radiation or matter. Stability conditions for matter solutions are obtained and we show how it is possible for there to exist an extended period of `matter domination' opening up the possibility that we can generate cosmological structures, and recover a consistent cosmology even in the presence of a large cosmological constant.

We propose a completely new parametrisation of the dark energy equation of state, which uses the dark energy density, $\Omega_e$ as a cosmic clock. We expand the equation of state in a series of orthogonal polynomials, with $\Omega_e$ as the expansion parameter and determine the expansion coefficients by fitting to SNIa and $H(z)$ data. Assuming that $\Omega_e$ is a monotonic function of time, we show that our parametrisation performs better than the popular Chevallier--Polarski--Linder (CPL) and Gerke and Efstathiou (GE) parametrisations, and we demonstrate that it is robust to the choice of prior. Expanding in orthogonal polynomials allows us to relate models of dark energy directly to our parametrisation, which we illustrate by placing constraints on the expansion coefficients extracted from two popular quintessence models. Finally, we comment on how this parametrisation could be modified to accommodate high redshift data, where any non--monotonicity of $\Omega_e$ would need to be accounted for.

Many models have been studied that contain more than one species of dark matter and some of these couple the Cold Dark Matter (CDM) to a light scalar field. In doing this we introduce additional long range forces, which in turn can significantly affect our estimates of cosmological parameters if not properly accounted for. It is, therefore, important to study these models and their resulting cosmological implications. We present a model in which a fraction of the total cold dark matter density is coupled to a scalar field. We study the background and perturbation evolution and calculate the resulting Cosmic Microwave Background anisotropy spectra. The greater the fraction of dark matter coupled to the scalar field and the stronger the coupling strength, the greater the deviation of the background evolution from LCDM. Previous work, with a single coupled dark matter species, has found an upper limit on the coupling strength of order O(0.1). We find that with a coupling of this magnitude more than half the dark matter can be coupled to a scalar field without producing any significant deviations from LCDM.

We find a covariant completion of the flat-space multi-galileon theory, preserving second-order field equations. We then generalise this to arrive at an enlarged class of second order theories describing multiple scalars and a single tensor, and conjecture that these are a multi-scalar version of Horndeski's most general scalar-tensor theory.

The theory of a single massive graviton has a cutoff much below its Planck scale, because the extra modes from the graviton multiplet involve higher derivative self-interactions, controlled by a scale convoluted from the small graviton mass. Generically, these correct the propagator by environmental effects. The resulting effective cutoff depends on the environmental parameters and the graviton mass. Requiring the theory to be perturbative down to ${\cal O}(1) mm$, we derive bounds on the graviton mass, corresponding to $\gtrsim {\cal O}(1) meV$ for the generic case, and somewhat weaker bounds in cases of fine-tuning. In all cases the mass is required to be much too large for the theory to conform with GR at cosmological distances. Similar results also hold in quartic and quintic Galileon theory.

We continue to introduce bi-galileon theory, the generalisation of the single galileon model introduced by Nicolis et al. The theory contains two coupled scalar fields and is described by a Lagrangian that is invariant under Galilean shifts in those fields. This paper is the second of two, and focuses on the phenomenology of the theory. We are particularly interesting in models that admit solutions that are asymptotically self accelerating or asymptotically self tuning. In contrast to the single galileon theories, we find examples of self accelerating models that are simultaneously free from ghosts, tachyons and tadpoles, able to pass solar system constraints through Vainshtein screening, and do not suffer from problems with superluminality, Cerenkov emission or strong coupling. We also find self tuning models and discuss how Weinberg's no go theorem is evaded by breaking Poincar\'e invariance in the scalar sector. Whereas the galileon description is valid all the way down to solar system scales for the self-accelerating models, unfortunately the same cannot be said for self tuning models owing to the scalars backreacting strongly on to the geometry.

We present a novel idea for screening the vacuum energy contribution to the overall value of the cosmological constant, thereby enabling us to choose the bare value of the vacuum curvature empirically, without any need to worry about the zero-point energy contributions of each particle. The trick is to couple matter to a metric that is really a composite of other fields, with the property that the square-root of its determinant is the integrand of a topological invariant, and/or a total derivative. This ensures that the vacuum energy contribution to the Lagrangian is non-dynamical. We then give an explicit example of a theory with this property that is free from Ostrogradski ghosts, and is consistent with solar system physics and cosmological tests.

We show how to generalize the classical duals found by Gabadadze {\it et al} to a very large class of self-interacting theories. This enables one to adopt a perturbative description beyond the scale at which classical perturbation theory breaks down in the original theory. This is particularly relevant if we want to test modified gravity scenarios that exhibit Vainshtein screening on solar system scales. We recognise the duals as being related to the Legendre transform of the original Lagrangian, and present a practical method for finding the dual in general; our methods can also be applied to self-interacting theories with a hierarchy of strong coupling scales, and with multiple fields. We find the classical dual of the full quintic galileon theory as an example.

In this review we present a thoroughly comprehensive survey of recent work on modified theories of gravity and their cosmological consequences. Amongst other things, we cover General Relativity, Scalar-Tensor, Einstein-Aether, and Bimetric theories, as well as TeVeS, f(R), general higher-order theories, Horava-Lifschitz gravity, Galileons, Ghost Condensates, and models of extra dimensions including Kaluza-Klein, Randall-Sundrum, DGP, and higher co-dimension braneworlds. We also review attempts to construct a Parameterised Post-Friedmannian formalism, that can be used to constrain deviations from General Relativity in cosmology, and that is suitable for comparison with data on the largest scales. These subjects have been intensively studied over the past decade, largely motivated by rapid progress in the field of observational cosmology that now allows, for the first time, precision tests of fundamental physics on the scale of the observable Universe. The purpose of this review is to provide a reference tool for researchers and students in cosmology and gravitational physics, as well as a self-contained, comprehensive and up-to-date introduction to the subject as a whole.

We have recently proposed a special class of scalar tensor theories known as the Fab Four. These arose from attempts to analyse the cosmological constant problem within the context of Horndeski's most general scalar tensor theory. The Fab Four together give rise to a model of self-tuning, with the relevant solutions evading Weinberg's no-go theorem by relaxing the condition of Poincare invariance in the scalar sector. The Fab Four are made up of four geometric terms in the action with each term containing a free potential function of the scalar field. In this paper we rigorously derive this model from the general model of Horndeski, proving that the Fab Four represents the only classical scalar tensor theory of this type that has any hope of tackling the cosmological constant problem. We present the full equations of motion for this theory, and give an heuristic argument to suggest that one might be able to keep radiative corrections under control. We also give the Fab Four in terms of the potentials presented in Deffayet et al's version of Horndeski.

We study the fields of spherically symmetric thin shell sources, a.k.a. Dyson spheres, in a {\it fully nonlinear covariant} theory of gravity with the simplest galileon field. We integrate exactly all the field equations once, reducing them to first order nonlinear equations. For the simplest galileon, static solutions come on {\it six} distinct branches. On one, a Dyson sphere surrounds itself with a galileon hair, which far away looks like a hair of any Brans-Dicke field. The hair changes below the Vainshtein scale, where the extra galileon terms dominate the minimal gradients of the field. Their hair looks more like a fuzz, because the galileon terms are suppressed by the derivative of the volume determinant. It shuts off the `hair bunching' over the `angular' 2-sphere. Hence the fuzz remains dilute even close to the source. This is really why the Vainshtein's suppression of the modifications of gravity works close to the source. On the other five branches, the static solutions are all {\it singular} far from the source, and shuttered off from asymptotic infinity. One of them, however, is really the self-accelerating branch, and the singularity is removed by turning on time dependence. We give examples of regulated solutions, where the Dyson sphere explodes outward, and its self-accelerating side is nonsingular. These constructions may open channels for nonperturbative transitions between branches, which need to be addressed further to determine phenomenological viability of multi-branch gravities.

Starting from the most general scalar-tensor theory with second order field equations in four dimensions, we establish the unique action that will allow for the existence of a consistent self-tuning mechanism on FLRW backgrounds, and show how it can be understood as a combination of just four base Lagrangians with an intriguing geometric structure dependent on the Ricci scalar, the Einstein tensor, the double dual of the Riemann tensor and the Gauss-Bonnet combination. Spacetime curvature can be screened from the net cosmological constant at any given moment because we allow the scalar field to break Poincar\'e invariance on the self-tuning vacua, thereby evading the Weinberg no-go theorem. We show how the four arbitrary functions of the scalar field combine in an elegant way opening up the possibility of obtaining non-trivial cosmological solutions.

A sizable fraction of the total energy density of the universe may be in heavy particles with a net dark $U(1)'$ charge comparable to its mass. When the charges have the same sign the cancellation between their gravitational and gauge forces may lead to a mismatch between different measures of masses in the universe. Measuring galactic masses by orbits of normal matter, such as galaxy rotation curves or lensing, will give the total mass, while the flows of dark matter agglomerates may yield smaller values if the gauge repulsion is not accounted for. If distant galaxies which house light beacons like SNe Ia contain such dark particles, the observations of their cosmic recession may mistake the weaker forces for an extra `antigravity', and infer an effective dark energy equation of state smaller than the real one. In some cases, including that of a cosmological constant, these effects can mimic $w<-1$. They can also lead to a {\it local} variation of galaxy-galaxy forces, yielding a larger `Hubble Flow' in those regions of space that could be taken for a dynamical dark energy, or superhorizon effects.

We study the spectrum of gravitational perturbations around a vacuum de Sitter brane in a 5D asymmetric braneworld model, with induced curvature on the brane. This generalises the stealth acceleration model proposed by Charmousis, Gregory and Padilla (CGP) which realises the Cardassian cosmology in which power law cosmic acceleration can be driven by ordinary matter. Whenever the bulk has infinite volume we find that there is always a perturbative ghost propagating on the de Sitter brane, in contrast to the Minkowski brane case analysed by CGP. We discuss the implication of this ghost for the stealth acceleration model, and identify a limiting case where the ghost decouples as the de Sitter curvature vanishes.

We consider a magnetic Bianchi I braneworld, embedded in between two Schwarzschild-AdS spacetimes, boosted equal amounts in opposite directions and compare them to the analagous solution in four-dimensional General Relativity. The efficient dissipation of anisotropy on the brane is explicitly demonstrated, a process we dub braneworld isotropization. From the bulk point of view, we attribute this to anisotropic energy being carried into the bulk by hot gravitons leaving the brane. From the brane point of view this can be interpreted in terms of the production of particles in the dual CFT. We explain how this result enables us to gain a better understanding of the behaviour of anisotropic branes already studied in the literature. We also show how there is evidence of particles being over-produced, and comment on how this may ultimately provide a possible observational signature of braneworlds.

We revisit the problem of defining non-minimal gravity in the first order formalism. Specializing to scalar-tensor theories, which may be disguised as `higher-derivative' models with the gravitational Lagrangians that depend only on the Ricci scalar, we show how to recast these theories as Palatini-like gravities. The correct formulation utilizes the Lagrange multiplier method, which preserves the canonical structure of the theory, and yields the conventional metric scalar-tensor gravity. We explain the discrepancies between the na\"ive Palatini and the Lagrange multiplier approach, showing that the na\"ive Palatini approach really swaps the theory for another. The differences disappear only in the limit of ordinary General Relativity, where an accidental redundancy ensures that the na\"ive Palatini works there. We outline the correct decoupling limits and the strong coupling regimes. As a corollary we find that the so-called `Modified Source Gravity' models suffer from strong coupling problems at very low scales, and hence cannot be a realistic approximation of our universe. We also comment on a method to decouple the extra scalar using the chameleon mechanism.

We examine brane induced gravity on codimension-1 branes, a.k.a DGP gravity, as a theory of five-dimensional gravity containing a certain class four-dimensional branes. From this perspective, the model suffers from a number of pathologies which went unnoticed before. By generalizing the 5D geometry from Minkowski to Schwarzschild, we find that when the bulk mass is large enough, the brane hits a pressure singularity at finite radius. Further, on the self-accelerating branch, the five-dimensional energy is unbounded from below, implying that the self-accelerating backgrounds are unstable. Even in an empty Minkowski bulk, standard Euclidean techniques suggest that the spontaneous nucleation of self-accelerating branes is unsuppressed. If so, quantum effects will strongly modify any classical intuition about the theory. We also note that unless considered as Z_2-orbifold boundaries, self-accelerating branes correspond to `wormhole' configurations, which introduces the usual problematic issues associated with wormholes. Altogether these pathologies present a serious challenge that any proposed UV completion of the DGP model must overcome.

We show how to construct consistent braneworld models which exhibit late time acceleration. Unlike self-acceleration, which has a de Sitter vacuum state, our models have the standard Minkowski vacuum and accelerate only in the presence of matter, which we dub ``stealth-acceleration''. We use an effective action for the brane which includes an induced gravity term, and allow for an asymmetric set-up. We study the linear stability of flat brane vacua and find the regions of parameter space where the set-up is stable. The 4-dimensional graviton is only quasi-localised in this set-up and as a result gravity is modified at late times. One of the two regions is strongly coupled and the scalar mode is eaten up by an extra symmetry that arises in this limit. Having filtered the well-defined theories we then focus on their cosmology. When the graviton is quasi-localised we find two main examples of acceleration. In each case, we provide an illustrative model and compare it to LambdaCDM.

We systematically explore the spectrum of gravitational perturbations in codimension-1 DGP braneworlds, and find a 4D ghost on the self-accelerating branch of solutions. The ghost appears for any value of the brane tension, although depending on the sign of the tension it is either the helicity-0 component of the lightest localized massive tensor of mass $0<m^2 < 2H^2$ for positive tension, the scalar `radion' for negative tension, or their admixture for vanishing tension. Because the ghost is gravitationally coupled to the brane-localized matter, the self-accelerating solutions are not a reliable benchmark for cosmic acceleration driven by gravity modified in the IR. In contrast, the normal branch of solutions is ghost-free, and so these solutions are perturbatively safe at large distance scales. We further find that when the $\mathbb{Z}_2$ orbifold symmetry is broken, new tachyonic instabilities, which are much milder than the ghosts, appear on the self-accelerating branch. Finally, using exact gravitational shock waves we analyze what happens if we relax boundary conditions at infinity. We find that non-normalizable bulk modes, if interpreted as 4D phenomena, may open the door to new ghost-like excitations.

We consider a single 3-brane sitting in between two different five dimensional spacetimes. On each side of the brane, the bulk is a solution to Gauss-Bonnet gravity, although the bare cosmological constant, funda mental Planck scale, and Gauss-Bonnet coupling can differ. This asymmetry leads to weighted junction conditions across the brane and interesting brane cosmology. We focus on two special cases: a generalized Randall-Sundrum model without any Gauss-Bonnet terms, and a stringy model, without any bare cosmological constants, and positive Gauss-Bonnet coupling. Even though we assume there is no vacuum energy on the brane, we find late time de Sitter cosmologies can occur. Remarkably, in certain parameter regions, this acceleration is preceded by a period of matter/radiation domination, with $H^2 \propto \rho$, all the way back to nucleosynthesis.

We consider a single Minkowski brane sandwiched in between two copies of anti-de Sitter space. We allow the bulk Planck mass and cosmological constant to differ on either side of the brane. Linearised perturbations about this background reveal that gravity can be modified in the infra-red. At intermediate scales, the braneworld propagator mimics four-dimensional GR in that it has the correct momentum dependance. However it has the wrong tensor structure. Beyond a source dependant scale, we show that quadratic brane bending contributions become important, and conspire to correct the tensor structure of the propagator. We argue that even higher order terms can consistently be ignored up to very high energies, and suggest that there is no problem with strong coupling. We also consider scalar and vector perturbations in the bulk, checking for scalar ghosts.

We consider a generalisation of the DGP model, by adding a second brane with localised curvature, and allowing for a bulk cosmological constant and brane tensions. We study radion and graviton fluctuations in detail, enabling us to check for ghosts and tachyons. By tuning our parameters accordingly, we find bigravity models that are free from ghosts and tachyons. These models will lead to large distance modifications of gravity that could be observable in the near future.