The discovery of cosmic acceleration has triggered a consistent body of theoretical work aimed at modeling its phenomenology and understanding its fundamental physical nature. In recent years, a powerful formalism that accomplishes both these goals has been developed, the so-called effective field theory of dark energy. It can capture the behavior of a wide class of modified gravity theories and classify them according to the imprints they leave on the smooth background expansion history of the Universe and on the evolution of linear perturbations. The effective field theory of dark energy is based on a Lagrangian description of cosmological perturbations which depends on a number of functions of time, some of which are non-minimal couplings representing genuine deviations from General Relativity. Such a formalism is thus particularly convenient to fit and interpret the wealth of new data that will be provided by future galaxy surveys. Despite its recent appearance, this formalism has already allowed a systematic investigation of what lies beyond the General Relativity landscape and provided a conspicuous amount of theoretical predictions and observational results. In this review, we report on these achievements.

The key probes of the growth of large-scale structure are its rate $f$ and amplitude $\sigma_8$. Redshift space distortions in the galaxy power spectrum allow us to measure only the combination $f\sigma_8$, which can be used to constrain the standard cosmological model or alternatives. By using measurements of the galaxy-galaxy lensing cross-correlation spectrum or of the galaxy bispectrum, it is possible to break the $f\sigma_8$ degeneracy and obtain separate estimates of $f$ and $\sigma_8$ from the same galaxy sample. Currently there are only a handful of such separate measurements, but even this allows for improved constraints on cosmological models. We consider what can be achieved by a future nominal sample that delivers a $\sim 1\%$ constraint on $f$ and $\sigma_8$ separately, compared to the case with a similar precision on the combination $f\sigma_8$. For the six cosmological parameters of $\Lambda$CDM, we find improvements of $\sim\! 5$--$50\%$ on their constraints. For modified gravity models in the Horndeski class, the improvements on these standard parameters are $\sim\! 0$--$15\%$. However, the precision on the sum of neutrino masses improves by 65\% and there is a significant increase in the precision on the background and perturbation Horndeski parameters.

We use growth of structure data to constrain the effective field theory of dark energy. Considering as case study Horndeski theories with the speed of gravitational waves equal to that of light, we show how constraints on the free parameters and the large-scale structure phenomenological functions can be improved by two ingredients: firstly by complementing the set of redshift-space distortions data with the three recent measurements of the growth rate $f$ and the amplitude of matter fluctuations $\sigma_8$ from the VIPERS and SDSS collaborations; secondly by applying a local Solar System bound on the variation of the Newton constant. This analysis allows us to conclude that: $i)$ despite firmly restricting the predictions of weaker gravity, the inclusion of the Solar System bound does not prevent suppressed growth relative to the standard model $\Lambda$CDM at low redshifts; $ii)$ the same bound in conjunction with the growth of structure data strongly restricts the redshift evolution of the gravitational slip parameter to be close to unity and the present value is constrained to one at the $10^{-3}$ level; $iii)$ the growth of structure data favours a fifth force contribution to the effective gravitational coupling at low redshifts and at more than one sigma at present time.

We summarise the effective field theory of dark energy construction to explore observable predictions of linear Horndeski theories. Based on \cite{Perenon:2016blf}, we review the diagnostic of these theories on the correlation of the large-scale structure phenomenological functions: the effective Newton constant, the light deflection parameter and the growth function of matter perturbations. We take this opportunity to discuss the evolution of the bounds the propagation speed of gravitational waves has undergone and use the most restrictive one to update the diagnostic.

We study the effects of Horndeski models of dark energy on the observables of the large-scale structure in the late time universe. A novel classification into {\it Late dark energy}, {\it Early dark energy} and {\it Early modified gravity} scenarios is proposed, according to whether such models predict deviations from the standard paradigm persistent at early time in the matter domination epoch. We discuss the physical imprints left by each specific class of models on the effective Newton constant $\mu$, the gravitational slip parameter $\eta$, the light deflection parameter $\Sigma$ and the growth function $f\sigma_8$ and demonstrate that a convenient way to dress a complete portrait of the viability of the Horndeski accelerating mechanism is via two, redshift-dependent, diagnostics: the $\mu(z)-\Sigma(z)$ and the $f\sigma_8(z)-\Sigma(z)$ planes. If future, model-independent, measurements point to either $\Sigma-1<0$ at redshift zero or $\mu-1<0$ with $\Sigma-1>0$ at high redshifts or $\mu-1>0$ with $\Sigma-1<0$ at high redshifts, Horndeski theories are effectively ruled out. If $f\sigma_8$ is measured to be larger than expected in a $\Lambda$CDM model at $z>1.5$ then Early dark energy models are definitely ruled out. On the opposite case, Late dark energy models are rejected by data if $\Sigma<1$, while, if $\Sigma>1$, only Early modifications of gravity provide a viable framework to interpret data.

We present a systematic study of modified gravity (MG) models containing a single scalar field non-minimally coupled to the metric. Despite a large parameter space, exploiting the effective field theory of dark energy (EFT of DE) formulation and imposing simple physical constraints such as stability conditions and (sub-)luminal propagation of perturbations, we arrive at a number of generic predictions about the large scale structures.

We present a systematic exploration of dark energy and modified gravity models containing a single scalar field non-minimally coupled to the metric. Even though the parameter space is large, by exploiting an effective field theory (EFT) formulation and by imposing simple physical constraints such as stability conditions and (sub-)luminal propagation of perturbations, we arrive at a number of generic predictions. (1) The linear growth rate of matter density fluctuations is generally suppressed compared to $\Lambda$CDM at intermediate redshifts ($0.5 \lesssim z \lesssim 1$), despite the introduction of an attractive long-range scalar force. This is due to the fact that, in self-accelerating models, the background gravitational coupling weakens at intermediate redshifts, over-compensating the effect of the attractive scalar force. (2) At higher redshifts, the opposite happens; we identify a period of super-growth when the linear growth rate is larger than that predicted by $\Lambda$CDM. (3) The gravitational slip parameter $\eta$ - the ratio of the space part of the metric perturbation to the time part - is bounded from above. For Brans-Dicke-type theories $\eta$ is at most unity. For more general theories, $\eta$ can exceed unity at intermediate redshifts, but not more than about $1.5$ if, at the same time, the linear growth rate is to be compatible with current observational constraints. We caution against phenomenological parametrization of data that do not correspond to predictions from viable physical theories. We advocate the EFT approach as a way to constrain new physics from future large-scale-structure data.