Improving initial conditions for cosmological N-body simulations

Lehman H. Garrison, Daniel J. Eisenstein, Douglas Ferrer, Marc V. Metchnik, Philip A Pinto

Research output: Contribution to journalArticle

18 Citations (Scopus)

Abstract

In cosmological N-body simulations, the representation of dark matter as discrete 'macroparticles' suppresses the growth of structure, such that simulations no longer reproduce linear theory on small scales near kNyquist. Marcos et al. demonstrate that this is due to sparse sampling of modes near kNyquist and that the often-assumed continuum growing modes are not proper growing modes of the particle system. We develop initial conditions (ICs) that respect the particle linear theory growing modes and then rescale the mode amplitudes to account for growth suppression. These ICs also allow us to take advantage of our very accurate N-body code ABACUS to implement second-order Lagrangian perturbation theory (2LPT) in configuration space. The combination of 2LPT and rescaling improves the accuracy of the late-time power spectra, halo mass functions, and halo clustering. In particular, we achieve 1 per cent accuracy in the power spectrum down to kNyquist, versus kNyquist/4 without rescaling or kNyquist/13 without 2LPT, relative to an oversampled reference simulation.We anticipate that our 2LPT will be useful for large simulations where fast Fourier transforms are expensive and that rescaling will be useful for suites of medium-resolution simulations used in cosmic emulators and galaxy survey mock catalogues. Code to generate ICs is available at https://github.com/lgarrison/zeldovich-PLT.

Original languageEnglish (US)
Pages (from-to)4125-4145
Number of pages21
JournalMonthly Notices of the Royal Astronomical Society
Volume461
Issue number4
DOIs
StatePublished - Oct 1 2016

Fingerprint

simulation
power spectra
halos
Fourier transform
catalogs
dark matter
perturbation theory
sampling
perturbation
retarding
galaxies
continuums
configurations
code
particle

Keywords

  • Galaxies: haloes
  • Large-scale structure of Universe
  • Methods: numerical

ASJC Scopus subject areas

  • Astronomy and Astrophysics
  • Space and Planetary Science

Cite this

Improving initial conditions for cosmological N-body simulations. / Garrison, Lehman H.; Eisenstein, Daniel J.; Ferrer, Douglas; Metchnik, Marc V.; Pinto, Philip A.

In: Monthly Notices of the Royal Astronomical Society, Vol. 461, No. 4, 01.10.2016, p. 4125-4145.

Research output: Contribution to journalArticle

Garrison, Lehman H. ; Eisenstein, Daniel J. ; Ferrer, Douglas ; Metchnik, Marc V. ; Pinto, Philip A. / Improving initial conditions for cosmological N-body simulations. In: Monthly Notices of the Royal Astronomical Society. 2016 ; Vol. 461, No. 4. pp. 4125-4145.
@article{029b566023c2474fbc66cd3fd13b46a8,
title = "Improving initial conditions for cosmological N-body simulations",
abstract = "In cosmological N-body simulations, the representation of dark matter as discrete 'macroparticles' suppresses the growth of structure, such that simulations no longer reproduce linear theory on small scales near kNyquist. Marcos et al. demonstrate that this is due to sparse sampling of modes near kNyquist and that the often-assumed continuum growing modes are not proper growing modes of the particle system. We develop initial conditions (ICs) that respect the particle linear theory growing modes and then rescale the mode amplitudes to account for growth suppression. These ICs also allow us to take advantage of our very accurate N-body code ABACUS to implement second-order Lagrangian perturbation theory (2LPT) in configuration space. The combination of 2LPT and rescaling improves the accuracy of the late-time power spectra, halo mass functions, and halo clustering. In particular, we achieve 1 per cent accuracy in the power spectrum down to kNyquist, versus kNyquist/4 without rescaling or kNyquist/13 without 2LPT, relative to an oversampled reference simulation.We anticipate that our 2LPT will be useful for large simulations where fast Fourier transforms are expensive and that rescaling will be useful for suites of medium-resolution simulations used in cosmic emulators and galaxy survey mock catalogues. Code to generate ICs is available at https://github.com/lgarrison/zeldovich-PLT.",
keywords = "Galaxies: haloes, Large-scale structure of Universe, Methods: numerical",
author = "Garrison, {Lehman H.} and Eisenstein, {Daniel J.} and Douglas Ferrer and Metchnik, {Marc V.} and Pinto, {Philip A}",
year = "2016",
month = "10",
day = "1",
doi = "10.1093/mnras/stw1594",
language = "English (US)",
volume = "461",
pages = "4125--4145",
journal = "Monthly Notices of the Royal Astronomical Society",
issn = "0035-8711",
publisher = "Oxford University Press",
number = "4",

}

