Abstract: Type Ia Supernovae (SNeIa) used as standardizable candles have been instrumental in the discovery of cosmic acceleration, usually attributed to some form of dark energy (DE). Recent studies have raised the issue of whether intrinsic SNeIa luminosities might evolve with redshift. While the evidence for cosmic acceleration is robust to this possible systematic, the question remains of how much the latter can affect the inferred properties of the DE component responsible for cosmic acceleration. This is the question we address in this work. We use SNeIa distance moduli measurements from the Pantheon and JLA samples. We consider models where the DE equation of state is a free parameter, either constant or time-varying, as well as models where DE and dark matter interact, and finally a model-agnostic parametrization of effects due to modified gravity (MG). When SNeIa data are combined with Cosmic Microwave Background (CMB) temperature and polarization anisotropy measurements, we find strong degeneracies between parameters governing the SNeIa systematics, the DE parameters, and the Hubble constant H-0. These degeneracies significantly broaden the DE parameter uncertainties, in some cases leading to O(sigma) shifts in the central values. However, including low-redshift Baryon Acoustic Oscillation and Cosmic Chronometer measurements, as well as CMB lensing measurements, considerably improves the previous constraints, and the only remaining effect of the examined systematic is a less than or similar to 40% broadening of the uncertainties on the DE parameters. The constraints we derive on the MG parameters are instead basically unaffected by the systematic in question. We therefore confirm the overall soundness of dark energy properties.

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.