ICE-COLA: Towards fast and accurate synthetic galaxy catalogues optimizing a quasi $N$-body method

Abstract

Next generation galaxy surveys demand the development of massive ensembles of galaxy mocks to model the observables and their covariances, what is computationally prohibitive using $N$-body simulations. COLA is a novel method designed to make this feasible by following an approximate dynamics but with up to 3 orders of magnitude speed-ups when compared to an exact $N$-body. In this paper we investigate the optimization of the code parameters in the compromise between computational cost and recovered accuracy in observables such as two-point clustering and halo abundance. We benchmark those observables with a state-of-the-art $N$-body run, the MICE Grand Challenge simulation (MICE-GC). We find that using 40 time steps linearly spaced since $z_i\sim 20$, and a force mesh resolution three times finer than that of the number of particles, yields a matter power spectrum within $1\%$ for $k\lesssim 1h{\mathrm{Mpc}}^{-1}$ and a halo mass function within $5\%$ of those in the $N$-body. In turn the halo bias is accurate within $2\%$ for $k\lesssim 0.7h{\mathrm{Mpc}}^{-1}$ whereas, in redshift space, the halo monopole and quadrupole are within $4\%$ for $k\lesssim 0.4h{\mathrm{Mpc}}^{-1}$. These results hold for a broad range in redshift ($0<z<1$) and for all halo mass bins investigated ($M>10^{12.5}{h}^{-1}{\mathrm{M}}_{\odot }$). To bring accuracy in clustering to one percent level we study various methods that re-calibrate halo masses and/or velocities. We thus propose an optimized choice of COLA code parameters as a powerful tool to optimally exploit future galaxy surveys.

16 pages, 13 figures; matches the version accepted by MNRAS; some results have been extended to higher redshifts; we added Appendix B

References in corpus (16)

• The Clustering of the SDSS DR7 Main Galaxy Sample I: A 4 per cent Distance Measure at z=0.15
• Transients from Initial Conditions in Cosmological Simulations
• The Clustering of the SDSS Main Galaxy Sample II: Mock galaxy catalogues and a measurement of the growth of structure from Redshift Space Distortions at $z=0.15$
• FastPM: a new scheme for fast simulations of dark matter and halos
• An algorithm to build mock galaxy catalogues using MICE simulations
• N-body simulations of gravitational dynamics
• Simulations of Baryon Acoustic Oscillations II: Covariance matrix of the matter power spectrum
• Dark Sky Simulations: Early Data Release
• nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics
• Fast and accurate mock catalogue generation for low-mass galaxies
• Constraining the halo bispectrum in real and redshift space from perturbation theory and non-linear stochastic bias
• Bayesian analysis of the dynamic cosmic web in the SDSS galaxy survey
• Matter Power Spectrum Covariance Matrix from the DEUS-PUR ΛCDM simulations: Mass Resolution and non-Gaussian Errors
• Comparing halo bias from abundance and clustering
• sCOLA: The N-body COLA Method Extended to the Spatial Domain
• DEUS Full Observable ΛCDM Universe Simulation: the numerical challenge

Cited by in corpus (81)

