N-body Simulations for f(R) Gravity using a Self-adaptive Particle-Mesh Code

TL;DR

This paper presents high-resolution N-body simulations for f(R) gravity using a self-adaptive particle-mesh code, confirming previous results and extending the resolution to better understand the Chameleon mechanism's impact on dark matter clustering.

Contribution

The authors develop a self-adaptive particle-mesh code for f(R) gravity simulations, fully incorporating the Chameleon mechanism and extending the resolution of matter power spectrum measurements.

Findings

Confirmed previous simulation results on matter power spectrum and halo properties.

Extended the resolution of the matter power spectrum up to k~20 h/Mpc.

Analyzed the effects of the Chameleon mechanism on non-linear dark matter clustering.

Abstract

We perform high resolution N-body simulations for f(R) gravity based on a self-adaptive particle- mesh code MLAPM. The Chameleon mechanism that recovers General Relativity on small scales is fully taken into account by self-consistently solving the non-linear equation for the scalar field. We independently confirm the previous simulation results, including the matter power spectrum, halo mass function and density profiles, obtained by Oyaizu et al. (Phys.Rev.D 78, 123524, 2008) and Schmidt et al. (Phys.Rev.D 79, 083518, 2009), and extend the resolution up to k~20 h/Mpc for the measurement of the matter power spectrum. Based on our simulation results, we discuss how the Chameleon mechanism affects the clustering of dark matter and halos on full non-linear scales.