Joint lensing and dynamical mass profile determinations of galaxy clusters are an excellent tool to constrain modification of gravity at cosmological scales. However, search for tiny departures from General Relativity calls for an accurate control of the systematics affecting the method. In this analysis we concentrate on the systematics in the reconstruction of mass profiles from the dynamics of cluster member galaxies, while assuming that lensing provides unbiased mass profile reconstructions. In particular, in the case study of linear $f(R)$ gravity, we aim at veryfying whether in realistic simulations of cluster formation a spurious detection of departure from GR can be detected due to violation of the main assumptions (e.g. dynamical equilibrium and spherical symmetry) on which the method is based. We aim at identifying and calibrating the impact of those systematics by analyzing a set of Dark Matter halos taken from $\Lambda$CDM N-body cosmological simulations performed with the GADGET-3 code. [...] If no selection criteria are applied, $\sim 60\%$ of clusters in a $\Lambda$CDM Universe (where GR is assumed) produce a spurious detection of modified gravity. We find that the probability of finding cluster in agreement with GR predictions $P_{GR}$ mainly depends on the properties of the halo's projected phase-space and on shape orientation of the cluster along the line-of-sight projection. We define two observational criteria which correlate with the probability to find clusters in agreement with GR predictions and which can be used to select [...] those objects that are more suitable for the application of the proposed method. In particular, we find that according to these criteria the percentage of spurious detection can be lowered down to $\sim 20\%$ in the best case. Our results are relevant in view of data that will be available with the next generation surveys.

The gravitational slip parameter is an important discriminator between large classes of gravity theories at cosmological and astrophysical scales. In this work we use a combination of simulated information of galaxy cluster mass profiles, inferred by Strong+Weak lensing analyses and by the study of the dynamics of the cluster member galaxies, to reconstruct the gravitational slip parameter $\eta$ and predict the accuracy with which it can be constrained with current and future galaxy cluster surveys. Performing a full-likelihood statistical analysis, we show that galaxy cluster observations can constrain $\eta$ down to the percent level already with a few tens of clusters. We discuss the significance of possible systematics, and show that the cluster masses and numbers of galaxy members used to reconstruct the dynamics mass profile have a mild effect on the predicted constraints.

We perform a maximum likelihood kinematic analysis of the two dynamically relaxed galaxy clusters MACS J1206.2-0847 at $z=0.44$ and RXC J2248.7-4431 at $z=0.35$ to determine the total mass profile in modified gravity models, using a modified version of the MAMPOSSt code of Mamon, Biviano and Bou\'e. Our work is based on the kinematic and lensing mass profiles derived using the data from the Cluster Lensing And Supernova survey with Hubble (hereafter CLASH) and the spectroscopic follow-up with the Very Large Telescope (hereafter CLASH-VLT). We assume a spherical Navarro-Frenk-White (NFW hereafter) profile in order to obtain a constraint on the fifth force interaction range $\lambda$ for models in which the dependence of this parameter on the enviroment is negligible at the scale considered (i.e. $\lambda=const$) and fixing the fifth force strength to the value predicted in $f(R)$ gravity. We then use information from lensing analysis to put a prior on the other NFW free parameters. In the case of MACSJ 1206 the joint kinematic+lensing analysis leads to an upper limit on the effective interaction range $\lambda\le1.61\,\mbox{Mpc}$ at $\Delta\chi^{2}=2.71$ on the marginalized distribution. For RXJ 2248 instead a possible tension with the $\Lambda$CDM model appears when adding lensing information, with a lower limit $\lambda\ge0.14\,\mbox{Mpc}$ at $\Delta\chi^{2}=2.71$. This is consequence of the slight difference between the lensing and kinematic data, appearing in GR for this cluster, that could be explained in terms of modifications of gravity. We discuss the impact of systematics and the limits of our analysis as well as future improvements of the results obtained.[...]

We use high-precision kinematic and lensing measurements of the total mass profile of the dynamically relaxed galaxy cluster MACS J1206.2-0847 at $z=0.44$ to estimate the value of the ratio $\eta=\Psi/\Phi$ between the two scalar potentials in the linear perturbed Friedmann-Lemaitre-Robertson-Walker metric.[...] Complementary kinematic and lensing mass profiles were derived from exhaustive analyses using the data from the Cluster Lensing And Supernova survey with Hubble (CLASH) and the spectroscopic follow-up with the Very Large Telescope (CLASH-VLT). Whereas the kinematic mass profile tracks only the time-time part of the perturbed metric (i.e. only $\Phi$), the lensing mass profile reflects the contribution of both time-time and space-space components (i.e. the sum $\Phi+\Psi$). We thus express $\eta$ as a function of the mass profiles and perform our analysis over the radial range $0.5\,Mpc\le r\le r_{200}=1.96\,Mpc$. Using a spherical Navarro-Frenk-White mass profile, which well fits the data, we obtain $\eta(r_{200})=1.01\,_{-0.28}^{+0.31}$ at the 68\% C.L. We discuss the effect of assuming different functional forms for mass profiles and of the orbit anisotropy in the kinematic reconstruction. Interpreting this result within the well-studied $f(R)$ modified gravity model, the constraint on $\eta$ translates into an upper bound to the interaction length (inverse of the scalaron mass) smaller than 2 Mpc. This tight constraint on the $f(R)$ interaction range is however substantially relaxed when systematic uncertainties in the analysis are considered. Our analysis highlights the potential of this method to detect deviations from general relativity, while calling for the need of further high-quality data on the total mass distribution of clusters and improved control on systematic effects.