No Arabic abstract
We present an emulator for the two-point clustering of biased tracers in real space. We construct this emulator using neural networks calibrated with more than $400$ cosmological models in a 8-dimensional cosmological parameter space that includes massive neutrinos an dynamical dark energy. The properties of biased tracers are described via a Lagrangian perturbative bias expansion which is advected to Eulerian space using the displacement field of numerical simulations. The cosmology-dependence is captured thanks to a cosmology-rescaling algorithm. We show that our emulator is capable of describing the power spectrum of galaxy formation simulations for a sample mimicking that of a typical Emission-Line survey at $z sim 1$ with an accuracy of $1-2%$ up to nonlinear scales $k sim 0.7 h mathrm{Mpc}^{-1}$.
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 < 5 h/Mpc, 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.
We present the one-loop perturbation theory for the power spectrum of the marked density field of matter and biased tracers in real- and redshift-space. The statistic has been shown to yield impressive constraints on cosmological parameters; to exploit this, we require an accurate and computationally inexpensive theoretical model. Comparison with $N$-body simulations demonstrates that linear theory fails on all scales, but inclusion of one-loop Effective Field Theory terms gives a substantial improvement, with $sim 5%$ accuracy at $z = 1$. The expansion is less convergent in redshift-space (achieving $sim 10%$ accuracy), but there are significant improvements for biased tracers due to the freedom in the bias coefficients. The large-scale theory contains non-negligible contributions from all perturbative orders; we suggest a reorganization of the theory that contains all terms relevant on large-scales, discussing both its explicit form at one-loop and structure at infinite-loop. This motivates a low-$k$ correction term, leading to a model that is sub-percent accurate on large scales, albeit with the inclusion of two (three) free coefficients in real- (redshift-)space. We further consider the effects of massive neutrinos, showing that beyond-EdS corrections to the perturbative kernels are negligible in practice. It remains to see whether the purported gains in cosmological parameters remain valid for biased tracers and can be captured by the theoretical model.
We present the BACCO project, a simulation framework specially designed to provide highly-accurate predictions for the distribution of mass, galaxies, and gas as a function of cosmological parameters. In this paper, we describe our main suite of simulations (L $sim2$ Gpc and $4320^3$ particles) and present various validation tests. Using a cosmology-rescaling technique, we predict the nonlinear mass power spectrum over the redshift range $0<z<1.5$ and over scales $10^{-2} < k/(h Mpc^{-1} ) < 5$ for 800 points in an 8-dimensional cosmological parameter space. For an efficient interpolation of the results, we build an emulator and compare its predictions against several widely-used methods. Over the whole range of scales considered, we expect our predictions to be accurate at the 2% level for parameters in the minimal $Lambda$ CDM model and to 3% when extended to dynamical dark energy and massive neutrinos. We make our emulator publicly available under http://www.dipc.org/bacco
Upcoming galaxy redshift surveys promise to significantly improve current limits on primordial non-Gaussianity (PNG) through measurements of 2- and 3-point correlation functions in Fourier space. However, realizing the full potential of this dataset is contingent upon having both accurate theoretical models and optimized analysis methods. Focusing on the local model of PNG, parameterized by $f_{rm NL}$, we perform a Monte-Carlo Markov Chain analysis to confront perturbation theory predictions of the halo power spectrum and bispectrum in real space against a suite of N-body simulations. We model the halo bispectrum at tree-level, including all contributions linear and quadratic in $f_{rm NL}$, and the halo power spectrum at 1-loop, including tree-level terms up to quadratic order in $f_{rm NL}$ and all loops induced by local PNG linear in $f_{rm NL}$. Keeping the cosmological parameters fixed, we examine the effect of informative priors on the linear non-Gaussian bias parameter on the statistical inference of $f_{rm NL}$. A conservative analysisof the combined power spectrum and bispectrum, in which only loose priors are imposed and all parameters are marginalized over, can improve the constraint on $f_{rm NL}$ by more than a factor of 5 relative to the power spectrum-only measurement. Imposing a strong prior on $b_phi$, or assuming bias relations for both $b_phi$ and $b_{phidelta}$ (motivated by a universal mass function assumption), improves the constraints further by a factor of few. In this case, however, we find a significant systematic shift in the inferred value of $f_{rm NL}$ if the same range of wavenumber is used. Likewise, a Poisson noise assumption can lead to significant systematics, and it is thus essential to leave all the stochastic amplitudes free.
We study the effect of the Einstein - de Sitter (EdS) approximation on the one-loop power spectrum of galaxies in redshift space in the Effective Field Theory of Large-Scale Structure. The dark matter density perturbations and velocity divergence are treated with exact time dependence. Splitting the density perturbation into its different temporal evolutions naturally gives rise to an irreducible basis of biases. While, as in the EdS approximation, at each time this basis spans a seven-dimensional space, this space is a slightly different one, and the difference is captured by a single calculable time- and $vec k$-dependent function. We then compute the redshift-space galaxy one-loop power spectrum with the EdS approximation ($P^{text{EdS-approx}}$) and without ($P^{text{Exact}}$). For the monopole we find $P_{text{0}}^{text{Exact}}/P_{text{0}}^{text{EdS-approx}}sim 1.003$ and for the quadrupole $P_{text{2}}^{text{Exact}}/P_{text{2}}^{text{EdS-approx}}sim 1.007$ at $z=0.57$, and sharply increasing at lower redshifts. Finally, we show that a substantial fraction of the effect remains even after allowing the bias coefficients to shift within a physically allowed range. This suggests that the EdS approximation can only fit the data to a level of precision that is roughly comparable to the precision of the next generation of cosmological surveys. Furthermore, we find that implementing the exact time dependence formalism is not demanding and is easily applicable to data. Both of these points motivate a direct study of this effect on the cosmological parameters.