TY - JOUR

T1 - Improving initial conditions for cosmological N-body simulations

AU - Garrison, Lehman H.

AU - Eisenstein, Daniel J.

AU - Ferrer, Douglas

AU - Metchnik, Marc V.

AU - Pinto, Philip A

PY - 2016/10/1

Y1 - 2016/10/1

N2 - In cosmological N-body simulations, the representation of dark matter as discrete 'macroparticles' suppresses the growth of structure, such that simulations no longer reproduce linear theory on small scales near kNyquist. Marcos et al. demonstrate that this is due to sparse sampling of modes near kNyquist and that the often-assumed continuum growing modes are not proper growing modes of the particle system. We develop initial conditions (ICs) that respect the particle linear theory growing modes and then rescale the mode amplitudes to account for growth suppression. These ICs also allow us to take advantage of our very accurate N-body code ABACUS to implement second-order Lagrangian perturbation theory (2LPT) in configuration space. The combination of 2LPT and rescaling improves the accuracy of the late-time power spectra, halo mass functions, and halo clustering. In particular, we achieve 1 per cent accuracy in the power spectrum down to kNyquist, versus kNyquist/4 without rescaling or kNyquist/13 without 2LPT, relative to an oversampled reference simulation.We anticipate that our 2LPT will be useful for large simulations where fast Fourier transforms are expensive and that rescaling will be useful for suites of medium-resolution simulations used in cosmic emulators and galaxy survey mock catalogues. Code to generate ICs is available at https://github.com/lgarrison/zeldovich-PLT.

AB - In cosmological N-body simulations, the representation of dark matter as discrete 'macroparticles' suppresses the growth of structure, such that simulations no longer reproduce linear theory on small scales near kNyquist. Marcos et al. demonstrate that this is due to sparse sampling of modes near kNyquist and that the often-assumed continuum growing modes are not proper growing modes of the particle system. We develop initial conditions (ICs) that respect the particle linear theory growing modes and then rescale the mode amplitudes to account for growth suppression. These ICs also allow us to take advantage of our very accurate N-body code ABACUS to implement second-order Lagrangian perturbation theory (2LPT) in configuration space. The combination of 2LPT and rescaling improves the accuracy of the late-time power spectra, halo mass functions, and halo clustering. In particular, we achieve 1 per cent accuracy in the power spectrum down to kNyquist, versus kNyquist/4 without rescaling or kNyquist/13 without 2LPT, relative to an oversampled reference simulation.We anticipate that our 2LPT will be useful for large simulations where fast Fourier transforms are expensive and that rescaling will be useful for suites of medium-resolution simulations used in cosmic emulators and galaxy survey mock catalogues. Code to generate ICs is available at https://github.com/lgarrison/zeldovich-PLT.

KW - Galaxies: haloes

KW - Large-scale structure of Universe

KW - Methods: numerical

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

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

U2 - 10.1093/mnras/stw1594

DO - 10.1093/mnras/stw1594

M3 - Article

AN - SCOPUS:84988624568

VL - 461

SP - 4125

EP - 4145

JO - Monthly Notices of the Royal Astronomical Society

JF - Monthly Notices of the Royal Astronomical Society

SN - 0035-8711

IS - 4

ER -