51�������

Abstract:

To maximize the amount of information extracted from cosmological datasets, simulations that accurately represent these observations are necessary. However, traditional simulations that evolve particles under gravity by estimating particle-particle interactions (𝑁-body simulations) are computationally expensive and prohibitive to scale to the large volumes and resolutions necessary for the upcoming datasets. Moreover, modeling the distribution of galaxies typically involves identifying virialized dark matter halos, which is also a time- and memory-consuming process for large 𝑁-body simulations, further exacerbating the computational cost. In this study, we introduce CHARM, a novel method for creating mock halo catalogs by matching the spatial, mass, and velocity statistics of halos directly from the large-scale distribution of the dark matter density field. We develop multistage neural spline flow-based networks to learn this mapping at redshift 𝑧 =0.5 directly with computationally cheaper low-resolution particle mesh simulations instead of relying on the high-resolution 𝑁-body simulations. We show that the mock halo catalogs and painted galaxy catalogs have the same statistical properties as obtained from 𝑁-body simulations in both real space and redshift space. Finally, we use these mock catalogs for cosmological inference using redshift-space galaxy power spectrum, bispectrum, and wavelet-based statistics using simulation-based inference, performing the first inference with accelerated forward model simulations and finding unbiased cosmological constraints with well-calibrated posteriors.

Abstract:

Context. Baryonic physics has a considerable impact on the distribution of matter in our Universe on scales probed by current and future cosmological surveys, acting as a key systematic in such analyses. Aims. We seek simple symbolic parametrisations for the impact of baryonic physics on the matter power spectrum for a range of physically motivated models, as a function of wavenumber, redshift, cosmology, and parameters controlling the baryonic feedback. Methods. We used symbolic regression to construct analytic approximations for the ratio of the matter power spectrum in the presence of baryons to that without such effects. We obtained separate functions of each of four distinct sub-grid prescriptions of baryonic physics from the CAMELS suite of hydrodynamical simulations (Astrid, IllustrisTNG, SIMBA, and Swift-EAGLE) as well as for a baryonification algorithm. We also provide functions that describe the uncertainty on these predictions, due to both the stochastic nature of baryonic physics and the errors on our fits. Results. The error on our approximations to the hydrodynamical simulations is comparable to the sample variance estimated through varying initial conditions, and our baryonification expression has a root mean squared error of better than one percent, although this increases on small scales. These errors are comparable to those of previous numerical emulators for these models. Our expressions are enforced to have the physically correct behaviour on large scales and at high redshift. Due to their analytic form, we are able to directly interpret the impact of varying cosmology and feedback parameters, and we can identify parameters that have little to no effect. Conlcusions. Each function is based on a different implementation of baryonic physics, and can therefore be used to discriminate between these models when applied to real data. We provide a publicly available code for all symbolic approximations found.

Abstract:

<jats:p><jats:italic>Context.</jats:italic> Current and future large-scale structure surveys aim to constrain the neutrino mass and the equation of state of dark energy. To do this efficiently, rapid yet accurate evaluation of the matter power spectrum in the presence of these effects is essential.</jats:p> <jats:p><jats:italic>Aims.</jats:italic> We aim to construct accurate and interpretable symbolic approximations of the linear and nonlinear matter power spectra as a function of cosmological parameters in extended ΛCDM models that contain massive neutrinos and nonconstant equations of state for dark energy. This constitutes an extension of the S<jats:sc>YREN-HALOFIT</jats:sc> emulators to incorporate these two effects, which we call S<jats:sc>YREN-NEW</jats:sc> (SYmbolic-Regression-ENhanced power spectrum emulator with NEutrinos and <jats:italic>W</jats:italic><jats:sub>0</jats:sub>−<jats:italic>w</jats:italic><jats:sub><jats:italic>a</jats:italic></jats:sub>). We also wish to obtain a simple approximation of the derived parameter, <jats:italic>σ</jats:italic><jats:sub>8</jats:sub>, as a function of the cosmological parameters for these models.</jats:p> <jats:p><jats:italic>Methods.</jats:italic> We utilizedd symbolic regression to efficiently search through candidate analytic expressions to approximate the various quantities of interest. Our results for the linear power spectrum are designed to emulate C<jats:sc>LASS</jats:sc>, whereas for the nonlinear case we aim to match the results of E<jats:sc>UCLIDEMULATOR</jats:sc>2. We compared our results to existing emulators and <jats:italic>N</jats:italic>-body simulations.</jats:p> <jats:p><jats:italic>Results.</jats:italic> Our analytic emulators for <jats:italic>σ</jats:italic><jats:sub>8</jats:sub>, and the linear and nonlinear power spectra achieve root mean squared errors of 0.1%, 0.3%, and 1.3%, respectively, across a wide range of cosmological parameters, redshifts and wavenumbers. The error on the nonlinear power spectrum is reduced by approximately a factor of 2 when considering observationally plausible dark energy models and neutrino masses. We verify that emulator-related discrepancies are subdominant compared to observational errors and other modeling uncertainties when computing shear power spectra for LSST-like surveys. Our expressions have similar accuracy to existing (numerical) emulators, but are at least an order of magnitude faster, both on a CPU and a GPU.</jats:p> <jats:p><jats:italic>Conclusions.</jats:italic> Our work greatly improves the accuracy, speed, and applicability range of current symbolic approximations of the linear and nonlinear matter power spectra. These now cover the same range of cosmological models as many numerical emulators with similar accuracy, but are much faster and more interpretable. We provide publicly available code for all symbolic approximations found.</jats:p>

Abstract:

<jats:p><jats:italic>Context.N</jats:italic>-body simulations are computationally expensive and machine learning (ML) based emulation techniques have thus emerged as a way to increase their speed. Surrogate models are indeed fast, however, they are limited in terms of their trustworthiness due to potentially substantial emulation errors that current approaches are not equipped to correct.</jats:p> <jats:p><jats:italic>Aims.</jats:italic> To alleviate this problem, we have introduced COmoving Computer Acceleration (COCA), a hybrid framework interfacing ML algorithm with an <jats:italic>N</jats:italic>-body simulator. The correct physical equations of motion are solved in an emulated frame of reference, so that any emulation error is corrected by design. Thus, we are able to find a solution for the perturbation of particle trajectories around the ML solution. This approach is computationally cheaper than obtaining the full solution and it is guaranteed to converge to the truth as the number of force evaluations is increased.</jats:p> <jats:p><jats:italic>Methods.</jats:italic> Even though it is applicable to any ML algorithm and <jats:italic>N</jats:italic>-body simulator, we assessed this approach in the particular case of particle-mesh (PM) cosmological simulations in a frame of reference predicted by a convolutional neural network. In such cases, the time dependence is encoded as an additional input parameter to the network.</jats:p> <jats:p><jats:italic>Results.</jats:italic> We find that COCA efficiently reduces emulation errors in particle trajectories, requiring far fewer force evaluations than running the corresponding simulation without ML. As a consequence, we were able to obtain accurate final density and velocity fields for a reduced computational budget. We demonstrate that this method exhibits robustness when applied to examples outside the range of the training data. When compared to the direct emulation of the Lagrangian displacement field using the same training resources, COCA’s ability to correct emulation errors results in more accurate predictions.</jats:p> <jats:p><jats:italic>Conclusions.</jats:italic> Therefore, COCA makes <jats:italic>N</jats:italic>-body simulations cheaper by skipping unnecessary force evaluations, while still solving the correct equations of motion and correcting for emulation errors made by ML.</jats:p>