MNRAS 000,1{14 (2020) Preprint 1 December 2020 Compiled using MNRAS LATEX style file v3.0 The BACCO Simulation Project: A baryonification emulator with Neural Networks Giovanni Aric`o,1? Raul E. Angulo,1;2y Sergio Contreras,1 Lurdes Ondaro-Mallea,1 Marcos Pellejero-Iba~nez,1 & Matteo Zennaro1 1Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal, 4, 20018, Donostia-San Sebasti´an,Guipuzkoa, Spain. 2IKERBASQUE, Basque Foundation for Science, 48013, Bilbao, Spain. Accepted XXX. Received YYY; in original form ZZZ ABSTRACT We present a neural-network emulator for baryonic effects in the non-linear matter power spectrum. We calibrate this emulator using more than 50; 000 measurements in a 15-dimensional parameters space, varying cosmology and baryonic physics. Baryonic physics is described through a baryonification algorithm, that has been shown to accurately capture the relevant effects on the power spectrum and bispectrum in state- of-the-art hydrodynamical simulations. Cosmological parameters are sampled using a cosmology-rescaling approach including massive neutrinos and dynamical dark energy. The specific quantity we emulate is the ratio between matter power spectrum with baryons and gravity-only, and we estimate the overall precision of the emulator to be 1 − 2%, at all scales 0:01 < k < 5h Mpc−1, and redshifts 0 < z < 1:5. We also obtain an accuracy of 1 − 2%, when testing the emulator against a collection of 74 different cosmological hydrodynamical simulations and their respective gravity-only counterparts. We show also that only one baryonic parameter, namely Mc, which set the gas fraction retained per halo mass, is enough to have accurate and realistic predictions of the baryonic feedback at a given epoch. Our emulator will become publicly available in http://www.dipc.org/bacco. Key words: large-scale structure { numerical methods { cosmological parameters 1 INTRODUCTION a Euclid-like survey (Laureijs et al. 2011), would result in a ≈ 5σ bias on cosmological parameters constraints (Sem- Gravity is the dominant force that shapes the structure of boloni et al. 2011, 2013; Schneider et al. 2020). the Universe on very large scales. On small scales, however, An alternative to deal with this uncertainty is to restrict hydrodynamical forces and astrophysical processes such as data analyses to scales large enough to be reasonably un- cooling, star formation, supernovae and AGN feedback, be- affected by baryons. However, this would imply discarding come important and cosmic structure is the result of their a huge amount of cosmological information, which could be joint co-evolution and interaction with gravity. Details of crucial to distinguish between, for instance, competing theo- such processes are not fully understood nor sufficiently con- ries for cosmic acceleration or the nature of the dark matter strained by current observations to formulate a predictive arXiv:2011.15018v1 [astro-ph.CO] 30 Nov 2020 particle. theory of structure formation on all scales. Arguably, the most complete and complex method to evolve In fact, our limited knowledge about the baryonic physics simultaneously gravity and baryon physics is currently given is one of the main uncertainties in modelling, and cosmo- by magneto-hydrodynamical simulations. These implement logical interpretation, of ongoing weak gravitational lens- a large number of different physical processes, and recently ing surveys, which map the distribution of density fluctu- the optimisations of the algorithms and the increasing power ations in the Universe. For upcoming surveys, sources of of super computers have opened-up to simulations of rel- statistical noise will be dramatically reduced and the rela- atively large cosmological boxes (Schaye et al. 2010; van tive importance of the baryonic uncertainties will increase. Daalen et al. 2011; Le Brun et al. 2014; Vogelsberger et al. For instance, ignoring the impact of baryonic processes in 2013; Dubois et al. 2014; Schaye et al. 2015; McCarthy et al. 2017; Springel et al. 2018). Unfortunately, these simulations are still too computationally expensive to be directly used ? E-mail:giovanni [email protected] in data analysis. Thus, over the last years, a number of y E-mail:[email protected] © 2020 The Authors 2 Aric`oet al. different approaches have been proposed for fast modelling Cosmology Ω Ω h ns of baryons. A non-complete list is: extensions to the halo cdm b model (e.g. Semboloni et al. 2011; Mohammed et al. 2014; Nenya 0.265 0.050 0.60 1.01 Fedeli 2014; Mead et al. 2015; Debackere et al. 2019; Mead Narya 0.310 0.050 0.70 1.01 et al. 2020); Principal Component Analyses (Eifler et al. Vilya 0.210 0.060 0.65 0.92 2015; Huang et al. 2019); machine learning-based meth- TheOne 0.259 0.048 0.68 0.96 ods (Tr¨osteret al. 2019; Villaescusa-Navarro et al. 2020); Effective Field Theory of Large-Scale Structure (Bragan¸ca Table 1. Cosmological parameters set used in the BACCO simu- lation project. All the cosmologies assume a flat geometry, mass- et al. 2020). Among these, a promising method is the so- less neutrinos (M = 0), a dark energy equation of state with called baryon correction model, or baryonification (Schnei- ν w0 = −1 and wa = 0, an amplitude of matter fluctuations der & Teyssier 2015; Schneider et al. 2019; Aric`oet al. σ8 = 0:9, and optical depth at recombination τ = 0:0952. 2020b). The main idea of baryonification is to modify the outputs of gravity-only simulations according to physically motivated recipes for the spatial distribution of baryons. As but does not explicitly incorporate the effect of cosmology a result, the three-dimensional total matter density field can on the properties of the underlying mass field. be predicted as a function of cosmology and baryonic physics Here, we exploit the framework described in Aric`oet al. (Aric`oet al. 2020a). Some main advantages of this approach (2020b,a) and in Angulo et al.(2020), to build an emulator are that its free parameters have a clear physical interpre- which encodes all the dependencies on 7 baryon parame- tation, it does not rely on the correctness of any partic- ters and also on 8 cosmological parameters, including mas- ular hydrodynamical simulation, it can be linked directly sive neutrinos and dynamical dark energy. We train a feed- with observations, and that it predicts the full density field forward Neural Network with tens of thousands of different and not only specific correlation functions or power spectra. combinations of cosmology and baryonification parameters The aforementioned advantages come at the price of a sig- to learn the changes on the matter power spectrum induced nificant computational cost in comparison to e.g. the halo by baryonic physics. model and other analytic approaches. The computational We will show that our emulation is accurate at a ∼ 1 − 2% cost of baryonification, in fact, roughly scales linearly with level over our whole 15-dimensional parameter space, which the number of haloes in a simulation and thus can be rather covers scales 0:01 ≤ k ≤ 5h Mpc−1 and redshifts 0 ≤ z ≤ expensive for large simulations. For instance, to apply the 1:5. We will then validate the performance of our emula- heavily optimised baryonification algorithm of Aric`oet al. tor against a library of 74 hydro-dynamical simulations and (2020b,a) on top of a 512 h−1Mpc-box simulation with 15363 their respective gravity-only counterparts. In particular, we particles takes approximately 10 CPU-minutes. Although make use of BAHAMAS, Cosmo-OWLS, OWLS, EAGLE, this is a negligible time for many applications, for inten- Illustris, Illustris TNG, and Horizon simulations. We will sive parameter-space sampling where hundreds of thousands show that we can accurately reproduce all of their power of evaluations might be required, the baryonification could spectra, and subsequently will explore the minimal baryonic make data analysis unfeasible. To overcome this problem, parameterisation that could be needed in future weak lens- in this paper we build and validate an emulator for the ef- ing data analyses. This paper is structured as follows: in §2 fects of baryonification on the density power spectrum. This we describe the numerical techniques we employ; in §3 we emulator is able to provide very accurate predictions in a present our baryonification emulator. In §4 we explore the fraction of second. An additional advantage of the emulator baryonic effects dependencies on cosmology, and in §5 we over the full algorithm is its extreme flexibility: the emulator discuss the model constraints obtained by fitting measure- does not need to directly handle the outputs of N-body sim- ment from hydrodynamical simulations. We recap our main ulations, and thus can be easily incorporated to any lensing findings and conclude in §6. analysis pipeline. Recently, emulators have been increasingly popular in cosmology, as they provide a mean to exploit the statistical power of large numerical simulations with negli- gible computational memory and time requirements. In the 2 NUMERICAL METHODS last years, emulators of different observables have been built, In this section we describe the main numerical methods we e.g. matter power spectrum (e.g. Heitmann et al. 2014; Kn- will employ throughout this paper. We start by describing abenhans et al. 2019; Winther et al. 2019; Angulo et al. our cosmological N-body simulations (§2.1). We continue by 2020; Euclid Collaboration et al. 2020), mass function (e.g. recapping the algorithms with which we will model bary- McClintock et al. 2019; Bocquet et al. 2020), and includ- onic effects on these simulations (§2.2) and modify their ing galaxy correlation function and Lyman−α Forest (Zhai cosmological parameters (§2.3). We then briefly discuss our et al. 2019; Bird et al. 2019). Alternatively, emulators of the power spectrum measurements (§2.4) and finish by perform- Likelihood function have been proposed for cosmology infer- ing tests to quantify the accuracy of our predictions (§2.5).