Using some of the latest cosmological data sets publicly available, we derive the strongest bounds in the literature on the sum of the three active neutrino masses, M-nu, within the assumption of a background flat Lambda CDM cosmology. In the most conservative scheme, combining Planck cosmic microwave background temperature anisotropies and baryon acoustic oscillations (BAO) data, as well as the up-to-date constraint on the optical depth to reionization (tau), the tightest 95% confidence level upper bound we find is M-nu < 0.151 eV. The addition of Planck high-l polarization data, which, however, might still be contaminated by systematics, further tightens the bound to M-nu < 0.118 eV. A proper model comparison treatment shows that the two aforementioned combinations disfavor the inverted hierarchy at similar to 64% C.L. and similar to 71% C.L., respectively. In addition, we compare the constraining power of measurements of the full- shape galaxy power spectrum versus the BAO signature, from the BOSS survey. Even though the latest BOSS full-shape measurements cover a larger volume and benefit from smaller error bars compared to previous similar measurements, the analysis method commonly adopted results in their constraining power still being less powerful than that of the extracted BAO signal. Our work uses only cosmological data; imposing the constraint M-nu > 0.06 eV from oscillations data would raise the quoted upper bounds by O(0.1 sigma) and would not affect our conclusions.

One of the most powerful cosmological data sets when it comes to constraining neutrino masses is represented by galaxy power spectrum measurements, P-gg(k). The constraining power of P-gg(k) is however severely limited by uncertainties in the modeling of the scale-dependent galaxy bias b(k). In this work we present a new proof-of-principle for a method to constrain b(k) by using the cross-correlation between the cosmic microwave background (CMB) lensing signal and galaxy maps (C-l(kappa g)) using a simple but theoretically well-motivated parametrization for b(k). We apply the method using C-l(kappa g) measured by cross-correlating Planck lensing maps and the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 (DR11) CMASS galaxy sample, and P-gg(k) measured from the BOSS DR12 CMASS sample. We detect a nonzero scale-dependence at moderate significance, which suggests that a proper modeling of b(k) is necessary in order to reduce the impact of nonlinearities and minimize the corresponding systematics. The accomplished increase in constraining power of P-gg(k) is demonstrated by determining a 95% confidence level upper bound on the sum of the three active neutrino masses M-nu of M-nu < 0.19 eV. This limit represents a significant improvement over previous bounds with comparable data sets. Our method will prove especially powerful and important as future large-scale structure surveys will overlap more significantly with the CMB lensing kernel providing a large cross-correlation signal.

It is a well known fact that galaxies are biased tracers of the distribution of matter in the Universe. The galaxy bias is usually factored as a function of redshift and scale, and approximated as being scale-independent on large, linear scales. In cosmologies with massive neutrinos, the galaxy bias defined with respect to the total matter field (cold dark matter, baryons, and non-relativistic neutrinos) also depends on the sum of the neutrino masses M-nu, and becomes scale-dependent even on large scales. This effect has been usually neglected given the sensitivity of current surveys. However, it becomes a severe systematic for future surveys aiming to provide the first detection of non-zero M-nu. The effect can be corrected for by defining the bias with respect to the density field of cold dark matter and baryons, rather than the total matter field. In this work, we provide a simple prescription for correctly mitigating the neutrino-induced scale-dependent bias effect in a practical way. We clarify a number of subtleties regarding how to properly implement this correction in the presence of redshift-space distortions and non-linear evolution of perturbations. We perform a Markov Chain Monte Carlo analysis on simulated galaxy clustering data that match the expected sensitivity of the Euclid survey. We find that the neutrino-induced scale-dependent bias can lead to important shifts in both the inferred mean value of M-nu, as well as its uncertainty, and provide an analytical explanation for the magnitude of the shifts. We show how these shifts propagate to the inferred values of other cosmological parameters correlated with M-nu, such as the cold dark matter physical density Omega(cdm)h(2) and the scalar spectral index n(s). In conclusion, we find that correctly accounting for the neutrino-induced scale-dependent bias will be of crucial importance for future galaxy clustering analyses. We encourage the cosmology community to correctly account for this effect using the simple prescription we present in our work. The tools necessary to easily correct for the neutrino-induced scale-dependent bias will be made publicly available in an upcoming release of the Boltzmann solver CLASS.

