Neutrino cosmology and Planck

Relic neutrinos play an important role in the evolution of the Universe, modifying some of the cosmological observables. We summarize the main aspects of cosmological neutrinos and describe how the precision of present cosmological data can be used to learn about neutrino properties. In particular, we discuss how cosmology provides information on the absolute scale of neutrino masses, complementary to beta decay and neutrinoless double-beta decay experiments. We explain why the combination of Planck temperature data with measurements of the baryon acoustic oscillation angular scale provides a strong bound on the sum of neutrino masses, 0.23 eV at the 95% confidence level, while the lensing potential spectrum and the cluster mass function measured by Planck are compatible with larger values. We also review the constraints from current data on other neutrino properties. Finally, we describe the very good perspectives from future cosmological measurements, which are expected to be sensitive to neutrino masses close the minimum values guaranteed by flavour oscillations.


Introduction
The role of neutrinos in cosmology is one of the best examples of the very close ties that have developed between nuclear physics, particle physics, astrophysics and cosmology. Here we focus on the most interesting aspects related to the case of massive (and light) relic neutrinos, but many others that were left out can be found in specialised books [1] and reviews [2,3,4,5,6,7].
We begin with a description of the properties and evolution of the background of relic neutrinos that fills the Universe. The largest part of this paper is devoted to the impact of massive neutrinos on cosmological observables, that can be used to extract bounds on neutrino masses from present data, with emphasis on the results of the Planck satellite. Next, we review the implication of current cosmological data on other neutrino properties (relic density, leptonic asymmetry, extra sterile neutrino species). Finally, we discuss the sensitivities of future cosmological experiments to neutrino masses.

The cosmic neutrino background
The existence of a relic sea of neutrinos is a generic feature of the standard hot big bang model, in number only slightly below that of relic photons that constitute the cosmic microwave background (CMB). This cosmic neutrino background (CNB) has not been detected yet, but its presence is indirectly established by the accurate agreement between the calculated and observed primordial abundances of light elements, as well as from the analysis of the power spectrum of CMB anisotropies and other cosmological observables.
Produced at large temperatures by frequent weak interactions, flavour neutrinos (ν e,µ,τ ) were kept in equilibrium until these processes became ineffective in the course of the expansion of the early Universe. While coupled to the rest of the primeval plasma (relativistic particles such as electrons, positrons and photons), neutrinos had a momentum spectrum with an equilibrium Fermi-Dirac form with temperature T , which is just one example of the general case of particles in equilibrium (fermions or bosons, relativistic or non-relativistic), as shown e.g. in [8]. In the previous equation we have included a neutrino chemical potential µ ν that would exist in the presence of a neutrino-antineutrino asymmetry, but it has been shown that even if it exists its contribution can not be very relevant [9]. As the Universe cools, the weak interaction rate falls below the expansion rate and neutrinos decouple from the rest of the plasma. An estimate of the decoupling temperature T dec can be found by equating the thermally averaged value of the weak interaction rate Γ ν = σ ν n ν , where σ ν ∝ G 2 F is the cross section of the electron-neutrino processes with G F the Fermi constant and n ν is the neutrino number density, with the expansion rate given by the Hubble parameter H = (8πρ/3M 2 P ) 1/2 . Here ρ ∝ T 4 is the total energy density, dominated by radiation, and M P = 1/G 1/2 is the Planck mass. If we approximate the numerical factors to unity, with Γ ν ≈ G 2 F T 5 and H ≈ T 2 /M P , we obtain the rough estimate T dec ≈ 1 MeV. More accurate calculations give slightly higher values of T dec which are flavour dependent because electron neutrinos and antineutrinos are in closer contact with e ± , as shown e.g. in [2].
Although neutrino decoupling is not described by a unique T dec , it can be approximated as an instantaneous process. The standard picture of instantaneous neutrino decoupling is very simple (see e.g. [8,10]) and reasonably accurate. In this approximation, the spectrum in eq. (1) is preserved after decoupling, because both neutrino momenta and temperature redshift identically with the expansion of the Universe. In other words, the number density of non-interacting neutrinos remains constant in a comoving volume since decoupling. We will see later that neutrinos cannot possess masses much larger than 1 eV, so they were ultra-relativistic at decoupling. This is the reason why the momentum distribution in eq. (1) does not depend on the neutrino masses, even after decoupling, i.e. there is no neutrino energy in the exponential of f eq (p).
When calculating quantities related to relic neutrinos, one must consider the various possible degrees of freedom per flavour. If neutrinos are massless or Majorana particles, there are two degrees of freedom for each flavour, one for neutrinos (one negative helicity state) and one for antineutrinos (one positive helicity state). Instead, for Dirac neutrinos there are in principle twice more degrees of freedom, corresponding to the two helicity states. However, the extra degrees of freedom should be included in the computation only if they are populated and brought into equilibrium before neutrinos decouple. In practice, the Dirac neutrinos with the "wrong-helicity" states do not interact with the plasma at MeV temperatures and have a vanishingly small density with respect to the usual left-handed neutrinos (unless neutrinos have masses close to the keV range, as explained in section 6.4 of [2], but this possibility is excluded). Thus the relic density of active neutrinos does not depend on their nature, either Dirac or Majorana particles.
Shortly after neutrino decoupling the temperature drops below the electron mass, favouring e ± annihilations into photons. If one assumes that this entropy transfer did not affect the neutrinos because they were already completely decoupled, it is easy to calculate the change in the photon temperature before any e ± annihilation and after the e ± pairs disappear by assuming entropy conservation of the electromagnetic plasma. The result, T after γ /T before γ = (11/4) 1/3 1.40102, is also the ratio between the temperatures of relic photons and neutrinos T γ /T ν . During the process of e ± annihilations T γ decreases with the expansion less than the inverse of the scale factor a. Instead the temperature of the decoupled neutrinos always falls as 1/a. It turns out that the processes of neutrino decoupling and e ± annihilations are sufficiently close in time so that some relic interactions between e ± and neutrinos exist. These relic processes are more efficient for larger neutrino energies, leading to nonthermal distortions in the neutrino spectra at the percent level and a slightly smaller increase of the comoving photon temperature, as noted in a series of works listed in [1,2]. These changes modify the contribution of relativistic relic neutrinos to the total energy density which is taken into account using N eff 3.046 [11], as defined later in eq. (6). In practice, these distortions only have small consequences on the evolution of cosmological perturbations, and for many purposes they can be safely neglected. Any quantity related to relic neutrinos can be calculated after decoupling with the spectrum in eq. (1) and T ν . For instance, the number density per flavour is which leads to a present value of 113 neutrinos and antineutrinos of each flavour per cm 3 . Instead, the energy density for massive neutrinos should in principle be calculated numerically, with two well-defined analytical limits, Let us now discuss the evolution of the CNB after decoupling in the expanding Universe, which is described by the Friedmann-Robertson-Walker metric [10] where we assumed negligible spatial curvature. Here a(t) is the scale factor usually normalized to unity now (a(t 0 ) = 1) and related to the redshift z as a = 1/(1 + z).
General relativity tells us the relation between the metric and the matter and energy via the Einstein equations, whose time-time component is the Friedmann equation that gives the Hubble rate in terms of the total energy density ρ. At any time, the critical density ρ c is defined as ρ c = 3H 2 /8πG, with a current value ρ 0 c = 1.8788 × 10 −29 h 2 g cm −3 , where h ≡ H 0 /(100 km s −1 Mpc −1 ). The different contributions to the total energy density are ρ = ρ γ +ρ cdm +ρ b +ρ ν +ρ Λ , and the evolution of each component is given by the energy conservation law in an expanding Universeρ = −3H(ρ+p), where p is the pressure. Thus the homogeneous density of photons ρ γ scales like a −4 , that of non-relativistic matter (ρ cdm for cold dark matter and ρ b for baryons) like a −3 , and the cosmological constant density ρ Λ is of course time-independent. Instead, the energy density of neutrinos contributes to the radiation density at early times but they behave as matter after the non-relativistic transition.
The evolution of all densities is depicted in Fig. 1, starting at MeV temperatures until now. The density fractions Ω i ≡ ρ i /ρ c are shown in this figure, where it is easy to see which of the Universe components dominantes, fixing its expansion rate: first radiation in the form of photons and neutrinos (Radiation Domination), then matter which can be CDM, baryons and massive neutrinos at late times (Matter Domination) and finally the cosmological constant density takes over at low z (typically z < 0.5).
Massive neutrinos are the only known particles that present a late transition from radiation to matter, when their density is clearly enhanced (upper solid lines in Fig.  1). Obviously the contribution of massive neutrinos to the energy density in the nonrelativistic limit is a function of the mass (or the sum of all masses for which m i T ν ), and the present value Ω ν could be of order unity for eV masses (see section 4).

