
We study the perturbations to General Relativistic black holes (i.e. those
without scalar hair) in Horndeski scalartensor gravity. First, we derive the
equations of odd and even parity perturbations of both the metric and scalar
field in the case of a Schwarzschild black hole, and show that the
gravitational waves emitted from such a system contain a mixture of
quasinormal mode frequencies from the usual General Relativistic spectrum and
those from the new scalar field spectrum, with the new scalar spectrum
characterised by just two free parameters. We then specialise to the subfamily
of Horndeski theories in which gravitational waves propagate at the speed of
light $c$ on cosmological backgrounds; the scalar quasinormal mode spectrum of
such theories is characterised by just a single parameter $\mu$ acting as an
effective mass of the scalar field. Analytical expressions for the quasinormal
mode frequencies of the scalar spectrum in this subfamily of theories are
provided for both static and slowly rotating black holes. In both regimes
comparisons to quasinormal modes calculated numerically show good agreement
with those calculated analytically in this work.

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

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

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

Future surveys will access large volumes of space and hence very long
wavelength fluctuations of the matter density and gravitational field. It has
been argued that the set of secondary effects that affect the galaxy
distribution, relativistic in nature, will bring new, complementary
cosmological constraints. We study this claim in detail by focusing on a subset
of widearea future surveys: Stage4 cosmic microwave background experiments
and photometric redshift surveys. In particular, we look at the magnification
lensing contribution to galaxy clustering and general relativistic corrections
to all observables. We quantify the amount of information encoded in these
effects in terms of the tightening of the final cosmological constraints as
well as the potential bias in inferred parameters associated with neglecting
them. We do so for a wide range of cosmological parameters, covering neutrino
masses, standard darkenergy parametrizations and scalartensor gravity
theories. Our results show that, while the effect of lensing magnification to
number counts does not contain a significant amount of information when galaxy
clustering is combined with cosmic shear measurements, this contribution does
play a significant role in biasing estimates on a host of parameter families if
unaccounted for. Since the amplitude of the magnification term is controlled by
the slope of the source number counts with apparent magnitude, $s(z)$, we also
estimate the accuracy to which this quantity must be known to avoid systematic
parameter biases, finding that future surveys will need to determine $s(z)$ to
the $\sim$510\% level. On the contrary, largescale generalrelativistic
corrections are irrelevant both in terms of information content and parameter
bias for most cosmological parameters, but significant for the level of
primordial nonGaussianity.

We consider the cosmological dynamics of a scalar field in a potential with
multiple troughs and peaks. We show that the dynamics of the scalar field will
evolve from light dark matterlike behaviour (such as that of a light axion) to
a combination of heavy dark matterlike and dark energylike behaviour. We
discuss the phenomenology of such a model, explaining how it might lead to a
small cosmological constant, as well as how it can decouple the dark sector
densities between the time of recombination and today, for both the homogeneous
background and perturbations. The final form of the dark matter is axionlike,
but with abundance and primordial isocurvature modes taking very different
values from traditional, axionic, dark matter.

Weyl invariant theories of scalars and gravity can generate all mass scales
spontaneously, initiated by a dynamical process of "inertial spontaneous
symmetry breaking" that does not involve a potential. This is dictated by the
structure of the Weyl current, $K_\mu$, and a cosmological phase during which
the universe expands and the EinsteinHilbert effective action is formed.
Maintaining exact Weyl invariance in the renormalised quantum theory is
straightforward when renormalisation conditions are referred back to the VEV's
of fields in the action of the theory, which implies a conserved Weyl current.
We do not require scale invariant regulators. We illustrate the computation of
a Weyl invariant ColemanWeinberg potential.

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

Tests of gravity at the galaxy scale are in their infancy. As a first step to
systematically uncovering the gravitational significance of galaxies, we map
three fundamental gravitational variables  the Newtonian potential,
acceleration and curvature  over the galaxy environments of the local
universe to a distance of approximately 200 Mpc. Our method combines the
contributions from galaxies in an allsky redshift survey, halos from an Nbody
simulation hosting lowluminosity objects, and linear and quasilinear modes of
the density field. We use the ranges of these variables to determine the extent
to which galaxies expand the scope of generic tests of gravity and are capable
of constraining specific classes of model for which they have special
significance. Finally, we investigate the improvements afforded by upcoming
galaxy surveys.

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

We use KaluzaKlein reduction to construct flatspace massive spin2
Lagrangians based on a kinetic term that has local Weyl and
transversediffeomorphism gauge symmetries in the massless limit. This yields
Lagrangians describing a free massive spin 2 which differ from the usual
FierzPauli theory, but are physically equivalent to it. These Lagrangians
require the use of auxiliary fields, which appear naturally from the
higherdimensional construction. We discuss how these Lagrangians are related
to each other and to the FierzPauli theory through St\"uckelberg
transformations, gauge fixings, and field replacements, and we use this to
generalize the construction to spin2 fields on Einstein manifolds.

