REVIEW 9 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 9 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.
-
Post-collapse Lagrangian perturbation theory in three dimensions
The authors introduce post-collapse Lagrangian perturbation theory (PCPT) that asymptotically captures early post-shell-crossing dynamics in 3D by using high-order LPT until crossing and then 1D pancake results for gr...
-
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.
-
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.
-
Parity Violation in Galaxy Shapes: Primordial Non-Gaussianity
The parity-odd intrinsic alignment power spectrum probes the collapsed limit of the parity-odd primordial trispectrum and can tighten constraints on parity-violating PNG when bias parameters are calibrated from N-body...
-
$\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...
-
Constraints on Neutrino Physics from DESI DR2 BAO and DR1 Full Shape
DESI DR2 BAO and full-shape data plus CMB yield ∑m_ν < 0.0642 eV (95% CL) under ΛCDM, in 3σ tension with oscillation lower limits, relaxed to <0.163 eV in w0waCDM.
Discussion (0). Sign in to comment.