We pioneer the black hole mass gap as a powerful new tool for constraining new particles. A new particle that couples to the Standard Model---such as an axion---acts as an additional source of loss in the cores of population-III stars, suppressing mass lost due to winds and quenching the pair-instability. This results in heavier astrophysical black holes. As an example, using stellar simulations we show that the solar axion explanation of the recent XENON1T excess implies astrophysical black holes of ~ 56 MS, squarely within the black hole mass gap predicted by the Standard Model.

The local distance ladder measurement of the Hubble constant requires a connection between geometric distances at low redshift and Type Ia supernovae in the Hubble flow, which may be achieved through either the Cepheid period--luminosity relation or the luminosity of the Tip of the Red Giant Branch (TRGB) feature of the Hertzsprung--Russell diagram. Any potential solution to the Hubble tension that works by altering the distance ladder must produce consistency of both the Cepheid and TRGB $H_0$ calibrations with the CMB. In this paper we extend our models of screened fifth forces (Desmond et al 2019) to cover the TRGB framework. A fifth force lowers TRGB luminosity, so a reduction in inferred $H_0$ requires that the stars that calibrate the luminosity---currently in the LMC---are on average less screened than those that calibrate the supernova magnitude. We show that even under pessimistic assumptions for the extinction to the LMC, full consistency with Planck can be achieved for a fifth force strength in unscreened RGB stars $\sim$0.2 that of Newtonian gravity. This is allowed by the comparison of Cepheid and TRGB distance measurements to nearby galaxies. Our results indicate that the framework of Desmond et al (2019) is more versatile than initially demonstrated, capable of ameliorating the Hubble tension on a second front.

The Hubble tension can be significantly eased if there is an early component of dark energy that becomes active around the time of matter-radiation equality. Early dark energy models suffer from a coincidence problem -- the physics of matter-radiation equality and early dark energy are completely disconnected, so some degree of fine-tuning is needed in order for them to occur nearly simultaneously. In this paper we propose a natural explanation for this coincidence. If the early dark energy scalar couples to neutrinos then it receives a large injection of energy around the time that neutrinos become non-relativistic. This is precisely when their temperature is of order their mass, which, coincidentally, occurs around the time of matter-radiation equality. Neutrino decoupling therefore provides a natural trigger for early dark energy by displacing the field from its minimum just before matter-radiation equality. We discuss various theoretical aspects of this proposal, potential observational signatures, and future directions for its study.

The discrepancy between the values of the Hubble constant $H_0$ derived from the local distance ladder and the cosmic microwave background provides a tantalising hint of new physics. We explore a potential resolution involving screened fifth forces in the local Universe, which alter the Cepheid calibration of supernova distances. In particular, if the Cepheids with direct distance measurements from parallax or water masers are screened but a significant fraction of those in other galaxies are not, neglecting the difference between their underlying period--luminosity relations biases the local $H_0$ measurement high. This difference derives from a reduction in the Cepheid pulsation period and possible increase in luminosity under a fifth force. We quantify the internal and environmental gravitational properties of the Riess et al. distance ladder galaxies to assess their degrees of screening under a range of phenomenological models, and propagate this information into the $H_0$ posterior as a function of fifth force strength. We consider well-studied screening models in scalar--tensor gravity theories such as chameleon, K-mouflage and Vainshtein, along with a recently-proposed mechanism based on baryon--dark matter interactions in which screening is governed by local dark matter density. We find that a fifth force strength $\sim5-30\%$ that of gravity can alleviate (though not resolve) the $H_0$ tension in some scenarios, around the sensitivity level at which tests based on other distance ladder data can constrain this strength. Although our analysis is exploratory and based on screening models not necessarily realised in full theories, our results demonstrate that new physics-based local resolutions of the $H_0$ tension are possible, supplementing those already known in the pre-recombination era.

We derive the details of a new screening mechanism where the interactions of baryons and dark matter can be screened according to the local dark matter density. In this mechanism, the value of Newton's constant is dark matter density-dependent, allowing for the possibility that astrophysical phenomena are very different in galaxies less dense than the Milky Way. The parameterized post Newtonian parameter $\gamma$, which quantifies the difference between kinematical and lensing probes, also depends on dark matter density. We calculate the effects of varying $G$ on various stages of stellar evolution, focusing on observables that impact cosmology: the Cepheid period--luminosity relation and the supernova Ia magnitude--redshift relation. Other potential tests of the model are also investigated including main-sequence, post-main sequence, and low mass dwarf stars. Finally, we discuss how extragalactic tests of $\gamma$ could provide complementary constraints.

