REVIEW 6 cited by
Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
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
read the original abstract
Inaccuracies in the initial conditions for cosmological N-body simulations could easily be the largest source of systematic error in predicting the non-linear large-scale structure. From the theory side, initial conditions are usually provided by using low-order truncations of the displacement field from Lagrangian perturbation theory, with the first and second-order approximations being the most common ones. Here we investigate the improvement brought by using initial conditions based on third-order Lagrangian perturbation theory (3LPT). We show that with 3LPT, truncation errors are vastly suppressed, thereby opening the portal to initializing simulations accurately as late as z=12 (for the resolution we consider). We analyse the competing effects of perturbative truncation and particle discreteness on various summary statistics. Discreteness errors are essentially decaying modes and thus get strongly amplified for earlier initialization times. We show that late starting times with 3LPT provide the most accurate configuration, which we find to coincide with the continuum fluid limit within 1 per cent for the power- and bispectrum at z=0 up to the particle Nyquist wave number of our simulations (k ~ 3h/Mpc). In conclusion, to suppress non-fluid artefacts, we recommend initializing simulations as late as possible with 3LPT. We make our 3LPT initial condition generator publicly available.
Forward citations
Cited by 6 Pith papers
-
Ab Initio Cosmological Simulations: From Inflation to Present-Day Structure Formation
Nonlinear axion-U(1) inflation, simulated on a lattice and fed into N-body simulations, boosts small-scale structure and high-redshift halo counts by up to ~200% relative to standard single-field initial conditions.
-
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.
-
gevolution 2.0: GPU-accelerated relativistic N-body simulations for cosmology
gevolution 2.0 accelerates relativistic N-body cosmological simulations on Nvidia GPUs, achieving a ~4.5x node-hour speed-up over its CPU-only predecessor.
-
DISCO-DJ II: a differentiable particle-mesh code for cosmology
A GPU-accelerated, differentiable particle-mesh N-body code achieves per-cent-level power-spectrum accuracy with few time steps and recovers sigma_8 plus initial conditions from a noisy mock field.
-
$\texttt{GENGARS}$: Accurate non-Gaussian initial conditions with arbitrary bispectrum for N-body simulations
GENGARS uses a Schwinger-parameterized reduced bispectrum kernel to generate N-body initial conditions for arbitrary separable PNG shapes, reducing spurious power-spectrum contributions relative to 2LPT-PNG.
-
Cosmological constraining power of the redshifts, heights, and angular clustering of weak gravitational lensing peaks
In an idealized Euclid-like simulation forecast, the combined redshift distribution, height distribution, and angular clustering of weak-lensing peaks constrain cosmology better than shear two-point correlation functi...
Discussion (0). Sign in to comment.