## Abstract

The standard procedure to generate initial conditions in numerical simulations of structure formations is to use the Zel'dovich approximation (ZA). Although the ZA correctly reproduces the linear growing modes of density and velocity perturbations, non-linear growth is inaccurately represented, particularly for velocity perturbations because of the ZA failure to conserve momentum. This implies that it takes time for the actual dynamics to establish the correct statistical properties of density and velocity fields. We extend the standard formulation of non-linear perturbation theory (PT) to include transients as non-linear excitations of decaying modes caused by the initial conditions. These new non-linear solutions interpolate between the initial conditions and the late-time solutions given by the exact non-linear dynamics. To quantify the magnitude of transients, we focus on higher order statistics of the density contrast δ and velocity divergence Θ, characterized by the S_{p} and T_{p} parameters. These describe the non-Gaussianity of the probability distribution through its connected moments 〈δ^{p}〉_{c} ≡ S_{p} 〈δ^{2}〉^{p-1}, 〈Θ^{p}〉_{c}& ≡ T_{p} 〈Θ^{2}〉^{p-1}. We calculate S_{p}(a) and T_{p}(a) to leading order in PT with top-hat smoothing as a function of the scale factor a. We find that the time-scale of transients is determined, at a given order p, by the effective spectral index n_{eff}. The skewness factor S_{3}(T_{3}) attains 10 per cent accuracy only after a ≈ 6 (a ≈ 15) for n_{eff} ≈ 0, whereas higher (lower) n_{eff} demands more (less) expansion away from the initial conditions. These requirements become much more stringent as p increases, always showing slower decay of transients for T_{p} than S_{p}. For models with density parameter Ω ≠ 1, the conditions above apply to the linear growth factor; thus an Ω = 0.3 open model requires roughly a factor of 2 larger expansion than a critical density model to reduce transients by the same amount. The predicted transients in S_{p} are in good agreement with numerical simulations. More accurate initial conditions can be achieved by using second-order Lagrangian PT (2LPT), which reproduces growing modes up to second order and thus eliminates transients in the skewness parameters. We show that for p > 3 this scheme can reduce the required expansion by more than an order of magnitude compared to the ZA. Setting up 2LPT initial conditions requires only minimal, inexpensive changes to ZA codes. We suggest simple steps for its implementation.

Original language | English (US) |
---|---|

Pages (from-to) | 1097-1118 |

Number of pages | 22 |

Journal | Monthly Notices of the Royal Astronomical Society |

Volume | 299 |

Issue number | 4 |

DOIs | |

State | Published - Oct 1 1998 |

## Keywords

- Galaxies: clusters: general
- Large-scale structure of Universe
- Methods: numerical
- Methods: statistical

## ASJC Scopus subject areas

- Astronomy and Astrophysics
- Space and Planetary Science