
A Hamiltonian density bounded from below implies that the lowestenergy state
is stable. We point out, contrary to common lore, that an unbounded Hamiltonian
density does not necessarily imply an instability: this is a
coordinatedependent statement. We discuss this fact in detail, starting from
kessence and extending our discussion to general field theories. We give the
correct stability criterion, using the relative orientation of the causal cones
for all propagating degrees of freedom. We then apply this criterion to an
exact Schwarzschildde Sitter solution of a beyondHorndeski theory, while
taking into account the recent experimental constraint regarding the speed of
gravitational waves. We extract the spin2 and spin0 causal cones by analyzing
respectively all the oddparity and the $\ell = 0$ evenparity modes. Contrary
to a claim in the literature, we prove that this solution does not exhibit any
kinetic instability for a given range of parameters defining the theory.

Both kessence and the pressureless perfect fluid develop caustic
singularities at finite time. We further explore the connection between the two
and show that they belong to the same class of models, which admits the caustic
free completion by means of the canonical complex scalar field. Specifically,
the free massive/selfinteracting complex scalar reproduces dynamics of
pressureless perfect fluid/shiftsymmetric kessence under certain initial
conditions in the limit of large mass/sharp selfinteracting potential. We
elucidate a mechanism of resolving caustic singularities in the complete
picture. The collapse time is promoted to complex number. Hence, the
singularity is not developed in real time. The same conclusion holds for a
collection of collisionless particles modelled by means of the Schroedinger
equation, or ultralight axions (generically, coherent oscillations of bosons
in the BoseEinstein condensate state).

We consider a generic scalartensor theory involving a shiftsymmetric scalar
field and minimally coupled matter fields. We prove that the Noether current
associated with shiftsymmetry vanishes in regular, spherically symmetric and
static spacetimes. We use this fact to prove the absence of scalar hair for
spherically symmetric and static stars in Horndeski and beyond theories. We
carefully detail the validity of this nohair theorem.

We study a massive vector tensor theory that acquires mass via a standard
Proca term but also via a higher order term containing an explicit coupling to
curvature. We find static solutions that are asymptotically flat, adS or
Lifshitz black holes. Since the higher order term regularizes the effect of the
Proca mass term, generically solutions are asymptotically regular for arbitrary
couplings. This is true in particular for asymptotically flat black holes. For
a particular coupling we find particle like solitons that have a regular and
non trivial geometry everywhere. In all adS solutions the Proca mass term plays
the role of an effective cosmological constant distinctly different from the
bare cosmological constant.

Theories of gravity in the beyond Horndeski class encompass a wide range of
scalartensor theories that will be tested on cosmological scales over the
coming decade. In this work, we investigate the possibility of testing them in
the strongfield regime by looking at the properties of compact
objectsneutron, hyperon, and quark starsembedded in an asymptotically de
Sitter spacetime, 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.

We find spherically symmetric and static black holes in shiftsymmetric
Horndeski and beyond Horndeski theories. They are asymptotically flat and
sourced by a non trivial static scalar field. The first class of solutions is
constructed in such a way that the Noether current associated with shift
symmetry vanishes, while the scalar field cannot be trivial. This in certain
cases leads to hairy black hole solutions (for the quartic Horndeski
Lagrangian), and in others to singular solutions (for a GaussBonnet term).
Additionally, we find the general spherically symmetric and static solutions
for a pure quartic Lagrangian, the metric of which is Schwarzschild. We show
that under two requirements on the theory in question, any vacuum GR solution
is also solution to the quartic theory. As an example, we show that a Kerr
black hole with a nontrivial scalar field is an exact solution to these
theories.

Motivated by the projectable HoravaLifshitz model/mimetic matter scenario,
we consider a particular modification of standard gravity, which manifests as
an imperfect low pressure fluid. While practically indistinguishable from a
collection of nonrelativistic weakly interacting particles on cosmological
scales, it leaves drastically different signatures in the Solar system. The
main effect stems from gravitational focusing of the flow of Imperfect Dark
Matter passing near the Sun. This entails strong amplification of Imperfect
Dark Matter energy density compared to its average value in the surrounding
halo. The enhancement is many orders of magnitude larger than in the case of
Cold Dark Matter, provoking deviations of the metric in the second order in the
Newtonian potential. Effects of gravitational focusing are prominent enough to
substantially affect the planetary dynamics. Using the existing bound on the
PPN parameter $\beta_{PPN}$, we deduce a stringent constraint on the unique
constant of the model.

