• A Hamiltonian density bounded from below implies that the lowest-energy state is stable. We point out, contrary to common lore, that an unbounded Hamiltonian density does not necessarily imply an instability: this is a coordinate-dependent statement. We discuss this fact in detail, starting from k-essence 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 Schwarzschild-de Sitter solution of a beyond-Horndeski theory, while taking into account the recent experimental constraint regarding the speed of gravitational waves. We extract the spin-2 and spin-0 causal cones by analyzing respectively all the odd-parity and the $\ell = 0$ even-parity 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 k-essence 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/self-interacting complex scalar reproduces dynamics of pressureless perfect fluid/shift-symmetric k-essence under certain initial conditions in the limit of large mass/sharp self-interacting 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 ultra-light axions (generically, coherent oscillations of bosons in the Bose--Einstein condensate state).
  • We consider a generic scalar-tensor theory involving a shift-symmetric scalar field and minimally coupled matter fields. We prove that the Noether current associated with shift-symmetry 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 no-hair 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 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.
  • We find spherically symmetric and static black holes in shift-symmetric 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 Gauss-Bonnet 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 non-trivial scalar field is an exact solution to these theories.
  • Motivated by the projectable Horava--Lifshitz 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 non-relativistic 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 self-tuning and the local predictions (inside the Solar system) of the most general shift-symmetric beyond Horndeski theory. We first show that the cosmological self-tuning 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 Schwarzschild-de Sitter solution around a static and spherically symmetric source. As a by-product of our study, we identify a general subclass of beyond Horndeski theory for which regular self-tuning black hole solutions exist, in presence of a time-dependent scalar field. We discuss possible future development of the present work.
  • 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.
  • We consider generic linear perturbations of a nonbidiagonal class of static black-hole 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 black-hole 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 spin-2 particle. Massive spin-2 fields and their gravitational interactions are uniquely described by ghost-free 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 spin-2 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 spin-2 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 spin-2 mass to be in the 1 to 100 TeV range.
  • We review black hole and star solutions for Horndeski theory. For non-shift symmetric theories, black holes involve a Kaluza-Klein 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 Gauss-Bonnet 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 3-dimensional solutions that are akin to the BTZ solutions but with a non-trivial scalar field that modifies the effective cosmological constant. We then study the 4-dimensional 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 de-Sitter 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 ghost-free extension of General Relativity involving a massless and a massive spin-2 field, automatically contains a perfect Dark Matter candidate. We demonstrate that the massive spin-2 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 freeze-in and suggest possible particle physics and gravitational signatures of our bimetric Dark Matter model.
  • We study propagation of waves and appearance of caustics in k-essence and galileon theories. First we show that previously known solutions for travelling waves in k-essence and galileon models correspond to very specific fine-tuned initial conditions. On the contrary, as we demonstrate by the method of characteristics, generic initial conditions leads to a wave in k-essence which ends up with formation of caustics. Finally, we find that any wave solution in pure k-essence is also a solution for a galileon theory with the same k-essence term. Thus in the Horndeski theory with a k-essence 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 ghost-free 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 bi-Schwarzschild black hole exhibits an unstable spherically symmetric mode, while the bi-Kerr 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 gauge-scalar 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 Reissner-Nordstrom 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 Eddington-Finkelstein-like coordinates: the physical metric has the original Kerr line element, while the fiducial metric is flat, but written in a rotating Eddington-Finkelstein form. For the bi-gravity 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 Kerr-Newman solution in General Relativity. We also discuss further possible ways to generalise the solutions.
  • We present a class of exact scalar-tensor black holes for a shift-symmetric part of the Horndeski action. The action includes a higher order scalar tensor interaction term. We find that for a static and spherically symmetric space-time, the scalar field, if time dependent, can be non-trivial and regular thus circumventing in an interesting way no-hair arguments for gallileons. Furthermore, within this class we find a stealth Schwarzschild and a partially self-tuned de-Sitter 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 Reissner-Nordstrom-de Sitter metric written in the Eddington-Finkelstein 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 Eddington-Finkelstein coordinates. We analyse the stability of these solutions against radial perturbations. First we recover the previously obtained result on the instability of the bidiagonal bi-Schwarzschild solutions. In the non-bidiagonal 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 bi-gravity theory, free from the Boulware-Deser 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 k-mouflage, allows to hide via non linear effects - typically for source distances smaller than a so-called 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 Fierz-Pauli 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 self-contained 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 Boulware-Deser ghost.
  • We show that linear perturbations around the simplest black hole solution of massive bi-gravity theories, the bi-Schwarzschild solution, exhibit an unstable mode featuring the Gregory-Laflamme instability of higher-dimensional black strings. The result is obtained for the massive gravity theory which is free from the Boulware-Deser 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 matter-scalar coupling, due to the Galileon-graviton kinetic braiding, therefore the solution for the Galileon field is non trivial even if the bare matter-scalar 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 well-posedness of the Cauchy problem and the positivity of energy for scalar and tensor perturbations, by diagonalizing the kinetic terms of the spin-2 and spin-0 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.