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

Understand

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 $\Lambda$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 $k\sim2.5$ h/Mpc.