Screened modified gravity theories evade the solar system tests that have proved prohibitive for classical alternative gravity theories. In many cases, they do not fit into the PPN formalism. The environmental dependence of the screening has motivated a concerted effort to find new and novel probes of gravity using objects that are well-studied but have hitherto not been used to test gravity. Astrophysical objects---stars, galaxies, clusters---have proved competitive tools for this purpose since they occupy the partially-screened regime between solar system and the Hubble flow. In this article we review the current astrophysical tests of screened modified gravity theories.

We introduce The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists. This review focuses on tests of gravity and is intended to be of use primarily to observers, but also to theorists with interest in the development of experimental tests. It is twinned with a separate review on tests of dark matter self-interactions (Adhikari et al., in prep.). Our focus is on astrophysical probes of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. These are complementary to both strong-field tests and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified gravity models lead to characteristic signals specifically on astrophysical scales. The constraining power of these signals is not limited by cosmic variance, but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, clusters and large scale structure. In this review we lay the groundwork for a thorough exploration of the astrophysical regime with an eye to using the current and next generation of observations for tests of gravity. We begin by setting the scene for how theories beyond General Relativity are expected to behave, focusing primarily on screened fifth forces. We describe the analytic and numerical techniques for exploring the pertinent astrophysical systems, as well as the signatures of modified gravity. With these in hand we present a range of observational tests, and discuss prospects for future measurements and theoretical developments.

Despite two decades of tremendous experimental and theoretical progress, the riddle of the accelerated expansion of the Universe remains to be solved. On the experimental side, our understanding of the possibilities and limitations of the major dark energy probes has evolved; here we summarize the major probes and their crucial challenges. On the theoretical side, the taxonomy of explanations for the accelerated expansion rate is better understood, providing clear guidance to the relevant observables. We argue that: i) improving statistical precision and systematic control by taking more data, supporting research efforts to address crucial challenges for each probe, using complementary methods, and relying on cross-correlations is well motivated; ii) blinding of analyses is difficult but ever more important; iii) studies of dark energy and modified gravity are related; and iv) it is crucial that R&D for a vibrant dark energy program in the 2030s be started now by supporting studies and technical R&D that will allow embryonic proposals to mature. Understanding dark energy, arguably the biggest unsolved mystery in both fundamental particle physics and cosmology, will remain one of the focal points of cosmology in the forthcoming decade.

Scalar-tensor theories of gravity where a new scalar degree of freedom couples to the Gauss-Bonnet invariant can exhibit the phenomenon of spontaneous black hole scalarization. These theories admit both the classic black holes predicted by general relativity as well as novel hairy black hole solutions. The stability of the hairy black holes is strongly dependent on the precise form of the scalar-gravity coupling. A radial stability investigation revealed that all scalarized black hole solutions are unstable when the coupling between the scalar field and the Gauss-Bonnet invariant is quadratic in the scalar, whereas stable solutions exist for exponential couplings. Here we elucidate this behavior. We demonstrate that, while the quadratic term controls the onset of the tachyonic instability that gives rise to the black hole hair, the higher-order coupling terms control the nonlinearities that quench that instability, and hence also control the stability of the hairy black hole solutions.

It has recently been shown that nontrivial couplings between a scalar and the Gauss-Bonnet invariant can give rise to black hole spontaneous scalarization. Theories that exhibit this phenomenon are among the leading candidates for testing gravity with upcoming black hole observations. All models considered so far have focused on specific forms for the coupling, neglecting scalar self-interactions. In this work, we take the first steps towards placing this phenomenon on a more robust theoretical footing by considering the leading-order scalar self-interactions as well as the scalar-Gauss-Bonnet coupling. Our approach is consistent with the principles of effective field theory and yields the simplest and most natural model. We find that a mass term for the scalar alters the threshold for the onset of scalarization, and we study the mass range over which scalarized black hole solutions exist. We also demonstrate that the quartic self-coupling is sufficient to produce scalarized solutions that are stable against radial perturbations, without the need to resort to higher-order terms in the Gauss-Bonnet coupling function. Our model therefore represents a canonical model that can be studied further, with the ultimate aim of developing falsifiable tests of black hole scalarization.

We study the phenomenology of a recent string construction with a quantum mechanically stable dark energy. A mild supersymmetry protects the vacuum energy but also allows $O(10 - 100)$ TeV scale superpartner masses. The construction is holographic in the sense that the 4D spacetime is generated from "pixels" originating from five-branes wrapped over metastable five-cycles of the compactification. The cosmological constant scales as $\Lambda \sim 1/N$ in the pixel number. An instability in the construction leads to cosmic expansion. This also causes more five-branes to wind up in the geometry, leading to a slowly decreasing cosmological constant which we interpret as an epoch of inflation followed by (pre-)heating when a rare event occurs in which the number of pixels increases by an order one fraction. The sudden appearance of radiation triggers an exponential increase in the number of pixels. Dark energy has a time varying equation of state with $w_a=-3\Omega_{m,0}(1+w_0)/2$, which is compatible with current bounds, and could be constrained further by future data releases. The pixelated nature of the Universe also implies a large-$l$ cutoff on the angular power spectrum of cosmological observables with $l_{\rm max} \sim O(N)$. We also use this pixel description to study the thermodynamics of de Sitter space, finding rough agreement with effective field theory considerations.

