No Arabic abstract
Interpreting observations of the Lyman-$alpha$ forest flux power spectrum requires interpolation between a small number of expensive simulations. We present a Gaussian process emulator modelling the 1D flux power spectrum as a function of the amplitude and slope of the small-scale linear matter power spectrum, and the state of the intergalactic medium at the epoch of interest ($2 < z < 4$). This parameterisation enables the prediction of the flux power spectrum in extended cosmological models that are not explicitly included in the training set, eliminating the need to construct bespoke emulators for a number of extensions to $Lambda$CDM. Our emulator is appropriate for cosmologies in which the linear matter power spectrum is described to percent level accuracy by just an amplitude and slope across the epoch of interest, and in the regime probed by eBOSS/DESI data. We demonstrate this for massive neutrino cosmologies, where the emulator is able to predict the flux power spectrum in a $Sigma m_ u=0.3$ eV neutrino cosmology to sub-percent accuracy, without including massive neutrinos in the training simulations. Further parameters would be required to describe models with sharp features in the linear power, such as warm or light axion dark matter. This work will facilitate the combination of upcoming DESI data with observations of the cosmic microwave background, to obtain constraints on neutrino mass and other extensions to $Lambda$CDM cosmology.
We present a method that extends the capabilities of the PINpointing Orbit-Crossing Collapsed HIerarchical Objects (PINOCCHIO) code, allowing it to generate accurate dark matter halo mock catalogues in cosmological models where the linear growth factor and the growth rate depend on scale. Such cosmologies comprise, among others, models with massive neutrinos and some classes of modified gravity theories. We validate the code by comparing the halo properties from PINOCCHIO against N-body simulations, focusing on cosmologies with massive neutrinos: $ uLambda$CDM. We analyse the halo mass function, halo two-point correlation function, halo power spectrum and the moments of the halo density field, showing that PINOCCHIO reproduces the results from simulations with the same level of precision as the original code ($sim5-10%$). We demonstrate that the abundance of halos in cosmologies with massless and massive neutrinos from PINOCCHIO matches very well the outcome of simulations, and point out that PINOCCHIO can reproduce the $Omega_ u-sigma_8$ degeneracy that affects the halo mass function. We show that the clustering properties of the halos from PINOCCHIO matches accurately those from simulations both in real and redshift-space, in the latter case up to $k=0.3~h~{rm Mpc}^{-1}$. We finally point out that the first moments of the halo density field from simulations are precisely reproduced by PINOCCHIO. We emphasize that the computational time required by PINOCCHIO to generate mock halo catalogues is orders of magnitude lower than the one needed for N-body simulations. This makes this tool ideal for applications like covariance matrix studies within the standard $Lambda$CDM model but also in cosmologies with massive neutrinos or some modified gravity theories.
Fuzzy Dark Matter (FDM) represents an alternative and intriguing description of the standard Cold Dark Matter (CDM) fluid, able to explain the lack of direct detection of dark matter particles in the GeV sector and to alleviate small scales tensions in the cosmic large-scale structure formation. Cosmological simulations of FDM models in the literature were performed either with very expensive high-resolution grid-based simulations of individual haloes or through N-body simulations encompassing larger cosmic volumes but resorting on significant approximations in the FDM non-linear dynamics to reduce their computational cost. With the use of the new N-body cosmological hydrodynamical code AX-GADGET, we are now able not only to overcome such numerical problems, but also to combine a fully consistent treatment of FDM dynamics with the presence of gas particles and baryonic physical processes, in order to quantify the FDM impact on specific astrophysical observables. In particular, in this paper we perform and analyse several hydrodynamical simulations in order to constrain the FDM mass by quantifying the impact of FDM on Lyman-$alpha$ forest observations, as obtained for the first time in the literature in a N-body setup without approximating the FDM dynamics. We also study the statistical properties of haloes, exploiting the large available sample, to extract information on how FDM affects the abundance, the shape, and density profiles of dark matter haloes.
We propose a new method for fitting the full-shape of the Lyman-$alpha$ (Ly$alpha$) forest three-dimensional (3D) correlation function in order to measure the Alcock-Paczynski (AP) effect. Our method preserves the robustness of baryon acoustic oscillations (BAO) analyses, while also providing extra cosmological information from a broader range of scales. We compute idealized forecasts for the Dark Energy Spectroscopic Instrument (DESI) using the Ly$alpha$ auto-correlation and its cross-correlation with quasars, and show how this type of analysis improves cosmological constraints. The DESI Ly$alpha$ BAO analysis is expected to measure $H(z_mathrm{eff})r_mathrm{d}$ and $D_mathrm{M}(z_mathrm{eff})/r_mathrm{d}$ with a precision of $sim0.9%$ each, where $H$ is the Hubble parameter, $r_mathrm{d}$ is the comoving BAO scale, $D_mathrm{M}$ is the comoving angular diameter distance and the effective redshift of the measurement is $z_mathrm{eff}simeq2.3$. By fitting the AP parameter from the full shape of the two correlations, we show that we can obtain a precision of $sim0.5-0.6%$ on each of $H(z_mathrm{eff})r_mathrm{d}$ and $D_mathrm{M}(z_mathrm{eff})/r_mathrm{d}$. Furthermore, we show that a joint full-shape analysis of the Ly$alpha$ auto-correlation and its cross-correlation with quasars can measure the linear growth rate times the amplitude of matter fluctuations in spheres of $8;h^{-1}$Mpc, $fsigma_8(z_mathrm{eff})$. Such an analysis could provide the first ever measurement of $fsigma_8(z_mathrm{eff})$ at redshift $z_mathrm{eff}>2$. By combining this with the quasar auto-correlation in a joint analysis of the three high-redshift two-point correlation functions, we show that DESI could be able to measure $fsigma_8(z_mathrm{eff}simeq2.3)$ with a precision of $5-12%$, depending on the smallest scale fitted.
Cosmological constraints are usually derived under the assumption of a $6$ parameters $Lambda$-CDM theoretical framework or simple one-parameter extensions. In this paper we present, for the first time, cosmological constraints in a significantly extended scenario, varying up to $12$ cosmological parameters simultaneously, including the sum of neutrino masses, the neutrino effective number, the dark energy equation of state, the gravitational waves background and the running of the spectral index of primordial perturbations. Using the latest Planck 2015 data release (with polarization) we found no significant indication for extensions to the standard $Lambda$-CDM scenario, with the notable exception of the angular power spectrum lensing amplitude, $A_{rm lens}$ that is larger than the expected value at more than two standard deviations even when combining the Planck data with BAO and supernovae type Ia external datasets. In our extended cosmological framework, we find that a combined Planck+BAO analysis constrains the value of the r.m.s. density fluctuation parameter to $sigma_8=0.781_{-0.063}^{+0.065}$ at $95 %$ c.l., helping to relieve the possible tensions with the CFHTlenS cosmic shear survey. We also find a lower value for the reionization optical depth $tau=0.058_{-0.043}^{+0.040}$ at $95$ % c.l. respect to the one derived under the assumption of $Lambda$-CDM. The scalar spectral index $n_S$ is now compatible with a Harrison-Zeldovich spectrum to within $2.5$ standard deviations. Combining the Planck dataset with the HST prior on the Hubble constant provides a value for the equation of state $w < -1$ at more than two standard deviations while the neutrino effective number is fully compatible with the expectations of the standard three neutrino framework.
Cosmological hydrodynamic simulations can accurately predict the properties of the intergalactic medium (IGM), but only under the condition of retaining high spatial resolution necessary to resolve density fluctuations in the IGM. This resolution constraint prohibits simulating large volumes, such as those probed by BOSS and future surveys, like DESI and 4MOST. To overcome this limitation, we present Iteratively Matched Statistics (IMS), a novel method to accurately model the Lyman-alpha forest with collisionless N-body simulations, where the relevant density fluctuations are unresolved. We use a small-box, high-resolution hydrodynamic simulation to obtain the probability distribution function (PDF) and the power spectrum of the real-space Lyman-alpha forest flux. These two statistics are iteratively mapped onto a pseudo-flux field of an N-body simulation, which we construct from the matter density. We demonstrate that our method can perfectly reproduce line-of-sight observables, such as the PDF and power spectrum, and accurately reproduce the 3D flux power spectrum (5-20%). We quantify the performance of the commonly used Gaussian smoothing technique and show that it has significantly lower accuracy (20-80%), especially for N-body simulations with achievable mean inter-particle separations in large-volume simulations. In addition, we show that IMS produces reasonable and smooth spectra, making it a powerful tool for modeling the IGM in large cosmological volumes and for producing realistic mock skies for Lyman-alpha forest surveys.