Transients from initial conditions based on Lagrangian perturbation theory in N-body simulations

Takayuki Tatekawa, Shuntaro Mizuno

Research output: Contribution to journalArticle

15 Citations (Scopus)

Abstract

We explore the initial conditions for cosmological N-body simulations suitable for calculating the skewness and kurtosis of the density field. In general, the initial conditions based on the perturbation theory provide incorrect second-order and higher-order growth. These errors implied by the use of the perturbation theory to set up the initial conditions in N-body simulations are called transients. Unless these transients are completely suppressed compared with the dominant growing mode, we cannot reproduce the correct evolution of cumulants with orders higher than two, even though there is no problem with the numerical scheme. We investigate the impact of transients on the observable statistical quantities by performing N-body simulations with initial conditions based on Lagrangian perturbation theory. We show that the effects of transients on the kurtosis from the initial conditions, based on second-order Lagrangian perturbation theory (2LPT), have almost disappeared by z∼5, as long as the initial conditions are set at z>30. This means that for practical purposes, the initial conditions based on 2LPT are accurate enough for numerical calculations of skewness and kurtosis.

Original languageEnglish
Article number014
JournalJournal of Cosmology and Astroparticle Physics
Volume2007
Issue number12
DOIs
Publication statusPublished - 2007 Dec 1
Externally publishedYes

Fingerprint

perturbation theory
kurtosis
simulation
skewness

Keywords

  • Cosmological simulations
  • N-body simulations

ASJC Scopus subject areas

  • Physics and Astronomy(all)

Cite this

Transients from initial conditions based on Lagrangian perturbation theory in N-body simulations. / Tatekawa, Takayuki; Mizuno, Shuntaro.

In: Journal of Cosmology and Astroparticle Physics, Vol. 2007, No. 12, 014, 01.12.2007.

Research output: Contribution to journalArticle

@article{3aefe5e15f2649998c9ca7436c296b22,
title = "Transients from initial conditions based on Lagrangian perturbation theory in N-body simulations",
abstract = "We explore the initial conditions for cosmological N-body simulations suitable for calculating the skewness and kurtosis of the density field. In general, the initial conditions based on the perturbation theory provide incorrect second-order and higher-order growth. These errors implied by the use of the perturbation theory to set up the initial conditions in N-body simulations are called transients. Unless these transients are completely suppressed compared with the dominant growing mode, we cannot reproduce the correct evolution of cumulants with orders higher than two, even though there is no problem with the numerical scheme. We investigate the impact of transients on the observable statistical quantities by performing N-body simulations with initial conditions based on Lagrangian perturbation theory. We show that the effects of transients on the kurtosis from the initial conditions, based on second-order Lagrangian perturbation theory (2LPT), have almost disappeared by z∼5, as long as the initial conditions are set at z>30. This means that for practical purposes, the initial conditions based on 2LPT are accurate enough for numerical calculations of skewness and kurtosis.",
keywords = "Cosmological simulations, N-body simulations",
author = "Takayuki Tatekawa and Shuntaro Mizuno",
year = "2007",
month = "12",
day = "1",
doi = "10.1088/1475-7516/2007/12/014",
language = "English",
volume = "2007",
journal = "Journal of Cosmology and Astroparticle Physics",
issn = "1475-7516",
publisher = "IOP Publishing Ltd.",
number = "12",

}

TY - JOUR

T1 - Transients from initial conditions based on Lagrangian perturbation theory in N-body simulations

AU - Tatekawa, Takayuki

AU - Mizuno, Shuntaro

PY - 2007/12/1

Y1 - 2007/12/1

N2 - We explore the initial conditions for cosmological N-body simulations suitable for calculating the skewness and kurtosis of the density field. In general, the initial conditions based on the perturbation theory provide incorrect second-order and higher-order growth. These errors implied by the use of the perturbation theory to set up the initial conditions in N-body simulations are called transients. Unless these transients are completely suppressed compared with the dominant growing mode, we cannot reproduce the correct evolution of cumulants with orders higher than two, even though there is no problem with the numerical scheme. We investigate the impact of transients on the observable statistical quantities by performing N-body simulations with initial conditions based on Lagrangian perturbation theory. We show that the effects of transients on the kurtosis from the initial conditions, based on second-order Lagrangian perturbation theory (2LPT), have almost disappeared by z∼5, as long as the initial conditions are set at z>30. This means that for practical purposes, the initial conditions based on 2LPT are accurate enough for numerical calculations of skewness and kurtosis.

AB - We explore the initial conditions for cosmological N-body simulations suitable for calculating the skewness and kurtosis of the density field. In general, the initial conditions based on the perturbation theory provide incorrect second-order and higher-order growth. These errors implied by the use of the perturbation theory to set up the initial conditions in N-body simulations are called transients. Unless these transients are completely suppressed compared with the dominant growing mode, we cannot reproduce the correct evolution of cumulants with orders higher than two, even though there is no problem with the numerical scheme. We investigate the impact of transients on the observable statistical quantities by performing N-body simulations with initial conditions based on Lagrangian perturbation theory. We show that the effects of transients on the kurtosis from the initial conditions, based on second-order Lagrangian perturbation theory (2LPT), have almost disappeared by z∼5, as long as the initial conditions are set at z>30. This means that for practical purposes, the initial conditions based on 2LPT are accurate enough for numerical calculations of skewness and kurtosis.

KW - Cosmological simulations

KW - N-body simulations

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

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

U2 - 10.1088/1475-7516/2007/12/014

DO - 10.1088/1475-7516/2007/12/014

M3 - Article

AN - SCOPUS:38049049284

VL - 2007

JO - Journal of Cosmology and Astroparticle Physics

JF - Journal of Cosmology and Astroparticle Physics

SN - 1475-7516

IS - 12

M1 - 014

ER -