Sterile neutrinos with a mass in the eV range have been invoked as a possible explanation of a variety of short baseline (SBL) neutrino oscillation anomalies. However, if one considers neutrino oscillations between active and sterile neutrinos, such neutrinos would have been fully thermalised in the early universe, and would be therefore in strong conflict with cosmological bounds. In this study we first update cosmological bounds on the mass and energy density of eV-scale sterile neutrinos. We then perform an updated study of a previously proposed model in which the sterile neutrino couples to a new light pseudoscalar degree of freedom. Consistently with previous analyses, we find that the model provides a good fit to all cosmological data and allows the high value of $H_0$ measured in the local universe to be consistent with measurements of the cosmic microwave background. However, new high $\ell$ polarisation data constrain the sterile neutrino mass to be less than approximately 1~eV in this scenario. Finally, we combine the cosmological bounds on the pseudoscalar model with a Bayesian inference analysis of SBL data and conclude that only a sterile mass in a narrow range around 1~eV remains consistent with both cosmology and SBL data.

We study the impact of assumptions made about the neutrino mass ordering on cosmological parameter estimation with the purpose of understanding whether in the future it will be possible to infer the specific neutrino mass distribution from cosmological data. We find that although the commonly used assumption of a degenerate neutrino hierarchy is manifestly wrong and leads to changes in cosmological observables such as the cosmic microwave background and large scale structure compared to the correct (normal or inverted) neutrino hierarchy, the induced changes are so small that even with extremely optimistic assumptions about future data they will remain undetectable. We are thus able to conclude that while cosmology can probe the neutrino contribution to the cosmic energy density extremely precisely (and hence provide a detection of a non-zero total neutrino mass at high significance), it will not be possible to directly measure the individual neutrino masses.

We consider an Early Dark Energy (EDE) cosmological model, and perform an analysis which takes into account both background and perturbation effects via the parameters $c^{2}_{\rm eff}$ and $c^{2}_{\rm vis}$, representing effective sound speed and viscosity, respectively. By using the latest available data we derive constraints on the amount of dark energy at early times and the present value of the equation of state. Our focus is on the effect that early dark energy has on the Cosmic Microwave Background (CMB) data, including polarization and lensing, in a generalized parameter space including a varying total neutrino mass, and tensor to scalar ratio, besides the 6 standard parameters of the minimal cosmological model. We find that the inclusion of Baryonic Acoustic Oscillations (BAO) data and CMB lensing significantly improves the constraints on the EDE parameters, while other high redshift data like the Quasar Hubble diagram and the Lyman-$\alpha$ forest BAO have instead a negligible impact. We find $\Omega_{\rm eDE} < 0.00391 $ and $w_{0} < -0.9500 $ at the $95 \%$ C.L. for EDE accounting for its clustering through the inclusion of perturbation dynamics. This limit becomes weaker $\Omega_{\rm eDE} < 0.00337$ if perturbations are neglected. The constraints on the EDE parameters are remarkably stable even when $\Sigma m_{\nu}$, and $r$ parameters are varied, with weak degeneracies between $\Omega_{\rm eDE}$ and $r$ or $\Sigma m_{\nu}$. Perturbation parameters are not constrained with current data sets, and tensions between the CMB derived $H_0$ and $\sigma_8$ values and those measured with local probes are not eased. This work demonstrates the capability of CMB probes to constrain the total amount of EDE well below the percent level.

We present a tomographic cosmic shear analysis of the Kilo-Degree Survey (KiDS) combined with the VISTA Kilo-Degree Infrared Galaxy Survey (VIKING). This is the first time that a full optical to near-infrared data set has been used for a wide-field cosmological weak lensing experiment. This unprecedented data, spanning $450~$deg$^2$, allows us to improve significantly the estimation of photometric redshifts, such that we are able to include robustly higher-redshift sources for the lensing measurement, and - most importantly - solidify our knowledge of the redshift distributions of the sources. Based on a flat $\Lambda$CDM model we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.737_{-0.036}^{+0.040}$ in a blind analysis from cosmic shear alone. The tension between KiDS cosmic shear and the Planck-Legacy CMB measurements remains in this systematically more robust analysis, with $S_8$ differing by $2.3\sigma$. This result is insensitive to changes in the priors on nuisance parameters for intrinsic alignment, baryon feedback, and neutrino mass. KiDS shear measurements are calibrated with a new, more realistic set of image simulations and no significant B-modes are detected in the survey, indicating that systematic errors are under control. When calibrating our redshift distributions by assuming the 30-band COSMOS-2015 photometric redshifts are correct (following the Dark Energy Survey and the Hyper Suprime-Cam Survey), we find the tension with Planck is alleviated. The robust determination of source redshift distributions remains one of the most challenging aspects for future cosmic shear surveys.

