We use two new hydrodynamical simulations of Lambda cold dark matter (LambdaCDM) and f(R) gravity to test the methodology used by Wilcox et al. (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 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 LambdaCDM simulation we obtain |fR0| < 8.3 × 10-5 (95 per cent CL), which is in good agreement with the W15 measurement of |fR0| < 6 × 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.
Wilcox et al. (Wed,) studied this question.