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: Recent Cosmic Microwave Background (CMB) results from the Planck satellite, combined with previous CMB data and Hubble constant measurements from the Hubble Space Telescope, provide a constraint on the effective number of relativistic degrees of freedom 3.62(-0.48)(+0.50) at 95% CL. New Planck data provide a unique opportunity to place limits on models containing relativistic species at the decoupling epoch. We present here the bounds on sterile neutrino models combining Planck data with galaxy clustering information. Assuming N-eff active plus sterile massive neutrino species, in the case of a Planck+WP+HighL+HST analysis we find m(nu,sterile)(eff) < 0.36 eV and 3.14 < N-eff < 4.15 at 95% CL, while using Planck+WP+HighL data in combination with the full shape of the galaxy power spectrum from the Baryon Oscillation Spectroscopic Survey BOSS Data Relase 9 measurements, we find that 3.30 < N-eff < 4.43 and m(nu,sterile)(eff) < 0.33 eV both at 95% CL with the three active neutrinos having the minimum mass allowed in the normal hierarchy scheme, i.e. Sigma m(nu) similar to 0.06 eV. These values compromise the viability of the (3 + 2) massive sterile neutrino models for the parameter region indicated by global fits of neutrino oscillation data. Within the (3 + 1) massive sterile neutrino scenario, we find m(nu,sterile)(eff) < 0.34 eV at 95% CL. While the existence of one extra sterile massive neutrino state is compatible with current oscillation data, the values for the sterile neutrino mass preferred by oscillation analyses are significantly higher than the current cosmological bound. We review as well the bounds on extended dark sectors with additional light species based on the latest Planck CMB observations.

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.