We present constraints on Horndeski gravity from a combined analysis of cosmic shear, galaxy-galaxy lensing and galaxy clustering from $450\,\mathrm{deg}^2$ of the Kilo-Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey. The Horndeski class of dark energy/modified gravity models includes the majority of universally coupled extensions to $\Lambda$CDM with one scalar field in addition to the metric. We study the functions of time that fully describe the evolution of linear perturbations in Horndeski gravity. Our results are compatible throughout with a $\Lambda$CDM model. By imposing gravitational wave constraints, we fix the tensor speed excess to zero and consider a subset of models including e.g. quintessence and $f(R)$ theories. Assuming proportionality of the Horndeski functions $\alpha_B$ and $\alpha_M$ (kinetic braiding and the Planck mass run rate, respectively) to the dark energy density fraction $\Omega_{\mathrm{DE}}(a) = 1 - \Omega_{\mathrm{m}}(a)$, we find for the proportionality coefficients $\hat{\alpha}_B = 0.20_{-0.33}^{+0.20} \,$ and $\, \hat{\alpha}_M = 0.25_{-0.29}^{+0.19}$. Our value of $S_8 \equiv \sigma_8 \sqrt{\Omega_{\mathrm{m}}/0.3}$ is in better agreement with the $Planck$ estimate when measured in the enlarged Horndeski parameter space than in a pure $\Lambda$CDM scenario. In our joint three-probe analysis we report a downward shift of the $S_8$ best fit value from the $Planck$ measurement of $\Delta S_8 = 0.016_{-0.046}^{+0.048}$ in Horndeski gravity, compared to $\Delta S_8 = 0.059_{-0.039}^{+0.040}$ in $\Lambda$CDM. Our constraints are robust to the modelling uncertainty of the non-linear matter power spectrum in Horndeski gravity. Our likelihood code for multi-probe analysis in both $\Lambda$CDM and Horndeski gravity is publicly available at http://github.com/alessiospuriomancini/KiDSHorndeski .

Several interesting Dark Matter (DM) models invoke a dark sector leading to two types of relic particles, possibly interacting with each other: non-relativistic DM, and relativistic Dark Radiation (DR). These models have interesting consequences for cosmological observables, and could in principle solve problems like the small-scale cold DM crisis, Hubble tension, and/or low $\sigma_8$ value. Their cosmological behaviour is captured by the ETHOS parametrisation, which includes a DR-DM scattering rate scaling like a power-law of the temperature, $T^n$. Scenarios with $n=0$, $2$, or $4$ can easily be realised in concrete dark sector set-ups. Here we update constraints on these three scenarios using recent CMB, BAO, and high-resolution Lyman-$\alpha$ data. We introduce a new Lyman-$\alpha$ likelihood that is applicable to a wide range of cosmological models with a suppression of the matter power spectrum on small scales. For $n=2$ and $4$, we find that Lyman-$\alpha$ data strengthen the CMB+BAO bounds on the DM-DR interaction rate by many orders of magnitude. However, models offering a possible solution to the missing satellite problem are still compatible with our new bounds. For $n=0$, high-resolution Lyman-$\alpha$ data bring no stronger constraints on the interaction rate than CMB+BAO data, except for extremely small values of the DR density. Using CMB+BAO data and a theory-motivated prior on the minimal density of DR, we find that the $n=0$ model can reduce the Hubble tension from $4.1\sigma$ to $2.7\sigma$, while simultaneously accommodating smaller values of the $\sigma_8$ and $S_8$ parameters hinted by cosmic shear data.

We forecast the sensitivity of thirty-five different combinations of future Cosmic Microwave Background and Large Scale Structure data sets to cosmological parameters and to the total neutrino mass. We work under conservative assumptions accounting for uncertainties in the modelling of systematics. In particular, for galaxy redshift surveys, we remove the information coming from non-linear scales. We use Bayesian parameter extraction from mock likelihoods to avoid Fisher matrix uncertainties. Our grid of results allows for a direct comparison between the sensitivity of different data sets. We find that future surveys will measure the neutrino mass with high significance and will not be substantially affected by potential parameter degeneracies between neutrino masses, the density of relativistic relics, and a possible time-varying equation of state of Dark Energy.

