paper

Efficient simulations of large scale structure in modified gravity cosmologies with comoving Lagrangian acceleration

arXiv:1612.06469 · doi:10.1103/PhysRevD.95.103515

Abstract

We implement an adaptation of the COLA approach, a hybrid scheme that combines Lagrangian perturbation theory with an N-body approach, to model non-linear collapse in chameleon and symmetron modified gravity models. Gravitational screening is modeled effectively through the attachment of a suppression factor to the linearized Klein-Gordon equations. The adapted COLA approach is benchmarked, with respect to an N-body code both for the CDM scenario and for the modified gravity theories. It is found to perform well in the estimation of the dark matter power spectra, with consistency of 1 % to h/Mpc. Redshift space distortions are shown to be effectively modeled through a Lorentzian parameterization with a velocity dispersion fit to the data. We find that COLA performs less well in predicting the halo mass functions, but has consistency, within uncertainties of our simulations, in the relative changes to the mass function induced by the modified gravity models relative to CDM. The results demonstrate that COLA, proposed to enable accurate and efficient, non-linear predictions for CDM, can be effectively applied to a wider set of cosmological scenarios, with intriguing properties, for which clustering behavior needs to be understood for upcoming surveys such as LSST, DESI, Euclid and WFIRST.

15 pages, 9 figures. Updated to match version published in PRD

References in corpus (13)

Cited by in corpus (48)