Simulation tests of galaxy cluster constraints on chameleon gravity
  We use two new hydrodynamical simulations of $\Lambda$CDM and $f(R)$ gravity to test the methodology used by Wilcox et al. 2015 (W15) in constraining the effects of a fifth force on the profiles of clusters of galaxies. We construct realistic simulated stacked weak lensing and X-ray surface brightness cluster profiles from these cosmological simulations, and then use these data projected along various lines-of-sight to test the spherical symmetry of our stacking procedure. We also test the applicability of the NFW profile to model weak lensing profiles of clusters in $f(R)$ gravity. Finally, we test the validity of the analytical model developed in W15 against the simulated profiles. Overall, we find our methodology is robust and broadly agrees with these simulated data. We also apply our full Markov Chain Monte Carlo (MCMC) analysis from W15 to our simulated X-ray and lensing profiles, providing consistent constraints on the modified gravity parameters as obtained from the real cluster data, e.g. for our $\Lambda$CDM simulation we obtain $|f_{\rm{R}0}| < 8.3 \times 10^{-5}$ (95% CL), which is in good agreement with the W15 measurement of $|f_{\rm{R}0}| < 6 \times 10^{-5}$. Overall, these tests confirm the power of our methodology which can now be applied to larger cluster samples available with the next generation surveys. 
Harry Wilcox, Robert C. Nichol, Gong-bo Zhao, David Bacon, Kazuya   Koyama, A. Kathy Romer
astro-ph.CO