We investigate the existence and behavior of oscillons in theories in which higher derivative terms are present in the Lagrangian, such as galileons. Such theories have emerged in a broad range of settings, from higher-dimensional models, to massive gravity, to models for late-time cosmological acceleration. By focusing on the simplest example---massive galileon effective field theories---we demonstrate that higher derivative terms can lead to the existence of completely new oscillons (quasi-breathers). We illustrate our techniques in the artificially simple case of 1 + 1 dimensions, and then present the complete analysis valid in 2 + 1 and 3 + 1 dimensions, exploring precisely how these new solutions are supported entirely by the non-linearities of the quartic galileon. These objects have the novel peculiarity that they are of the differentiability class $C^1$.

Chameleon theories of gravity predict that the gaseous component of isolated dwarf galaxies rotates with a faster velocity than the stellar component. In this paper, we exploit this effect to obtain new constraints on the model parameters using the measured rotation curves of six low surface brightness galaxies. For $f(R)$ theories, we rule out values of $f_{R0}>10^{-6}$. For more general theories, we find that the constraints from Cepheid variable stars are currently more competitive than the bounds we obtain here but we are able to rule out self-screening parameters $\chi_c>10^{-6}$ for fifth-force strengths (coupling of the scalar to matter) as low as $0.05$ the Newtonian force. This region of parameter space has hitherto been inaccessible to astrophysical probes. We discuss the future prospects for improving these bounds.

The splashback radius is a physical scale in dark matter halos that is set by the gravitational dynamics of recently accreted shells. We use analytical models and N-body simulations to study the dependence of splashback on dark energy and screened modified gravity theories. In modified gravity models, the transition from screened to unscreened regions typically occurs in the cluster outskirts, suggesting potentially observable signatures in the splashback feature. We investigate the location of splashback in both chameleon and Vainshtein screened models and find significant differences compared with $\Lambda$CDM predictions. We also find an interesting interplay between dynamical friction and modified gravity, providing a distinctive signature for modified gravity models in the behavior of the splashback feature as a function of galaxy luminosity.

We introduce a novel method to circumvent Weinberg's no-go theorem for self-tuning the cosmological vacuum energy: a Lorentz-violating finite-temperature superfluid can counter the effects of an arbitrarily large cosmological constant. Fluctuations of the superfluid result in the graviton acquiring a Lorentz-violating mass and we identify a unique class of theories that are pathology free, phenomenologically viable, and do not suffer from instantaneous modes. This new and hitherto unidentified phase of massive gravity propagates the same degrees of freedom as general relativity with an additional Lorentz-violating scalar that is introduced by higher-derivative operators in a UV insensitive manner. The superfluid is therefore a consistent infrared modification of gravity. We demonstrate how the superfluid can degravitate a cosmological constant and discuss its phenomenology.

We identify a class of scalar-tensor theories with coupling between the scalar and the Gauss-Bonnet invariant that exhibit spontaneous scalarization for both black holes and compact stars. In particular, these theories formally admit all of the stationary solutions of general relativity, but these are not dynamically preferred if certain conditions are satisfied. Remarkably, black holes exhibit scalarization if their mass lies within one of many narrow bands. We find evidence that scalarization can occur in neutron stars as well.

Future space-based tests of relativistic gravitation-laser ranging to Phobos, accelerometers in orbit, and optical networks surrounding Earth-will constrain the theory of gravity with unprecedented precision by testing the inverse-square law, the strong and weak equivalence principles, and the deflection and time-delay of light by massive bodies. In this paper, we estimate the bounds that could be obtained on alternative gravity theories that use screening mechanisms to suppress deviations from general relativity in the solar system: chameleon, symmetron, and galileon models. We find that space-based tests of the parameterized post-Newtonian parameter $\gamma$ will constrain chameleon and symmetron theories to new levels in the solar system, and that tests of the inverse-square law using laser ranging to Phobos will provide the most stringent constraints on galileon theories to date. We end by discussing the potential for constraining these theories using upcoming tests of the weak equivalence principle, and conclude that further theoretical modeling is required in order to fully utilize the data.