We explore cosmological constraints on the sum of the three active neutrino masses M-v in the context of dynamical dark energy (DDE) models with equation of state (EoS) parametrized as a function of redshift z by w(z) = w(0) + w(a)z/ (1 + z), and satisfying w(z) >= -1 for all z. We make use of cosmic microwave background data from the Planck satellite, baryon acoustic oscillation measurements, and supernovae la luminosity distance measurements, and perform a Bayesian analysis. We show that, within these models, the bounds on M-v do not degrade with respect to those obtained in the Lambda CDM case; in fact, the bounds arc slightly tighter, despite the enlarged parameter space. We explain our results based on the observation that, for fixed choices of w(0), w(a) such that w(z) >= -1 (but not w = -1 for all z), the upper limit on M-v is tighter than the Lambda CDM limit because of the well-known degeneracy between w and M-v. The Bayesian analysis we have carried out then integrates over the possible values of w(0)-w(a) such that w(z) >= -1, all of which correspond to tighter limits on M-v than the Lambda CDM limit. We find a 95% credible interval (C.I.) upper bound of M-v < 0.13 eV. This bound can be compared with the 95% C.I. upper bounds of M-v < 0.16 eV, obtained within the Lambda CDM model, and M-v < 0.41 eV, obtained in a DDE model with arbitrary EoS (which allows values of w < -1). Contrary to the results derived for DDE models with arbitrary EoS, we find that a dark energy component with w(z) >= -1 is unable to alleviate the tension between high-redshift observables and direct measurements of the Hubble constant H o . Finally, in light of the results of this analysis, we also discuss the implications for DDE models of a possible determination of the neutrino mass ordering by laboratory searches.

We study the impact of assumptions about neutrino properties on the estimation of inflationary parameters from cosmological data, with a specific focus on the allowed contours in the n(s)/r plane, where n(s) is the scalar spectral index and r is the tensor-to-scalar ratio. We study the following neutrino properties: (i) the total neutrino mass M-i = Sigma(i)m(i) (where the index i = 1, 2, 3 runs over the three neutrino mass eigenstates); (ii) the number of relativistic degrees of freedom N-eff at the time of recombination; and (iii) the neutrino hierarchy. Whereas previous literature assumed three degenerate neutrino masses or two massless neutrino species (approximations that clearly do not match neutrino oscillation data), we study the cases of normal and inverted hierarchy. Our basic result is that these three neutrino properties induce < 1 sigma shift of the probability contours in the n(s)/r plane with both current or upcoming data. We find that the choice of neutrino hierarchy (normal, inverted, or degenerate) has a negligible impact. However, the minimal cutoff on the total neutrino mass M-v,M-min = 0 that accompanies previous works using the degenerate hierarchy does introduce biases in the n(s)/r plane and should be replaced by M-v,M-min = 0.059 eV as required by oscillation data. Using current cosmic microwave background (CMB) data from Planck and Bicep/Keck, marginalizing over the total neutrino mass M-v and over r can lead to a shift in the mean value of ns of similar to 0.3 sigma toward lower values. However, once baryon acoustic oscillation measurements are included, the standard contours in the n(s)/r plane are basically reproduced. Larger shifts of the contours in the n(s)/r plane (up to 0.8 sigma) arise for nonstandard values of N-eff. We also provide forecasts for the future CMB experiments Cosmic Origins Explorer (COrE, satellite) and Stage-IV (ground-based) and show that the incomplete knowledge of neutrino properties, taken into account by a marginalization over M-v, could induce a shift of similar to 0.4 sigma toward lower values in the determination of ns (or a similar to 0.8 sigma shift if one marginalizes over N-eff). Comparison to specific inflationary models is shown. Imperfect knowledge of neutrino properties must be taken into account properly, given the desired precision in determining whether or not inflationary models match the future data.