We study both the cosmological selftuning and the local predictions (inside
the Solar system) of the most general shiftsymmetric beyond Horndeski theory.
We first show that the cosmological selftuning is generic in this class of
theories: By adjusting a mass parameter entering the action, a large bare
cosmological constant can be effectively reduced to a small observed one.
Requiring then that the metric should be close enough to the Schwarzschild
solution in the Solar system, to pass the experimental tests of general
relativity, and taking into account the renormalization of Newton's constant,
we select a subclass of models which presents all desired properties: It is
able to screen a big vacuum energy density, while predicting an exact
Schwarzschildde Sitter solution around a static and spherically symmetric
source. As a byproduct of our study, we identify a general subclass of beyond
Horndeski theory for which regular selftuning black hole solutions exist, in
presence of a timedependent scalar field. We discuss possible future
development of the present work.

This work studies relativistic stars in beyond Horndeski scalartensor
theories that exhibit a breaking of the Vainshtein mechanism inside matter,
focusing on a model based on the quartic beyond Horndeski Lagrangian. We
selfconsistently derive the scalar field profile for static spherically
symmetric objects in asymptotically de Sitter spacetime and show that the
Vainshtein breaking branch of the solutions is the physical branch thereby
resolving several ambiguities with nonrelativistic frameworks. The geometry
outside the star is shown to be exactly Schwarzschildde Sitter and therefore
the PPN parameter $\beta_{\rm PPN}=1$, confirming that the external screening
works at the postNewtonian level. The TolmanOppenheimerVolkoff (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 reexamine white dwarf stars and show that
postNewtonian corrections are important in beyond Horndeski theories and
therefore the bounds coming from previous analyses should be revisited.

We consider generic linear perturbations of a nonbidiagonal class of static
blackhole solutions in massive (bi)gravity. We show that the quasinormal
spectrum of these solutions coincides with that of a Schwarzschild black hole
in general relativity, thus proving that these solutions are mode stable. This
is in contrast to the case of bidiagonal blackhole solutions which are
affected by a radial instability. On the other hand, the full set of
perturbation equations is generically richer than that of a Schwarzschild black
hole in general relativity, and this affects the linear response of the black
hole to external perturbations. Finally, we argue that the generalization of
these solutions to the spinning case does not suffer from the superradiant
instability, despite the fact that the theory describes a massive graviton.

We provide further details on a recent proposal addressing the nature of the
dark sectors in cosmology and demonstrate that all current observations related
to Dark Matter can be explained by the presence of a heavy spin2 particle.
Massive spin2 fields and their gravitational interactions are uniquely
described by ghostfree bimetric theory, which is a minimal and natural
extension of General Relativity. In this setup, the largeness of the physical
Planck mass is naturally related to extremely weak couplings of the heavy
spin2 field to baryonic matter and therefore explains the absence of signals
in experiments dedicated to Dark Matter searches. It also ensures the
phenomenological viability of our model as we confirm by comparing it with
cosmological and local tests of gravity. At the same time, the spin2 field
possesses standard gravitational interactions and it decays universally into
all Standard Model fields but not into massless gravitons. Matching the
measured DM abundance together with the requirement of stability constrains the
spin2 mass to be in the 1 to 100 TeV range.

We review black hole and star solutions for Horndeski theory. For nonshift
symmetric theories, black holes involve a KaluzaKlein reduction of higher
dimensional Lovelock solutions. On the other hand, for shift symmetric theories
of Horndeski and beyond Horndeski, black holes involve two classes of
solutions: those that include, at the level of the action, a linear coupling to
the GaussBonnet term and those that involve time dependence in the galileon
field. We analyze the latter class in detail for a specific subclass of
Horndeski theory, discussing the general solution of a static and spherically
symmetric spacetime. We then discuss stability issues, slowly rotating
solutions as well as black holes coupled to matter. The latter case involves a
conformally coupled scalar field as well as an electromagnetic field and the
(primary) hair black holes thus obtained. We review and discuss the recent
results on neutron stars in Horndeski theories.

We find and study the properties of black hole solutions for a subclass of
Horndeski theory including the cubic Galileon term. The theory under study has
shift symmetry but not reflection symmetry for the scalar field. The Galileon
is assumed to have linear time dependence characterized by a velocity
parameter. We give analytic 3dimensional solutions that are akin to the BTZ
solutions but with a nontrivial scalar field that modifies the effective
cosmological constant. We then study the 4dimensional asymptotically flat and
de Sitter solutions. The latter present three different branches according to
their effective cosmological constant. For two of these branches, we find
families of black hole solutions, parametrized by the velocity of the scalar
field. These spherically symmetric solutions, obtained numerically, are
different from GR solutions close to the black hole event horizon, while they
have the same deSitter asymptotic behavior. The velocity parameter represents
black hole primary hair.

Observational evidence for the existence of Dark Matter is limited to its
gravitational effects. The extensive program for dedicated searches has yielded
null results so far, challenging the most popular models. Here we propose that
this is the case because the very existence of cold Dark Matter is a
manifestation of gravity itself. The consistent bimetric theory of gravity, the
only known ghostfree extension of General Relativity involving a massless and
a massive spin2 field, automatically contains a perfect Dark Matter candidate.
We demonstrate that the massive spin2 particle can be heavy, stable on
cosmological scales, and that it interacts with matter only through a
gravitational type of coupling. Remarkably, these features persist in the same
region of parameter space where bimetric theory satisfies the current gravity
tests. We show that the observed Dark Matter abundance can be generated via
freezein and suggest possible particle physics and gravitational signatures of
our bimetric Dark Matter model.

We study propagation of waves and appearance of caustics in kessence and
galileon theories. First we show that previously known solutions for travelling
waves in kessence and galileon models correspond to very specific finetuned
initial conditions. On the contrary, as we demonstrate by the method of
characteristics, generic initial conditions leads to a wave in kessence which
ends up with formation of caustics. Finally, we find that any wave solution in
pure kessence is also a solution for a galileon theory with the same kessence
term. Thus in the Horndeski theory with a kessence term formation of caustics
is generic. We discuss physical consequences of the caustics formation and
possible ways to cure the problem.

We review the black hole solutions of the ghostfree massive gravity theory
and its bimetric extension and outline the main results on the stability of
these solutions against small perturbations. Massive (bi)gravity accommodates
exact black hole solutions, analogous to those of General Relativity. In
addition to these solutions, hairy black holes  solutions with no
correspondent in General Relativity  have been found numerically, whose
existence is a natural consequence of the absence of Birkhoff's theorem in
these theories. The existence of extra propagating degrees of freedom, makes
the stability properties of these black holes richer and more complex than
those of General Relativity. In particular, the biSchwarzschild black hole
exhibits an unstable spherically symmetric mode, while the biKerr geometry is
also generically unstable, both against the spherical mode and against
superradiant instabilities. If astrophysical black holes are described by these
solutions, the superradiant instability of the Kerr solution imposes stringent
bounds on the graviton mass.

We consider an Abelian gauge field coupled to a particular truncation of
Horndeski theory. The Galileon field has translation symmetry and couples non
minimally both to the metric and the gauge field. When the gaugescalar
coupling is zero the gauge field reduces to a standard Maxwell field. By taking
into account the symmetries of the action, we construct charged black hole
solutions. Allowing the scalar field to softly break symmetries of spacetime we
construct black holes where the scalar field is regular on the black hole event
horizon. Some of these solutions can be interpreted as the equivalent of
ReissnerNordstrom black holes of scalar tensor theories with a non trivial
scalar field. A self tuning black hole solution found previously is extended to
the presence of dyonic charge without affecting whatsoever the self tuning of a
large positive cosmological constant. Finally, for a general shift invariant
scalar tensor theory we demonstrate that the scalar field Ansatz and method we
employ are mathematically compatible with the field equations. This opens up
the possibility for novel searches of hairy black holes in a far more general
setting of Horndeski theory.

We present a solution for rotating black holes in massive gravity. We first
give a solution of massive gravity with one dynamical metric. Both metrics of
this solution are expressed in the advanced EddingtonFinkelsteinlike
coordinates: the physical metric has the original Kerr line element, while the
fiducial metric is flat, but written in a rotating EddingtonFinkelstein form.
For the bigravity theory we give an analogue of this solution: the two metrics
have the original Kerr form, but, in general, different black hole masses. The
generalisation of the solution to include the electric charge is also given, it
is an analogue of the KerrNewman solution in General Relativity. We also
discuss further possible ways to generalise the solutions.

We present a class of exact scalartensor black holes for a shiftsymmetric
part of the Horndeski action. The action includes a higher order scalar tensor
interaction term. We find that for a static and spherically symmetric
spacetime, the scalar field, if time dependent, can be nontrivial and regular
thus circumventing in an interesting way nohair arguments for gallileons.
Furthermore, within this class we find a stealth Schwarzschild and a partially
selftuned deSitter Schwarzschild black hole, both exhibiting a non trivial
and regular space and time dependent scalar. In the latter solution the bulk
vacuum energy is screened from a necessarily smaller geometric effective de
Sitter vacuum via an integration constant associated to the time dependent
scalar field.

We present a new class of solutions describing charged black holes in massive
(bi)gravity. For a generic choice of the parameters of the massive gravity
action, the solution is the ReissnerNordstromde Sitter metric written in the
EddingtonFinkelstein coordinates for both metrics. We also study a special
case of the parameters, for which the space of solutions contains an extra
symmetry.

We consider the analytic solutions of massive (bi)gravity which can be
written in a simple form using advanced EddingtonFinkelstein coordinates. We
analyse the stability of these solutions against radial perturbations. First we
recover the previously obtained result on the instability of the bidiagonal
biSchwarzschild solutions. In the nonbidiagonal case (which contains, in
particular, the Schwarzschild solution with Minkowski fiducial metric) we show
that generically there are physical spherically symmetric perturbations, but no
unstable modes.

We study static spherically symmetric solutions of massive bigravity theory,
free from the BoulwareDeser ghost. We show the recovery of General Relativity
via the Vainshtein mechanism, in the weak limit of the physical metric. We find
a single polynomial equation determining the behavior of the solution for
distances smaller than the inverse graviton mass. This equation is generically
of the seventh order, while for a specific choice of the parameters of the
theory it can be reduced to lower orders. The solution is analytic in different
regimes: for distances below the Vainshtein radius (where General Relativity is
recovered), and in the opposite regime, beyond the Vainshtein radius, where the
solution approaches the flat metric.

We introduce the Vainshtein mechanism which plays a crucial role in massive
gravities, as well as in related theories such as Galileons and their
extensions. This mechanism, also known as kmouflage, allows to hide via non
linear effects  typically for source distances smaller than a socalled
Vainshtein radius which depends on the source and on the theory considered 
some degrees of freedom whose effects are then only left important at large
distances, e.g. for cosmology. It is introduced here in non linear FierzPauli
theories (massive gravities), including the dRGT theories, in their decoupling
limits, as well as in other models such as DGP model or generalized Galileons.
This presentation is selfcontained and before discussing the Vainshtein
mechanism we introduce some useful results and concepts concerning massive
gravity, such as the vDVZ discontinuity, the decoupling limits or the
BoulwareDeser ghost.

We show that linear perturbations around the simplest black hole solution of
massive bigravity theories, the biSchwarzschild solution, exhibit an unstable
mode featuring the GregoryLaflamme instability of higherdimensional black
strings. The result is obtained for the massive gravity theory which is free
from the BoulwareDeser ghost, as well as for its extension with two dynamical
metrics. These results may indicate that static black holes in massive gravity
do not exist. For graviton mass of order of Hubble scale, however, the
instability timescale is of order of the Hubble time.

We study spherically symmetric solutions of the cubic covariant Galileon
model in curved spacetime in presence of a matter source, in the test scalar
field approximation. We show that a cosmological time evolution of the Galileon
field gives rise to an induced matterscalar coupling, due to the
Galileongraviton kinetic braiding, therefore the solution for the Galileon
field is non trivial even if the bare matterscalar coupling constant is set to
zero. The local solution crucially depends on the asymptotic boundary
conditions, and in particular, Minkowski and de Sitter asymptotics correspond
to different branches of the solution. We study the stability of these
solutions, namely, the wellposedness of the Cauchy problem and the positivity
of energy for scalar and tensor perturbations, by diagonalizing the kinetic
terms of the spin2 and spin0 degrees of freedom. In addition, we find that in
presence of a cosmological time evolution of the Galileon field, its kinetic
mixing with the graviton leads to a friction force, resulting to efficient
damping of scalar perturbations within matter.