Theories of modified gravity where light scalars with non-trivial self-interactions and non-minimal couplings to matter-chameleon and symmetron theories-dynamically suppress deviations from general relativity in the solar system. On other scales, the environmental nature of the screening means that such scalars may be relevant. The highly-nonlinear nature of screening mechanisms means that they evade classical fifth-force searches, and there has been an intense effort towards designing new and novel tests to probe them, both in the laboratory and using astrophysical objects, and by reinterpreting existing datasets. The results of these searches are often presented using different parametrizations, which can make it difficult to compare constraints coming from different probes. The purpose of this review is to summarize the present state-of-the-art searches for screened scalars coupled to matter, and to translate the current bounds into a single parametrization to survey the state of the models. Presently, commonly studied chameleon models are well-constrained but less commonly studied models have large regions of parameter space that are still viable. Symmetron models are constrained well by astrophysical and laboratory tests, but there is a desert separating the two scales where the model is unconstrained. The coupling of chameleons to photons is tightly constrained but the symmetron coupling has yet to be explored. We also summarize the current bounds on $f(R)$ models that exhibit the chameleon mechanism (Hu \& Sawicki models). The simplest of these are well constrained by astrophysical probes, but there are currently few reported bounds for theories with higher powers of $R$. The review ends by discussing the future prospects for constraining screened modified gravity models further using upcoming and planned experiments.

The LIGO/VIRGO collaboration has recently announced the detection of gravitational waves from a neutron star-neutron star merger (GW170817) and the simultaneous measurement of an optical counterpart (the gamma-ray burst GRB 170817A). The close arrival time of the gravitational and electromagnetic waves limits the difference in speed of photons and gravitons to be less than about one part in $10^{15}$. This has three important implications for cosmological scalar-tensor gravity theories that are often touted as dark energy candidates and alternatives to $\Lambda$CDM. First, for the most general scalar-tensor theories---beyond Horndeski models---three of the five parameters appearing in the effective theory of dark energy can now be severely constrained on astrophysical scales; we present the results of combining the new gravity wave results with galaxy cluster observations. Second, the combination with the lack of strong equivalence principle violations exhibited by the supermassive black hole in M87, constrains the quartic galileon model to be cosmologically irrelevant. Finally, we derive a new bound on the disformal coupling to photons that implies that such couplings are irrelevant for the cosmic evolution of the field.

Generic scalar-tensor theories of gravity predict deviations from Newtonian physics inside astrophysical bodies. In this paper, we point out that low mass stellar objects, red and brown dwarf stars, are excellent probes of these theories. We calculate two important and potentially observable quantities: the radius of brown dwarfs and the minimum mass for hydrogen burning in red dwarfs. The brown dwarf radius can differ significantly from the GR prediction and upcoming surveys that probe the mass-radius relation for stars with masses $<\mathcal{O}(0.1M_\odot)$ have the potential to place new constraints. The minimum mass for hydrogen burning can be larger than several presently observed Red Dwarf stars. This places a new and extremely stringent constraint on the parameters that appear in the effective field theory of dark energy and rules out several well-studied dark energy models.

Ghost-free theories beyond the Horndeski class exhibit a partial breaking of the Vainshtein mechanism inside non-relativistic sources of finite extent. We exploit this breaking to identify new and novel astrophysical probes of these theories. Non-relativistic objects feel a gravitational force that is weaker than that predicted by general relativity. The new equation of hydrostatic equilibrium equation is derived and solved to predict the modified behaviour of stars. It is found that main-sequence stars are dimmer and cooler than their general relativity counterparts but the red giant phase is largely indistinguishable. The rotation curves and lensing potential of Milky Way-like galaxies are calculated. The circular velocities are smaller than predicted by general relativity at fixed radius and the lensing mass is larger than the dynamical mass. We discuss potential astrophysical probes of these theories and identify strong lensing as a particularly promising candidate.

The most general scalar-tensor theories of gravity predict a weakening of the gravitational force inside astrophysical bodies. There is a minimum mass for hydrogen burning in stars that is set by the interplay of plasma physics and the theory of gravity. We calculate this for alternative theories of gravity, and find that it is always significantly larger than the general relativity prediction. The observation of several low mass Red Dwarf stars therefore rules out a large class of scalar-tensor gravity theories, and places strong constraints on the cosmological parameters appearing in the effective field theory of dark energy.

We study the quasi-normal modes of asymptotically anti-de Sitter black holes in a class of shift-symmetric Horndeski theories where a gravitational scalar is derivatively coupled to the Einstein tensor. The space-time differs from exact Schwarzschild-anti-de Sitter, resulting in a different effective potential for the quasi-normal modes and a different spectrum. We numerically compute this spectrum for a massless test scalar coupled both minimally to the metric, and non-minimally to the gravitational scalar. We find interesting differences from the Schwarzschild-anti-de Sitter black hole found in general relativity.

