Abstract: Interacting dark energy models may play a crucial role in explaining several important observational issues in modern cosmology and also may provide a solution to current cosmological tensions. Since the phenomenology of the dark sector could be extremely rich, one should not restrict the interacting models to have a coupling parameter which is constant in cosmic time, rather allow for its dynamical behaviour, as it is common practice in the literature when dealing with other dark energy properties, as the dark energy equation of state. We present here a compendium of the current cosmological constraints on a large variety of interacting models, investigating scenarios where the coupling parameter of the interaction function and the dark energy equation of state can be either constant or dynamical. For the most general schemes, in which both the coupling parameter of the interaction function and the dark energy equation of state are dynamical, we find 95% CL evidence for a dark energy component at early times and slightly milder evidence for a dynamical dark coupling for the most complete observational data set exploited here, which includes CMB, BAO and Supernova Ia measurements. Interestingly, there are some cases where a dark energy component different from the cosmological constant case at early times together with a coupling different from zero today, can alleviate both the H-0 and S-8 tension for the full dataset combination considered here. Due to the energy exchange among the dark sectors, the current values of the matter energy density and of the clustering parameter sigma(8) are shifted from their ACDM-like values. This fact makes future surveys, especially those focused on weak lensing measurements, unique tools to test the nature and the couplings of the dark energy sector. (c) 2023 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons .org /licenses /by /4 .0/).

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: 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.