Perfectly parallel cosmological simulations using spatial comoving Lagrangian acceleration

Abstract

Existing cosmological simulation methods lack a high degree of parallelism due to the long-range nature of the gravitational force, which limits the size of simulations that can be run at high resolution. To solve this problem, we propose a new, perfectly parallel approach to simulate cosmic structure formation, which is based on the spatial COmoving Lagrangian Acceleration (sCOLA) framework. Building upon a hybrid analytical and numerical description of particles' trajectories, our algorithm allows for an efficient tiling of a cosmological volume, where the dynamics within each tile is computed independently. As a consequence, the degree of parallelism is equal to the number of tiles. We optimised the accuracy of sCOLA through the use of a buffer region around tiles and of appropriate Dirichlet boundary conditions around sCOLA boxes. As a result, we show that cosmological simulations at the degree of accuracy required for the analysis of the next generation of surveys can be run in drastically reduced wall-clock times and with very low memory requirements. The perfect scalability of our algorithm unlocks profoundly new possibilities for computing larger cosmological simulations at high resolution, taking advantage of a variety of hardware architectures.

arXiv comment24 pages, 8 figures, 2 tables. Typos corrected with respect to A&A published version. The code is publicly available at https://bitbucket.org/florent-leclercq/simbelmyne/

References in corpus (43)

• The NumPy array: a structure for efficient numerical computation
• Planck 2015 results. XIII. Cosmological parameters
• Advanced LIGO
• Baryonic Features in the Matter Transfer Function
• Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES
• Large-Scale Structure of the Universe and Cosmological Perturbation Theory
• The Simons Observatory: Science goals and forecasts
• The DESI Experiment Part I: Science,Targeting, and Survey Design
• Euclid Definition Study Report
• Power Spectra for Cold Dark Matter and its Variants
• Enzo: An Adaptive Mesh Refinement Code for Astrophysics
• Multi-scale initial conditions for cosmological simulations
• The effects of galaxy formation on the matter power spectrum: A challenge for precision cosmology
• eROSITA Science Book: Mapping the Structure of the Energetic Universe
• Solving Large Scale Structure in Ten Easy Steps with COLA
• Cosmic Dawn (CoDa): the First Radiation-Hydrodynamics Simulation of Reionization and Galaxy Formation in the Local Universe
• Quantifying baryon effects on the matter power spectrum and the weak lensing shear correlation
• The impact of baryons on the matter power spectrum from the Horizon-AGN cosmological hydrodynamical simulation
• Calibrating the Nonlinear Matter Power Spectrum: Requirements for Future Weak Lensing Surveys
• A simple multigrid scheme for solving the Poisson equation with arbitrary domain boundaries
• Massive optimal data compression and density estimation for scalable, likelihood-free inference in cosmology
• Neutrino masses and cosmological parameters from a Euclid-like survey: Markov Chain Monte Carlo forecasts including theoretical errors
• SPTpol: an instrument for CMB polarization measurements with the South Pole Telescope
• The Properties of Cosmic Velocity Fields
• A High-Fidelity Realization of the Euclid Code Comparison $N$-body Simulation with Abacus
• ICE-COLA: Towards fast and accurate synthetic galaxy catalogues optimizing a quasi $N$-body method
• Fast and accurate mock catalogue generation for low-mass galaxies
• Bayesian optimisation for likelihood-free cosmological inference
• Bayesian analysis of the dynamic cosmic web in the SDSS galaxy survey
• Particle-Mesh code for cosmological simulations
• EMMA: an AMR cosmological simulation code with radiative transfer
• Nonlinear evolution of density perturbations using approximate constancy of gravitational potential
• Can Standard Cosmological Models Explain the Observed Abell Cluster Bulk Flow?
• sCOLA: The N-body COLA Method Extended to the Spatial Domain
• The phase-space structure of nearby dark matter as constrained by the SDSS
• DEUS Full Observable ΛCDM Universe Simulation: the numerical challenge
• Estimating CDM Particle Trajectories in the Mildly Non-Linear Regime of Structure Formation. Implications for the Density Field in Real and Redshift Space
• CUBE: An Information-optimized parallel Cosmological $N$-body Algorithm
• Primordial power spectrum and cosmology from black-box galaxy surveys
• CUBE -- Towards an Optimal Scaling of Cosmological N-body Simulations
• SWIFT: Fast algorithms for multi-resolution SPH on multi-core architectures
• SWIFT: task-based hydrodynamics and gravity for cosmological simulations
• Bayesian large-scale structure inference and cosmic web analysis

Cited by in corpus (12)

• Large-scale dark matter simulations
• CARPool: fast, accurate computation of large-scale structure statistics by pairing costly and cheap cosmological simulations
• Bayesian forward modelling of cosmic shear data
• Measuring the Tidal Response of Structure Formation: Anisotropic Separate Universe Simulations using TreePM
• Nonlinear 3D Cosmic Web Simulation with Heavy-Tailed Generative Adversarial Networks
• Sesame: A power spectrum emulator pipeline for beyond-$\Lambda$CDM models
• Bayesian Control Variates for optimal covariance estimation with pairs of simulations and surrogates
• $\text{PineTree}$: A generative, fast, and differentiable halo model for wide-field galaxy surveys
• COmoving Computer Acceleration (COCA): $N$-body simulations in an emulated frame of reference
• Dynamic Zoom Simulations: a fast, adaptive algorithm for simulating lightcones
• Towards an optimal extraction of cosmological parameters from galaxy cluster surveys using convolutional neural networks
• The Manticore Project II: Bayesian digital twins of cosmic structure across the SDSS and BOSS volumes