We discuss the groundbreaking science that will be possible with a wide area
survey, using the MeerKAT telescope, known as MeerKLASS (MeerKAT Large Area
Synoptic Survey). The current specifications of MeerKAT make it a great fit for
science applications that require large survey speeds but not necessarily high
angular resolutions. In particular, for cosmology, a large survey over $\sim
4,000 \, {\rm deg}^2$ for $\sim 4,000$ hours will potentially provide the first
ever measurements of the baryon acoustic oscillations using the 21cm intensity
mapping technique, with enough accuracy to impose constraints on the nature of
dark energy. The combination with multiwavelength data will give unique
additional information, such as exquisite constraints on primordial
nonGaussianity using the multitracer technique, as well as a better handle on
foregrounds and systematics. Such a wide survey with MeerKAT is also a great
match for HI galaxy studies, providing unrivalled statistics in the preSKA era
for galaxies resolved in the HI emission line beyond local structures at z >
0.01. It will also produce a large continuum galaxy sample down to a depth of
about 5\,$\mu$Jy in Lband, which is quite unique over such large areas and
will allow studies of the largescale structure of the Universe out to high
redshifts, complementing the galaxy HI survey to form a transformational
multiwavelength approach to study galaxy dynamics and evolution. Finally, the
same survey will supply unique information for a range of other science
applications, including a large statistical investigation of galaxy clusters as
well as produce a rotation measure map across a huge swathe of the sky. The
MeerKLASS survey will be a crucial step on the road to using SKA1MID for
cosmological applications and other commensal surveys, as described in the top
priority SKA key science projects (abridged).

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

We present the public version of hi_class (www.hiclasscode.net), an
extension of the Boltzmann code CLASS to a broad ensemble of modifications to
general relativity. In particular, hi_class can calculate predictions for
models based on Horndeski's theory, which is the most general scalartensor
theory described by secondorder equations of motion and encompasses any
perfectfluid dark energy, quintessence, BransDicke, $f(R)$ and covariant
Galileon models. hi_class has been thoroughly tested and can be readily used to
understand the impact of alternative theories of gravity on linear structure
formation as well as for cosmological parameter extraction.

Imaging surveys of galaxies will have a high number density and angular
resolution yet a poor redshift precision. Intensity maps of neutral hydrogen
(HI) will have accurate redshift resolution yet will not resolve individual
sources. Using this complementarity, we show how the clustering redshifts
approach, proposed for spectroscopic surveys can also be used in combination
with intensity mapping observations to calibrate the redshift distribution of
galaxies in an imaging survey and, as a result, reduce uncertainties in
photometric redshift measurements. We show how the intensity mapping surveys to
be carried out with the MeerKAT, HIRAX and SKA instruments can improve
photometric redshift uncertainties to well below the requirements of DES and
LSST. The effectiveness of this method as a function of instrumental
parameters, foreground subtraction and other potential systematic errors is
discussed in detail.

The next generation of surveys will greatly improve our knowledge of
cosmological gravity. In this paper we focus on how Stage IV photometric
redshift surveys, including weak lensing and multiple tracers of the matter
distribution and radio experiments combined with measurements of the cosmic
microwave background will lead to precision constraints on deviations from
General Relativity. We use a broad subclass of Horndeski scalartensor theories
to forecast the accuracy with which we will be able to determine these
deviations and their degeneracies with other cosmological parameters. Our
analysis includes relativistic effects, does not rely on the quasistatic
evolution and makes conservative assumptions about the effect of screening on
small scales. We define a figure of merit for cosmological tests of gravity and
show how the combination of different types of surveys, probing different
length scales and redshifts, can be used to pin down constraints on the
gravitational physics to better than a few percent, roughly an order of
magnitude better than present probes. Future cosmological experiments will be
able to constrain the BransDicke parameter at a level comparable to Solar
System and astrophysical tests

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

Scalar fields, $\phi_i$ can be coupled nonminimally to curvature and satisfy
the general criteria: (i) the theory has no mass input parameters, including
the Planck mass; (ii) the $\phi_i$ have arbitrary values and gradients, but
undergo a general expansion and relaxation to constant values that satisfy a
nontrivial constraint, $K(\phi_i) =$ constant; (iii) this constraint breaks
scale symmetry spontaneously, and the Planck mass is dynamically generated;
(iv) there can be adequate inflation associated with slow roll in a scale
invariant potential subject to the constraint; (v) the final vacuum can have a
small to vanishing cosmological constant (vi) large hierarchies in vacuum
expectation values can naturally form; (vii) there is a harmless dilaton which
naturally eludes the usual constraints on massless scalars. These models are
governed by a global Weyl scale symmetry and its conserved current, $K_\mu$ .
At the quantum level the Weyl scale symmetry can be maintained by an invariant
specification of renormalized quantities.

We revisit the possibility that the Planck mass is spontaneously generated in
scale invariant scalartensor theories of gravity, typically leading to a
"dilaton." The fifth force, arising from the dilaton, is severely constrained
by astrophysical measurements. We explore the possibility that nature is
fundamentally Weylscale invariant and argue that, as a consequence, the fifth
force effects are dramatically suppressed and such models are viable. We
discuss possible obstructions to maintaining scale invariance and how these
might be resolved.