Extra radiation and the effective number of neutrinos
Together with photons, in the standard case neutrinos fix the expansion rate while the Universe is dominated by radiation. Their contribution to the total radiation content can be parametrized in terms of the effective number of neutrinos N eff , defined as where we have normalized ρ r to the photon energy density because its value today is known from the measurement of the CMB temperature. This equation is valid when neutrino decoupling is complete and holds as long as all neutrinos are relativistic. We know that the number of light neutrinos sensitive to weak interactions (flavour or active neutrinos) equals three from the analysis of the invisible Z-boson width at LEP, N ν = 2.9840±0.0082 [12], and we saw in a previous section from the analysis of neutrino decoupling that they contribute as N eff 3.046. Any departure of N eff from this last value would be due to non-standard neutrino features or to the contribution of other relativistic relics. For instance, the energy density of a hypothetical scalar particle φ in equilibrium with the same temperature as neutrinos would be ρ φ = (π/30) T 4 ν , leading to a departure of N eff from the standard value of 4/7. A detailed discussion of cosmological scenarios where N eff is not fixed to three can be found in [1,2,13].
Relativistic particles, such as neutrinos, fix the expansion rate during Big Bang Nucleosynthesis (BBN) ‡ which in turn fixes the produced abundances of light elements, and in particular that of 4 He. Thus the value of N eff is constrained at BBN from the comparison of theoretical predictions and experimental data on the primordial abundances of light elements [14,15]. This is why BBN gave the first allowed range of the number of neutrino species before accelerators. In a recent analysis of the BBN constraints on N eff [16] (see the references therein for other analyses), the authors discuss a new and more conservative approach, motivated by growing concerns on the reliability of astrophysical determinations of primordial 4 He. According to [16], BBN limits the extra radiation to ∆N eff ≤ 1 at 95% C.L.
In addition, a value of N eff different from the standard one would affect the CMB observables, as will be explained in section 7.1. We will see that the Planck satellite measurement is in very good agreement with both the standard prediction of N eff 3.046 and BBN results, in spite of a marginal preference for extra relativistic degrees of freedom (exacerbated if astrophysical measurements of H 0 are included).