Theoretical uncertainties on non-linear scales are among the main obstacles to exploit the sensitivity of forthcoming galaxy and hydrogen surveys like Euclid or the Square Kilometre Array (SKA). Here, we devise a new method to model the theoretical error that goes beyond the usual cut-off on small scales. The advantage of this more efficient implementation of the non-linear uncertainties is tested through a Markov-Chain-Monte-Carlo (MCMC) forecast of the sensitivity of Euclid and SKA to the parameters of the standard $\Lambda$CDM model, including massive neutrinos with total mass $M_\nu$, and to 3 extended scenarios, including 1) additional relativistic degrees of freedom ($\Lambda$CDM + $M_\nu$ + $N_\mathrm{eff}$), 2) a deviation from the cosmological constant ($\Lambda$CDM + $M_\nu$ + $w_0$), and 3) a time-varying dark energy equation of state parameter ($\Lambda$CDM + $M_\nu$ + $\left(w_0,w_a \right)$). We compare the sensitivity of 14 different combinations of cosmological probes and experimental configurations. For Euclid combined with Planck, assuming a plain cosmological constant, our method gives robust predictions for a high sensitivity to the primordial spectral index $n_{\rm s}$ ($\sigma(n_s)=0.00085$), the Hubble constant $H_0$ ($\sigma(H_0)=0.141 \, {\rm km/s/Mpc}$), the total neutrino mass $M_\nu$ ($\sigma(M_\nu)=0.020 \, {\rm eV}$). Assuming dynamical dark energy we get $\sigma(M_\nu)=0.030 \, {\rm eV}$ for the mass and $(\sigma(w_0), \sigma(w_a)) = (0.0214, 0.071)$ for the equation of state parameters. The predicted sensitivity to $M_\nu$ is mostly stable against the extensions of the cosmological model considered here. Interestingly, a significant improvement of the constraints on the extended model parameters is also obtained when combining Euclid with a low redshift HI intensity mapping survey by SKA1, demonstrating the importance of the synergy of Euclid and SKA.

Halos and galaxies are tracers of the underlying dark matter structures. While their bias is well understood in the case of a simple Universe composed dominantly of dark matter, the relation becomes more complex in the presence of massive neutrinos. Indeed massive neutrinos introduce rich dynamics in the process of structure formation leading to scale-dependent bias. We study this process from the perspective of general relativity employing a simple spherical collapse model. We find a characteristic signature at the neutrino free-streaming scale in addition to a large-scale feature from general relativity. The scale-dependent halo bias opposes the suppression in the matter distribution due to neutrino free-streaming and leads to corrections of a few percent in the halo power spectrum. It is not only sensitive to the sum of the neutrino-masses, but respond to the individual masses. Accurate models for the neutrino bias are a crucial ingredient for the future data analysis and play an important role in constraining the neutrino masses.

In cosmologies with massive neutrinos, the galaxy bias defined with respect to the total matter field (cold dark matter, baryons, and non-relativistic neutrinos) depends on the sum of the neutrino masses $M_{\nu}$, and becomes scale-dependent even on large scales. This effect has been usually neglected given the sensitivity of current surveys, but becomes a severe systematic for future surveys aiming to provide the first detection of non-zero $M_{\nu}$. The effect can be corrected for by defining the bias with respect to the density field of cold dark matter and baryons instead of the total matter field. In this work, we provide a simple prescription for correctly mitigating the neutrino-induced scale-dependent bias effect in a practical way. We clarify a number of subtleties regarding how to properly implement this correction in the presence of redshift-space distortions and non-linear evolution of perturbations. We perform a MCMC analysis on simulated galaxy clustering data that match the expected sensitivity of the \textit{Euclid} survey. We find that the neutrino-induced scale-dependent bias can lead to important shifts in both the inferred mean value of $M_{\nu}$, as well as its uncertainty. We show how these shifts propagate to other cosmological parameters correlated with $M_{\nu}$, such as the cold dark matter physical density $\Omega_{cdm} h^2$ and the scalar spectral index $n_s$. In conclusion, we find that correctly accounting for the neutrino-induced scale-dependent bias will be of crucial importance for future galaxy clustering analyses. We encourage the cosmology community to correctly account for this effect using the simple prescription we present in our work. The tools necessary to easily correct for the neutrino-induced scale-dependent bias will be made publicly available in an upcoming release of the Boltzmann solver \texttt{CLASS}.

We present a Bayesian analysis of the combination of current neutrino oscillation, neutrinoless double beta decay and CMB observations. Our major goal is to carefully investigate the possibility to single out one neutrino mass ordering, Normal Ordering or Inverted Ordering, with current data. Two possible parametrizations (three neutrino masses versus the lightest neutrino mass plus the two oscillation mass splittings) and priors (linear versus logarithmic) are examined. We find that the preference for NO is only driven by neutrino oscillation data. Moreover, the values of the Bayes factor indicate that the evidence for NO is strong only when the scan is performed over the three neutrino masses with logarithmic priors; for every other combination of parameterization and prior, the preference for NO is only weak. As a by-product of our Bayesian analyses, we are able to a) compare the Bayesian bounds on the neutrino mixing parameters to those obtained by means of frequentist approaches, finding a very good agreement; b) determine that the lightest neutrino mass plus the two mass splittings parametrization, motivated by the physical observables, is strongly preferred over the three neutrino mass eigenstates scan and c) find that there is a weak-to-moderate preference for logarithmic priors. These results establish the optimal strategy to successfully explore the neutrino parameter space, based on the use of the oscillation mass splittings and a logarithmic prior on the lightest neutrino mass. We also show that the limits on the total neutrino mass $\sum m_\nu$ can change dramatically when moving from one prior to the other. These results have profound implications for future studies on the neutrino mass ordering, as they crucially state the need for self-consistent analyses which explore the best parametrization and priors, without combining results that involve different assumptions.

