The aim of this paper is to constrain modified gravity with redshift space distortion observations and supernovae measurements. Compared with a standard LCDM analysis, we include three additional free parameters, namely the initial conditions of the matter perturbations,the overall perturbation normalization, and a scale-dependent Y parameter modifying the Poisson equation, in an attempt to perform a more model-independent analysis. First, we constrain the Poisson parameter Y (also called Geff ) by using currently available fsigma_8 data and the recent SN catalog JLA. We find that the inclusion of the additional free parameters makes the constraints significantly weaker than when fixing them to the standard cosmological value. Second, we constrain Y by using forecast growth-rate data for Euclid and SKA missions. Here again we point out the weakening of the constraints when the additional parameters are included. Finally, we adopt as modified gravity Poisson parameter the specific Horndeski form, and use scale-dependent forecasts to build an exclusion plot for the Yukawa potential akin to the ones realized in laboratory experiments, both for the Euclid and the SKA surveys.

Most cosmological constraints on modified gravity are obtained assuming that the cosmic evolution was standard $\Lambda$CDM in the past and that the present matter density and power spectrum normalization are the same as in a $\Lambda$CDM model. Here we examine how the constraints change when these assumptions are lifted. We focus in particular on the parameter $Y$ (also called $G_{\mathrm{eff}}$) that quantifies the deviation from the Poisson equation. This parameter can be estimated by comparing with the model-independent growth rate quantity $f\sigma_{8}(z)$ obtained through redshift distortions. We reduce the model dependency in evaluating $Y$ by marginalizing over $\sigma_{8}$ and over the initial conditions, and by absorbing the degenerate parameter $\Omega_{m,0}$ into $Y$. We use all currently available values of $f\sigma_{8}(z)$. We find that the combination $\hat{Y}=Y\Omega_{m,0}$, assumed constant in the observed redshift range, can be constrained only very weakly by current data, $\hat{Y}=0.28_{-0.23}^{+0.35}$ at 68\% c.l. We also forecast the precision of a future estimation of $\hat{Y}$ in a Euclid-like redshift survey. We find that the future constraints will reduce substantially the uncertainty, $\hat{Y}=0.30_{-0.09}^{+0.08}$ , at 68\% c.l., but the relative error on $\hat{Y}$ around the fiducial remains quite high, of the order of 30\%. The main reason for these weak constraints is that $\hat{Y}$ is strongly degenerate with the initial conditions, so that large or small values of $\hat{Y}$ are compensated by choosing non-standard initial values of the derivative of the matter density contrast. Finally, we produce a forecast of a cosmological exclusion plot on the Yukawa strength and range parameters, which complements similar plots on laboratory scales but explores scales and epochs reachable only with large-scale galaxy surveys. (abridged)

We study the gravitational clustering of spherically symmetric overdensities and the statistics of the resulting dark matter halos in the "symmetron model", in which a new long range force is mediated by a $Z_2$ symmetric scalar field. Depending on the initial radius of the overdensity, we identify two distinct regimes: for small initial radii the symmetron mediated force affects the spherical collapse at all redshifts; for initial radii larger than some critical size this force vanishes before collapse because of the symmetron screening mechanism. In both cases overdensities collapse earlier than in the $\Lambda$CDM and statistically tend to form more massive dark matter halos. Regarding the halo-mass function of these objects, we observe order one departures from standard $\Lambda$CDM predictions. The formalism developed here can be easily applied to other models where fifth-forces participate to the dynamics of the gravitational collapse.