No Arabic abstract
We estimate the velocity field in a large set of $N$-body simulations including massive neutrino particles, and measure the auto-power spectrum of the velocity divergence field as well as the cross-power spectrum between the cold dark matter density and the velocity divergence. We perform these measurements at four different redshifts and within four different cosmological scenarios, covering a wide range in neutrino masses. We find that the nonlinear correction to the velocity power spectra largely depend on the degree of nonlinear evolution with no specific dependence on the value of neutrino mass. We provide a fitting formula, based on the value of the r.m.s. of the matter fluctuations in spheres of $8h^{-1}$Mpc, describing the nonlinear corrections with 3% accuracy on scales below $k=0.7; h$ Mpc$^{-1}$.
The set-up of the initial conditions in cosmological N-body simulations is usually implemented by rescaling the desired low-redshift linear power spectrum to the required starting redshift consistently with the Newtonian evolution of the simulation. The implementation of this practical solution requires more care in the context of massive neutrino cosmologies, mainly because of the non-trivial scale-dependence of the linear growth that characterises these models. In this work we consider a simple two-fluid, Newtonian approximation for cold dark matter and massive neutrinos perturbations that can reproduce the cold matter linear evolution predicted by Boltzmann codes such as CAMB or CLASS with a 0.1% accuracy or below for all redshift relevant to nonlinear structure formation. We use this description, in the first place, to quantify the systematic errors induced by several approximations often assumed in numerical simulations, including the typical set-up of the initial conditions for massive neutrino cosmologies adopted in previous works. We then take advantage of the flexibility of this approach to rescale the late-time linear power spectra to the simulation initial redshift, in order to be as consistent as possible with the dynamics of the N-body code and the approximations it assumes. We implement our method in a public code providing the initial displacements and velocities for cold dark matter and neutrino particles that will allow accurate, i.e. one-percent level, numerical simulations for this cosmological scenario.
Cosmic voids are a promising environment to characterize neutrino-induced effects on the large-scale distribution of matter in the universe. We perform a comprehensive numerical study of the statistical properties of voids, identified both in the matter and galaxy distributions, in massive and massless neutrino cosmologies. The matter density field is obtained by running several independent $N$-body simulations with cold dark matter and neutrino particles, while the galaxy catalogs are modeled by populating the dark matter halos in simulations via a halo occupation distribution (HOD) model to reproduce the clustering properties observed by the Sloan Digital Sky Survey (SDSS) II Data Release 7. We focus on the impact of massive neutrinos on the following void statistical properties: number density, ellipticities, two-point statistics, density and velocity profiles. Considering the matter density field, we find that voids in massive neutrino cosmologies are less evolved than those in the corresponding massless neutrinos case: there is a larger number of small voids and a smaller number of large ones, their profiles are less evacuated, and they present a lower wall at the edge. Moreover, the degeneracy between $sigma_8$ and $Omega_{ u}$ is broken when looking at void properties. In terms of the galaxy density field, we find that differences among cosmologies are difficult to detect because of the small number of galaxy voids in the simulations. Differences are instead present when looking at the matter density and velocity profiles around these voids.
We present a new method for fitting peculiar velocity models to complete flux limited magnitude-redshifts catalogues, using the luminosity function of the sources as a distance indicator.The method is characterised by its robustness. In particular, no assumptions are made concerning the spatial distribution of sources and their luminosity function. Moreover, selection effects in redshift are allowed. Furthermore the inclusion of additional observables correlated with the absolute magnitude -- such as for example rotation velocity information as described by the Tully-Fisher relation -- is straightforward. As an illustration of the method, the predicted IRAS peculiar velocity model characterised by the density parameter beta is tested on two samples. The application of our method to the Tully-Fisher MarkIII MAT sample leads to a value of beta=0.6 pm 0.125, fully consistent with the results obtained previously by the VELMOD and ITF methods on similar datasets. Unlike these methods however, we make a very conservative use of the Tully-Fisher information. Specifically, we require to assume neither the linearity of the Tully-Fisher relation nor a gaussian distribution of its residuals. Moreover, the robustness of the method implies that no Malmquist corrections are required. A second application is carried out, using the fluxes of the IRAS 1.2 Jy sample as the distance indicator. In this case the effective depth of the volume in which the velocity model is compared to the data is almost twice the effective depth of the MarkIII MAT sample. The results suggest that the predicted IRAS velocity model, while successful in reproducing locally the cosmic flow, fails to describe the kinematics on larger scales.
Halos and galaxies are tracers of the underlying dark matter structures. While their bias is well understood in the case of a simple Universe composed dominantly of dark matter, the relation becomes more complex in the presence of massive neutrinos. Indeed massive neutrinos introduce rich dynamics in the process of structure formation leading to scale-dependent bias. We study this process from the perspective of general relativity employing a simple spherical collapse model. We find a characteristic signature at the neutrino free-streaming scale in addition to a large-scale feature from general relativity. The scale-dependent halo bias opposes the suppression in the matter distribution due to neutrino free-streaming and leads to corrections of a few percent in the halo power spectrum. It is not only sensitive to the sum of the neutrino-masses, but respond to the individual masses. Accurate models for the neutrino bias are a crucial ingredient for the future data analysis and play an important role in constraining the neutrino masses.
We present a convolutional neural network to classify distinct cosmological scenarios based on the statistically similar weak-lensing maps they generate. Modified gravity (MG) models that include massive neutrinos can mimic the standard concordance model ($Lambda$CDM) in terms of Gaussian weak-lensing observables. An inability to distinguish viable models that are based on different physics potentially limits a deeper understanding of the fundamental nature of cosmic acceleration. For a fixed redshift of sources, we demonstrate that a machine learning network trained on simulated convergence maps can discriminate between such models better than conventional higher-order statistics. Results improve further when multiple source redshifts are combined. To accelerate training, we implement a novel data compression strategy that incorporates our prior knowledge of the morphology of typical convergence map features. Our method fully distinguishes $Lambda$CDM from its most similar MG model on noise-free data, and it correctly identifies among the MG models with at least 80% accuracy when using the full redshift information. Adding noise lowers the correct classification rate of all models, but the neural network still significantly outperforms the peak statistics used in a previous analysis.