Abstract
We implement an adaptation of the cola approach, a hybrid scheme that combines Lagrangian perturbation theory with an N-body approach, to model nonlinear 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 cold dark matter () 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 . Redshift space distortions are shown to be effectively modeled through a Lorentzian parametrization 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 . The results demonstrate that cola, proposed to enable accurate and efficient, nonlinear predictions for , 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.
2 More- Received 27 December 2016
DOI:https://doi.org/10.1103/PhysRevD.95.103515
© 2017 American Physical Society