One of the open questions in modern cosmology is the small scale crisis of the cold dark matter paradigm. Increasing attention has recently been devoted to self-interacting dark matter models as a possible answer. However, solving the so-called "missing satellites" problem requires in addition the presence of an extra relativistic particle (dubbed dark radiation) scattering with dark matter in the early universe. Here we investigate the impact of different theoretical models devising dark matter dark radiation interactions on large scale cosmological observables. We use cosmic microwave background data to put constraints on the dark radiation component and its coupling to dark matter. We find that the values of the coupling allowed by the data imply a cut-off scale of the halo mass function consistent with the one required to match the observations of satellites in the Milky Way.

The interplay between cosmology and earth based experiments is crucial in order to pin down neutrino physics. Indeed cosmology can provide very tight, yet model dependent, constraints on some neutrino properties. Here we focus on the neutrino mass sum, reviewing the up to date current bounds and showing the results of our forecast of the sensitivity of future experiments. Finally, we discuss the case for sterile neutrinos, explaining how non standard sterile neutrino self-interactions can reconcile the oscillation anomalies with cosmology.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

Future Cosmic Microwave Background experiments together with upcoming galaxy and 21-cm surveys will provide extremely accurate measurements of different cosmological observables located at different epochs of the cosmic history. The new data will be able to constrain the neutrino mass sum with the best precision ever. In order to exploit the complementarity of the different redshift probes, a deep understanding of the physical effects driving the impact of massive neutrinos on CMB and large scale structures is required. The goal of this work is to describe these effects, assuming a summed neutrino mass close to its minimum allowed value. We find that parameter degeneracies can be removed by appropriate combinations, leading to robust and model independent constraints. A joint forecast of the sensitivity of Euclid and DESI surveys together with a CORE-like CMB experiment leads to a $1\sigma$ uncertainty of $14$~meV on the summed neutrino mass. However this particular combination gives rise to a peculiar degeneracy between $M_\nu$ and the optical depth at reionization. Independent constraints from 21-cm surveys can break this degeneracy and decrease the $1\sigma$ uncertainty down to $12$~meV.

The Short BaseLine (SBL) neutrino oscillation anomalies hint at the presence of a sterile neutrino with a mass of around 1 eV. However, such a neutrino is incompatible with cosmological data, in particular observations of the Cosmic Microwave Background (CMB) anisotropies. However, this conclusion can change by invoking new physics. One possibility is to introduce a secret interaction in the sterile neutrino sector mediated by a light pseudoscalar. In this pseudoscalar model, CMB data prefer a sterile neutrino mass that is fully compatible with the mass ranges suggested by SBL anomalies. In addition, this model predicts a value of the Hubble parameter which is completely consistent with local measurements.

We study in detail how neutrino perturbations can be followed in linear theory by using only terms up to $l=2$ in the Boltzmann hierarchy. We provide a new approximation to the third moment and demonstrate that the neutrino power spectrum can be calculated to a precision of better than $\sim$ 5% for masses up to $\sim$ 1 eV. The matter and CMB power spectra can be calculated far more precisely and typically at least a factor of a few better than with existing approximations. We then proceed to study how the neutrino power spectrum can be reliably calculated even in the presence of non-linear gravitational clustering by using the full non-linear gravitational potential derived from semi-analytic methods based on $N$-body simulations in the Boltzmann evolution hierarchy. Our results agree extremely well with results derived from $N$-body simulations.

Sterile neutrinos in the electronvolt mass range are hinted at by a number of terrestrial neutrino experiments. However, such neutrinos are highly incompatible with data from the Cosmic Microwave Background and large scale structure. This paper discusses how charging sterile neutrinos under a new pseudoscalar interaction can reconcile eV sterile neutrinos with terrestrial neutrino data. We show that this model can reconcile eV sterile neutrinos in cosmology, providing a fit to all available data which is way better than the standard $\Lambda$CDM model with one additional fully thermalized sterile neutrino. In particular it also prefers a value of the Hubble parameter much closer to the locally measured value.

