We use cookies to ensure that we give you the best experience on our website. By continuing to browse this repository, you give consent for essential cookies to be used. You can read more about our Privacy and Cookie Policy.

Durham Research Online
You are in:

Speeding up N-body simulations of modified gravity : chameleon screening models.

Bose, Sownak and Li, Baojiu and Barreira, Alexandre and He, Jian-hua and Hellwing, Wojciech A. and Koyama, Kazuya and Llinares, Claudio and Zhao, Gong-Bo (2017) 'Speeding up N-body simulations of modified gravity : chameleon screening models.', Journal of cosmology and astroparticle physics., 2017 (02). 050.


We describe and demonstrate the potential of a new and very efficient method for simulating certain classes of modified gravity theories, such as the widely studied f(R) gravity models. High resolution simulations for such models are currently very slow due to the highly nonlinear partial differential equation that needs to be solved exactly to predict the modified gravitational force. This nonlinearity is partly inherent, but is also exacerbated by the specific numerical algorithm used, which employs a variable redefinition to prevent numerical instabilities. The standard Newton-Gauss-Seidel iterative method used to tackle this problem has a poor convergence rate. Our new method not only avoids this, but also allows the discretised equation to be written in a form that is analytically solvable. We show that this new method greatly improves the performance and efficiency of f(R) simulations. For example, a test simulation with 5123 particles in a box of size 512 Mpc/h is now 5 times faster than before, while a Millennium-resolution simulation for f(R) gravity is estimated to be more than 20 times faster than with the old method. Our new implementation will be particularly useful for running very high resolution, large-sized simulations which, to date, are only possible for the standard model, and also makes it feasible to run large numbers of lower resolution simulations for covariance analyses. We hope that the method will bring us to a new era for precision cosmological tests of gravity.

Item Type:Article
Full text:(AM) Accepted Manuscript
Download PDF
Publisher Web site:
Publisher statement:This is an author-created, un-copyedited version of an article published in Journal of Cosmology and Astroparticle Physics. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at
Date accepted:17 February 2017
Date deposited:03 March 2017
Date of first online publication:28 February 2017
Date first made open access:28 February 2018

Save or Share this output

Look up in GoogleScholar