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

Understand

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.