Matter power spectrum and the challenge of percent accuracy

Aurel Schneider, Romain Teyssier, Doug Potter, Joachim Stadel, Julian Onions, Darren S. Reed, Robert E. Smith, Volker Springel, Frazer R. Pearce, Roman Scoccimarro

    Research output: Contribution to journalArticle

    Abstract

    Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day N-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used N-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at k1 h Mpc-1 and to within three percent at k10 h Mpc-1. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for k 2 h Mpc-1. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than L=0.5 h-1Gpc to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of Mp=109 h-1Mo is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

    Original languageEnglish (US)
    Article number047
    JournalJournal of Cosmology and Astroparticle Physics
    Volume2016
    Issue number4
    DOIs
    StatePublished - 2016

    Fingerprint

    power spectra
    boxes
    simulation
    particle mass
    mesh
    synchronism
    coverings
    perturbation theory
    galaxies

    Keywords

    • cosmological simulations
    • power spectrum

    ASJC Scopus subject areas

    • Astronomy and Astrophysics

    Cite this

    Schneider, A., Teyssier, R., Potter, D., Stadel, J., Onions, J., Reed, D. S., ... Scoccimarro, R. (2016). Matter power spectrum and the challenge of percent accuracy. Journal of Cosmology and Astroparticle Physics, 2016(4), [047]. https://doi.org/10.1088/1475-7516/2016/04/047

    Matter power spectrum and the challenge of percent accuracy. / Schneider, Aurel; Teyssier, Romain; Potter, Doug; Stadel, Joachim; Onions, Julian; Reed, Darren S.; Smith, Robert E.; Springel, Volker; Pearce, Frazer R.; Scoccimarro, Roman.

    In: Journal of Cosmology and Astroparticle Physics, Vol. 2016, No. 4, 047, 2016.

    Research output: Contribution to journalArticle

    Schneider, A, Teyssier, R, Potter, D, Stadel, J, Onions, J, Reed, DS, Smith, RE, Springel, V, Pearce, FR & Scoccimarro, R 2016, 'Matter power spectrum and the challenge of percent accuracy', Journal of Cosmology and Astroparticle Physics, vol. 2016, no. 4, 047. https://doi.org/10.1088/1475-7516/2016/04/047
    Schneider, Aurel ; Teyssier, Romain ; Potter, Doug ; Stadel, Joachim ; Onions, Julian ; Reed, Darren S. ; Smith, Robert E. ; Springel, Volker ; Pearce, Frazer R. ; Scoccimarro, Roman. / Matter power spectrum and the challenge of percent accuracy. In: Journal of Cosmology and Astroparticle Physics. 2016 ; Vol. 2016, No. 4.
    @article{7f4bb20dd4424ad3941256ec8275b455,
    title = "Matter power spectrum and the challenge of percent accuracy",
    abstract = "Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day N-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used N-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at k1 h Mpc-1 and to within three percent at k10 h Mpc-1. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for k 2 h Mpc-1. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than L=0.5 h-1Gpc to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of Mp=109 h-1Mo is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.",
    keywords = "cosmological simulations, power spectrum",
    author = "Aurel Schneider and Romain Teyssier and Doug Potter and Joachim Stadel and Julian Onions and Reed, {Darren S.} and Smith, {Robert E.} and Volker Springel and Pearce, {Frazer R.} and Roman Scoccimarro",
    year = "2016",
    doi = "10.1088/1475-7516/2016/04/047",
    language = "English (US)",
    volume = "2016",
    journal = "Journal of Cosmology and Astroparticle Physics",
    issn = "1475-7516",
    publisher = "IOP Publishing Ltd.",
    number = "4",

    }

    TY - JOUR

    T1 - Matter power spectrum and the challenge of percent accuracy

    AU - Schneider, Aurel

    AU - Teyssier, Romain

    AU - Potter, Doug

    AU - Stadel, Joachim

    AU - Onions, Julian

    AU - Reed, Darren S.

    AU - Smith, Robert E.

    AU - Springel, Volker

    AU - Pearce, Frazer R.

    AU - Scoccimarro, Roman

    PY - 2016

    Y1 - 2016

    N2 - Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day N-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used N-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at k1 h Mpc-1 and to within three percent at k10 h Mpc-1. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for k 2 h Mpc-1. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than L=0.5 h-1Gpc to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of Mp=109 h-1Mo is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

    AB - Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day N-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used N-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at k1 h Mpc-1 and to within three percent at k10 h Mpc-1. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for k 2 h Mpc-1. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than L=0.5 h-1Gpc to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of Mp=109 h-1Mo is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

    KW - cosmological simulations

    KW - power spectrum

    UR - http://www.scopus.com/inward/record.url?scp=85012878162&partnerID=8YFLogxK

    UR - http://www.scopus.com/inward/citedby.url?scp=85012878162&partnerID=8YFLogxK

    U2 - 10.1088/1475-7516/2016/04/047

    DO - 10.1088/1475-7516/2016/04/047

    M3 - Article

    AN - SCOPUS:85012878162

    VL - 2016

    JO - Journal of Cosmology and Astroparticle Physics

    JF - Journal of Cosmology and Astroparticle Physics

    SN - 1475-7516

    IS - 4

    M1 - 047

    ER -