Short baseline neutrino oscillation experiments have shown hints of the existence of additional sterile neutrinos in the eV mass range. Such sterile neutrinos are incompatible with cosmology because they suppress structure formation unless they can be prevented from thermalising in the early Universe or removed by subsequent decay or annihilation. Here we present a novel scenario in which both sterile neutrinos and dark matter are coupled to a new, light pseudoscalar. This can prevent thermalisation of sterile neutrinos and make dark matter sufficiently self-interacting to have an impact on galactic dynamics and possibly resolve some of the known problems with the standard cold dark matter scenario. Even more importantly it leads to a strongly self-interacting plasma of sterile neutrinos and pseudoscalars at late times and provides an excellent fit to CMB data. The usual cosmological neutrino mass problem is avoided by sterile neutrino annihilation to pseudoscalars. The preferred value of $H_0$ is substantially higher than in standard $\Lambda$CDM and in much better agreement with local measurements.

We study the potential of a future, large-volume photometric survey to constrain the axion mass $m_a$ in the hot dark matter limit. Future surveys such as Euclid will have significantly more constraining power than current observations for hot dark matter. Nonetheless, the lowest accessible axion masses are limited by the fact that axions lighter than $\sim 0.15$ eV decouple before the QCD epoch, assumed here to occur at a temperature $T_{\rm QCD} \sim 170$ MeV; this leaves an axion population of such low density that its late-time cosmological impact is negligible. For larger axion masses, $m_a \gtrsim 0.15$ eV, where axions remain in equilibrium until after the QCD phase transition, we find that a Euclid-like survey combined with Planck CMB data can detect $m_a$ at very high significance. Our conclusions are robust against assumptions about prior knowledge of the neutrino mass. Given that the proposed IAXO solar axion search is sensitive to $m_a\lesssim 0.2$ eV, the axion mass range probed by cosmology is nicely complementary.

Despite the great progress of current cosmological measurements, the nature of the dominant component of the universe, coined dark energy, is still an open question. Early Dark Energy is a possible candidate which may also alleviate some fine tuning issues of the standard paradigm. Using the latest available cosmological data, we find that the 95% CL upper bound on the early dark energy density parameter is $\Omega_{\textrm{eDE}}$. On the other hand, the dark energy component may be a stressed and inhomogeneous fluid. If this is the case, the effective sound speed and the viscosity parameters are unconstrained by current data. Future omniscope-like $21$cm surveys, combined with present CMB data, could be able to distinguish between standard quintessence scenarios from other possible models with $2\sigma$ significance, assuming a non-negligible early dark energy contribution. The precision achieved on the $\Omega_{\textrm{eDE}}$ parameter from these $21$ cm probes could be below $\mathcal{O} (10\%)$.

We provide updated bounds on non-standard neutrino interactions based on data from the Planck satellite as well as auxiliary cosmological measurements. Two types of models are studied - A Fermi-like 4-point interaction and an interaction mediated by a light pseudoscalar - and we show that these two models are representative of models in which neutrinos either decouple or recouple in the early Universe. Current cosmological data constrain the effective 4-point coupling to be $G_X \leq \left(0.06 \, {\rm GeV}\right)^{-2}$, corresponding to $G_X \leq 2.5 \times 10^7 G_F$. For non-standard pseudoscalar interactions we set a limit on the diagonal elements of the dimensionless coupling matrix, $g_{ij}$, of $g_{ii} \leq 1.2 \times 10^{-7}$. For the off-diagonal elements which induce neutrino decay the bound is significantly stronger, corresponding to $g_{ij} \leq 2.3 \times 10^{-11}(m/0.05 \, {\rm eV})^{-2}$, or a lifetime constraint of $\tau \geq 1.2 \times 10^{9} \, {\rm s} \, (m/0.05 \, {\rm eV})^{3} \,$. This is currently the strongest known bound on this particular type of neutrino decay. We finally note that extremely strong neutrino self-interactions which completely suppress anisotropic stress over all of cosmic history are very highly disfavored by current data ($\Delta \chi^2 \sim 10^4$).

The recent discovery of B-modes in the polarization pattern of the Cosmic Microwave Background by the BICEP2 experiment has important implications for neutrino physics. We revisit cosmological bounds on light sterile neutrinos and show that they are compatible with all current cosmological data provided that the mass is relatively low. Using CMB data, including BICEP-2, we find an upper bound of $m_s < 0.85$ eV ($2\sigma$ Confidence Level). This bound is strengthened to 0.48 eV when HST measurements of $H_0$ are included. However, the inclusion of SZ cluster data from the Planck mission and weak gravitational measurements from the CFHTLenS project favours a non-zero sterile neutrino mass of $0.44^{+0.11}_{-0.16}$ eV. Short baseline neutrino oscillations, on the other hand, indicate a new mass state around 1.2 eV. This mass is highly incompatible with cosmological data if the sterile neutrino is fully thermalised ($\Delta \chi^2>10$). However, if the sterile neutrino only partly thermalises it can be compatible with all current data, both cosmological and terrestrial.