We put forward a new proposal for generating the baryon asymmetry of the universe by making use of the dynamics of a $\mathrm{U}(1)$ scalar field coupled to dark matter. High dark matter densities cause the $\mathrm{U}(1)$ symmetry to break spontaneously so that the field acquires a large vacuum expectation value. The symmetry is restored when the density redshifts below a critical value, resulting in the coherent oscillation of the scalar field. A net $B-L$ number can be generated either via baryon number-conserving couplings to the standard model or through small symmetry-violating operators and the subsequent decay of the scalar condensate.

Scalar-tensor theories of gravity generally violate the strong equivalence principle, namely compact objects have a suppressed coupling to the scalar force, causing them to fall slower. A black hole is the extreme example where such a coupling vanishes, i.e. black hole has no scalar hair. Following earlier work, we explore observational scenarios for detecting strong equivalence principle violation, focusing on galileon gravity as an example. For galaxies in-falling towards galaxy clusters, the supermassive black hole can be offset from the galaxy center away from the direction of the cluster. Hence, well resolved images of galaxies around nearby clusters can be used to identify the displaced black hole via the star cluster bound to it. We show that this signal is accessible with imaging surveys, both ongoing ones such as the Dark Energy Survey, and future ground and space based surveys. Already, the observation of the central black hole in M~87 places new constraints on the galileon parameters, which we present here. $\mathcal{O}(1)$ matter couplings are disfavored for a large region of the parameter space. We also find a novel phenomenon whereby the black hole can escape the galaxy completely in less than one billion years.

Lorentz-violating theories of gravity typically contain constrained vector fields. We show that the lowest-order coupling of such vectors to $\mathrm{U}(1)$-symmetric scalars can naturally give rise to baryogenesis in a manner akin to the Affleck-Dine mechanism. We calculate the cosmology of this new mechanism, demonstrating that a net $B-L$ can be generated in the early Universe, and that the resulting baryon-to-photon ratio matches that which is presently observed. We discuss constraints on the model using solar system and astrophysical tests of Lorentz violation in the gravity sector. Generic Lorentz-violating theories can give rise to the observed matter-antimatter asymmetry without violating any current bounds.

Theories of gravity in the beyond Horndeski class recover the predictions of general relativity in the solar system whilst admitting novel cosmologies, including late-time de Sitter solutions in the absence of a cosmological constant. Deviations from Newton's law are predicted inside astrophysical bodies, which allow for falsifiable, smoking-gun tests of the theory. In this work we study the pulsations of stars by deriving and solving the wave equation governing linear adiabatic oscillations to find the modified period of pulsation. Using both semi-analytic and numerical models, we perform a preliminary survey of the stellar zoo in an attempt to identify the best candidate objects for testing the theory. Brown dwarfs and Cepheid stars are found to be particularly sensitive objects and we discuss the possibility of using both to test the theory.

Theories of gravity in the beyond Horndeski class encompass a wide range of scalar-tensor theories that will be tested on cosmological scales over the coming decade. In this work, we investigate the possibility of testing them in the strong-field regime by looking at the properties of compact objects-neutron, hyperon, and quark stars-embedded in an asymptotically de Sitter space-time, for a specific subclass of theories. We extend previous works to include slow rotation and find a relation between the dimensionless moment of intertia, ($\bar{I}=Ic^2/G_{\rm N} M^3$), and the compactness, $\cal{C}=G_{\rm N} M/Rc^2$ (an $\bar{I}$-$\cal{C}$ relation), independent of the equation of state, that is reminiscent of but distinct from the general relativity prediction. Several of our equations of state contain hyperons and free quarks, allowing us to revisit the hyperon puzzle. We find that the maximum mass of hyperon stars can be larger than $2M_\odot$ for small values of the beyond Horndeski parameter, thus providing a resolution of the hyperon puzzle based on modified gravity. Moreover, stable quark stars exist when hyperonic stars are unstable, which means that the phase transition from hyperon to quark stars is predicted just as in general relativity, albeit with larger quark star masses. Two important and potentially observable consequences of some of the theories we consider are the existence of neutron stars in a range of masses significantly higher than in GR, and $\bar{I}$-$\mathcal{C}$ relations that differ from their GR counterparts. In the former case, we find objects that, if observed, could not be accounted for in GR because they violate the usual GR causality condition. We end by discussing several difficult technical issues that remain to be addressed in order to reach more realistic predictions that may be tested using gravitational wave searches or neutron star observations.

