After a decade and a half of research motivated by the accelerating universe, theory and experiment have a reached a certain level of maturity. The development of theoretical models beyond Λ or smooth dark energy, often called modified gravity, has led to broader insights into a path forward, and a host of observational and experimental tests have been developed. In this review we present the current state of the field and describe a framework for anticipating developments in the next decade. We identify the guiding principles for rigorous and consistent modifications of the standard model, and discuss the prospects for empirical tests.We begin by reviewing recent attempts to consistently modify Einstein gravity in the infrared, focusing on the notion that additional degrees of freedom introduced by the modification must "screen" themselves from local tests of gravity. We categorize screening mechanisms into three broad classes: mechanisms which become active in regions of high Newtonian potential, those in which first derivatives of the field become important, and those for which second derivatives of the field are important. Examples of the first class, such as f (R) gravity, employ the familiar chameleon or symmetron mechanisms, whereas examples of the last class are galileon and massive gravity theories, employing the Vainshtein mechanism. In each case, we describe the theories as effective theories and discuss prospects for completion in a more fundamental theory. We describe experimental tests of each class of theories, summarizing laboratory and solar system tests and describing in some detail astrophysical and cosmological tests. Finally, we discuss prospects for future tests which will be sensitive to different signatures of new physics in the gravitational sector.The review is structured so that those parts that are more relevant to theorists vs. observers/experimentalists are clearly indicated, in the hope that this will serve as a useful reference for both audiences, as well as helping those interested in bridging the gap between them.

Understanding the reason for the observed accelerated expansion of the Universe represents one of the fundamental open questions in physics. In cosmology, a classification has emerged among physical models for the acceleration, distinguishing between Dark Energy and Modified Gravity. In this review, we give a brief overview of models in both categories as well as their phenomenology and characteristic observable signatures in cosmology. We also introduce a rigorous distinction between Dark Energy and Modified Gravity based on the strong and weak equivalence principles.

We show that the galileons can be thought of as Wess-Zumino terms for the spontaneous breaking of space-time symmetries. Wess-Zumino terms are terms which are not captured by the coset construction for phenomenological Lagrangians with broken symmetries. Rather they are, in d space-time dimensions, d-form potentials for (d+1)-forms which are non-trivial co-cycles in Lie algebra cohomology of the full symmetry group relative to the unbroken symmetry group. We introduce the galileon algebras and construct the non-trivial (d + 1)form co-cycles, showing that the presence of galileons and multi-galileons in all dimensions is counted by the dimensions of particular Lie algebra cohomology groups. We also discuss the DBI and conformal galileons from this point of view, showing that they are not Wess-Zumino terms, with one exception in each case.

A key insight of the bootstrap approach to cosmological correlations is the fact that all correlators of slow-roll inflation can be reduced to a unique building block — the four-point function of conformally coupled scalars, arising from the exchange of a massive scalar. Correlators corresponding to the exchange of particles with spin are then obtained by applying a spin-raising operator to the scalar-exchange solution. Similarly, the correlators of massless external fields can be derived by acting with a suitable weight-raising operator. In this paper, we present a systematic and highly streamlined derivation of these operators (and their generalizations) using tools of conformal field theory. Our results greatly simplify the theoretical foundations of the cosmological bootstrap program.

We show that there is a special choice of parameters for which the galileon theory is invariant under an enhanced shift symmetry whose non-linear part is quadratic in the coordinates. This symmetry fixes the theory to be equivalent to one with only even powers of the field, with no free coefficients, and accounts for the improved soft limit behavior observed in the quartic galileon S-matrix.

We extend the cosmological bootstrap to correlators involving massless spinning particles, focusing on spin-1 and spin-2. In de Sitter space, these correlators are constrained both by symmetries and by locality. In particular, the de Sitter isometries become conformal symmetries on the future boundary of the spacetime, which are reflected in a set of Ward identities that the boundary correlators must satisfy. We solve these Ward identities by acting with weight-shifting operators on scalar seed solutions. Using this weight-shifting approach, we derive three- and four-point correlators of massless spin-1 and spin-2 fields with conformally coupled scalars. Four-point functions arising from tree-level exchange are singular in particular kinematic configurations, and the coefficients of these singularities satisfy certain factorization properties. We show that in many cases these factorization limits fix the structure of the correlators uniquely, without having to solve the conformal Ward identities. The additional constraint of locality for massless spinning particles manifests itself as current conservation on the boundary. We find that the four-point functions only satisfy current conservation if the s, t, and u-channels are related to each other, leading to nontrivial constraints on the couplings between the conserved currents and other operators in the theory. For spin-1 currents this implies charge conservation, while for spin-2 currents we recover the equivalence principle from a purely boundary perspective. For multiple spin-1 fields, we recover the structure of Yang--Mills theory. Finally, we apply our methods to slow-roll inflation and derive a few phenomenologically relevant scalar-tensor three-point functions.

We derive the quadratic action for the physical degrees of freedom of massless spin-0, spin-1, and spin-2 perturbations on a Schwarzschild-(A)dS background in arbitrary dimensions. We then use these results to compute the static response of asymptotically flat Schwarzschild black holes to external fields. Our analysis reproduces known facts about black hole Love numbers — in particular that they vanish for all types of perturbation in four spacetime dimensions — but also leads to new results. For instance, we find that neutral Schwarzschild black holes polarize in the presence of an electromagnetic background in any number of spacetime dimensions except four. Moreover, we calculate for the first time black hole Love numbers for vector-type gravitational perturbations in higher dimensions and find that they generically do not vanish. Along the way, we shed some light on an apparent discrepancy between previous results in the literature, and clarify some aspects of the matching between perturbative calculations of static response on a Schwarzschild background and the point-particle effective theory.

We show that the Dirac-Born-Infeld conformal galileons, derived from the world-volume theory of a 3-brane moving in an anti-de Sitter bulk, admit a background, stable under quantum corrections, which violates the null energy condition. The perturbations around this background are stable and propagate subluminally. Unlike other known examples of null energy condition violation, such as ghost condensation and conformal galileons, this theory also admits a stable, Poincaré-invariant vacuum. The 2 → 2 amplitude satisfies standard analyticity conditions. The full S matrix is likely not analytic, however, since perturbations around deformations of the Poincaré invariant vacuum propagate superluminally.