We discuss models involving two scalar fields coupled to classical gravity
that satisfy the general criteria: (i) the theory has no mass input parameters,
(ii) classical scale symmetry is broken only through $\frac{1}{12}\varsigma
\phi^2 R$ couplings where $\varsigma$ departs from the special conformal value
of $1$; (iii) the Planck mass is dynamically generated by the vacuum
expectations values (VEVs) of the scalars (iv) there is a stage of viable
inflation associated with slow roll in the twoscalar potential; (v) the final
vacuum has a small to vanishing cosmological constant and an hierarchically
small ratio of the VEVs and the ratio of the scalar masses to the Planck scale.
This assumes the paradigm of classical scale symmetry as a custodial symmetry
of large hierarchies.

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

Euclid is a European Space Agency medium class mission selected for launch in
2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to
understand the origin of the accelerated expansion of the universe. Euclid will
explore the expansion history of the universe and the evolution of cosmic
structures by measuring shapes and redshifts of galaxies as well as the
distribution of clusters of galaxies over a large fraction of the sky. Although
the main driver for Euclid is the nature of dark energy, Euclid science covers
a vast range of topics, from cosmology to galaxy evolution to planetary
research. In this review we focus on cosmology and fundamental physics, with a
strong emphasis on science beyond the current standard models. We discuss five
broad topics: dark energy and modified gravity, dark matter, initial
conditions, basic assumptions and questions of methodology in the data
analysis. This review has been planned and carried out within Euclid's Theory
Working Group and is meant to provide a guide to the scientific themes that
will underlie the activity of the group during the preparation of the Euclid
mission.

Future groundbased CMB experiments will generate competitive largescale
structure datasets by precisely characterizing CMB secondary anisotropies over
a large fraction of the sky. We describe a method for constraining the growth
rate of structure to sub1% precision out to $z\approx 1$, using a combination
of galaxy cluster peculiar velocities measured using the kinetic
SunyaevZel'dovich (kSZ) effect, and the velocity field reconstructed from
galaxy redshift surveys. We consider only thermal SZselected cluster samples,
which will consist of $\mathcal{O}(10^410^5)$ sources for Stage 3 and 4 CMB
experiments respectively. Three different methods for separating the kSZ effect
from the primary CMB are compared, including a novel blind "constrained
realization" method that improves signaltonoise by a factor of $\sim 2$ over
a commonlyused aperture photometry technique. Measurements of the integrated
tSZ $y$parameter are used to break the kSZ velocityoptical depth degeneracy,
and the effects of including CMB polarization and SZ profile uncertainties are
also considered. A combination of future Stage 4 experiments should be able to
measure the product of the growth and expansion rates, $\alpha\equiv f H$, to
better than 1% in bins of $\Delta z = 0.1$ out to $z \approx 1$  competitive
with contemporary redshiftspace distortion constraints from galaxy surveys.

We measure the redshiftspace correlation function from a spectroscopic
sample of 2783 emission line galaxies from the FastSound survey. The survey,
which uses the Subaru Telescope and covers the redshift ranges of
$1.19<z<1.55$, is the first cosmological study at such high redshifts. We
detect clear anisotropy due to redshiftspace distortions (RSD) both in the
correlation function as a function of separations parallel and perpendicular to
the line of sight and its quadrupole moment. RSD has been extensively used to
test general relativity on cosmological scales at $z<1$. Adopting a LCDM
cosmology with the fixed expansion history and no velocity dispersion
$\sigma_{\rm v}=0$, and using the RSD measurements on scales above 8Mpc/h, we
obtain the first constraint on the growth rate at the redshift,
$f(z)\sigma_8(z)=0.482\pm 0.116$ at $z\sim 1.4$ after marginalizing over the
galaxy bias parameter $b(z)\sigma_8(z)$. This corresponds to $4.2\sigma$
detection of RSD. Our constraint is consistent with the prediction of general
relativity $f\sigma_8\sim 0.392$ within the $1\sigma$ confidence level. When
we allow $\sigma_{\rm v}$ to vary and marginalize it over, the growth rate
constraint becomes $f\sigma_8=0.494^{+0.126}_{0.120}$. We also demonstrate
that by combining with the lowz constraints on $f\sigma_8$, highz galaxy
surveys like the FastSound can be useful to distinguish modified gravity models
without relying on CMB anisotropy experiments.

We present a complete derivation of the observationally motivated definition
of the modified gravity statistic $E_G$. Using this expression, we investigate
how variations to theory and survey parameters may introduce uncertainty in the
general relativistic prediction of $E_G$. We forecast errors on $E_G$ for
measurements using two combinations of upcoming surveys, and find that
theoretical uncertainties may dominate for a futuristic measurement. Finally,
we compute predictions of $E_G$ under modifications to general relativity in
the quasistatic regime, and comment on the pros and cons of using $E_G$ to test
gravity with future surveys.