Abstract: Up-to-date cosmological data analyses have shown that (sigma) a closed universe is preferred by the Planck data at more than 99% CL, and (b) interacting scenarios offer a very compelling solution to the Hubble constant tension. In light of these two recent appealing scenarios, we consider here an interacting dark matter-dark energy model with a non-zero spatial curvature component and a freely varying dark energy equation of state in both the quintessential and phantom regimes. When considering Cosmic Microwave Background data only, a phantom and closed universe can perfectly alleviate the Hubble tension, without the necessity of a coupling among the dark sectors. Accounting for other possible cosmological observations compromises the viability of this very attractive scenario as a global solution to current cosmological tensions, either by spoiling its effectiveness concerning the H-0 problem, as in the case of Supernovae Ia data, or by introducing a strong disagreement in the preferred value of the spatial curvature, as in the case of Baryon Acoustic Oscillations.

Abstract: We address the issue of setting up generic non-Gaussian initial conditions for N-body simulations. We consider inflationary-motivated primordial non-Gaussianity where the perturbations in the Bardeen potential are given by a dominant Gaussian part plus a non-Gaussian part specified by its bispectrum. The approach we explore here is suitable for any bispectrum, i.e. it does not have to be of the so-called separable or factorizable form. The procedure of generating a non-Gaussian field with a given bispectrum (and a given power spectrum for the Gaussian component) is not univocal, and care must be taken so that higher-order corrections do not leave a too large signature on the power spectrum. This is so far a limiting factor of our approach. We then run N-body simulations for the most popular inflationary-motivated non-Gaussian shapes. The halo mass function and the non-linear power spectrum agree with theoretical analytical approximations proposed in the literature, even if they were so far developed and tested only for a particular shape (the local one). We plan to make the simulations outputs available to the community via the non-Gaussian simulations comparison project web site http://icc.ub.edu/similar to liciaverde/NGSCP.html.

Abstract: We employ state-of-the art cosmological observables including supernova surveys and BAO information to provide constraints on the mass and mixing angle of a non-resonantly produced sterile neutrino species, showing that cosmology can effectively rule out sterile neutrinos which decay between BBN and the present day. The decoupling of an additional heavy neutrino species can modify the time dependence of the Universe's expansion between BBN and recombination and, in extreme cases, lead to an additional matter-dominated period; while this could naively lead to a younger Universe with a larger Hubble parameter, it could later be compensated by the extra radiation expected in the form of neutrinos from sterile decay. However, recombination-era observables including the Cosmic Microwave Background (CMB), the shift parameter R-CMB and the sound horizon r(s) from Baryon Acoustic Oscillations (BAO) severely constrain this scenario. We self-consistently include the full time-evolution of the coupled sterile neutrino and standard model sectors in an MCMC, showing that if decay occurs after BBN, the sterile neutrino is essentially bounded by the constraint sin(2) theta less than or similar to 0.026(m(s)/eV)(-2).

Abstract: We study whether non-linear gravitational effects of relic neutrinos on the development of clustering and large-scale structure may be observable by weak gravitational lensing. We compute the density profile of relic massive neutrinos in a spherical model of a cluster of galaxies, for several neutrino mass schemes and cluster masses. Relic neutrinos add a small perturbation to the mass profile, making it more extended in the outer parts. In principle, this non-linear neutrino perturbation is detectable in an all-sky weak lensing survey such as EUCLID by averaging the shear profile of a large fraction of the visible massive clusters in the universe, or from its signature in the general weak lensing power spectrum or its cross-spectrum with galaxies. However, correctly modeling the distribution of mass in baryons and cold dark matter and suppressing any systematic errors to the accuracy required for detecting this neutrino perturbation is severely challenging.

Abstract: We investigate the non-linear evolution of the relic cosmic neutrino background by running large box-size, high resolution N-body simulations which incorporate cold dark matter (CDM) and neutrinos as independent particle species. Our set of simulations explore the properties of neutrinos in a reference Lambda CDM model with total neutrino masses between 0.05-0.60 eV in cold dark matter haloes of mass 10(11) – 10(15) h(-1) M-circle dot, over a redshift range z = 0 – 2. We compute the halo mass function and show that it is reasonably well fitted by the Sheth-Tormen formula, once the neutrino contribution to the total matter is removed. More importantly, we focus on the CDM and neutrino properties of the density and peculiar velocity fields in the cosmological volume, inside and in the outskirts of virialized haloes. The dynamical state of the neutrino particles depends strongly on their momentum: whereas neutrinos in the low velocity tail behave similarly to CDM particles, neutrinos in the high velocity tail are not affected by the clustering of the underlying CDM component. We find that the neutrino (linear) unperturbed momentum distribution is modified and mass and redshift dependent deviations from the expected Fermi-Dirac distribution are in place both in the cosmological volume and inside haloes. The neutrino density profiles around virialized haloes have been carefully investigated and a simple fitting formula is provided. The neutrino profile, unlike the cold dark matter one, is found to be cored with core size and central density that depend on the neutrino mass, redshift and mass of the halo, for halos of masses larger than similar to 10(13.5) h(-1) M-circle dot. For lower masses the neutrino profile is best fitted by a simple power-law relation in the range probed by the simulations. The results we obtain are numerically converged in terms of neutrino profiles at the 10% level for scales above similar to 200 h(-1) kpc at z = 0, and are stable with respect to box-size and starting redshift of the simulation. Our findings are particularly important in view of upcoming large-scale structure surveys, like Euclid, that are expected to probe the non-linear regime at the percent level with lensing and clustering observations.