Abstract: We examine interactions between dark matter and dark energy in light of the latest cosmological observations, focusing on a specific model with coupling proportional to the dark energy density. Our data includes Cosmic Microwave Background (CMB) measurements from the Planck 2018 legacy data release, late-time measurements of the expansion history from Baryon Acoustic Oscillations (BAO) and Supernovae Type Ia (SNeIa), galaxy clustering and cosmic shear measurements from the Dark Energy Survey Year 1 results, and the 2019 local distance ladder measurement of the Hubble constant H-0 from the Hubble Space Telescope. Considering Planck data both in combination with BAO or SNeIa data reduces the H-0 tension to a level which could possibly be compatible with a statistical fluctuation. The very same model also significantly reduces the Omega(m) – sigma(8) tension between CMB and cosmic shear measurements. Interactions between the dark sectors of our Universe remain therefore a promising joint solution to these persisting cosmological tensions.

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: One of the targets of future cosmic microwave background (CMB) and baryon acoustic oscillation measurements is to improve the current accuracy in the neutrino sector and reach a much better sensitivity on extra dark radiation in the early Universe. In this paper, we study how these improvements can be translated into constraining power for well-motivated extensions of the standard model of elementary particles that involve axions thermalized before the quantum chromodynamics (QCD) phase transition by scatterings with gluons. Assuming a fiducial Lambda cold dark matter cosmological model, we simulate future data for Stage-IV CMB-like and Dark Energy Spectroscopic Instrument (DESI)-like surveys and analyse a mixed scenario of axion and neutrino hot dark matter. We further account also for the effects of these QCD axions on the light element abundances predicted by big bang nucleosynthesis. The most constraining forecasted limits on the hot relic masses are m(a) less than or similar to 0.92 eV and n-ary sumation m(nu) less than or similar to 0.12 eV at 95 per cent Confidence Level, showing that future cosmic observations can substantially improve the current bounds, supporting multimessenger analyses of axion, neutrino, and primordial light element properties.

Abstract: We present up-to-date cosmological bounds on the sum of active neutrino masses as well as on extended cosmological scenarios with additional thermal relics, as thermal axions or sterile neutrino species. Our analyses consider all the current available cosmological data in the beginning of year 2014, including the very recent and most precise baryon acoustic oscillation measurements from the Baryon Oscillation Spectroscopic Survey. In the minimal three-active-neutrino scenario, we find Sigma m(nu) < 0.22 eV at 95% C.L. from the combination of cosmic microwave background (CMB), baryon acoustic oscillation, and Hubble Space Telescope measurements of the Hubble constant. A nonzero value for the sum of the three active neutrino masses of similar to 0.3 eV is significantly favored at more than three standard deviations when adding the constraints on s 8 and Om from the Planck cluster catalog on galaxy number counts. This preference for nonzero thermal relic masses disappears almost completely in both the thermal axion and massive sterile neutrino schemes. Extra light species contribute to the effective number of relativistic degrees of freedom, parametrized via N-eff. We found that when the recent detection of B mode polarization from the BICEP2 experiment is considered, an analysis of the combined CMB data in the framework of LCDM + r models gives N-eff = 3.90 +/- 0.42, suggesting the presence of an extra relativistic relic at more than 95% C.L. from CMB-only data.