PySCo-EFT and ECOSMOG-EFT: a tandem of N-body simulation codes for the Effective Field Theory of Dark Energy

TL;DR

This paper introduces two new N-body simulation codes, PySCo-EFT and ECOSMOG-EFT, designed to accurately model the non-linear effects of the Effective Field Theory of Dark Energy on cosmic structure formation.

Contribution

The paper presents the development and validation of two novel simulation codes for EFTofDE cosmologies, incorporating non-linear screening mechanisms and achieving high accuracy.

Findings

Codes agree within 0.5% with linear theory on large scales.

Agreement between the two codes is within 2% on non-linear scales.

Numerical effects are limited to below 2% at high wavenumbers.

Abstract

Modified gravity theories constitute viable alternatives to the standard cosmological model for explaining the observed late-time accelerated expansion of the Universe. The Effective Field Theory of Dark Energy (EFTofDE) is an efficient framework to describe a wide range of such theories with a limited number of parameters. To robustly constrain them by comparison with clustering and weak lensing data from upcoming large-scale structure surveys, high-resolution cosmological N-body simulations are required to obtain accurate predictions for the matter distribution on non-linear scales. We introduce two new N-body simulation codes for EFTofDE cosmologies: PySCo-EFT, a Python-based particle mesh code, and ECOSMOG-EFT, a RAMSES-based code with adaptive mesh refinement. We consider Horndeski models with a luminal gravitational wave speed. We use iterative solvers and multigrid schemes to solve for the additional scalar field equation in both codes, incorporating the non-linear Vainshtein screening mechanism. We present validation and convergence tests of the codes. We obtain a sub-0.5 percent agreement with linear theory on large scales and a similar agreement between the two codes on non-linear scales. The dominant numerical effects on the matter-power-spectrum boost are mass resolution, finite-volume effects, refinement threshold, and starting redshift, but they are limited to below 2% at the largest wavenumbers (k=10 h/Mpc) for the range of tested values. We investigate the impact of the EFTofDE parameters on the matter-power-spectrum ratios between EFTofDE and $\Lambda$CDM cases. Depending on the EFTofDE parameters, the screening plays a negligible or dominant role compared to the linearised field equations. Our codes provide tools for generating fast and accurate predictions of the impact of the EFTofDE on the clustering of matter, incorporating non-linear screening.