Euclid preparation: II. The EUCLIDEMULATOR - a tool to compute the cosmology dependence of the nonlinear matter power spectrum

Knabenhans, M., Stadel, J., Marelli, S., Potter, D., Teyssier, R., Legrand, L., Schneider, A., Sudret, B., Blot, L., Awan, S., Burigana, C., Carvalho, C. S., Kurki-Suonio, H., Sirri, G. & Collaboration, E., 1-Apr-2019, In : Monthly Notices of the Royal Astronomical Society. 484, p. 5509-5529 21 p.

Research output: Contribution to journalArticleAcademicpeer-review

  • Mischa Knabenhans
  • Joachim Stadel
  • Stefano Marelli
  • Doug Potter
  • Romain Teyssier
  • Laurent Legrand
  • Aurel Schneider
  • Bruno Sudret
  • Linda Blot
  • Saeeda Awan
  • Carlo Burigana
  • Carla Sofia Carvalho
  • Hannu Kurki-Suonio
  • Gabriele Sirri
  • Euclid Collaboration
We present a new power spectrum emulator named EuclidEmulator that estimates the nonlinear correction to the linear dark matter power spectrum depending on the six cosmological parameters ωb, ωm, ns, h, w0, and σ8. It is constructed using the uncertainty quantification software UQLab using a spectral decomposition method called polynomial chaos expansion. All steps in its construction have been tested and optimized: the large high-resolution N-body simulations carried out with PKDGRAV3 were validated using a simulation from the Euclid Flagship campaign and demonstrated to have converged up to wavenumbers k≈ 5 h Mpc^{-1} for redshifts z ≤ 5. The emulator is based on 100 input cosmologies simulated in boxes of (1250 Mpc/h)3 using 20483 particles. We show that by creating mock emulators it is possible to successfully predict and optimize the performance of the final emulator prior to performing any N-body simulations. The absolute accuracy of the final nonlinear power spectrum is as good as one obtained with N-body simulations, conservatively, {̃ } 1 per cent for k≲ 1 h Mpc^{-1} and z ≲ 1. This enables efficient forward modelling in the nonlinear regime, allowing for estimation of cosmological parameters using Markov Chain Monte Carlo methods. EuclidEmulator has been compared to HALOFIT, CosmicEmu, and NGenHalofit, and shown to be more accurate than these other approaches. This work paves a new way for optimal construction of future emulators that also consider other cosmological observables, use higher resolution input simulations, and investigate higher dimensional cosmological parameter spaces.

Original languageEnglish
Pages (from-to)5509-5529
Number of pages21
JournalMonthly Notices of the Royal Astronomical Society
Publication statusPublished - 1-Apr-2019


  • methods: numerical, methods: statistical, cosmological parameters, large-scale structure of Universe, BARYON PHYSICS, SIMULATIONS, STATISTICS

Download statistics

No data available

ID: 118173244