REVIEW 6 cited by
Transients from Initial Conditions: A Perturbative Analysis
Not yet reviewed by Pith; the record is open.
This paper has not been read by Pith yet. Machine review is queued; the pith claim, tier, and objections will appear here once it completes.
SPECIMEN: schema-true, not a live event
T0 review · schema-true
One-sentence machine reading of the paper's core claim.
pith:XXXXXXXX · record.json · timestamp
Signed reviews
read the original abstract
The standard procedure to generate initial conditions (IC) in numerical simulations 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 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 perturbation theory (PT) to include transients as non-linear excitations of decaying modes caused by the IC. We focus on higher-order statistics of the density contrast and velocity divergence, characterized by the S_p and T_p parameters. We find that the time-scale of transients is determined, at a given order p, by the spectral index n. The skewness factor S_3 (T_3) attains 10% accuracy only after a=6 (a=15) for n=0, whereas higher (lower) n demands more (less) expansion away from the IC. These requirements become much more stringent as p increases. An Omega=0.3 model requires a factor of two larger expansion than an Omega=1 model to reduce transients by the same amount. The predicted transients in S_p are in good agreement with numerical simulations. More accurate IC can be achieved by using 2nd order Lagrangian PT (2LPT), which reproduces growing modes up to 2nd order and thus eliminates transients in the skewness. We show that for p>3 this reduces the required expansion by more than an order of magnitude compared to the ZA. Setting up 2LPT IC only requires minimal, inexpensive changes to ZA codes. We suggest simple steps for its implementation.
Forward citations
Cited by 6 Pith papers
-
Generating Moving Field Initial Conditions with Spatially Varying Boost
A 'spatially varying boost' algorithm assigns arbitrary, position-dependent bulk velocities to field initial data by composing local Lorentz boosts, demonstrated on solitons, Proca fields, and spin-1 wave dark matter.
-
Perturbation-theory informed integrators for cosmological simulations
LPT-matched integrators for cosmological simulations outperform FastPM with O(1-100) timesteps while convergence is limited to order 3/2 post-shell-crossing due to acceleration field irregularity.
-
Renormalized Perturbation Theory at Field-level: the LSS bootstrap in GridSPT
A renormalized field-level perturbation theory is shown to recover the LSS bootstrap parameter consistently across different grid cutoffs, validated at third and fifth order against N-body simulations.
-
Shapes and orientations of massive halos in the statistically anisotropic universe
Quadrupolar statistical anisotropy in initial conditions leaves cluster shapes unchanged but aligns the major axes of massive halos, more strongly for higher mass.
-
Biased tracers, Hybrid Effective Field Theory and Modified Gravity
Extends HEFT to f(R) gravity for loop-corrected biased power spectra and outlines emulator extensions from LambdaCDM codes.
-
Sahyadri: A simulation suite for the cosmology dependence of the Cosmic Web
Sahyadri is a new N-body simulation suite with 25x better halo mass resolution than AbacusSummit and seed-matched parameter variations, showing coherent Omega_m signals in power spectra, mass functions, and cosmic-web...
Discussion (0). Continue with ORCID to comment.