• FastPM: a new scheme for fast simulations of dark matter and halos
• Large-scale dark matter simulations
• The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: 1000 multi-tracer mock catalogues with redshift evolution and systematics for galaxies and quasars of the final data release
• Galaxy Power Spectrum Multipoles Covariance in Perturbation Theory
• The Websky Extragalactic CMB Simulations
• Dark Energy Survey Year 3 Results: A 2.7% measurement of Baryon Acoustic Oscillation distance scale at redshift 0.835
• Cosmological constraints from BOSS with analytic covariance matrices
• The mass-Peak Patch algorithm for fast generation of deep all-sky dark matter halo catalogues and its N-Body validation
• COLA with scale-dependent growth: applications to screened modified gravity models
• Comparing approximate methods for mock catalogues and covariance matrices II: Power spectrum multipoles
• The Effect of Fiber Collisions on the Galaxy Power Spectrum Multipole
• Comparing approximate methods for mock catalogues and covariance matrices I: correlation function
• Comparing approximate methods for mock catalogues and covariance matrices III: Bispectrum
• Dark matter statistics for large galaxy catalogs: power spectra and covariance matrices
• Testing the theory of gravity with DESI: estimators, predictions and simulation requirements
• Improving fast generation of halo catalogs with higher-order Lagrangian perturbation theory
• Accurate cosmic shear errors: do we need ensembles of simulations?
• The BACCO simulation project: biased tracers in real space
• Theoretical Systematics of Future Baryon Acoustic Oscillation Surveys
• CARPool: fast, accurate computation of large-scale structure statistics by pairing costly and cheap cosmological simulations
• A volumetric deep Convolutional Neural Network for simulation of mock dark matter halo catalogues
• FlowPM: Distributed TensorFlow Implementation of the FastPM Cosmological N-body Solver
• Field Level Neural Network Emulator for Cosmological N-body Simulations
• Dark Energy Survey: A 2.1% measurement of the angular Baryonic Acoustic Oscillation scale at redshift $z_{\mathrm{eff}}$=0.85 from the final dataset
• Correlations between Triaxial Shapes and Formation History of Dark Matter Haloes
• BAM: Bias Assignment Method to generate mock catalogs
• Fast full N-body simulations of generic modified gravity: conformal coupling models
• A gradient based method for modeling baryons and matter in halos of fast simulations
• CoLoRe: fast cosmological realisations over large volumes with multiple tracers
• ICE-COLA: fast simulations for weak lensing observables
• Non-linear contributions to angular power spectra
• Fast full $N$-body simulations of generic modified gravity: derivative coupling models
• Beyond two-point statistics: using the minimum spanning tree as a tool for cosmology
• Nonlinear 3D Cosmic Web Simulation with Heavy-Tailed Generative Adversarial Networks
• $\text{Hi-COLA}$: Fast, approximate simulations of structure formation in Horndeski gravity
• Euclid: Modelling massive neutrinos in cosmology -- a code comparison
• Suppressing cosmic variance with paired-and-fixed cosmological simulations: average properties and covariances of dark matter clustering statistics
• Enabling matter power spectrum emulation in beyond-$\Lambda$CDM cosmologies with COLA
• A Bayesian method for combining theoretical and simulated covariance matrices for large-scale structure surveys
• The bias of dark matter tracers: assessing the accuracy of mapping techniques
• Modelling galaxy clustering in redshift space with a Lagrangian bias formalism and $N$-body simulations
• Euclid preparation. XLI. Galaxy power spectrum modelling in real space
• Does jackknife scale really matter for accurate large-scale structure covariances?
• Testing approximate predictions of displacements of cosmological dark matter halos
• Fast generation of mock galaxy catalogues in modified gravity models with COLA
• The cosmic web connection to the dark matter halo distribution through gravity
• Dark Energy Survey Year 3 Results: Measurement of the Baryon Acoustic Oscillations with Three-dimensional Clustering
• Towards Accurate Field-Level Inference of Massive Cosmic Structures
• Efficient parallel algorithm for estimating higher-order polyspectra
• Accurate halo-galaxy mocks from automatic bias estimation and particle mesh gravity solvers
• Estimating Cosmological Parameters from the Dark Matter Distribution
• Clustering with general photo-$z$ uncertainties: Application to Baryon Acoustic Oscillations
• Fast generation of weak lensing maps by the inverse-Gaussianization method
• Sesame: A power spectrum emulator pipeline for beyond-$\Lambda$CDM models
• Perfectly parallel cosmological simulations using spatial comoving Lagrangian acceleration
• Revisiting Vainshtein Screening for fast N-body simulations
• Bayesian Control Variates for optimal covariance estimation with pairs of simulations and surrogates
• Assessing non-linear models for galaxy clustering III: Theoretical accuracy for Stage IV surveys
• $\text{PineTree}$: A generative, fast, and differentiable halo model for wide-field galaxy surveys
• Primordial non-Gaussianity with Angular correlation function: Integral constraint and validation for DES
• Fast numerical method to generate halo catalogs in modified gravity (part I): second-order Lagrangian Perturbation Theory
• BullFrog: Multi-step perturbation theory as a time integrator for cosmological simulations
• A field-level emulator for modified gravity
• The halo 3-point correlation function: a methodological analysis
• Accurate relativistic observables from post-processing light cone catalogues
• Perturbation-theory informed integrators for cosmological simulations
• COmoving Computer Acceleration (COCA): $N$-body simulations in an emulated frame of reference
• Cosmological inference including massive neutrinos from the matter power spectrum: biases induced by uncertainties in the covariance matrix
• Covariance matrices for variance-suppressed simulations
• Fast simulation mapping: from standard to modified gravity cosmologies using the bias assignment method
• Including relativistic and primordial Non-Gaussianity contributions in cosmological simulations by modifying the initial condition
• Optimizing Redshift Distribution Inference through Joint Self-Calibration and Clustering-Redshift Synergy
• Approximate methods for the generation of dark matter halo catalogs in the age of precision cosmology
• K-mouflage at high $k$: extending the reach of $\text{Hi-COLA}$
• Dark Energy Survey: DESI-Independent Angular BAO Measurement
• Improving the Accuracy of Halo Mass Based Statistics For Fast Approximate N-body Simulations
• Fast Generation of Weak Lensing Maps in Modified Gravity with COLA
• Dark Energy Survey Year 3 Results: Galaxy mock catalogs for BAO analysis
• Hybrid $P_{\ell }(k)$: general, unified, non-linear matter power spectrum in redshift space
• The Manticore Project II: Bayesian digital twins of cosmic structure across the SDSS and BOSS volumes
• Generative modeling of convergence maps based on predicted one-point statistics