Massive neutrinos as Dark Matter
Nowadays the existence of Dark Matter (DM), the dominant non-baryonic component of the matter density in the Universe, is well established. A priori, massive neutrinos are excellent DM candidates, in particular because we are certain that they exist, in contrast with other candidate particles. Their energy density in units of the critical value is Here i m i includes all masses of the neutrino states which are non-relativistic today. It is also useful to define the neutrino density fraction with respect to the matter density In order to find the contribution of relic neutrinos to the present values of Ω ν or f ν , we should consider which neutrino masses are allowed by non-cosmological data. Oscillation experiments measure ∆m 2 21 = m 2 2 − m 2 1 and ∆m 2 31 = m 2 3 − m 2 1 , the relevant differences of squared neutrino masses for solar and atmospheric neutrinos, respectively. As a reference, we take the 3σ ranges of mixing parameters from [17] (see also [18,19] Here s 2 ij = sin 2 θ ij , where θ ij (ij = 12, 23 or 13) are the three mixing angles. Unfortunately oscillation experiments are insensitive to the absolute scale of neutrino masses, because the knowledge of ∆m 2 21 > 0 and |∆m 2 31 | leads to the two possible ‡ In addition, BBN is the last cosmological epoch sensitive to neutrino flavour, because electron neutrinos and antineutrinos play a direct role in the weak processes. schemes shown in fig. 1 of [4], but leaves one neutrino mass unconstrained. These two schemes are known as normal (NH, ∆m 2 31 > 0) and inverted (IH, ∆m 2 31 < 0) mass hierarchies. In the above equation the values in parentheses correspond to the IH, otherwise the allowed regions are the same for both hierarchies. For small values of the lightest neutrino mass m 0 , i.e. m 1 (m 3 ) for NH (IH), the mass states follow a hierarchical scenario, while for masses much larger than the differences all neutrinos share in practice the same mass (degenerate region). In general, the relation between the individual masses and the total neutrino mass can be found numerically.
There are two types of laboratory experiments searching for the absolute scale of neutrino masses, a crucial piece of information for constructing models of neutrino masses and mixings. The neutrinoless double beta decay (Z, A) → (Z + 2, A) + 2e − (in short 0ν2β) is a rare nuclear processes where lepton number is violated and whose observation would mean that neutrinos are Majorana particles. If the 0ν2β process is mediated by a light neutrino, the results from neutrinoless double beta decay experiments are converted into an upper bound or a measurement of the effective mass m ββ = |c 2 12 c 2 13 m 1 + s 2 12 c 2 13 m 2 e iφ 2 + s 2 13 m 3 e iφ 3 | , where φ 1,2 are two Majorana phases that appear in lepton-number-violating processes. An important issue for 0ν2β results is related to the uncertainties on the corresponding nuclear matrix elements. For more details and the current experimental results, see [20]. Beta decay experiments, which involve only the kinematics of electrons, are in principle the best strategy for measuring directly the neutrino mass [21]. The current limits from tritium beta decay apply only to the range of degenerate neutrino masses, is the relevant parameter for beta decay experiments. The bound at 95% CL is m 0 < 2.05 − 2.3 eV from the Troitsk and Mainz experiments, respectively. This value is expected to be improved by the KATRIN project to reach a discovery potential for 0.3 − 0.35 eV masses (or a sensitivity of 0.2 eV at 90% CL). Taking into account this upper bound and the minimal values of the total neutrino mass in the normal (inverted) hierarchy, the sum of neutrino masses M ν ≡ i m i is restricted to the approximate range 0.06 (0.1) eV M ν 6 eV (12) As we discuss in the next sections, cosmology is at first order sensitive to the total neutrino mass if all states have the same number density, providing information on m 0 but blind to neutrino mixing angles or possible CP violating phases. Thus cosmological results are complementary to terrestrial experiments. The interested reader can find the allowed regions in the parameter space defined by any pair of parameters (M ν , m ββ , m β ) in [22,23,24]. The two cases involving M ν are shown in Fig. 2. Now we can find the possible present values of Ω ν in agreement the approximate bounds of eq. (12). Note that even if the three neutrinos are non-degenerate in mass, eq. (7) can be safely applied, because we know from neutrino oscillation data that at least two of the neutrino states are non-relativistic today, because both |∆m 2 31 | 1/2 0.05 eV and (∆m 2 21 ) 1/2 0.009 eV are larger than the temperature T ν 1.96 K 1.7 × 10 −4 eV. If the third neutrino state is very light and still relativistic, its relative contribution to Ω ν is negligible and eq. (7) remains an excellent approximation of the total density. One finds that Ω ν is restricted to the approximate range 0.0013 (0.0022) Ω ν 0.13 (13) where we already included that h ≈ 0.7. This applies only to the standard case of three light active neutrinos, while in general a cosmological upper bound on Ω ν has been used since the 1970s to constrain the possible values of neutrino masses. For instance, if we demand that neutrinos should not be heavy enough to overclose the Universe (Ω ν < 1), we obtain an upper bound M ν 45 eV (again fixing h = 0.7). Moreover, from present analyses of cosmological data we know that the approximate contribution of matter is Ω m 0.3, the neutrino masses should obey the stronger bound M ν 15 eV. Thus with this simple argument one obtains a bound which is roughly only a factor 2 worse than the bound from tritium beta decay, but of course with the caveats that apply to any cosmological analysis. In the three-neutrino case, these bounds should be understood in terms of m 0 = M ν /3. Dark matter particles with a large velocity dispersion such as that of neutrinos are called hot dark matter (HDM). The role of neutrinos as HDM particles has been widely discussed since the 1970s (see e.g. the historical review [25]). It was realized in the mid-1980s that HDM affects the evolution of cosmological perturbations in a particular way: it erases the density contrasts on wavelengths smaller than a mass-dependent freestreaming scale. In a universe dominated by HDM, this suppression contradicts various observations. For instance, large objects such as superclusters of galaxies form first, while smaller structures like clusters and galaxies form via a fragmentation process. This top-down scenario is at odds with the fact that galaxies seem older than clusters.
Given the failure of HDM-dominated scenarios, the attention then turned to cold dark matter (CDM) candidates, i.e. particles which were non-relativistic at the epoch when the universe became matter-dominated, which provided a better agreement with observations. Still in the mid-1990s it appeared that a small mixture of HDM in a universe dominated by CDM fitted better the observational data on density fluctuations at small scales than a pure CDM model. However, within the presently favoured ΛCDM model dominated at late times by a cosmological constant (or some form of dark energy) there is no need for a significant contribution of HDM. Instead, one can use the available cosmological data to find how large the neutrino contribution can be, as we will see later.
Before concluding this section, we would like to mention the case of a sterile neutrino with a mass of the order of a few keV's and a very small mixing with the flavour neutrinos. Such "heavy" neutrinos could be produced by active-sterile oscillations but not fully thermalized, so that they could replace the usual CDM component. But due to their large thermal velocity (slightly smaller than that of active neutrinos), they would behave as Warm Dark Matter and erase small-scale cosmological structures. Their mass can be bounded from below using Lyman-α forest data from quasar spectra, and from above using X-ray observations. The viability of this scenario was reviewed in [26].

Effects of neutrino masses on cosmology
Here we briefly describe the main cosmological observables, and their sensitivity to neutrino masses. A more detailed discussion of the effects of massive neutrinos on the evolution of cosmological perturbations can be found in sections 5.3.3 and 6.1.4 of [1].

Brief description of cosmological observables
Although there exist many different types of cosmological measurements, here we will restrict the discussion to those that are at present the more important for obtaining an upper bound or eventually a measurement of neutrino masses.
First of all, we have the CMB temperature anisotropy power spectrum, defined as the angular two-point correlation function of CMB maps δT /T (n) (n being a direction in the sky). This function is usually expanded in Legendre multipoles where P l (x) are the Legendre polynomials. For Gaussian fluctuations, all the information is encoded in the multipoles C l which probe correlations on angular scales θ = π/l. There exists interesting complementary information to the temperature power spectrum if the CMB polarization is measured, and currently we have some less precise data on the temperature × E-polarization (TE) correlation function and the E-polarization self-correlation spectrum (EE). The current Large Scale Structure (LSS) of the Universe is probed by the matter power spectrum at a given time or redhsift z. It is defined as the two-point correlation function of non-relativistic matter fluctuations in Fourier space where δ m = δρ m /ρ m . Usually P (k) refers to the matter power spectrum evaluated today (at z = 0). In the case of several fluids (e.g. CDM, baryons and non-relativistic neutrinos), the total matter perturbation can be expanded as δ m = iρ i δ i / iρ i . Because the energy density is related to the mass density of non-relativistic matter through E = mc 2 , δ m represents indifferently the energy or mass power spectrum. The shape of the matter power spectrum is affected in a scale-dependent way by the freestreaming caused by small neutrino masses of O(eV) and thus it is the key observable for constraining m ν with cosmological methods.

Neutrino free-streaming
After thermal decoupling, relic neutrinos constitute a collisionless fluid, where the individual particles free-stream with a characteristic velocity that, in average, is the thermal velocity v th . It is possible to define a horizon as the typical distance on which particles travel between time t i and t. When the Universe was dominated by radiation or matter t t i , this horizon is, as usual, asymptotically equal to v th /H, up to a numerical factor of order one. Similar to the definition of the Jeans length (see section 4.4 in [4]), we can define the neutrino free-streaming wavenumber As long as neutrinos are relativistic, they travel at the speed of light and their freestreaming length is simply equal to the Hubble radius. When they become nonrelativistic, their thermal velocity decays like where we used for the present neutrino temperature T 0 ν (4/11) 1/3 T 0 γ and T 0 γ 2.726 K. This gives for the free-streaming wavenumber during matter or Λ domination where Ω Λ and Ω m are the cosmological constant and matter density fractions, respectively, evaluated today. After the non-relativistic transition and during matter domination, the free-streaming length continues to increase, but only like (aH) −1 ∝ t 1/3 , i.e. more slowly than the scale factor a ∝ t 2/3 . Therefore, the comoving free- for neutrinos becoming non-relativistic during matter domination, the comoving freestreaming wavenumber passes through a minimum k nr at the time of the transition, i.e. when m = p = 3.15T ν and a 0 /a = (1 + z) = 2.0 × 10 3 (m/1 eV). This minimum value is found to be The physical effect of free-streaming is to damp small-scale neutrino density fluctuations: neutrinos cannot be confined into (or kept outside of) regions smaller than the freestreaming length, because their velocity is greater than the escape velocity from gravitational potential wells on those scales. Instead, on scales much larger than the free-streaming scale, the neutrino velocity can be effectively considered as vanishing, and after the non-relativistic transition the neutrino perturbations behave like CDM perturbations. In particular, modes with k < k nr are never affected by free-streaming and evolve like being in a pure ΛCDM model.

Impact of massive neutrinos on the matter power spectrum
The small initial cosmological perturbations evolve within the linear regime at early times. During matter domination, the smallest cosmological scales start evolving nonlinearily, leading to the formation of the structures we see today. In the recent universe, the largest observable scales are still related to the linear evolution, while other scales can only be understood using non-linear N-body simulations. We will not review here all the details of this complicated evolution (see [1,6,4] and references therein), but we will emphasize the main effects caused by massive neutrinos on linear scales in the framework of the standard cosmological scenario: a Λ Mixed Dark Matter (ΛMDM) model, where Mixed refers to the inclusion of some HDM component. On large scales (i.e. on wave-numbers smaller than the value k nr defined in the previous subsection), neutrino free-streaming can be ignored, and neutrino perturbations are indistinguishable from CDM perturbations. On those scales, the matter power spectrum P (k, z) can be shown to depend only on the matter density fraction today (including neutrinos), Ω m , and on the primordial perturbation spectrum. If the neutrino mass is varied with Ω m fixed, the large -scale power spectrum remains invariant.
However, the small-scale matter power spectrum P (k ≥ k nr ) is reduced in presence of massive neutrinos for at least two reasons: by the absence of neutrino perturbations in the total matter power spectrum, and by a slower growth rate of CDM/baryon perturbations at late times. The third effect has the largest amplitude. At low redhsift z 0, the step-like suppression of P (k) starts at k ≥ k nr and saturates at k ∼ 1h/Mpc with a constant amplitude ∆P (k)/P (k) −8f ν . This result was obtained by fitting numerical simulations [27], but a more accurate approximation can be derived analytically [1,4]. The full matter power spectrum can be calculated at any time and scale by Boltzmann codes such as camb [28] or class [29], that solve numerically the evolution of the cosmological perturbations. The step-like suppression of the matter power spectrum induced by various values of f ν is shown in Fig. 3. On small scales and at late times, matter density perturbations enter into the regime of non-linear clustering. The neutrino mass impact on the non-linear matter power spectrum is now modeled with rather good precision, at least within one decade above k max in wave-number space [30,31,32]. It appears that the step-like suppression is enhanced by non-linear effects up to roughly ∆P (k)/P (k) −10f ν (at redshift zero and k ∼ 1 h Mpc −1 ), and is reduced above this scale. Hence, non-linear corrections render the neutrino mass effect even more characteristic than in the linear theory, and may help to increase the sensitivity of future experiments (see also [33,34,35] for the effect of neutrino masses on even smaller scales).
Until this point, we reduced the neutrino mass effect to that of f ν or M ν . In principle, the mass splitting between the three different families for a common total mass is visible in P (k). The time at which each species becomes non-relativistic depends on individual masses m i . Hence, both the scale of the step-like suppression induced by each neutrino and the amount of suppression in the small-scale power spectrum have a small dependence on individual masses. The differences between the power spectrum of various models with the same total mass and different mass splittings was computed numerically in [36] for the linear spectrum, and [32] for the non-linear spectrum. At the moment, it seems that even the most ambitious future surveys will not be able to distinguish these mass splitting effects with a reasonable significance [37,38].

Impact of massive neutrinos on the CMB anisotropy spectrum
For masses smaller than 0.6 eV, neutrinos are still relativistic at the time of photon decoupling, and their mass cannot impact the evolution of CMB perturbations. Therefore, the effect of the mass can only appear at two levels: that of the background evolution, and that of secondary anisotropies (related to the behaviour of photon perturbations after decoupling: Integrated Sachs-Wolfe (ISW) effect, weak lensing by large scale structure, etc.) Let us first review the background effects of massive neutrinos on the CMB. Because the temperature and polarization spectrum shape is the result of several intricate effects, one cannot discuss the neutrino mass impact without specifying which other parameters are kept fixed. Neutrinos with a mass in the range from 10 −3 eV to 1 eV should be counted as radiation at the time of equality, and as non-relativistic matter today: the total non-relativistic density, parametrized by ω m = Ω m h 2 , depends on the total neutrino mass M ν = i m i . Hence, when M ν is varied, there must be a variation either in the redshift of matter-to-radiation equality z eq , or in the matter density today ω m .
This can potentially impact the CMB in three ways. A shift in the redshift of equality affects the position and amplitude of the peaks. A change in the non-relativistic matter density at late times can impact both the angular diameter distance to the last scattering surface d A (z dec ), controlling the overall position of CMB spectrum features in multipole space, and the slope of the low-l tail of the CMB spectrum, due to the late Integrated Sachs-Wolfe (ISW) effect. Out of these three effects (changes in z eq , in d A and in the late ISW), only two can be cancelled by a simultaneous variation of the total neutrino mass and of other free parameters in the ΛMDM model. Hence, the CMB spectrum is sensitive to the background effect of the total neutrino mass. In practice however, the late ISW effect is difficult to measure due to cosmic variance and CMB data alone cannot provide a useful information on sub-eV neutrino masses. If one considers extensions of the ΛMDM, this becomes even more clear: by playing with the spatial curvature, one can neutralize all three effects simultaneously. But as soon as CMB data is used in combination with other background cosmology observations (constraining for instance the Hubble parameter, the cosmological constant value or the BAO scale), some bounds can be derived on M ν .
There exists another effect of massive neutrinos on the CMB at the level of secondary anisotropies: when neutrinos become non-relativistic, they reduce the time variation of the gravitational potential inside the Hubble radius. This affects the photon temperature through the early ISW effect, and leads to a depletion in the temperature spectrum of the order of (∆C l /C l ) ∼ −(M ν /0.3 eV)% on multipoles 20 < l < 500, with a dependence of the maximum depletion scale on individual masses m i [1,39]. This effect is roughly thirty times smaller than the depletion in the small-scale matter power spectrum, ∆P (k)/P (k) ∼ −(M ν /0.01 eV)%.
We show in Figure 4 the effect on the CMB temperature spectrum of increasing the neutrino mass while keeping z eq and d A fixed: the only observed differences are then for M ν =3x0.3eV, same z eq , l peak M ν =3x0.6eV, same z eq , l peak 2 < l < 50 (late ISW effect due to neutrino background evolution) and for 50 < l < 200 (early ISW effect due to neutrino perturbations). We conclude that the CMB alone is not a very powerful tool for constraining sub-eV neutrino masses, and should be used in combination with homogeneous cosmology constraints and/or measurements of the LSS power spectrum, for instance from galaxy clustering, galaxy lensing or CMB lensing.

Current bounds on neutrino masses
In this section, we review the neutrino mass bounds that can be derived from current cosmological data. Note that the confidence limits in the next subsections are all based on the Bayesian inference method, and are given at the 95% confidence level after marginalization over all free cosmological parameters. We refer the reader to section 5.1 of [4] for a detailed discussion on this statistical method.

CMB anisotropy bounds
The best measurement of CMB temperature anisotropies comes from the first release of the Planck satellite [40,41], which significantly improved over the nine-year data of WMAP [42] on large angular scales, and over ground-based or balloon-borne (ACT, SPT, . . . ) on small scales. Moreover, the Planck temperature data offers the advantage of covering all relevant angular scales with a single experiment, hence avoiding relative calibration issues between different data sets. The publication of polarisation data has been postponed to the next Planck release. In the meantime, Planck temperature data is usually analysed together with the l ≤ 23 multipoles of WMAP E-type polarisation (referred later as WP), in order to remove degeneracies between the optical depth to reionisation and other parameters. Moreover, Planck temperature data is used only up to the multipole l = 2500: for larger l's, it is better to use higher-resolution experiments (ACT, SPT), limited to l > 2500 to avoid any overlap, and called usually"highL". Note that including information on temperature at l > 2500 can be useful for separating the various foregrounds from the CMB signal, but it has a very minor impact on cosmological parameters, which are well constrained by the Planck+WP data alone [41].
Assuming the minimal ΛCDM model with six free parameters, and promoting the total neutrino mass M ν as a seventh free parameter, one gets M ν < 0.66 eV (95%; Planck+WP+highL) [41]. As explained in section 5.4, this bound comes mainly from two physical effects: the early-ISW-induced dip at intermediate l's, and the lensing effect causing a smoothing of the power spectrum at higher l's. To show that the latter effect is important, one can treat the amplitude of the lensing potential spectrum as a free parameter (while in reality, its value is fully predicted by the underlying cosmological model). This is equivalent to partially removing the piece of information coming from the lensing of the last scattering surface. Then, the mass bound degrades by 63%.
Assuming an extended underlying cosmological model, the mass bound gets weaker. This is true especially if one allows for non-zero spatial curvature in the universe: the bounds then degrades to M ν < 0.98 eV (95%; Planck+WP+highL). However, when the density of extra relativistic relics is promoted as a new free parameter, the bound on M ν does not change significantly, showing the current CMB data is able to resolve the degeneracy between M ν and N eff , that used to exist when using older data sets.

Adding information on the cosmological expansion
The measurement of CMB anisotropies can be complemented by some data on the cosmological expansion at low redshift (z < 2). The recent expansion history can be inferred from the luminosity of type-Ia supernovae, from the angular scale of Baryon Acoustic Oscillations (BAO) reconstructed from the galaxy power spectrum at various redshifts, or from direct measurements of the current expansion rate H 0 through nearby cepheids and supernovae.
These measurements are very useful probes of neutrino masses, not directly but indirectly. Indeed, we have seen in section 5.4 that if the neutrino mass is varied at the same time as parameters controlling the late cosmological evolution (like Ω Λ or Ω k in a non-flat universe), most characteristic quantites affecting the shape of the CMB spectrum can be kept constant, including the angular scale of the first acoustic peak. Then, the effect of M ν reduces to an irrelevant late ISW effect, plus the early ISW and lensing effects discussed above. But if instead these parameters are measured directly and constrained independently, the degeneracy is removed, and varying M ν will also change the scale of the first peak. In that way, the CMB is very sensitive to M ν , since the peak scale is the quantity best constrained by CMB observations. The problem is that as long as one assumes a minimal ΛCDM model, the combination of CMB data with either BAO, H 0 or supernovae data reveals small tensions (roughly at the 2-σ level). This could mean either that some data sets have slightly underestimated systematics, or that something is not correct in the assumed underlying cosmological model. The main tension is between Planck+WP and HST measurements of H 0 , giving respectively H 0 = 67.3 ± 2.54 km s −1 Mpc −1 (95%; Planck+WP) and H 0 = 73.8 ± 4.8 km s −1 Mpc −1 (95%; HST) [43]. BAO data are more consistent with Planck than with HST and exacerbate the tension, giving H 0 = 67.8 ± 1.5 km s −1 Mpc −1 (95%; Planck+WP+BAO). Note that direct H 0 measurements probe the local value of the expansion parameter, while BAO and CMB data measure the expansion rate averaged over very large scales. However, in the standard cosmological model, the local and global H 0 are expected to differ by a tiny amount, so the tension cannot be explained in that way [44]. Promoting the total neutrino mass as a seventh free parameter does not release the tension either, since for the ΛCDM+M ν model one finds H 0 = 67.7 ± 1.8 km s −1 Mpc −1 (95%; Planck+WP+BAO) [41].
The neutrino mass bound also depends on which data is used. For CMB+BAO, one gets M ν < 0.23 eV (95%; Planck+WP+highL+BAO), while for CMB+HST one gets a stronger bound M ν < 0.18 eV (95%; Planck+WP+highL+BAO), and for CMB + supernovae the result is slightly looser, M ν < 0.25 eV (95%; Planck+WP+highL+SNLS) [41] (some of these numbers are not given in the Planck paper, but in the publicly released grid of bounds available at the ESA website §). If one takes the conservative point of view that the BAO data is the less likely to be affected by systematics because it is a pure geometrical measurement of an angle in the sky, one should mainly remember the first of these bounds (0.23 eV). This is already a very spectacular result, only a factor four higher than the minimum allowed value M ν ∼ 0.06 eV.

Adding information on large scale structure
Large scale structure can be probed by Planck itself in several ways, using lensing extraction, foregrounds or spectral distorsions.
Among these probes, the most robust is expected to be lensing extraction, because it reconstructs the matter power spectrum mainly in the linear regime. Weak lensing of CMB photons by large scale structure leads to a distorsion of the CMB maps, and creates correlations in the observed map between multipoles on different scales (unlike in a pure gaussian map). The lensing extraction technique consists in using these correlations in order to reconstruct the lensing field. In the same way, one can compute the lensing power spectrum, and get indication on the matter power spectrum over a range of scales and redshifts (typically, 1 < z < 3) [45]. CMB lensing extraction was previously advocated as a way to improve the Planck error bar on the total neutrino mass by a factor two [46]. This cannot happen if the two data sets (anisotropy spectrum and extracted lensing spectrum) pull the results in opposite directions. Unfortunately, this is exactly what happens with the first Planck data release. We have seen that the Planck temperature spectrum puts strong limits on the mass. Instead, the extracted § http://www.sciops.esa.int/wikiSI/planckpla/ in section 6.9. lensing spectrum has a marginal preference for a non-zero M ν , because the step-like suppression induced by neutrino masses in the lensing spectrum gives a better fit. As a result, the bound gets looser when including lensing extraction: M ν < 0.84 eV (95%; Planck+WP+highL+lensing) or M ν < 0.25 eV (95%; Planck+WP+BAO+lensing). It will be interesting to see how these bounds will evolve in the next couple of years, with better data and better control over systematics in the lensing extraction process. The preference for a non-zero neutrino mass in the lensing data could then go away. Instead, there is still a possibility that lensing data correctly prefers M ν > 0, while some systematics in Planck data (for instance, in the low-l likelihood) artificially disfavour M ν = 0.
Spectral distorsions of CMB maps caused by the Sunyaev-Zel'dovitch effect allow to construct a map of galaxy clusters. The distortions provide information on the position and redshift of the clusters, and even allow to estimate their mass up to some bias factor. The data can be used to make a histogram of the number of clusters as a function of redshift, or the mass of clusters as a function of redshift. These histograms can be compared to the predictions of structure formation theory (in the mildly nonlinear regime). They give another handle on the underlying matter power spectrum, and hence on neutrino masses.
The situation with Planck SZ clusters is similar to the situation with lensing: the cluster data prefers a non-zero neutrino mass, unless the bias is assumed to take extreme value, that would normally be excluded by theoretical arguments. Ref. [47] gives a combined CMB + SZ cluster result preferring a non-zero neutrino mass at the 2.5σ level, M ν = 0.22 ± 0.18 eV (95%; Planck+WP+BAO+SZ). Ref. [48] shows how the bounds depend on several assumptions made in the SZ cluster analysis. A similar trend exists with essentially all other data on the cluster mass function. The conclusion is that something remains to be understood: either systematics are underestimated in measurements of the cluster mass function, or the Planck temperature data has some incorrect feature pushing down the neutrino masses (so far, the former assumption sounds more plausible).
The CMB data can be combined with many other probes of large scale structure (galaxy spectrum, flux spectrum of Lyman-α forests in quasars, weak lensing surveys). Until now, these other data sets tend to be less constraining than those discussed here .

Neutrino density, non-thermal distorsions and leptonic asymmetry
Previous estimates of the total neutrino mass rely on the standard neutrino decoupling model, leading to a neutrino density corresponding to N eff = 3.046 as long as neutrinos Except for a recent measurement of the amplitude of the matter power spectrum, based on the BOSS galaxy redshift survey, which leads to a preference for a non-zero neutrino mass roughly at the 3-σ level [49].
are relativistic. This model could be incorrect, or missing some physical ingredients. Hence neutrinos could in principle have a different density, coming from a different neutrino-to-photon temperature ratio, or from non-thermal distorsions. For instance, this would be the case if some exotic particles would decay after neutrino decoupling, either producing extra relic neutrinos, or reheating the thermal bath.
Another possible reason leading to a larger N eff would be the production of a very large neutrino-antineutrino asymmetry in the early universe. In that case, the phasespace distribution of neutrinos at CMB times depends on the initial leptonic asymmetry of each family, and on the efficiency of neutrino oscillations in the early universe; in any case, a leptonic asymmetry produced at early times would tend to enhance the total neutrino density at late times.
If neutrinos were massless, all the effects described above (over-or under-production of neutrinos, non-thermal distorsions, leptonic asymmetry) would entirely be described by the value of N eff , because the equations of evolution of cosmological perturbations in the neutrino sector could be integrated over momentum. When neutrino masses are taken into account, this is not true anymore. In that case, N eff explicitly refers to the radiation density at early times, before the non-relativistic transition of neutrinos. Still, as long as individual neutrino masses are not too large, it is useful to analyse cosmological models with free N eff and M ν set in first approximation to its minimal value (0.06 eV), or with N eff and M ν both promoted as free parameters, with an arbitrary mass splitting. To be very rigorous, each model would require an individual analysis, but in fairly good approximation, their properties are well accounted by the two parameters N eff and M ν . This is not true anymore when individual neutrino masses can be large (of the order of one or several electron-volts), as assumed in section 7.3 on massive sterile neutrinos.
Hence it is crucial to measure N eff in order to check whether we correctly understand neutrino cosmology. However, if a value of N eff larger than three was measured in CMB or LSS data, we still would not know if this come from physics in the neutrino sector, or from other relativistic relics; even assuming that such relics do not exist, we would not be able to discriminate between the different cases (shift in temperature, asymmetry or non-thermal distorsions). However, we could learn something more from other cosmological probes, such as the study of BBN, leptogenesis and baryogenesis, etc., or from laboratory experiments. For instance, we know from a joint analysis of BBN and neutrino oscillation data that in order to be compatible with measurements of primordial element abundances, the leptonic asymmetry cannot enhance the neutrino density (at CMB and current time) above N eff 3.1 [9].
The CMB is sensitive to N eff , first, through the time of equality: when this parameter increases while all other parameters are kept fixed, equality is postponed. However, this effect can be cancelled by increasing simultaneously the matter density, by exactly the same amount. Actually, one can also renormalise the cosmological constant in order to preserve all characteristic redshifts in the evolution of the universe (radiation/matter equality, matter/Λ equality). This transformation leaves the CMB almost invariant, but not quite. Indeed, a global increase of all densities goes with an increase of H 0 , that is visible on small angular scale: it enhances the Silk damping effect (i.e., the effect of diffusion just before photon decoupling, important for l > 1000). Generally speaking, a good way to describe the main effect of N eff on the CMB is to say that it changes the scale of Silk damping relative to the scale of the sound horizon at decoupling. Concretely, this will appear as a shift in the damping envelope of high peaks relative to the position of the first peak [1,50]. This effect is not the only one remaining when all densities are equally enhanced: more neutrinos also imply stronger gravitational interactions between photons and free-streaming species before decoupling. This is especially important for scales that just crossed the Hubble scale during radiation domination. The result is a small shift and damping of the acoustic peak ("baryon drag" effect [1,51]).
As a result of all these effects, for a minimal 7-parameter model (ΛCDM + N eff ), the CMB data alone gives N eff = 3.36 +0. 68 −0.64 (95%; Planck+WP+highL), well compatible with the standard prediction N eff = 3.046. An analysis with a free lensing amplitude A L of the CMB power spectrum is described in [52].
Like for neutrino masses, data on the cosmological expansion at late time allows to tighten the bounds on N eff , because it removes the partial degeneracy observed in the CMB when all densities are increased in the same proportions. Data on BAO, supernovae or direct measurements of H 0 tighten the constraints on H 0 and/or on the matter density, and forbid to increase N eff with fixed redshifts for the two equalities (radiation/matter equality, matter/Λ equality). It is worth noticing that unlike models with a free neutrino mass, models with a free N eff have enough freedom for relaxing the tensions between different data sets. For instance, in the ΛCDM + N eff model, the bound on the expansion rate is H 0 = 69.7 +5.8 −5.3 km s −1 Mpc −1 (95%; Planck+WP), well compatible with H 0 = 73.8 ± 4.8 km s −1 Mpc −1 (95%; HST) [43]. Even when including BAO data, one gets H 0 = 69.3 +3.5 −3.4 km s −1 Mpc −1 (95%; Planck+WP+highL+BAO), still compatible with HST measurements.
From the previous discussion on the effect of N eff on the CMB, one can easily infer that there is a positive correlation between measured values of N eff and H 0 . Hence, when Planck data is combined with HST data, which favours high values of the expansion rate, one gets more than 2σ evidence for enhanced radiation, N eff = 3.62 +0. 50 −0.48 (95%; Planck+WP+highL+HST). Instead, with BAO, the evidence disappears, N eff = 3.30 +0. 54 −0.51 (95%; Planck+WP+highL+BAO). Since models with a free N eff relax the tension between the different data sets, it makes sense to combine CMB data with BAO and HST at the same time, which gives N eff = 3.52 +0.48 −0.45 (95%; Planck+WP+highL+BAO+HST), slightly more than 2σ evidence for enhanced radiation (see also [53]).
The conclusion is that if HST data are robust, then one should consider seriously the possibility that N eff exceeds the standard value, because this is one of the simplest way to relax the constraint between HST and other data sets ¶. This excess could be caused by several effects (leptonic asymmetry, non-standard neutrino phase space distribution, or any type of relativistic relics). If instead we assume that HST results are biased by systematics and should not be included, then the evidence for N eff > 3.046 becomes very weak (although a high tensor-to-scalar ratio, like the one suggested by the very recent BICEP2 results [55], also tends to slightly favour an excess in N eff [56]).

Non-standard interactions
Non-standard neutrino interactions (beyond the weak force) could occur in many extensions of the standard model, and could affect cosmological observables in several different ways. For instance, if non-standard neutrino-electron interactions exist, it was shown that they could enhance N eff at most up to 3.1 [57]. Many other cases have been studied in the literature (see the list of references in section 5.3.4 of [1]), still not covering all possibilities. A few representative cases have been studied recently in [58], with two different assumptions concerning the type of self-interaction experienced in the neutrino sector. When the self-interaction is very efficient, neutrinos tend to behave like a relativistic fluid, instead of free-streaming particles with anisotropic stress. Strongly interacting neutrinos with vanishing anisotropic stress are ruled out by Planck data with huge significance. The best fit occurs for standard decoupled neutrinos, and the limits found by [58] on the self-interaction cross-section are very stringent.

Sterile neutrinos
There are lots of good reasons to postulate the existence of sterile neutrinos: to explain the origin of neutrino masses, the existence of neutrino oscillations, or the mechanism responsible for baryogenesis. Sterile neutrinos could also play the role of dark matter. However, the kind of sterile neutrinos that can be probed with CMB and LSS observations are light sterile neutrinos, with a mass in the eV range. These sterile neutrinos are not motivated by fundamental physics issues, but rather by a few anomalies in short baseline neutrino oscillation data (LSND, MiniBooNE, reactor experiments).
Sterile neutrinos can be generated through various mechanisms, including resonant or non-resonant oscillations with the flavour neutrinos produced in the thermal bath. Depending on the active-sterile neutrino mixing angle(s), sterile neutrinos could thermalise or not [59]. When sterile neutrinos are in the relativistic regime, they can be entirely described in terms of an enhanced N eff (close to four with one generation of thermalised sterile neutrinos). When they become non-relativistic, the parameterisation of the model is not trivial anymore: in principle, one should specify the phase-space distribution of the sterile neutrino, and vary as many mass parameters as there are neutrino species (3+1 in the simplest scenarios). Since sterile neutrinos have a significant mass, different models with the same (N eff , M ν ) parameters but with different mass splittings and/or a different phase space distribution for the relic sterile neutrino can have distinct signatures on the CMB.
The Planck collaboration has released results for the simple case of two active neutrinos with negligible mass, plus one active neutrino with m = 0.06 eV, plus finally one sterile neutrino with a mass m s and a phase-space distribution equal to the one of active neutrinos multiplied by a suppression factor χ s . This corresponds to the so-called Dodelson-Widrow (DW) scenario [60], in which sterile neutrinos are generated through non-resonant oscillations. This case can be remaped into the one of thermally distributed sterile neutrinos with a smaller temperature than active neutrinos. The analysis could be carried out in terms of two free parameters (N eff , m s ), but it is more convenient to use (N eff , ω s ), where ω s ≡ Ω s h 2 is the relic density of sterile neutrinos. Indeed, ω s is the parameter really probed by the CMB through the lensing and shift-in-the-peak effects. The Planck paper reports constraints on the quantity 94.1 ω s eV, called the effective sterile neutrino mass, because it would coincide with the real mass if the sterile neutrino distribution was the same as for active neutrinos in the instantaneous decoupling limit. Ref. [41] provides joint constraints on (N eff , 94.1 ω s eV) from CMB+BAO data, showing that a thermalised neutrino with the same temperature as active neutrinos could have a mass of at most half an electron-volt, while a DW neutrino with 1 eV mass should have χ s < 0.5. The results show no evidence at all for sterile neutrinos.
In the previous sections, we have seen that CMB anisotropy data alone prefers a vanishing neutrino mass, and is well compatible with the standard prediction N eff = 3.046. Instead, HST data pushes for extra radiation, while lensing extraction and SZ clusters push for a non-zero neutrino mass. When all these data are considered at the same time, one gets marginal evidence for both N eff > 3 and M ν > 0, which could be interpreted in terms of light sterile neutrinos. Refs. [56,61,62] performed a joint analysis of all these data sets for a model with three active neutrinos of negligible mass, and one massive DW sterile neutrino. They find marginal evidence for such a sterile neutrino, with a mass in the range m s = 0.59 ± 0.38 eV (95% C.L.) and with χ s = 0.61 ± 0.60 (95% C.L.). These results are in moderate tension with the 3+1 scenario that could explain the short baseline anomaly (see also [63,64]). Note that these results do not only apply to sterile neutrinos and can be easily transposed to the case of other light massive relics, such as thermalised axions (see e.g. [56,65,66]).

Sensitivity of future experiments to neutrino parameters
The next releases of the Planck CMB satellite will lead to better bounds on neutrino masses. The forecast presented in [46] predicts a neutrino mass sensitivity of σ(M ν ) ∼ 0.1 eV from Planck alone, using full temperature and polarisation data, and lensing extraction. However, if the small internal tension between temperature and lensing found in the first release survives in the next ones, the final error bar will remain larger.
Several galaxy surveys with better sensitivity and larger volume are about to release data or have been planned over the next decades, including the Baryon Oscillation Spectroscopic Survey (BOSS), the Dark Energy Survey (DES), the Large Synoptic Survey Telescope (LSST) or the Euclid satellite. Concerning cosmic shear surveys, spectacular improvements are expected from Pan-STARRS, or the DES, LSST and the Euclid surveys already mentioned above.
In a near future, the prediction of ref. [67] is that the combination of full Planck data with BAO scale information from the full BOSS survey could lower the error down to σ(M ν ) ∼ 0.06 eV. In addition, the authors of [68] find that adding Lyman alpha data from BOSS should lead to comparable sensitivities, and even better results might be expected from the addition of galaxy power spectrum data from the same survey.
In ref. [69] it was found that the measurement of the galaxy harmonic power spectrum in seven redshift bins by DES should lead to a sensitivity of σ(M ν ) ∼ 0.06 eV when combined with Planck data (without lensing extraction). Similar bounds were derived in [70] for another combination of comparable experiments. This shows that at the horizon of 2015, a total neutrino mass close to M ν 0.1 eV could be marginally detected at the 2-σ level by cosmological observations. Because this value coincides with the lowest possible total mass in the inverted hierarchy scenario, the latter could start to be marginally ruled out in case the data still prefers M ν = 0.
Many papers studied the sensitivity of Euclid to the total neutrino mass. The conservative analysis of [71], taking into account a theoretical error related to the difficulty to model non-linear effects on small scales, suggests an error of the order of σ(M ν ) ∼ 0.03 eV in combination with Planck data. Constraints based on the groundbased LSST should be slightly weaker [72]. Hence, in the early 2020's, we expect that a combination of cosmological data sets could detect the total neutrino mass of the normal hierarchy scenario, M ν 0.05 eV, at the 2-σ level. If the total mass is instead close to M ν 0.1 eV, it will be detected at the 4-σ level. However, in that case, available experiments would not have enough sensitivity to distinguish between an inverted and normal hierarchy scenario with the same M ν .
Even more progress could be provided by the promising technique of 21-cm surveys. Instead of mapping the distribution of hydrogen atoms trough the absorption rate of photons traveling from quasars, it should be possible to observe directly the photons emitted by these atoms at a wavelength λ 21 cm from the transition from one hyperfine level to the other. While travelling towards the observer, these photons are redshifted, and seen with a wavelength indicating the position of the emitting atoms in redshift space. Recent theoretical progresses in this field show that using this technique, future dedicated experiments should be able to map hydrogen and hence baryonic fluctuations at very high redshift (typically 6 < z < 12), and to probe the matter power spectrum deep inside the matter-dominated regime on linear scales [73,74]. This field is still in its infancy, and the forecasts presented so far have to be taken with care, due to the difficulty to make a realistic estimate of systematic errors in future data sets. A sensitivity of σ(M ν ) ∼ 0.075 eV for the combination of Planck with the Square Kilometer Array (SKA) project, or σ(M ν ) ∼ 0.0075 eV with the Fast Fourier Transform Telescope (FFTT), was found in [38]. However, the authors show that such impressive experiments would still fail in discriminating between the NH and IH scenario.
An eventual post-Planck CMB satellite or post-Euclid survey would also have a great potential. The forecast analysis in [75] shows that for a CMB satellite of next generation one could get σ(M ν ) ∼ 0.03 eV alone, thanks to a very precise reconstruction of the CMB lensing potential, while [76] discusses the potential of cluster surveys. Finally, the authors of [37] show how far the characteristics of an hypothetical galaxy or cosmic shear survey should be pushed in order to discriminate between two allowed NH and IH scenarios with the same total mass. In the future, the sensitivity of CMB and LSS experiments to N eff will also increase significantly. Ref. [46] find that using lensing extraction, the full Planck can lower the error bar down to σ(N eff ) ∼ 0.3. Refs. [77,78] find that the combination of Planck data with the Euclid galaxy survey or cosmic shear survey would give σ(N eff ) ∼ 0.1. However, Ref. [79] claims that the combination of Euclid's galaxy survey, cosmic shear survey and cluster mass function measurement would give σ(N eff ) ∼ 0.02.

Conclusions
Neutrinos, despite the weakness of their interactions and their small masses, can play an important role in cosmology that we have reviewed in this contribution. In addition, cosmological data can be used to constrain neutrino properties, providing information on these elusive particles that complements the efforts of laboratory experiments. In particular, the data on cosmological observables have been used to bound the radiation content of the Universe via the effective number of neutrinos, including a potential extra contribution from other relativistic particles.
But probably the most important contribution of cosmology to our knowledge of neutrino properties is the information that provides on the absolute scale of neutrino masses. We have seen that the analysis of cosmological data can lead to either a bound or a measurement of the sum of neutrino masses, an important result complementary to terrestrial experiments such as tritium beta decay and neutrinoless double beta decay experiments. In the next future, thanks to the data from new cosmological experiments we could even hope to test the minimal values of neutrino masses guaranteed by the present evidences for flavour neutrino oscillations. For this and many other reasons, we expect that neutrino cosmology will remain an active research field in the next years.