The $\text{Abacus}$ Cosmological $N$-body Code

Abstract

We present $\text{Abacus}$, a fast and accurate cosmological $N$-body code based on a new method for calculating the gravitational potential from a static multipole mesh. The method analytically separates the near- and far-field forces, reducing the former to direct $1/r^2$ summation and the latter to a discrete convolution over multipoles. The method achieves 70 million particle updates per second per node of the Summit supercomputer, while maintaining a median fractional force error of $10^{-5}$. We express the simulation time step as an event-driven "pipeline", incorporating asynchronous events such as completion of co-processor work, Input/Output, and network communication. $\text{Abacus}$ has been used to produce the largest suite of $N$-body simulations to date, the $\text{AbacusSummit}$ suite of 60 trillion particles (Maksimova et al., 2021), incorporating on-the-fly halo finding. $\text{Abacus}$ enables the production of mock catalogs of the volume and resolution required by the coming generation of cosmological surveys.

arXiv comment29 pages, 19 figures. Published in MNRAS

References in corpus (40)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• The NumPy array: a structure for efficient numerical computation
• The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
• The cosmological simulation code GADGET-2
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Athena: A New Code for Astrophysical MHD
• Solving Large Scale Structure in Ten Easy Steps with COLA
• The Uchuu Simulations: Data Release 1 and Dark Matter Halo Concentrations
• nbodykit: an open-source, massively parallel toolkit for large-scale structure
• AbacusSummit: A Massive Set of High-Accuracy, High-Resolution $N$-Body Simulations
• FastPM: a new scheme for fast simulations of dark matter and halos
• Corrfunc --- A Suite of Blazing Fast Correlation Functions on the CPU
• A new way of setting the phases for cosmological multi-scale Gaussian initial conditions
• Limitations to the "basic'' HOD model and beyond
• The Abacus Cosmos: A Suite of Cosmological N-body Simulations
• The Outer Rim Simulation: A Path to Many-Core Supercomputers
• A High-Fidelity Realization of the Euclid Code Comparison $N$-body Simulation with Abacus
• The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: exploring the Halo Occupation Distribution model for Emission Line Galaxies
• The host dark matter halos of [OII] emitters at 0.5< z< 1.5
• Dark matter statistics for large galaxy catalogs: power spectra and covariance matrices
• Improving Initial Conditions for Cosmological $N$-Body Simulations
• Evidence for galaxy assembly bias in BOSS CMASS redshift-space galaxy correlation function
• Emulating galaxy clustering and galaxy-galaxy lensing into the deeply nonlinear regime: methodology, information, and forecasts
• Linear perturbative theory of the discrete cosmological N-body problem
• An iterative reconstruction of cosmological initial density fields
• Cosmology with galaxy-galaxy lensing on non-perturbative scales: Emulation method and application to BOSS LOWZ
• Exploring the squeezed three-point galaxy correlation function with generalized halo occupation distribution models
• A Hybrid Deep Learning Approach to Cosmological Constraints From Galaxy Redshift Surveys
• Corrfunc: Blazing fast correlation functions with AVX512F SIMD Intrinsics
• Cosmology with Stacked Cluster Weak Lensing and Cluster-Galaxy Cross-Correlations
• The Effective Halo Model: Creating a Physical and Accurate Model of the Matter Power Spectrum and Cluster Counts
• Can Assembly Bias Explain the Lensing Amplitude of the BOSS CMASS Sample in a Planck Cosmology?
• CUBE: An Information-optimized parallel Cosmological $N$-body Algorithm
• Quantifying resolution in cosmological N-body simulations using self-similarity
• Covariance matrices for galaxy cluster weak lensing: from virial regime to uncorrelated large-scale structure
• Good and Proper: Self-similarity of N-body Simulations with Proper Force Softening
• Tests of Acoustic Scale Shifts in Halo-based Mock Galaxy Catalogues
• Generating Approximate Halo Catalogs for Blind Challenges in Precision Cosmology
• Decorrelating the errors of the galaxy correlation function with compact transformation matrices
• Checkpointing with cp: the POSIX Shared Memory System

Cited by in corpus (70)

• SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
• AbacusHOD: A highly efficient extended multi-tracer HOD framework and its application to BOSS and eBOSS data
• Consistent lensing and clustering in a low-$S_8$ Universe with BOSS, DES Year 3, HSC Year 1 and KiDS-1000
• \textsc{CompaSO}: A new halo finder for competitive assignment to spherical overdensities
• DESI 2024 II: Sample Definitions, Characteristics, and Two-point Clustering Statistics
• Stringent ${\sigma }_8$ constraints from small-scale galaxy clustering using a hybrid MCMC+emulator framework
• Going Beyond the Galaxy Power Spectrum: an Analysis of BOSS Data with Wavelet Scattering Transforms
• DESI luminous red galaxy samples for cross-correlations
• Euclid. V. The Flagship galaxy mock catalogue: a comprehensive simulation for the Euclid mission
• Aemulus $\nu$: Precise Predictions for Matter and Biased Tracer Power Spectra in the Presence of Neutrinos
• Precise Cosmological Constraints from BOSS Galaxy Clustering with a Simulation-Based Emulator of the Wavelet Scattering Transform
• Galaxy Clustering Analysis with SimBIG and the Wavelet Scattering Transform
• $\mathrm{S}\mathrm{IM}\mathrm{BIG}$: Mock Challenge for a Forward Modeling Approach to Galaxy Clustering
• The DESI $N$-body Simulation Project I: Testing the Robustness of Simulations for the DESI Dark Time Survey
• Reconstructing Cosmological Initial Conditions from Late-Time Structure with Convolutional Neural Networks
• Mitigation of DESI fiber assignment incompleteness effect on two-point clustering with small angular scale truncated estimators
• Semi-analytical covariance matrices for two-point correlation function for DESI 2024 data
• Precision Redshift-Space Galaxy Power Spectra using Zel'dovich Control Variates
• The DESI $N$-body Simulation Project -- II. Suppressing sample variance with fast simulations
• The halo light cone catalogues of \Abacus{AbacusSummit}
• CSST Cosmological Emulator I: Matter Power Spectrum Emulation with one percent accuracy to k = 10 h/Mpc
• Constructing high-fidelity halo merger trees in AbacusSummit
• A comparison between Shapefit compression and Full-Modelling method with PyBird for DESI 2024 and beyond
• HOD-Dependent Systematics in Emission Line Galaxies for the DESI 2024 BAO analysis
• Differentiable Cosmological Simulation with Adjoint Method
• Generating mock galaxy catalogues for flux-limited samples like the DESI Bright Galaxy Survey
• Mitigating the noise of DESI mocks using analytic control variates
• Velocity reconstruction in the era of DESI and Rubin (part I): Exploring spectroscopic, photometric & hybrid samples
• HOD-Dependent Systematics for Luminous Red Galaxies in the DESI 2024 BAO Analysis
• Full forward model of galaxy clustering statistics with AbacusSummit lightcones
• Detection of spatial clustering in the 1000 richest SDSS DR8 redMaPPer clusters with Nearest Neighbor distributions
• 2D k-th nearest neighbor statistics: a highly informative probe of galaxy clustering
• Planting a Lyman alpha forest on AbacusSummit
• Accuracy of power spectra in dissipationless cosmological simulations
• Cosmology from Planck CMB Lensing and DESI DR1 Quasar Tomography
• Model-Independent Measurement of the Matter-Radiation Equality Scale in DESI 2024
• Constraining galaxy-halo connection with high-order statistics
• Generalized framework for likelihood-based field-level inference of growth rate from velocity and density fields
• Magrathea-Pathfinder: A 3D adaptive-mesh code for geodesic ray tracing in $N$-body simulations
• A Comparison of Void-Finding Algorithms using Crossing Numbers
• DESIVAST: A Catalog of Low-Redshift Voids using Data from the DESI DR1 Bright Galaxy Survey
• Self-Similarity of $k$-Nearest Neighbor Distributions in Scale-Free Simulations
• Window convolution of the galaxy clustering bispectrum
• Unraveling emission line galaxy conformity at z~1 with DESI early data
• Perturbation-theory informed integrators for cosmological simulations
• An Analytically Tractable Marked Power Spectrum
• Exploring HOD-dependent systematics for the DESI 2024 Full-Shape galaxy clustering analysis
• A fully Lagrangian, non-parametric bias model for dark-matter halos
• Constraining $f_{NL}$ using the Large-Scale Modulation of Small-Scale Statistics
• Cosmology with voids from the Nancy Grace Roman Space Telescope
• Convergence of halo statistics: code comparison between Rockstar and CompaSO using scale-free simulations
• Field-Level Comparison and Robustness Analysis of Cosmological N-body Simulations
• ELG$\times$LRG distribution through dark matter halo dynamics
• Bounds on galaxy stochasticity from halo occupation distribution modeling
• PySCo: A fast Particle-Mesh $N$-body code for modified gravity simulations in Python
• Cosmological parameters from the joint analysis of Density Split and Second Order Statistics: an Emulator based on the Halo Occupation Distribution
• Constraining primordial non-Gaussianity from the large scale structure two-point and three-point correlation functions
• Studying baryon acoustic oscillations using photometric redshifts from the DESI Legacy Imaging survey DR9
• Non-parametric Lagrangian biasing from the insights of neural nets
• A theoretical approach to density-split clustering
• Constraining Galaxy-Halo Connection Using Machine Learning
• Simulation budgeting for hybrid effective field theories
• Void spin distribution as a powerful probe of ${\sigma }_8$
• The impact of the redshift-dependent selection effect of halos on the redshift-space power spectrum
• Cosmological Analysis with Calibrated Neural Quantile Estimation and Approximate Simulators
• Bridging Simulations and EFT: A Hybrid Model of the Lyman-Alpha Forest Field
• Constraining cosmological parameters using density split lensing and the conditional stellar mass function
• Emulating galaxy and peculiar velocity clustering on non-linear scales
• Growth rate measurements from a joint analysis of the large-scale galaxy clustering in Fourier and configuration space
• Clustering analysis of medium-band selected high-redshift galaxies