We use cosmological observations in the post-Planck era to derive limits on thermally produced cosmological axions. In the early universe such axions contribute to the radiation density and later to the hot dark matter fraction. We find an upper limit m_a < 0.67 eV at 95% C.L. after marginalising over the unknown neutrino masses, using CMB temperature and polarisation data from Planck and WMAP respectively, the halo matter power spectrum extracted from SDSS-DR7, and the local Hubble expansion rate H_0 released by the Carnegie Hubble Program based on a recalibration of the Hubble Space Telescope Key Project sample. Leaving out the local H_0 measurement relaxes the limit somewhat to 0.86 eV, while Planck+WMAP alone constrain the axion mass to 1.01 eV, the first time an upper limit on m_a has been obtained from CMB data alone. Our axion limit is therefore not very sensitive to the tension between the Planck-inferred H_0 and the locally measured value. This is in contrast with the upper limit on the neutrino mass sum, which we find here to range from 0.27 eV at 95% C.L. combining all of the aforementioned observations, to 0.84 eV from CMB data alone.

New measurements of the cosmic microwave background (CMB) by the Planck mission have greatly increased our knowledge about the Universe. Dark radiation, a weakly interacting component of radiation, is one of the important ingredients in our cosmological model which is testable by Planck and other observational probes. At the moment the possible existence of dark radiation is an unsolved question. For instance, the discrepancy between the value of the Hubble constant, H_0, inferred from the Planck data and local measurements of H_0 can to some extent be alleviated by enlarging the minimal LambdaCDM model to include additional relativistic degrees of freedom. From a fundamental physics point of view dark radiation is no less interesting. Indeed, it could well be one of the most accessible windows to physics beyond the standard model. An example of this is that sterile neutrinos, hinted at in terrestrial oscillation experiments, might also be a source of dark radiation, and cosmological observations can therefore be used to test specific particle physics models. Here we review the most recent cosmological results including a complete investigation of the dark radiation sector in order to provide an overview of models that are still compatible with new cosmological observations. Furthermore we update the cosmological constraints on neutrino physics and dark radiation properties focussing on tensions between data sets and degeneracies among parameters that can degrade our information or mimic the existence of extra species.

Cosmology and short baseline neutrino oscillation data both hint at the existence of light sterile neutrinos with masses in the 1 eV range. Here we perform a detailed analysis of the sterile neutrino scenario using both cosmological and SBL data. We have additionally considered the possibility that the extra neutrino degrees of freedom are not fully thermalised in the early universe. Even when analyzing only cosmological data we find a preference for the existence of massive sterile neutrinos in both (3+1) and (3+2) scenarios, and with the inclusion of SBL data the evidence is formally at the 3.3sigma level in the case of a (3+1) model. Interestingly, cosmological and SBL data both point to the same mass scale of approximately 1 eV. In the (3+1) framework WMAP9+SPT provide a value of the sterile mass eigenstate m_4 = (1.72 \pm 0.65) eV: this result is strenghtened by adding the prior from SBL posterior to m_4 = (1.27 \pm 0.12) eV (m_4 = (1.23 \pm 0.13) eV when SDSS is also considered in the cosmological analysis). In the (3+2) scheme, two additional, non--fully thermalized, neutrinos are compatible with the whole set of cosmological and SBL data, leading to mass values of m_4 = (0.95 \pm 0.30) eV and m_5 = (1.59 \pm 0.49) eV. The inclusion of Planck data does not change our considerations about the mass scale; concerning the extra neutrino degrees of freedom, invoking a partial thermalisation the 3+1 model is still consistent with the latest data.

Recent Cosmic Microwave Background measurements at high multipoles from the South Pole Telescope and from the Atacama Cosmology Telescope seem to disagree in their conclusions for the neutrino and dark radiation properties. In this paper we set new bounds on the dark radiation and neutrino properties in different cosmological scenarios combining the ACT and SPT data with the nine year data release of the Wilkinson Microwave Anisotropy Probe (WMAP9), Baryon Acoustic Oscillation data, Hubble Telescope measurements of the Hubble constant, and Supernovae Ia luminosity distance data. In the standard three massive neutrino case, the two high multipole probes give similar results if Baryon Acoustic Oscillation data are removed from the analyses and Hubble Telescope measurements are also exploited. A similar result is obtained within a standard cosmology with Neff massless neutrinos, although in this case the agreement between these two measurements is also improved when considering simultaneously Baryon Acoustic Oscillation data and Hubble Space Telescope measurements. In the Neff massive neutrino case the two high multipole probes give very different results regardless of the external data sets used in the combined analyses. In the case in which a dark radiation background with unknown clustering properties is also considered, SPT data seem to exclude the standard value for the dark radiation viscosity cvis=1/3 at the 2 sigma CL, finding evidence for massive neutrinos only when combining SPT data with BAO measurements.

