Abstract: We present robust, model-marginalized limits on both the total neutrino mass (E m1,) and abundances (Neff) to minimize the role of parametrizations, priors and models when extracting neutrino properties from cosmology. The cosmological observations we consider are cosmic microwave background temperature fluctuation and polarization measurements, supernovae Ia luminosity distances, baryon acoustic oscillation observations and determinations of the growth rate parameter from the Data Release 16 of the Sloan Digital Sky Survey IV. The degenerate neutrino mass spectrum (which implies the prior sigma m(1), > 0) is weakly or moderately preferred over the normal and inverted hierarchy possibilities, which imply the priors sigma m(1), > 0.06 and sigma m(1), > 0.1 eV respectively. Concerning the underlying cosmological model, the ACDM minimal scenario is almost always strongly preferred over the possible extensions explored here. The most constraining 95% CL bound on the total neutrino mass in the ACDM + sigma m(1), picture is sigma m(1), < 0.087 eV. The parameter N-eff is restricted to 3.08 +/- 0.17 (68% CL) in the ACDM + Neff model. These limits barely change when considering the ACDM + sigma m(1), + Neff scenario. Given the robustness and the strong constraining power of the cosmological measurements employed here, the model -marginalized posteriors obtained considering a large spectra of nonminimal cosmologies are very close to the previous bounds, obtained within the ACDM framework in the degenerate neutrino mass spectrum. Future cosmological measurements may improve the current Bayesian evidence favoring the degenerate neutrino mass spectra, challenging therefore the consistency between cosmological neutrino mass bounds and oscillation neutrino measurements, and potentially suggesting a more complicated cosmological model and/or neutrino sector.

Abstract: Cosmological bounds on neutrinos and additional hypothetical light thermal relics, such as QCD axions, are currently among the most restrictive ones. These limits mainly rely on cosmic microwave background temperature anisotropies. Nonetheless, one of the largest cosmological signatures of thermal relics is that on gravitational lensing, due to their free-streaming behavior before their nonrelativistic period. We investigate late-time only hot-relic mass constraints, primarily based on recently released lensing data from the Atacama Cosmology Telescope, both alone and in combination with lensing data from the Planck satellite. Additionally, we consider other local probes, such as baryon acoustic oscillations measurements, shear-shear, galaxy-galaxy, and galaxy-shear correlation functions from the dark energy survey, and distance moduli measurements from Type-Ia Supernovae. The tightest bounds we find are Sigma m(v) < 0.43 eV and m(a) < 1.1 eV, both at 95% CL Interestingly, these limits are still much stronger than those found on e.g., laboratory neutrino mass searches, reassessing the robustness of the extraction of thermal relic properties via cosmological observations. In addition, when considering lensing-only data, the significance of the Hubble constant tension is considerably reduced, while the clustering parameter sigma 8 controversy is completely absent.

Abstract: Axions currently provide the most compelling solution to the strong CP problem. These particles may be copiously produced in the early universe, including via thermal processes. Therefore, relic axions constitute a hot dark matter component and their masses are strongly degenerate with those of the three active neutrinos, as they leave identical signatures in the different cosmological observables. In addition, thermal axions, while still relativistic states, also contribute to the relativistic degrees of freedom, parameterized via N-eff. We present the cosmological bounds on the relic axion and neutrino masses, exploiting the full Planck mission data, which include polarization measurements. In the mixed hot dark matter scenario explored here, we find the tightest and more robust constraint to date on the sum of the three active neutrino masses, Sigma m nu < 0.136eV at 95% CL, as it is obtained in the very well-known linear perturbation regime. The Planck Sunyaev-Zeldovich cluster number count data further tightens this bound, providing a 95% CL upper limit of Sigma m nu < 0.126 eV in this very same mixed hot dark matter model, a value which is very close to the expectations in the inverted hierarchical neutrino mass scenario. Using this same combination of data sets we find the most stringent bound to date on the thermal axion mass, m(a) < 0.529 eV at 95% CL.