Solving Large Scale Structure in Ten Easy Steps with COLA
read the original abstract
We present the COmoving Lagrangian Acceleration (COLA) method: an N-body method for solving for Large Scale Structure (LSS) in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory (LPT). Unlike standard N-body methods, the COLA method can straightforwardly trade accuracy at small-scales in order to gain computational speed without sacrificing accuracy at large scales. This is especially useful for cheaply generating large ensembles of accurate mock halo catalogs required to study galaxy clustering and weak lensing, as those catalogs are essential for performing detailed error analysis for ongoing and future surveys of LSS. As an illustration, we ran a COLA-based N-body code on a box of size 100Mpc/h with particles of mass ~5*10^9Msolar/h. Running the code with only 10 timesteps was sufficient to obtain an accurate description of halo statistics down to halo masses of at least 10^11Msolar/h. This is only at a modest speed penalty when compared to mocks obtained with LPT. A standard detailed N-body run is orders of magnitude slower than our COLA-based code. The speed-up we obtain with COLA is due to the fact that we calculate the large-scale dynamics exactly using LPT, while letting the N-body code solve for the small scales, without requiring it to capture exactly the internal dynamics of halos. Achieving a similar level of accuracy in halo statistics without the COLA method requires at least 3 times more timesteps than when COLA is employed.
This paper has not been read by Pith yet.
Forward citations
Cited by 9 Pith papers
-
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.
-
Constraints on Horndeski Gravity with Phantom Crossing
ACG models embed the observationally preferred phantom-crossing dark energy behavior inside a consistent Horndeski Lagrangian and achieve data fits of similar quality to w0waCDM while being narrowed by perturbative probes.
-
Learning the Universe: Posterior Reliability of Neural Generative Models in High-Dimensional Field-Level Inference of Cosmic Initial Conditions
Generative models for cosmological field-level inference can reproduce posterior means and cross-correlations yet fail to capture correct uncertainty geometry when validated against HMC reference samples.
-
The Manticore Project II: Bayesian digital twins of cosmic structure across the SDSS and BOSS volumes
Manticore-Deep uses tiled Bayesian field-level inference on SDSS and BOSS data to produce posterior ensembles of 3D cosmic fields that are consistent with LCDM and validated by 7.4σ CMB lensing and 3.5σ kSZ detections.
-
Dominated-Convergence Failure in Cosmological Perturbation Theory and a Numerical Foundation for BBGKY+ZA
Identifies dominated-convergence obstruction as root cause of divergences in EPT and LPT, and supplies numerical building blocks for a BBGKY+ZA approach that evades the issue.
-
Modeling nonlinear scales for dynamical dark energy cosmologies with COLA
COLA-based hybrid emulator reproduces nonlinear power spectrum boosts in w0wa models to <2% error vs EuclidEmulator2 and produces <0.3σ shifts in LSST-like cosmic shear parameter constraints.
-
Impact of inhomogeneous curvature on growth rate measurements from magnitude fluctuations
Full-GR simulations find that inhomogeneous curvature produces only sub-dominant systematic offsets in growth-rate measurements from magnitude fluctuations at z ≲ 0.2 relative to current statistical errors.
-
The sound of dynamical dark energy and modified gravity
Using CMB, SN, BAO and shear data, the work shows dynamical dark energy in MG models correlates with deviations from GR below z=2 at >95% CL, a link that holds for varying sound speed but vanishes for a cosmological constant.
-
Large-scale structures of the Universe: physics, phenomenology, statistics
Lecture series on the physics, phenomenology, and statistics of large-scale cosmic structure evolution and non-Gaussian predictions.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.