Horndeski gravity was highly constrained from the recent gravitational wave observations by the LIGO Collaboration down to $|c_{g}/c-1|\gtrsim 10^{-15}$. In this Letter we study the tensorial perturbations in a flat cosmological background for an analogue version of Horndenki gravity which is based in Teleparallel Gravity constructed from a flat manifold with a nonvanishing torsion tensor. It is found that in this approach, one can construct a more general Horndeski theory satisfying $c_T=c_g/c=1$ without eliminating the coupling functions $G_5(\phi,X)$ and $G_4(\phi,X)$ that were highly constrained in standard Horndeski theory. Hence, in the Teleparallel approach one is able to restore these terms, creating an interesting way to revive Horndeski gravity.

Non-local theories of gravity have recently gained a lot of interest because they can suitably represent the behavior of gravitational interaction in the ultraviolet regime. Furthermore, at infrared scales, they give rise to notable cosmological effects which could be important to describe the dark energy behavior. In particular, exponential forms of the distortion function seem particularly useful for this purpose. Using Noether Symmetries, it can be shown that the only non-trivial form of the distortion function is the exponential one, which is working not only for cosmological mini-superspaces, but also in a spherically symmetric spacetime. Taking this result into account, we study the weak field approximation of this type of non-local gravity, and comparing with the orbits of S2 stars around the Galactic center (NTT/VLT data), we set constraints on the parameters of the theory. Non-local effects do not play a significant role on the orbits of S2 stars around Sgr A*, but give richer phenomenology at cosmological scales than the $\Lambda$CDM model. Also, we show that non-local gravity model gives better agreement between theory and astronomical observations than Keplerian orbits.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

Adopting Noether point symmetries, we classify and integrate dynamical systems coming from Horndeski cosmologies. The method is particularly effective both to select the form of Horndeski models and to derive exact cosmological solutions. Starting from the Lagrangians selected by the Noether symmetries, it is possible to derive several modified theories of gravity like $f(R)$ gravity, Brans-Dicke gravity, string inspired gravity and so on. In any case, exact solutions are found out.

The accelerating behavior of cosmic fluid opposes to the gravitational attraction, at present epoch, whereas standard gravity is dominant at small scales. As a consequence, there exists a \emph{point} where the effects are counterbalanced, dubbed \emph{turnaround radius}, $r_{\text{ta}}$. By construction, it provides a bound on maximum structure sizes of the observed universe. Once an upper bound on $r_{\text{ta}}$ is provided, i.e. $R_{\text{TA,max}}$, one can check whether cosmological models guarantee structure formation. Here, we focus on $f(R)$ gravity, without imposing \emph{a priori} the form of $f(R)$. We thus provide an analytic expression for the turnaround radius in the framework of $f(R)$ models. To figure this out, we compute the turnaround radius in two distinct cases: 1) under the hypothesis of static and spherically symmetric space-time, and 2) by using the cosmological perturbation theory. We thus find a criterion to enable large scale structures to be stable in $f(R)$ models, circumscribing the class of $f(R)$ theories as suitable alternative to dark energy. In particular, we get that for constant curvature, the viability condition becomes $R_{\text{dS}}f'(R_{\text{dS}}) \leq 5.48 \Lambda \Rightarrow f'(R_{\text{dS}}) \leq 1.37$, with $\Lambda$ and $R_{\text{dS}}$ respectively the observed cosmological constant and the Ricci curvature. This prescription rules out models which do not pass the aforementioned $R_{\text{TA,max}}$ limit.

We study a generalized nonlocal theory of gravity which, in specific limits, can become either the curvature non-local or teleparallel non-local theory. Using the Noether Symmetry Approach, we find that the coupling functions coming from the non-local terms are constrained to be either exponential or linear in form. It is well known that in some non-local theories, a certain kind of exponential non-local couplings are needed in order to achieve a renormalizable theory. In this paper, we explicitly show that this kind of coupling does not need to by introduced by hand, instead, it appears naturally from the symmetries of the Lagrangian in flat Friedmann-Robertson-Walker cosmology. Finally, we find de-Sitter and power law cosmological solutions for different nonlocal theories. The symmetries for the generalized non-local theory is also found and some cosmological solutions are also achieved under the full theory.

The maximum size of a cosmic structure is given by the maximum turnaround radius -- the scale where the attraction due to its mass is balanced by the repulsion due to dark energy. We derive generic formulae for the estimation of the maximum turnaround radius in any theory of gravity obeying the Einstein equivalence principle, in two situations: on a spherically symmetric spacetime and on a perturbed Friedman-Robertson-Walker spacetime. We show that the two formulae agree. As an application of our formula, we calculate the maximum turnaround radius in the case of the Brans-Dicke theory of gravity. We find that for this theory, such maximum sizes always lie above the \LCDM value, by a factor $1 + \frac{1}{3\omega}$, where $\omega\gg 1$ is the Brans-Dicke parameter, implying consistency of the theory with current data.

A generalized teleparallel cosmological model, $f(T_\mathcal{G},T)$, containing the torsion scalar $T$ and the teleparallel counterpart of the Gauss-Bonnet topological invariant $T_{\mathcal{G}}$, is studied in the framework of the Noether Symmetry Approach. As $f(\mathcal{G}, R)$ gravity, where $\mathcal{G}$ is the Gauss-Bonnet topological invariant and $R$ is the Ricci curvature scalar, exhausts all the curvature information that one can construct from the Riemann tensor, in the same way, $f(T_\mathcal{G},T)$ contains all the possible information directly related to the torsion tensor. In this paper, we discuss how the Noether Symmetry Approach allows to fix the form of the function $f(T_\mathcal{G},T)$ and to derive exact cosmological solutions.

The maximum size of a bound cosmic structure is computed perturbatively as a function of its mass in the framework of the cubic galileon, proposed recently to model the dark energy of our Universe. Comparison of our results with observations constrains the matter-galileon coupling of the model to $0.033\lesssim \alpha \lesssim 0.17$, thus improving previous bounds based solely on solar system physics.