This work studies relativistic stars in beyond Horndeski scalar-tensor theories that exhibit a breaking of the Vainshtein mechanism inside matter, focusing on a model based on the quartic beyond Horndeski Lagrangian. We self-consistently derive the scalar field profile for static spherically symmetric objects in asymptotically de Sitter space-time and show that the Vainshtein breaking branch of the solutions is the physical branch thereby resolving several ambiguities with non-relativistic frameworks. The geometry outside the star is shown to be exactly Schwarzschild-de Sitter and therefore the PPN parameter $\beta_{\rm PPN}=1$, confirming that the external screening works at the post-Newtonian level. The Tolman-Oppenheimer-Volkoff (TOV) equations are derived and a new lower bound on the Vainshtein breaking parameter $\Upsilon_1>-4/9$ is found by requiring the existence of static spherically symmetric stars. Focusing on the unconstrained case where $\Upsilon_1<0$, we numerically solve the TOV equations for polytropic and realistic equations of state and find stars with larger radii at fixed mass. Furthermore, the maximum mass can increase dramatically and stars with masses in excess of $3M_\odot$ can be found for relatively small values of the Vainshtein breaking parameter. We re-examine white dwarf stars and show that post-Newtonian corrections are important in beyond Horndeski theories and therefore the bounds coming from previous analyses should be revisited.

The Beyond Horndeski class of alternative gravity theories allow for Self-accelerating de-Sitter cosmologies with no need for a cosmological constant. This makes them viable alternatives to $\Lambda$CDM and so testing their small-scale predictions against General Relativity is of paramount importance. These theories generically predict deviations in both the Newtonian force law and the gravitational lensing of light inside extended objects. Therefore, by simultaneously fitting the X-ray and lensing profiles of galaxy clusters new constraints can be obtained. In this work, we apply this methodology to the stacked profiles of 58 high-redshift ($ 0.1<z<1.2$) clusters using X-ray surface brightness profiles from the XMM Cluster Survey and weak lensing profiles from CFHTLenS. By performing a multi-parameter Markov chain Monte Carlo analysis, we are able to place new constraints on the parameters governing deviations from Newton's law $\Upsilon_{1}=-0.11^{+0.93}_{-0.67}$ and light bending $\Upsilon_{2}=-0.22^{+1.22}_{-1.19}$. Both constraints are consistent with General Relativity, for which $\Upsilon_{1}=\Upsilon_{2}=0$. We present here the first observational constraints on $\Upsilon_{2}$, as well as the first extragalactic measurement of both parameters.

The chameleon model is a scalar field theory with a screening mechanism that explains how a cosmologically relevant light scalar can avoid the constraints of intra-solar-system searches for fifth-forces. The chameleon is a popular dark energy candidate and also arises in $f(R)$ theories of gravity. Whilst the chameleon is designed to avoid historical searches for fifth-forces it is not unobservable and much effort has gone into identifying the best observables and experiments to detect it. These results are not always presented for the same models or in the same language, a particular problem when comparing astrophysical and laboratory searches making it difficult to understand what regions of parameter space remain. Here we present combined constraints on the chameleon model from astrophysical and laboratory searches for the first time and identify the remaining windows of parameter space. We discuss the implications for cosmological chameleon searches and future small-scale probes.

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.

The Jordan frame action for general disformal theories is presented and studied for the first time, motivated by several unresolved mysteries that arise when working in the Einstein frame. We present the Friedmann equations and, specialising to exponential functions, study the late-time cosmology using dynamical systems methods and by finding approximate solutions. Our analysis reveals that either the disformal effects are irrelevant or the universe evolves towards a phantom phase where the equation of state of dark energy is $-3$. There is a marginal case where the asymptotic state of the universe depends on the model parameters and de-Sitter solutions can be obtained. Our findings indicate that the metric singularity found using the Einstein frame construction corresponds phantom behaviour in the Jordan frame and we argue that this is the case for general disformal theories.

Disformal theories of gravity are scalar-tensor theories where the scalar couples derivatively to matter via the Jordan frame metric. These models have recently attracted interest in the cosmological context since they admit accelerating solutions. We derive the solution for a static isolated mass in generic disformal gravity theories and transform it into the parameterised post-Newtonian form. This allows us to investigate constraints placed on such theories by local tests of gravity. The tightest constraints come from preferred-frame effects due to the motion of the Solar System with respect to the evolving cosmological background field. The constraints we obtain improve upon the previous solar system constraints by two orders of magnitude, and constrain the scale of the disformal coupling for generic models to $\mathcal{M} \gtrsim 100$ eV. These constraints render all disformal effects irrelevant for cosmology.

Chameleon and similar (symmetron and dilation) theories of gravity can exhibit new and interesting features on cosmological scales whilst screening the modifications on small scales thereby satisfying solar system tests of general relativity. This thesis explores the regime between these two scales: astrophysics. The majority of this thesis is focused on discerning new and novel astrophysical probes of chameleon gravity in the form of stellar structure and oscillation tests. These are used to place new constraints on the theory parameters and the implications of these are discussed, as are the future prospects for improving them using planned future surveys. The final two chapters review supersymmetric completions of these theories.