An extra dark radiation component can be present in the universe in the form of sterile neutrinos, axions or other very light degrees of freedom which may interact with the dark matter sector. We derive here the cosmological constraints on the dark radiation abundance, on its effective velocity and on its viscosity parameter from current data in dark radiation-dark matter coupled models. The cosmological bounds on the number of extra dark radiation species do not change significantly when considering interacting schemes. We also find that the constraints on the dark radiation effective velocity are degraded by an order of magnitude while the errors on the viscosity parameter are a factor of two larger when considering interacting scenarios. If future Cosmic Microwave Background data are analysed assuming a non interacting model but the dark radiation and the dark matter sectors interact in nature, the reconstructed values for the effective velocity and for the viscosity parameter will be shifted from their standard 1/3 expectation, namely ceff=0.34 (+0.006 -0.003) and cvis=0.29 (+0.002 -0.001) at 95% CL for the future COrE mission data.

Recent results from short--baseline neutrino oscillation experiments and Cosmic Microwave Background anisotropy measurements suggest the presence of additional sterile neutrinos. In this paper we properly combine these data sets to derive bounds on the sterile neutrino masses in the 3+1 and 3+2 frameworks, finding a potentially good agreement between the two datasets. However, when galaxy clustering is included in the analysis a tension between the oscillation and cosmological data is clearly present.

Recent analyses that include cosmic microwave background (CMB) anisotropy measurements from the Atacama Cosmology Telescope and the South Pole Telescope have hinted at the presence of a dark radiation component at more than two standard deviations. However, this result depends sensitively on the assumption of an HST prior on the Hubble constant, where $H_0=73.8\pm2.4$ km/s/Mpc at 68% c.l.. From a median statistics (MS) analysis of 537 non-CMB $H_0$ measurements from Huchra's compilation we derive $H_0=68 \pm2.8$ km/s/Mpc at 68% c.l., in good agreement with the results of a recent analysis of the full Huchra list of $H_0$ measurements. This result is also fully consistent with the value of $H_0=69.7\pm2.5$ km/s/Mpc at 68% c.l. obtained from CMB measurements under assumption of the standard $\Lambda$CDM model. We show that with the MS $H_0$ prior the evidence for dark radiation is weakened to $\sim 1.2$ standard deviations. Parametrizing the dark radiation component through the effective number of relativistic degrees of freedom $N_{eff}$, we find $N_{eff}=3.98\pm0.37$ at 68% c.l. with the HST prior and $N_{eff}=3.52\pm0.39$ at 68% c.l. with the MS prior. We also discuss the implications for current limits on neutrino masses and on primordial Helium abundances.

Cosmological data have provided new constraints on the number of neutrino species and the neutrino mass. However these constraints depend on assumptions related to the underlying cosmology. Since a correlation is expected between the number of effective neutrinos N_{eff}, the neutrino mass \sum m_\nu, and the curvature of the universe \Omega_k, it is useful to investigate the current constraints in the framework of a non-flat universe. In this paper we update the constraints on neutrino parameters by making use of the latest cosmic microwave background (CMB) data from the ACT and SPT experiments and consider the possibility of a universe with non-zero curvature. We first place new constraints on N_{eff} and \Omega_k, with N_{eff} = 4.03 +/- 0.45 and 10^3 \Omega_k = -4.46 +/- 5.24. Thus, even when \Omega_k is allowed to vary, N_{eff} = 3 is still disfavored with 95% confidence. We then investigate the correlation between neutrino mass and curvature that shifts the 95% upper limit of \sum m_\nu < 0.45 eV to \sum m_\nu < 0.95 eV. Thus, the impact of assuming flatness in neutrino cosmology is significant and an essential consideration with future experiments.

Recent cosmological data have provided evidence for a "dark" relativistic background at high statistical significance. Parameterized in terms of the number of relativistic degrees of freedom Neff, however, the current data seems to indicate a higher value than the one expected in the standard scenario based on three active neutrinos. This dark radiation component can be characterized not only by its abundance but also by its clustering properties, as its effective sound speed and its viscosity parameter. It is therefore crucial to study the correlations among the dark radiation properties and key cosmological parameters, as the dark energy equation of state or the running of the scalar spectral index, with current and future CMB data. We find that dark radiation with viscosity parameters different from their standard values may be misinterpreted as an evolving dark energy component or as a running spectral index in the power spectrum of primordial fluctuations.

