Fast generation of mock galaxy catalogues with COLA

Abstract

We investigate the feasibility of using COmoving Lagrangian Acceleration (COLA) technique to efficiently generate galaxy mock catalogues that can accurately reproduce the statistical properties of observed galaxies. Our proposed scheme combines the subhalo abundance matching (SHAM) procedure with COLA simulations, utilizing only three free parameters: the scatter magnitude (${\sigma }_{\mathrm{scat}}$) in SHAM, the initial redshift ($z_{\mathrm{init}}$) of the COLA simulation, and the time stride ($da$) used by COLA. In this proof-of-concept study, we focus on a subset of BOSS CMASS NGC galaxies within the redshift range $z\in [0.45,0.55]$. We perform $\mathtt{GADGET}$ simulation and low-resolution COLA simulations with various combinations of $(z_{\mathrm{init}},da)$, each using $1024^3$ particles in an $800h^{-1}\mathrm{Mpc}$ box. By minimizing the difference between COLA mock and CMASS NGC galaxies for the monopole of the two-point correlation function (2PCF), we obtain the optimal ${\sigma }_{\mathrm{scat}}$. We have found that by setting $z_{\mathrm{init}}=29$ and $da=1/30$, we achieve a good agreement between COLA mock and CMASS NGC galaxies within the range of 4 to $20h^{-1}\mathrm{Mpc}$, with a computational cost two orders of magnitude lower than that of the N-body code. Moreover, a detailed verification is performed by comparing various statistical properties, such as anisotropic 2PCF, three-point clustering, and power spectrum multipoles, which shows similar performance between GADGET mock and COLA mock catalogues with the CMASS NGC galaxies. Furthermore, we assess the robustness of the COLA mock catalogues across different cosmological models, demonstrating consistent results in the resulting 2PCFs. Our findings suggest that COLA simulations are a promising tool for efficiently generating mock catalogues for emulators and machine learning analyses in exploring the large-scale structure of the Universe.

24 pages, 14 figures, 4 tables

References in corpus (15)

• Array Programming with NumPy
• The NumPy array: a structure for efficient numerical computation
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation
• The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
• Modeling the color evolution of luminous red galaxies - improvements with empirical stellar spectra
• Semi-analytic forecasts for JWST -- VI. Simulated lightcones and galaxy clustering predictions
• Testing gravity on large scales by combining weak lensing with galaxy clustering using CFHTLenS and BOSS CMASS
• Cosmological constraints from the redshift dependence of the Alcock-Paczynski effect: application to the SDSS-III BOSS DR12 galaxies
• A Hybrid Deep Learning Approach to Cosmological Constraints From Galaxy Redshift Surveys
• ICE-COLA: fast simulations for weak lensing observables
• Semi-analytic forecasts for Roman -- the beginning of a new era of deep-wide galaxy surveys
• Graph Database Solution for Higher Order Spatial Statistics in the Era of Big Data
• Cosmological constraints from the redshift dependence of the Alcock-Paczynski test and volume effect: galaxy two-point correlation function
• Extracting cosmological parameters from N-body simulations using machine learning techniques
• Cosmic Velocity Field Reconstruction Using AI

Cited by in corpus (5)

• $\text{PineTree}$: A generative, fast, and differentiable halo model for wide-field galaxy surveys
• Improving SDSS cosmological constraints through $\beta$-skeleton weighted correlation functions
• 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
• Enhancing Cosmological Constraints by Two-dimensional $\beta$-cosmic-web Weighted Angular Correlation Functions