The late-time cosmological dynamics of disformal gravity are investigated using dynamical systems methods. It is shown that in the general case there are no stable attractors that screen fifth-forces locally and simultaneously describe a dark energy dominated universe. Viable scenarios have late-time properties that are independent of the disformal parameters and are identical to the equivalent conformal quintessence model. Our analysis reveals that configurations where the Jordan frame metric becomes singular are only reached in the infinite future, thus explaining the natural pathology resistance observed numerically by several previous works. The viability of models where this can happen is discussed in terms of both the cosmological dynamics and local phenomena. We identify a special parameter tuning such that there is a new fixed point that can match the presently observed dark energy density and equation of state. This model is unviable when the scalar couples to the visible sector but may provide a good candidate model for theories where only dark matter is disformally coupled.

This paper is concerned with theories of gravity that contain a scalar coupled both conformally and disformally to matter through the metric. By systematically deriving the non-relativistic limit, it is shown that no new non-linear screening mechanisms are present beyond the Vainshtein mechanism and chameleon-like screening. If one includes the cosmological expansion of the universe, disformal effects that are usually taken to be absent can be present in the solar system. When the conformal factor is absent, fifth-forces can be screened on all scales when the cosmological field is slowly-rolling. We investigate the cosmology of these models and use local tests of gravity to place new constraints on the disformal coupling and find $\mathcal{M}>\mathcal{O}(\textrm{eV})$, which is not competitive with laboratory tests. Finally, we discuss the future prospects for testing these theories and the implications for other theories of modified gravity. In particular, the Vainshtein radius of solar system objects can be altered from the static prediction when cosmological time-derivatives are non-negligible.

Modified theories of gravity have received a renewed interest due to their ability to account for the cosmic acceleration. In order to satisfy the solar system tests of gravity, these theories need to include a screening mechanism that hides the modifications on small scales. One popular and well-studied theory is chameleon gravity. Our own galaxy is necessarily screened, but less dense dwarf galaxies may be unscreened and their constituent stars can exhibit novel features. In particular, unscreened stars are brighter, hotter and more ephemeral than screened stars in our own galaxy. They also pulsate with a shorter period. In this essay, we exploit these new features to constrain chameleon gravity to levels three orders of magnitude lower the previous measurements. These constraints are currently the strongest in the literature.

We calculate the rate of energy loss from compact astrophysical objects due to a scalar field in screened modified gravity models of the chameleon, dilaton and symmetron types. The cosmological evolution of the field results in a time-variation of the scalar charge of screened objects implying the emission of scalar radiation. Focusing on binary objects, this leads to an additional decay in the orbital period complementing that due to the emission of gravitational waves. Using the Hulse-Taylor binary pulsar, the double pulsar PSR J0737-3039 and the pulsar-white dwarf system PSR J1738+033, we find a new observational bound on the time variation of the scalar charge of the earth in the Milky Way. We then translate this into a new bound on the range of the scalar interaction in the Milky Way. Ultimately, we find that pulsar tests are not competitive with current observational constraints.

Starting from the equations of modified gravity hydrodynamics, we derive the equ tions of motion governing linear, adiabatic, radial perturbations of stars in scalar-tensor theories. There are two new features: first, the eigenvalue equation for the period of stellar oscillations is modified such that the eigenfrequencies are always larger than predicted by General Relativity. Second, the General Relativity condition for stellar instability is altered so that the adiabatic index can fall below 4/3 before unstable modes appear. Stars are more stable in modified gravity theories. Specialising to the case of chameleon-like theories, we investigate these effects numerically using both polytropic Lane-Emden stars and models coming from modified gravity stellar structure simulations. The change in the oscillation period can be as large as 50% and the critical adiabatic index for instability falls by a composition dependent amount of order 10^(-1). By solving the new equation for Cepheid models, it is found that the change in the inferred distance using the period-luminosity relation can be up to three times larger than if one had only considered the modified equilibrium structure. We discuss the implications of these results for recent and up-coming astrophysical tests and estimate that previous methods can produce new constraints such that the modifications are screened in regions of Newtonian potential of order 10^(-8).

We investigate the cosmological dynamics of a class of supersymmetric chameleon models coupled to cold dark matter fermions. Supergravity corrections ensure that these models are efficiently screened in all astrophysical objects of interest, however this does not preclude the enhancement of gravity on linear cosmological scales. We analyse the background cosmology and solve the modified equations for the growth of cold dark matter density perturbations in closed form. Using this, we go on to derive the modified linear power spectrum which is characterised by two scales, the horizon size at matter-radiation equality and at the redshift when the chameleon reaches the minimum of its effective potential. The model includes a cosmological constant in the form of a Fayet-Illiopolous term, which emerges at late times due to the coupling of the chameleon to two charged scalars. We examine the conditions under which this leads to viable background cosmology and go on to analyse the deviations from the LCDM predictions in the linear regime. We find that for reasonable values in the model's parameter space there is generically a region where the model's cosmology is viable and current measurements can be reproduced. A small discrepancy of the matter power spectrum from its LCDM counterpart can be obtained in a smaller subset of the parameter space.

