Modeling nonlinear scales for dynamical dark energy cosmologies with COLA
Phys. Rev. D 113, 103513 – Published 7 May, 2026
DOI: https://doi.org/10.1103/1zgv-sx5s
Abstract
Upcoming galaxy surveys will bring a wealth of information about the clustering of matter at small scales, but modeling small-scale structure beyond cold dark matter () remains computationally challenging. While accurate -body emulators exist to model the matter power spectrum for and some limited extensions, it is unfeasible to generate -body simulation suites for all candidate models. Motivated by recent hints of an evolving dark energy equation of state from galaxy surveys, we assess the viability of employing the comoving Lagrangian acceleration (COLA) method to generate simulation suites assuming the dynamical dark energy model. Following up on our previous work, we combine COLA simulations with an existing high-precision emulator to extend its predictions into new regions of parameter space. We assess the precision of our emulator at the level of the matter power spectrum, finding that our emulator can reproduce the nonlinear boosts from euclidemulator2 at less than 2% error. Moreover, we perform an analysis of a simulated cosmic shear survey akin to future data from the first year of observations from the Legacy Survey of Space and Time (LSST), assessing the differences in parameter constraints between our COLA-based emulator and the benchmark emulator. We find our emulator to be in excellent agreement with the benchmark, achieving less than shifts in cosmological parameter constraints across multiple fiducial cosmologies. We further compare our emulator’s performance to a commonly used approach: assuming the boost can be employed for extended parameter spaces without modification. We find that our emulator yields a significantly smaller distribution, parameter constraint biases, and a more accurate figure of merit compared to this second approach. These results demonstrate that COLA emulators provide a computationally efficient and physically motivated path forward for modeling nonlinear structure in extended cosmologies, offering a practical alternative to full -body suites in the era of precision cosmology.