Abstract: We present model-marginalized limits on mixed hot dark matter scenarios, which consider both thermal neutrinos and thermal QCD axions. A novel aspect of our analyses is the inclusion of small-scale cosmic microwave background (CMB) observations from the Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT), together with those from the Planck satellite and baryon acoustic oscillation (BAO) data. After marginalizing over a number of well-motivated nonminimal background cosmologies, the tightest 95% Confidential Level (CL) upper bound we obtain is 0.21 eV, both for P m nu and ma, from the combination of ACT, Planck and BAO measurements. Restricting the analyses to the standard ?CDM picture, we find P m nu < 0.16 eV and ma < 0.18 eV, both at 95% CL Interestingly, the best background cosmology is never found within the minimal ?CDM plus hot relics, regardless of the datasets exploited in the analyses. The combination of Planck with either BAO, SPT or ACT prefers a universe with a nonzero value of the running in the primordial power spectrum with strong evidence. Small-scale CMB probes, both alone and combined with BAO, either prefer, with substantial evidence, nonflat universes (as in the case of SPT) or a model with a time varying dark energy component (as in the case of ACT).

Abstract: We revisit the joint constraints in the mixed hot dark matter scenario in which both thermally produced QCD axions and relic neutrinos are present. Upon recomputing the cosmological axion abundance via recent advances in the literature, we improve the state-of-the-art analyses and provide updated bounds on axion and neutrino masses. By avoiding approximate methods, such as the instantaneous decoupling approximation, and limitations due to the limited validity of the perturbative approach in QCD that forced to artificially divide the constraints from the axion-pion and the axion-gluon production channels, we find robust and self-consistent limits. We investigate the two most popular axion frameworks: KSVZ and DFSZ. From Big Bang Nucleosynthesis (BBN) light element abundances data we find for the KSVZ axion Delta N-eff < 0.31 and an axion mass bound m(a) < 0.53 eV (i.e., a bound on the axion decay constant f(a) > 1.07 x 10(7) GeV) both at 95% CL. These BBN bounds are improved to Delta N-eff < 0.14 and m(a) < 0.16 eV (f(a) > 3.56 x 10(7) GeV) if a prior on the baryon energy density from Cosmic Microwave Background (CMB) data is assumed. When instead considering cosmological observations from the CMB temperature, polarization and lensing from the Planck satellite combined with large scale structure data we find Delta N-eff < 0.23, m(a) < 0.28 eV (f(a) > 2.02 x 10(7) GeV) and Sigma m(nu) < 0.16 eV at 95% CL. This corresponds approximately to a factor of 5 improvement in the axion mass bound with respect to the existing limits. Very similar results are obtained for the DFSZ axion. We also forecast upcoming observations from future CMB and galaxy surveys, showing that they could reach percent level errors for m(a) similar to 1 eV.

Abstract: The simplest ΛCDM model provides a good fit to a large span of cosmological data but harbors large areas of phenomenology and ignorance. With the improvement of the number and the accuracy of observations, discrepancies among key cosmological parameters of the model have emerged. The most statistically significant tension is the 4 sigma to 6 sigma disagreement between predictions of the Hubble constant, H (0), made by the early time probes in concert with the 'vanilla' ΛCDM cosmological model, and a number of late time, model-independent determinations of H (0) from local measurements of distances and redshifts. The high precision and consistency of the data at both ends present strong challenges to the possible solution space and demands a hypothesis with enough rigor to explain multiple observations-whether these invoke new physics, unexpected large-scale structures or multiple, unrelated errors. A thorough review of the problem including a discussion of recent Hubble constant estimates and a summary of the proposed theoretical solutions is presented here. We include more than 1000 references, indicating that the interest in this area has grown considerably just during the last few years. We classify the many proposals to resolve the tension in these categories: early dark energy, late dark energy, dark energy models with 6 degrees of freedom and their extensions, models with extra relativistic degrees of freedom, models with extra interactions, unified cosmologies, modified gravity, inflationary models, modified recombination history, physics of the critical phenomena, and alternative proposals. Some are formally successful, improving the fit to the data in light of their additional degrees of freedom, restoring agreement within 1-2 sigma between Planck 2018, using the cosmic microwave background power spectra data, baryon acoustic oscillations, Pantheon SN data, and R20, the latest SH0ES Team Riess, et al (2021 Astrophys. J. 908 L6) measurement of the Hubble constant (H (0) = 73.2 +/- 1.3 km s(-1) Mpc(-1) at 68% confidence level). However, there are many more unsuccessful models which leave the discrepancy well above the 3 sigma disagreement level. In many cases, reduced tension comes not simply from a change in the value of H (0) but also due to an increase in its uncertainty due to degeneracy with additional physics, complicating the picture and pointing to the need for additional probes. While no specific proposal makes a strong case for being highly likely or far better than all others, solutions involving early or dynamical dark energy, neutrino interactions, interacting cosmologies, primordial magnetic fields, and modified gravity provide the best options until a better alternative comes along.

Abstract: Axions, if realized in nature, can be copiously produced in the early universe via thermal processes, contributing to the mass-energy density of thermal hot relics. In light of the most recent cosmological observations, we analyse two different thermal processes within a realistic mixed hot dark matter scenario which includes also massive neutrinos. Considering the axion-gluon thermalization channel, we derive our most constraining bounds on the hot relic masses m(a) < 7.46 eV and Sigma m(nu) < 0.114 eV both at 95 percent CL; while studying the axion-pion scattering, without assuming any specific model for the axion-pion interactions, and remaining in the range of validity of the chiral perturbation theory, our most constraining bounds are improved to m(a) < 0.91 eV and Sigma m(nu) < 0.105 eV, both at 95 percent CL. Interestingly, in both cases, the total neutrino mass lies very close to the inverted neutrino mass ordering prediction. If future terrestrial double beta decay and/or long-baseline neutrino experiments find that the nature mass ordering is the inverted one, this could rule out a wide region in the currently allowed thermal axion window. Our results therefore, strongly support multi messenger searches of axions and neutrino properties, together with joint analyses of their expected sensitivities.

Abstract: The concordance of the ACDM cosmological model in light of current observations has been the subject of an intense debate in recent months. The 2018 Planck Cosmic Microwave Background (CMB) temperature anisotropy power spectrum measurements appear at face value to favour a spatially closed Universe with curvature parameter Omega(K) < 0. This preference disappears if Baryon Acoustic Oscillation (BAO) measurements are combined with Planck data to break the geometrical degeneracy, although the reliability of this combination has been questioned due to the strong tension present between the two datasets when assuming a curved Universe. Here, we approach this issue from yet another point of view, using measurements of the full-shape (FS) galaxy power spectrum, P(k), from the Baryon Oscillation Spectroscopic Survey DR12 CMASS sample. By combining Planck data with FS measurements, we break the geometrical degeneracy and find Omega(K) = 0.0023 +/- 0.0028. This constrains the Universe to be spatially flat to sub-percent precision, in excellent agreement with results obtained using BAO measurements. However, as with BAO, the overall increase in the best-fit chi(2) suggests a similar level of tension between Planck and P(k) under the assumption of a curved Universe. While the debate on spatial curvature and the concordance between cosmological datasets remains open, our results provide new perspectives on the issue, highlighting the crucial role of FS measurements in the era of precision cosmology.