We present a framework for embedding scalar-tensor models of screened modifed gravity such as chameleons, symmetrons and environmental dilatons into global supersymmetry. This achieved by secluding the dark sector from both the observable and supersymmetry breaking sectors. We examine the resulting supersymmetric features in a model-independent manner and find that, when the theory follows from an underlying supergravity, the mediation of supersymmetry breaking to the dark sector induces a soft mass for the scalar of order the gravitino mass. This is enough to forbid the construction of supersymmetric symmetrons and ensures that when other screening mechanisms operate, no object in the universe is unscreened thereby precluding any observable signatures. In view of a possible origin of modifed gravity within fundamental physics, we find that no-scale models are the only ones that can circumvent these features. We also present a novel mechanism where the coupling of the scalar to two other scalars charged under U(1) can dynamically generate a small cosmological constant at late times in the form of a Fayet-Iliopoulos term.

We use distance measurements in the nearby universe to carry out new tests of gravity, surpassing other astrophysical tests by over two orders of magnitude for chameleon theories. The three nearby distance indicators -- cepheids, tip of the red giant branch (TRGB) stars, and water masers -- operate in gravitational fields of widely different strengths. This enables tests of scalar-tensor gravity theories because they are screened from enhanced forces to different extents. Inferred distances from cepheids and TRGB stars are altered (in opposite directions) over a range of chameleon gravity theory parameters well below the sensitivity of cosmological probes. Using published data we have compared cepheid and TRGB distances in a sample of unscreened dwarf galaxies within 10 Mpc. As a control sample we use a comparable set of screened galaxies. We find no evidence for the order unity force enhancements expected in these theories. Using a two-parameter description of the models (the coupling strength and background field value) we obtain constraints on chameleon and symmetron screening scenarios. In particular we show that f(R) models with background field values fR0 above 5e-7 are ruled out at the 95% confidence level. We also compare TRGB and maser distances to the galaxy NGC 4258 as a second test for larger field values. While there are several approximations and caveats in our study, our analysis demonstrates the power of gravity tests in the local universe. We discuss the prospects for additional, improved tests with future observations.

We investigate the effect of modifed gravity with screening mechanisms, such as the chameleon or symmetron models, upon the structure of main sequence stars. We find that unscreened stars can be significantly more luminous and ephemeral than their screened doppelgangers. By embedding these stars into dwarf galaxies, which can be unscreened for values of the parameters not yet ruled out observationally, we show that the cumulative effect of their increased luminosity can enhance the total galactic luminosity. We estimate this enhancement and find that it can be considerable given model parameters that are still under experimental scrutiny. By looking for systematic offsets between screened dwarf galaxies in clusters and unscreened galaxies in voids, these effects could form the basis of an independent observational test that can potentially lower the current experimental bounds on the model independent parameters of these theories by and order of magnitude or more.

The GalICS (Galaxies in Cosmological Simulations) semi-analytical model of hierar- chical galaxy formation is used to investigate the effects of different galactic properties, including star formation rate (SFR) and outflows, on the shape of the mass metallic- ity relation and to predict the relation for galaxies at redshift z=2.27 and z=3.54. Our version of GalICS has the chemical evolution implemented in great detail and is less heavily reliant on approximations such as instantaneous recycling. We vary the model parameters controlling both the efficiency and redshift dependence of the SFR as well as the efficiency of supernova feedback. We find that the factors controlling the SFR influence the relation significantly at all redshifts and require a strong redshift dependence, proportional to 1+z, in order to reproduce the observed relation at the low mass end. Indeed, at any redshift, the predicted relation flattens out at the high mass end resulting in a poorer agreement with observations in this regime. We also find that variation of the parameters associated with outflows has a minimal effect on the relation at high redshift but does serve to alter its shape in the more recent past. We thus conclude that the relation is one between SFR and mass and that outflows are only important in shaping the relation at late times. When the relation is stratified by SFR it is apparent that the predicted galaxies with increasing stellar masses have higher SFRs, supporting the view that galaxy downsizing is the origin of the relation. Attempting to reproduce the observed relation, we vary the parameters controlling the efficiency of star formation and its redshift dependence and compare the predicted relations with Erb et al. (2006) at z=2.27 and Maiolino et al. (2008) at z=3.54 in order to find the best-fitting parameters. (Abridged)