While the arcminute-scale Cosmic Microwave Background (CMB) anisotropies are due to secondary effects, point sources dominate the total anisotropy power spectrum. At high frequencies the point sources are primarily in the form of dusty, star-forming galaxies. Both Herschel and Planck have recently measured the anisotropy power spectrum of cosmic infrared background (CIB) generated by dusty, star-forming galaxies from degree to sub-arcminute angular scales, including the non-linear clustering of these galaxies at multipoles of 3000 to 6000 relevant to CMB secondary anisotropy studies. We scale the CIB angular power spectra to CMB frequencies and interpret the combined WMAP-7 year and arcminute-scale Atacama Cosmology Telescope (ACT) and South Pole Telescope (SPT) CMB power spectra measurements to constrain the Sunyaev-Zel'dovich (SZ) effects. Allowing the CIB clustering amplitude to vary, we constrain the amplitudes of thermal and kinetic SZ power spectra at 150 GHz.

The recent measurements of the Cosmic Microwave Background anisotropies at arcminute angular scales performed by the ACT and SPT experiments are probing the damping regime of CMB fluctuations. The analysis of these datasets unexpectedly suggests that the effective number of relativistic degrees of freedom is larger than the standard value of Neff = 3.04, and inconsistent with it at more than two standard deviations. In this paper we study the role of a mechanism that could affect the shape of the CMB angular fluctuations at those scales, namely a change in the recombination process through variations in the fine structure constant. We show that the new CMB data significantly improve the previous constraints on variations of {\alpha}, with {\alpha}/{\alpha}0 = 0.984 \pm 0.005, i.e. hinting also to a more than two standard deviation from the current, local, value {\alpha}0. A significant degeneracy is present between {\alpha} and Neff, and when variations in the latter are allowed the constraints on {\alpha} are relaxed and again consistent with the standard value. Deviations of either parameter from their standard values would imply the presence of new, currently unknown physics.

Cosmological measurements are affected by the energy density of both active and sterile massive neutrinos. We extend here a recent analysis of current cosmological data to non minimal cosmologies. Several possible scenarios are examined: a constant w \neq -1 dark energy equation of state, a non flat universe, a time varying dark energy component and coupled dark matter dark energy universes or modified gravity scenarios. When considering cosmological data only, (3+2) massive neutrino models with ~0.5 eV sterile species are allowed at 95% CL. This scenario has been shown to reconcile reactor, LSND and MiniBooNE positive signals with null results from other searches. Big Bang Nucleosynthesis bounds could compromise the viability of (3+2) models if the two sterile species are fully thermalized states at decoupling.

Sterile massive neutrinos are a natural extension of the Standard Model of elementary particles. The energy density of the extra sterile massive states affects cosmological measurements in an analogous way to that of active neutrino species. We perform here an analysis of current cosmological data and derive bounds on the masses of the active and the sterile neutrino states as well as on the number of sterile states. The so-called (3+2) models with three sub-eV active massive neutrinos plus two sub-eV massive sterile species is well within the 95% CL allowed regions when considering cosmological data only. If the two extra sterile states have thermal abundances at decoupling, Big Bang Nucleosynthesis bounds compromise the viability of (3+2) models. Forecasts from future cosmological data on the active and sterile neutrino parameters are also presented. Independent measurements of the neutrino mass from tritium beta decay experiments and of the Hubble constant could shed light on sub-eV massive sterile neutrino scenarios.

Combined analyses of recent cosmological data are showing interesting hints for the presence of an extra relativistic component, coined Dark Radiation. Here we perform a new search for Dark Radiation, parametrizing it with an effective number of relativistic degrees of freedom parameter, $\neff$. We show that the cosmological data we considered are clearly suggesting the presence for an extra relativistic component with $\neff=4.08_{-0.68}^{+0.71}$ at 95% c.l.. Performing an analysis on Dark Radiation sound speed $c_{\rm eff}$ and viscosity $c_{\rm vis}$ parameters, we found $\ceff=0.312\pm0.026$ and $\cvis=0.29_{-0.16}^{+0.21}$ at 95% c.l., consistent with the expectations of a relativistic free streaming component ($\ceff=\cvis$=1/3). Assuming the presence of 3 relativistic neutrinos we constrain the extra relativistic component with $\nnus=1.10_{-0.72}^{+0.79}$ and $\ceff=0.24_{-0.13}^{+0.08}$ at 95% c.l. while $\cvis$ results as unconstrained. Assuming a massive neutrino component we obtain further indications for Dark Radiation with $\nnus=1.12_{-0.74}^{+0.86}$ at 95% c.l. .

Current cosmic microwave background (CMB) bounds on the sum of the neutrino masses assume a sudden reionization scenario described by a single parameter that determines the onset of reionization. We investigate the bounds on the neutrino mass in a more general reionization scenario based on a principal component approach. We found the constraint on the sum of the neutrino masses from CMB data can be relaxed by a $\sim$40% in a generalized reionization scenario. Moreover, the amplitude of the r.m.s. mass fluctuations $\sigma_8$ is also considerably lower providing a better consistency with a low amplitude of the Sunyaev-Zel'dovich signal.