#### Abstract

Neutrinos can play an important role in the evolution of the universe, modifying some of the cosmological observables. In this contribution we summarize the main aspects of cosmological relic neutrinos, and we describe how the precision of present cosmological data can be used to learn about neutrino properties, in particular their mass, providing complementary information to beta decay and neutrinoless double-beta decay experiments. We show how the analysis of current cosmological observations, such as the anisotropies of the cosmic microwave background or the distribution of large-scale structure, provides an upper bound on the sum of neutrino masses of order 1 eV or less, with very good perspectives from future cosmological measurements which are expected to be sensitive to neutrino masses well into the sub-eV range.

#### 1. Introduction

The subject of this contribution is the role of neutrinos in cosmology, 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 [1, 2].

We begin with a description of the properties and evolution of the background of relic neutrinos that fills the universe. Then we review the possible effects of neutrino oscillations on cosmology. The topic of neutrinos and Big Bang Nucleosynthesis is reviewed in a different contribution to this special issue [3]. The largest part of this contribution is devoted to the impact of massive neutrinos on cosmological observables that can be used to extract bounds on neutrino masses from present data. Finally we discuss the sensitivities on neutrino masses from future cosmological experiments.

Note that massive neutrinos could also play a role in the generation of the baryon asymmetry of the universe from a previously created lepton asymmetry. In these leptogenesis scenarios, one can also obtain quite restrictive bounds on light neutrino masses, which are, however, strongly model dependent. We do not discuss this subject here, it is covered in other contribution to this special issue [4].

For further details, the reader is referred to recent reviews on neutrino cosmology such as [5–7] and in particular [8]. A more general review on the connection between particle physics and cosmology can be found in [9].

#### 2. 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. In this section we will summarize the evolution and main properties of the CNB.

##### 2.1. Relic Neutrino Production and Decoupling

Produced at large temperatures by frequent weak interactions, cosmic neutrinos of any flavour () 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 , which is just one example of the general case of particles in equilibrium (fermions or bosons, relativistic or nonrelativistic), as shown for example, in [10]. In the previous equation we have included a neutrino chemical potential that would exist in the presence of a neutrino-antineutrino asymmetry, but we will see later in Section 4.1 that even if it exists its contribution cannot be very relevant.

As the universe cools, the weak interaction rate falls below the expansion rate and one says that neutrinos decouple from the rest of the plasma. An estimate of the decoupling temperature can be found by equating the thermally averaged value of the weak interaction rate where is the cross section of the electron-neutrino processes with the Fermi constant and is the neutrino number density, with the expansion rate given by the Hubble parameter Here is the total energy density, dominated by relativistic particles, and is the Planck mass. If we approximate the numerical factors to unity, with and , we obtain the rough estimate MeV. More accurate calculations give slightly higher values of which are flavour dependent because electron neutrinos and antineutrinos are in closer contact with electrons and positrons, as shown for example, in [1].

Although neutrino decoupling is not described by a unique , it can be approximated as an instantaneous process. The standard picture of * instantaneous neutrino decoupling* is very simple (see e.g., [10] or [11]) and reasonably accurate. In this approximation, the spectrum in (2.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 noninteracting neutrinos remains constant in a comoving volume since the decoupling epoch. We will see later that active neutrinos cannot possess masses much larger than 1 eV, so they were ultrarelativistic at decoupling. This is the reason why the momentum distribution in (2.1) does not depend on the neutrino masses, even after decoupling, that is, there is no neutrino energy in the exponential of .

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 the time of neutrino decoupling. In practice, the Dirac neutrinos with the “wrong-helicity” states do not interact with the plasma at temperatures of the MeV order 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 [1], but such a large mass is excluded for active neutrinos). 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 annihilations that heat the 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 annihilation and after the electron-positron pairs disappear by assuming entropy conservation of the electromagnetic plasma. The result is which is also the ratio between the temperatures of relic photons and neutrinos . The evolution of this ratio during the process of annihilations is shown in Figure 1(a), while one can see in Figure 1(b) how in this epoch the photon temperature decreases with the expansion less than the inverse of the scale factor . Instead the temperature of the decoupled neutrinos always falls as .

**(a)**

**(b)**

It turns out that the standard picture of neutrino decoupling described above is slightly modified: the processes of neutrino decoupling and annihilations are sufficiently close in time so that some relic interactions between 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]. A proper calculation of the process of noninstantaneous neutrino decoupling demands solving the momentum-dependent Boltzmann equations for the neutrino spectra, a set of integrodifferential kinetic equations that are difficult to solve numerically. This problem was considered in [12] including the effect of flavour neutrino oscillations on the neutrino decoupling process. One finds an increase in the neutrino energy densities with respect to the instantaneous decoupling approximation ( and for 's and 's, resp.) and a value of the comoving photon temperature after annihilations which is a factor larger, instead of . These changes modify the contribution of relativistic relic neutrinos to the total energy density which is taken into account using , as defined later in (3.1). In practice, the distortions calculated in [12] 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 (2.1) and . For instance, the number density per flavour is fixed by the temperature, which leads to a present value of neutrinos and antineutrinos of each flavour per . Instead, the energy density for massive neutrinos should in principle be calculated numerically, with two well-defined analytical limits,

##### 2.2. Background Evolution

Let us discuss the evolution of the CNB after decoupling in the expanding universe, which is described by the Friedmann-Robertson-Walker metric [11] where we assumed negligible spatial curvature. Here is the scale factor usually normalized to unity now () and related to the redshift as . General relativity tells us the relation between the metric and the matter and energy in the universe 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 is defined as , and the current value of the Hubble parameter gives the critical density today where . The different contributions to the total energy density are and the evolution of each component is given by the energy conservation law in an expanding universe , where is the pressure. Thus the homogeneous density of photons scales like , that of nonrelativistic matter ( for cold dark matter and for baryons) like , 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 nonrelativistic transition.

The evolution of all densities is depicted in Figure 2, starting at MeV temperatures until now. The various density fractions are shown in this figure, where it is easy to see which of the universe components is dominant, fixing its expansion rate: first radiation in the form of photons and neutrinos (radiation domination or RD), then matter which can be CDM, baryons, and massive neutrinos at late times (matter domination or MD), and finally the cosmological constant density takes over at low redshift (typically ).

Massive neutrinos are the only particles that present the transition from radiation to matter, when their density is clearly enhanced (upper solid lines in Figure 2). 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 ), and the present value could be of order unity for eV masses (see Section 5).

Shortly after neutrino decoupling, the CNB plays an interesting role in Big Bang Nucleosynthesis (BBN), the period of the universe when the primordial abundances of light elements are created. This subject is described in the contribution [3] (for a recent review, see [13]). Here we just summarize the two main effects of relic neutrinos at BBN. The first one is that they contribute to the relativistic energy density of the universe (if ), thus fixing the expansion rate. This is why BBN gave the first allowed range of the number of neutrino species before accelerators (see the next section). On the other hand, BBN is the last cosmological epoch sensitive to neutrino flavour, because electron neutrinos and antineutrinos play a direct role in the weak processes. We will see some examples of BBN bounds on neutrinos (effective number or oscillations) in the following sections.

#### 3. Extra Radiation and the Effective Number of Neutrinos

Together with photons, in the standard case neutrinos fix the expansion rate during the cosmological era when the universe is dominated by radiation. Their contribution to the total radiation content can be parametrized in terms of the effective number of neutrinos , through the relation where we have normalized 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 -boson width at LEP, [14], and we saw in a previous section from the analysis of neutrino decoupling that they contribute as . Any departure of from this last value would be due to nonstandard 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 , leading to a departure of from the standard value of . A detailed discussion of cosmological scenarios where is not fixed to three can be found in [1, 2, 15].

The expansion rate during BBN fixes the produced abundances of light elements and in particular that of . Thus the value of can be constrained at the BBN epoch from the comparison of theoretical predictions and experimental data on the primordial abundances of light elements [3, 13]. In addition, a value of different from the standard one would modify the transition epoch from a radiation dominated to a matter dominated universe, which has some consequences on some cosmological observables such as the power spectrum of CMB anisotropies, leading to independent bounds on the radiation content. These are two complementary ways of constraining at very different epochs.

A recent analysis of the BBN constraints on can be found in [16] (see the references therein for a list of recent works). The authors have discussed a new and more conservative approach to derive BBN constraints on , motivated by growing concerns on the reliability of astrophysical determinations of primordial . According to [16], the extra radiation at the BBN epoch is limited to at 95% C.L. On the other hand, recent analyses of late cosmological observables seem to favor , with best-fit values of order , although with large errorbars as shown for instance in [17]. The considered data include CMB temperature anisotropies and polarization, combined with other data such as measurements of the present value of the Hubble parameter (), the power spectrum of Luminous Red Galaxies (LRG), or distance measurements from the baryon acoustic oscillations (BAO) in the distribution of galaxies. The allowed regions at 95% C.L. from [17] are (WMAP + BAO + ) and (WMAP + LRG + ). Other recent analyses of bounds from cosmological data can be found in [18–20]. These ranges are in reasonable agreement with the standard prediction of , although they show a marginal preference for extra relativistic degrees of freedom whose robustness is still unclear. Moreover, they show that there exists an allowed region of values that is common at early (BBN) and more recent epochs, although with large error bars. The upcoming CMB measurements by the PLANCK satellite will soon pin down the radiation content of the universe, clarifying whether one really needs new physics leading to relativistic degrees of freedom beyond the contribution of flavour neutrinos.

#### 4. Neutrino Oscillations in the Early Universe

Nowadays there exist compelling evidences for flavour neutrino oscillations from a variety of experimental data on solar, atmospheric, reactor, and accelerator neutrinos. These are very important results, because the existence of flavour change implies that neutrinos mix and have nonzero masses, which in turn requires particle physics beyond the standard model. Thus it is interesting to check whether neutrino oscillations can modify any of the cosmological observables. More on neutrino oscillations and their implications can be found in other contributions to this special issue or any of the existing reviews such as [21–23], to which we refer the reader for more details.

It turns out that in the standard cosmological picture all flavour neutrinos were produced with the same energy spectrum, as we saw in Section 2.1, so we do not expect any effect from the oscillations among these three states. This is true up to the small spectral distortion caused by the heating of neutrinos from annihilations [12], as described before. In this section we will briefly consider two cases where neutrino oscillations could have cosmological consequences: flavour oscillations with nonzero relic neutrino asymmetries and active-sterile neutrino oscillations.

##### 4.1. Active-Active Neutrino Oscillations: Relic Neutrino Asymmetries

A nonzero relic neutrino asymmetry exists when the number densities of neutrinos and antineutrinos of a given flavour are different and can be parameterized by the ratio . Such a putative asymmetry, that could have been produced by some mechanism well before the thermal decoupling epoch, is expected to be of the same order of the cosmological baryon number , that is, a few times , from the equilibration of lepton and baryon numbers by sphalerons in the very early universe. In such a case, cosmological neutrino asymmetries would be too small to have any observable consequence. However, values for these parameters which are orders of magnitude larger than are not excluded by observations. Actually, large are predicted in theoretical models where the generation of lepton asymmetry took place after the electroweak phase transition or the electroweak washing out of preexisting asymmetries is not effective.

Neutrino asymmetries can be quantified assuming that a given flavour is characterized by a Fermi-Dirac distribution as in (2.1) with a nonzero chemical potential or equivalently with the dimensionless degeneracy parameter (for antineutrinos, ). In such a case, sometimes one says that the relic neutrinos are degenerate (but not in the sense of equal masses). Degenerate electron neutrinos have a direct effect on BBN: any change in the spectra modifies the primordial neutron-to-proton ratio, which in this case is . Therefore, a positive decreases the primordial He mass fraction, while a negative increases it [24], leading to an allowed range compatible with and very restrictive for negative values. In addition a nonzero relic neutrino asymmetry always enhances the contribution of the CNB to the relativistic energy density, because for any one has a departure from the standard value of the effective number of neutrinos given, if the neutrino spectra follow an equilibrium form, by We have seen that this increased radiation modifies the outcome of BBN and that bounds on can be obtained. In addition, another consequence of the extra radiation density is that it postpones the epoch of matter-radiation equality, producing observable effects on the spectrum of CMB anisotropies and the distribution of cosmic large-scale structures (LSS). Both independent bounds on the radiation content can be translated into flavour-independent limits on .

Altogether these cosmological limits on the neutrino chemical potentials or relic neutrino asymmetries are not very restrictive, because at least for BBN their effect in the or sector can be compensated by a positive . For example, an analysis of the combined effect of a nonzero neutrino asymmetry on BBN and CMB/LSS yields the allowed regions [25] in agreement with similar but more updated bounds as cited in [8]. These limits allow for a very significant radiation contribution of degenerate neutrinos, leading many authors to discuss the implications of a large neutrino asymmetry in different physical situations (see e.g., [1]).

It is obvious that the limits in (4.3) would be modified if neutrino flavour oscillations were effective before BBN, changing the initial distribution of flavour asymmetries. Actually, it was shown in [26–28] that this is the case for the neutrino mixing parameters in the region favoured by present data. This result is obtained only after the proper inclusion of the refractive terms produced by the background neutrinos, which synchronize the oscillations of neutrinos with different momenta (which would evolve independently without them). If flavour equilibrium is reached before BBN and the momentum distributions of neutrinos keep a Fermi-Dirac form as in (2.1), the restrictive limits on in (4.1) apply to all flavours. The bounds on the common value of the neutrino degeneracy parameter would be at [29].

More recent analyses [30, 31] have shown that this conclusion does not always hold. For the present measured values of neutrino mixing parameters, the degree of flavour equilibration depends on the value of the mixing angle , which fixes the onset of flavour oscillations involving 's. This in turn determines whether neutrinos interact enough with electrons and positrons to transfer the excess of energy density due to the initial to the electromagnetic plasma.

In [31] the BBN bounds on the total neutrino asymmetry (cosmological lepton number) were found for a range of initial flavour neutrino asymmetries. An example of this analysis is shown in Figure 3, taken from [32]. From this plot one can easily see the effect of flavour oscillations on the BBN constraints on the total neutrino asymmetry. With no neutrino mixing the value of is severely constrained by data, while the asymmetry for other neutrino flavours could be much larger. On the other hand, flavour oscillations imply that an initially large can be compensated by an asymmetry in the other flavours with opposite sign. The most restrictive BBN bound on applies then to the total asymmetry, an effect that can be seen graphically in Figure 3 as a * rotation* of the allowed region from a quasi-horizontal one for zero mixing to an almost vertical region for , in particular for the inverted mass hierarchy. In all cases depicted in Figure 3 the allowed values of the asymmetries are fixed by both deuterium and , the latter imposing that the value of at BBN must be very close to zero. Finally, in [33] it has been shown that data still fix the bounds on when CMB results and other cosmological parameters such as neutrino masses are included in the analysis, and only for future CMB data the bounds on the asymmetry could be improved. At the same time the contribution of neutrino asymmetries to an enhancement of radiation is limited to for values of in the region allowed by oscillation data.

##### 4.2. Active-Sterile Neutrino Oscillations

In addition to the flavour or active neutrinos (three species as we saw from accelerator data), there could also exist extra massive neutrino states that are sterile, that is, singlets of the Standard Model gauge group and thus insensitive to weak interactions. Most of the current data on neutrino oscillations can be perfectly explained with only the three active species, but there exist a few experimental results, sometimes called anomalies, that cannot be explained in this framework. If neutrino oscillations are responsible for all the experimental data, a solution might require additional, sterile, neutrino species. These kinds of particles are predicted by many theoretical models beyond the SM, being neutral leptons insensitive to weak interactions whose only interaction is gravitational. Their masses are usually heavy, while lighter sterile neutrinos are rarer but possible. Here we will briefly summarize these anomalies observed in neutrino experiments, whereas an updated review of this subject is given in [34].

The long-standing evidence (more than ) for oscillations comes from the LSND experiment. Its results pointed out a , much larger than those required from a three-neutrino analysis, as shown later in (5.3). At the same time, the KARMEN experiment, very similar but not identical to LSND, provided no support for such evidence, while a recent data release of the MiniBooNE experiment, designed to check the LSND results with larger distance and energy but similar ratio , is still consistent with the LSND signal. In addition, an unexplained excess of electron-like events is observed in MiniBooNE at low energies. The simultaneous interpretation of LSND (antineutrino) and MiniBooNE (neutrino and antineutrino) results in terms of sterile neutrino oscillations requires CP-violation or some other exotic scenarios, as reviewed in [34].

A new anomaly supporting oscillations with sterile neutrinos appeared from a revaluation of reactor antineutrino fluxes, which found a increase relative to previous flux calculations. As a result, data from reactor neutrino experiments at very short distances can be interpreted as an apparent deficit of . This is known as the reactor antineutrino anomaly and is again compatible with sterile neutrinos having a . Finally, an independent experimental evidence for disappearance at very short baselines exists from the Gallium radioactive source experiments GALLEX and SAGE.

The existence of all these experimental hints for sterile neutrinos and a mass scale at the eV is intriguing, but so far a fully consistent picture has not emerged. Many analyses have been performed trying to explain all experimental data with 1 or 2 additional sterile neutrinos, known as the or schemes, with the corresponding additional mixing parameters. It seems that none of these schemes does describe well all data, as explained in detail in [34], but for the topics of this paper the potential existence of oscillations into sterile neutrinos would lead to important cosmological consequences, such as extra radiation from fully or partly thermalized sterile neutrinos or a larger Hot or Warm Dark Matter component. The required values of neutrino masses in these or scenarios are, as we will see in Section 7, in tension with the current cosmological bounds.

In any case, it is interesting to consider the main effects of additional sterile neutrino species in cosmology. Because sterile neutrinos are insensitive to weak interactions, they do not follow the behaviour of active neutrinos and are not expected to be present in the early universe at MeV temperatures. Even if they could interact through other kinds of interactions, significantly weaker than the standard weak ones, as predicted by extended particle physics models, they would have a thermal spectrum at very high temperatures, but their density would have been strongly diluted by many subsequent particle-antiparticle annihilations. Therefore, barring the nonthermal production from additional physics beyond the SM, the main way of obtaining a significant abundance of sterile neutrinos is through their mixing with the active ones.

In principle, the cosmological evolution of the active-sterile neutrino system should be found solving the corresponding Boltzmann kinetic equations for the density matrices. There exists a vast literature on this subject, where different analyses considered several approximations. For a list of references, see for instance [1, 2]. Although in general one should consider, at least, a mixing of three active neutrinos and one sterile species (with 4 masses, 6 mixing angles, and 3 CP-violating phases), let us first assume an admixture of one sterile state to electron neutrinos. In the early universe one expects that neutrino oscillations could be effective when the vacuum oscillation term becomes larger than the main matter potential term from charged leptons at a temperature . The evolution of the active-sterile neutrino system depends on the sign of . For negative values there could be resonant oscillations (or resonant production, RP). Instead, for one has the so-called nonresonant production of sterile neutrinos (NRP). Comparing the value of with the decoupling temperature of active neutrinos MeV, one obtains that for values: active-sterile oscillations are effective after neutrino decoupling. In such a case, the total comoving number density of active and sterile neutrinos will be constant because active neutrinos are no longer interacting with the rest of the primeval plasma. Correspondingly, distortions in the momentum spectra of neutrinos are expected which directly affect the production of at BBN when the active neutrinos are of the electron flavour (see e.g., [35], also for the case of nonzero initial abundance). Instead, for much larger values of oscillations are effective before neutrino decoupling, when weak interactions are frequent. In this case the actual growth of the sterile neutrino population depends on the interplay between oscillations and interactions, while the energy spectrum of the active ones will be kept in equilibrium. Their combined contribution to radiation can be as large as , depending on the specific value of the mixing angle and on the sign of . This “thermalization” of the sterile neutrinos is a well-known phenomenon that is very difficult to avoid unless the cosmological scenario is drastically modified.

In the NRP case, it was shown in [36] that the production probability of sterile neutrinos is where is the mixing angle in matter and the frequency of oscillations in the medium. Here is the production rate of active neutrinos in the plasma, and the averaging is made over the thermal background. If the oscillation frequency is very high one can substitute , and if is small one obtains for a small number density of and active neutrinos close to equilibrium. Therefore, a thermal or close to thermal population of sterile neutrinos is expected provided that such a production rate of is larger than the expansion rate of the universe, a condition that holds unless either or is very small.

For RP of sterile neutrinos at temperatures below the resonance propagates from small to large values of the neutrino comoving momentum, covering the whole momentum distribution while the active neutrinos are repopulated by interactions. The thermalization of is thus significantly enhanced, even for quite small values of the mixing angle.

In order to illustrate this discussion with an actual calculation of the active-sterile system with the kinetic equations in the two-flavour approximation, among the many published analyses we have chosen a recent one [37]. Their results for the final extra contribution of sterile neutrinos to radiation, in the case of zero initial lepton asymmetry, are shown as isocontours of in Figure 4 as a function of the mixing parameters (in units) and . In the NRP case (left panel) one can clearly see that the same final corresponds to constant values of . For RP 's are more efficiently brought into equilibrium, even for quite small values of the mixing angle. All calculations described so far correspond to the two-neutrino limit of one active and one sterile states, but a proper calculation should also include the unavoidable presence of mixing among the active neutrinos. A full four-flavour calculation has not been performed, but a few analyses did solve simplified kinetic equations taking into account active neutrino mixing, such as [38–40].

**(a)**

**(b)**

We note that for the active-sterile parameters needed to solve the oscillation anomalies described at the beginning of this section, the thermalization of sterile neutrinos is achieved, that is, . An example of the best-fit values of a particular calculation in the neutrino model is indicated in the plot. Therefore, it seems that such an extra radiation is guaranteed in these situations unless oscillations are suppressed, as in the case of a nonzero initial lepton asymmetry much larger than [42, 43].

#### 5. Massive Neutrinos as Dark Matter

Nowadays the existence of Dark Matter (DM), the dominant nonbaryonic 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. Together with CMB photons, relic neutrinos can be found anywhere in the universe with a number density given by the present value of (2.5) of 339 neutrinos and antineutrinos per cm^{3}, and their energy density in units of the critical value of the energy density (see (2.9)) is
Here includes all masses of the neutrino states which are nonrelativistic today. It is also useful to define the neutrino density fraction with respect to the total matter density

In order to check whether relic neutrinos can have a contribution of order unity to the present values of or , we should consider which neutrino masses are allowed by noncosmological data. Oscillation experiments measure the differences of squared neutrino masses and , the relevant ones for solar, and atmospheric neutrinos, respectively. As a reference, we take the following ranges of mixing parameters from [44] (see also [45, 46]): Here , where ( or ) are the three mixing angles. Unfortunately oscillation experiments are insensitive to the absolute scale of neutrino masses, because the knowledge of and leads to the two possible schemes shown in figure 1 of [8], but leaves one neutrino mass unconstrained. These two schemes are known as normal (NH) and inverted (IH) mass hierarchies, characterized by the sign of , positive and negative, respectively. In the above equation the values in parentheses correspond to the IH, otherwise the mixing parameters present the same allowed regions for both hierarchies. For small values of the lightest neutrino mass , that is, () 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 and then we say that they are degenerate. In general, the relation between the individual masses and the total neutrino mass can be found numerically, as shown in Figure 5.

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 (in short ) is a rare nuclear processes where lepton number is violated and whose observation would mean that neutrinos are Majorana particles. If the 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 where are the two Majorana phases that appear in lepton-number-violating processes. An important issue for results is related to the uncertainties on the corresponding nuclear matrix elements. For more details and the current experimental results, see [47].

Beta decay experiments, which involve only the kinematics of electrons, are in principle the best strategy for measuring directly the neutrino mass [48]. The current limits from tritium beta decay apply only to the range of degenerate neutrino masses, so that , where is the relevant parameter for beta decay experiments. The bound at CL is 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 eV at 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 is restricted to the approximate range

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 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 in [21, 49, 50]. The two cases involving are shown in Figure 6.

Now we can find the possible present values of in agreement the approximate bounds of (5.6). Note that even if the three neutrinos are nondegenerate in mass, (5.1) can be safely applied, because we know from neutrino oscillation data that at least two of the neutrino states are nonrelativistic today, because both eV and eV are larger than the temperature eV. If the third neutrino state is very light and still relativistic, its relative contribution to is negligible and (5.1) remains an excellent approximation of the total density. One finds that is restricted to the approximate range where we already included that . 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 (), we obtain an upper bound eV (again fixing ). Moreover, because from present analyses of cosmological data we know that the approximate contribution of matter is , the neutrino masses should obey the stronger bound eV. We see that 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 .

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, and the reader can find a historical review in [51]. 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 free-streaming scale. In a universe dominated by HDM, this suppression is in contradiction with 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, that is, particles which were nonrelativistic 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 play the role of dark matter and 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 is currently under careful examination (see e.g., [52]).

#### 6. Effects of Neutrino Masses on Cosmology

In this section we will briefly describe the main cosmological observables and the effects that neutrino masses cause on them. A more detailed discussion of the effects of massive neutrinos on the evolution of cosmological perturbations can be found in Sections 4.5 and 4.6 of [8].

##### 6.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 ( being a direction in the sky). This function is usually expanded in Legendre multipoles
where are the Legendre polynomials. For Gaussian fluctuations, all the information is encoded in the multipoles which probe correlations on angular scales . We have seen that each neutrino state can only have a mass of the order of eV, so that the transition of relic neutrinos to the nonrelativistic regime is expected to take place after the time of recombination between electrons and nucleons, that is, after photon decoupling. Because the shape of the CMB spectrum is related mainly to the physical evolution* before* recombination, it will only be marginally affected by the neutrino mass, except through a modified background evolution and some secondary anisotropy corrections. 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, observed with various techniques described in the next section (directly or indirectly, today or in the near past at redshift ). It is defined as the two-point correlation function of nonrelativistic matter fluctuations in Fourier space where . Usually refers to the matter power spectrum evaluated today (at ). In the case of several fluids (e.g., CDM, baryons, and nonrelativistic neutrinos), the total matter perturbation can be expanded as Because the energy density is related to the mass density of nonrelativistic matter through , represents indifferently the energy or mass power spectrum. The shape of the matter power spectrum is affected in a scale-dependent way by the free-streaming caused by small neutrino masses of (eV), and thus it is the key observable for constraining with cosmological methods.

##### 6.2. 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 . It is possible to define a horizon as the typical distance on which particles travel between time and . When the universe was dominated by radiation or matter , this horizon is, as usual, asymptotically equal to , up to a numerical factor of order one. Similar to the definition of the Jeans length (see section 4.4 in [8]), we can define the neutrino free-streaming wavenumber and length as As long as neutrinos are relativistic, they travel at the speed of light and their free-streaming 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 and K. This gives for the free-streaming wavelength and wavenumber during matter or domination where and are the cosmological constant and matter density fractions, respectively, evaluated today. After the nonrelativistic transition and during matter domination, the free-streaming length continues to increase, but only like , that is, more slowly than the scale factor . Therefore, the comoving free-streaming length actually decreases like . As a consequence, for neutrinos becoming nonrelativistic during matter domination, the comoving free-streaming wavenumber passes through a minimum at the time of the transition, that is, when and . 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 free-streaming 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 nonrelativistic transition the neutrino perturbations behave like CDM perturbations. In particular, modes with are never affected by free streaming and evolve like being in a pure CDM model.

##### 6.3. 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 nonlinearly, 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 nonlinear N-body simulations. We will not review here all the details of this complicated evolution (see [2, 7, 8] 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 (that is, on wave numbers smaller than the value 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 can be shown to depend only on the matter density fraction today (including neutrinos), , and on the primordial perturbation spectrum. If the neutrino mass is varied with fixed, the large-scale power spectrum remains invariant.

On small scales such that , the matter power spectrum is affected by neutrino masses for essentially the three following reasons.

(1) Massive neutrinos do not cluster on those scales. The matter power spectrum can be expanded as a function of the three nonrelativistic species,

On scales of interest and in the recent universe, baryon and CDM fluctuations are equal to each other, while . Hence, even if the evolution of was not affected by neutrino masses (which is not the case, see the remaining two points below), the power spectrum would be reduced by a factor with

(2) The redshift of radiation-to-matter equality or the baryon-to-CDM ratio might be slightly affected by neutrino masses, with a potential impact on the small-scale matter power spectrum. This depends very much on which other parameters are kept fixed when the neutrino mass varies. If neutrino masses are smaller than roughly 0.5 eV, they are still relativistic at the time of radiation-to-matter equality, and the redshift of equality depends on , not on . It is possible to increase and with fixed parameters , , (provided that the Hubble parameter also increases like the square root of ). In that case, there is no significant impact of massive neutrinos on the matter power spectrum through background effects, that is, through a change in the homogeneous cosmological evolution. However, there are some important perturbation effects that we will now summarize.

(3) The growth rate of cold dark matter perturbations is reduced through an absence of gravitational back-reaction effects from free-streaming neutrinos. This growth rate is set by an equation of the type where stands for the CDM relative density perturbation in Fourier space and for the metric perturbation playing the role of the Newtonian potential inside the Hubble radius. There is a similar equation for decoupled baryons, and very soon after baryon decoupling we can identify on scales of interest. On the right-hand side, we neglected time derivatives of metric fluctuations playing only a minor role. The right-hand side represents gravitational clustering and is given by the Poisson equation as a function of the total density fluctuation. The second term on the left-hand side represents Hubble friction, that is, the fact that the cosmological expansion enhances distances, reduces gravitational forces, and slows down gravitational clustering. The coefficient is given by the Friedmann equation as a function of the total background density. In a universe such that all species present in the Friedmann equation do cluster, as it is the case in a matter-dominated universe with and , the solution is simply given by : the so-called linear growth factor is proportional to the scale factor. But whenever one of the species contributing to the background expansion does not cluster efficiently, it can be neglected in the Poisson equation. In that case, the term on the right-hand side becomes smaller with respect to the Hubble friction term, and CDM (as well as baryons) clusters at a slower rate. This is the case in presence of massive neutrinos and for : the linear growth rate during matter domination is then equal to . During domination, this effect sums up with that of the cosmological constant (or of any nonclustering dark energy), which also tends to reduce the growth rate for the very same reason.

In summary, the small-scale matter power spectrum is reduced in presence of massive neutrinos for at least two reason: 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 , the step-like suppression of starts at and saturates at Mpc with a constant amplitude . This result was obtained by fitting numerical simulations [53], but a more accurate approximation can be derived analytically [2, 8]. As mentioned in the second item above, neutrino masses can have additional indirect effects through a change in the background evolution, depending on which cosmological parameters are kept fixed when varies.

When fitting data, one can use analytical approximations to the full MDM or MDM matter power spectrum, valid for arbitrary scales and redshifts, as listed in [8]. However, nowadays the analyses are performed using the matter power spectra calculated by Boltzmann codes such as CAMB [54] or CLASS [55] that solve numerically the evolution of the cosmological perturbations. The step-like suppression of the matter power spectrum induced by various values of is shown in Figure 7.

Is it possible to mimic the effect of massive neutrinos on the matter power spectrum with some combination of other cosmological parameters? If so, one would say that a parameter degeneracy exists, reducing the sensitivity to neutrino masses. This possibility depends on the intervals and in which can be accurately measured. Ideally, if we could have and , the effect of the neutrino mass would be nondegenerate, because of its very characteristic step-like effect. Moreover, because neutrinos render the linear growth factor scale dependent, with proportional to (resp., or ) for (resp., , the amplitude of the step-like suppression is redshift dependent. Using for example, weak lensing techniques or Lyman- forest data, one could get accurate measurements of the matter spectrum at several redshifts in the range . This will offer an opportunity to test the “time dependence” of the neutrino mass effect, and to distinguish it, for instance, from an equivalent step-like suppression in the primordial spectrum that would still imply a scale-independent growth factor.

The problem is that usually the matter power spectrum can only be accurately measured in the intermediate region where the mass effect is neither null nor maximal: in other words, many experiments only have access to the transition region in the step-like transfer function. In this region, the neutrino mass affects the slope of the matter power spectrum in a way which can be easily confused with the effect of other cosmological parameters. Because of these parameter degeneracies, the current LSS data alone cannot provide significant constraints on the neutrino mass, and it is necessary to combine them with other cosmological data, in particular the CMB anisotropy spectrum, which could lift most of the degeneracies. Still, for exotic models with for example, extrarelativistic degrees of freedom, a constant equation-of-state parameter of the dark energy different from , or a nonpower-law primordial spectrum, the neutrino mass bound can become significantly weaker.

The value of is not limited by observational sensitivities but by the range in which we trust predictions from linear theory. Beyond the scale of nonlinearity (of the order of Mpc at ), the data is only useful provided that one is able to make accurate predictions not only for the nonlinear power spectrum, but also for redshift-space distortions (coming from the fact that we observe the redshift of objects, not their distance away from us) and finally for the scale dependence of the light-to-mass bias (relating the power spectrum of a given category of observed compact objects to the underlying total matter power spectrum). Spectacular progresses are being carried on these three fronts, thanks to better N-body simulations and analytical techniques. Including massive neutrinos in such calculations is particularly difficult, but successful attempts were presented for example, in [56–58]. The neutrino mass impact on the nonlinear matter power spectrum is now modeled with rather good precision, at least within one decade above in wave-number space. It appears that the step-like suppression is enhanced by nonlinear effects up to roughly (at redshift zero and Mpc) and is reduced above this scale. Hence, nonlinear corrections render the neutrino mass effect even more characteristic than in the linear theory and may help to increase the sensitivity of future experiments.

Until this point, we reduced the neutrino mass effect to that of the parameter or . In principle, the mass splitting between the three different families for a common total mass is visible in the matter power spectrum. The time at which each species becomes nonrelativistic depends on individual masses . 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 were computed numerically in [59] for the linear spectrum and [58] for the nonlinear 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 [60, 61].

##### 6.4. Impact of Massive Neutrinos on the CMB Anisotropy Spectrum

For neutrino masses of the order of eV (about ) the three neutrino species are still relativistic at the time of photon decoupling, and the direct effect of free-streaming neutrinos on the evolution of the baryon-photon acoustic oscillations is the same in the CDM and MDM cases. Therefore, the effect of the mass can only appear at two levels that of the background evolution and that of secondary CMB anisotropies, related to the behavior of photon perturbations after decoupling. Both levels are potentially affected by the evolution of neutrinos after the time of their nonrelativistic transition. If neutrinos were heavier than a few eV, they would already be nonrelativistic at decoupling, and they could trigger more direct effects in the CMB, as described in [62]. However, we will see later that this situation is disfavoured by current upper bounds on the neutrino mass.

Let us first review the background effects of massive neutrinos on the CMB. Because the temperature and polarization spectrum shape are 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 eV to eV should be counted as radiation at the time of equality and as nonrelativistic matter today: the total nonrelativistic density, parametrized by , depends on the total neutrino mass . Hence, when is varied, there must be a variation either in the redshift of matter-to-radiation equality or in the matter density today .

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 nonrelativistic matter density at late times can impact both the angular diameter distance to the last scattering surface , controlling the overall position of CMB spectrum features in multipole space, and the slope of the low- tail of the CMB spectrum, due to the late Integrated Sachs-Wolfe (ISW) effect. Out of these three effects (changes in , in 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 .

There exists another effect of massive neutrinos on the CMB at the level of secondary anisotropies: when neutrinos become nonrelativistic, 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 on multipoles , with a dependence of the maximum depletion scale on individual masses . This effect is roughly ten times smaller than the depletion in the small-scale matter power spectrum, .

We show in Figure 8 the effect on the CMB temperature spectrum of increasing the neutrino mass while keeping and fixed: the only observed differences are then for (late ISW effect due to neutrino background evolution) and for (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 (see Section 8).

#### 7. Current Bounds on Neutrino Masses

In this section we review how the available cosmological data are used to get information on the absolute scale of neutrino masses, complementary to laboratory experiments. Note that the bounds in the next subsections are all based on the Bayesian inference method, and the upper bounds on the sum of neutrino masses are given at 95% C.L. after marginalization over all free cosmological parameters. We refer the reader to Section 5.1 of [8] for a detailed discussion on this statistical method, as well as for most of the references for the experimental data or parameter analysis.

##### 7.1. CMB Anisotropy Bounds

The experimental situation of the measurement of the CMB anisotropies is dominated by the seven-year release of WMAP data [17], which improved the already precise TT and TE angular power spectra of the previous releases and included a detection of the E-polarization self-correlation spectrum (EE). On similar or smaller angular scales than WMAP, we have results from experiments that are either ground based or balloon borne (ACT, SPT, etc.).

We saw in the previous section that the CMB spectrum has a small sensitivity to neutrino masses even when each mass is below 0.5 eV, and all nonrelativistic transitions take place after photon decoupling [63]. This sensitivity is due to background effects (mainly the late ISW effect if all other MDM parameters are left free) plus perturbation effects (mainly the early ISW effect). Therefore, it is possible to constrain neutrino masses using CMB experiments only, down to the level at which these small effects are masked by instrumental noise, by cosmic variance, or by parameter degeneracies in the case of some cosmological models beyond the minimal Mixed Dark matter framework. WMAP alone is able to set a limit eV depending on whether the dark energy component is a cosmological constant or not [17]. Because the neutrino mass effects in the CMB are visible mainly at , combining WMAP data with other CMB datasets (from e.g., ACT or SPT) does not improve this bound. On the contrary, adding more information to the background cosmological evolution is very useful, because it helps removing degeneracies between the various parameters: in that case, the CMB data can better probe neutrino masses through their background effects. For instance, the bound from WMAP combined with BAO scale measurements and a direct determination of by the Hubble space telescope Key Project [64] is significantly stronger: eV (95% CL) [17], while the combination of WMAP with a different determination at various redshift from early-type galaxy evolution gives eV (95% CL) [65]. These are important results, because they do not depend on the uncertainties from LSS data discussed next.

##### 7.2. Large-Scale Structure Observations

The matter power spectrum can be probed with various methods on different scales and redshifts. Let us review here the major techniques which have led so far to relevant neutrino mass bounds.

###### 7.2.1. Galaxy Power Spectrum

Galaxy maps (or, similarly, cluster maps) can be smoothed over small scales and Fourier transformed in order to provide a power spectrum. Relating such a spectrum to the total matter power spectrum is a tricky exercise, especially on small scales (corresponding to wave-numbers ), because it is difficult to make accurate predictions for nonlinear corrections, redshift-space distortions, and the light-to-mass bias. Current bounds are based on the analysis of linear scales only, for which the assumption of a scale-independent bias is well motivated. This bias is, however, left as a free parameter, in such way that galaxy spectrum data give indications on the shape but not on the global amplitude of . In the next subsection, we will report constraints from the halo power spectrum of Large Red Galaxies (denoted later as Gal-LRG), measured by [66], the spectrum of the MegaZ catalogue (Gal-MegaZ), used by [67], and the WiggleZ Dark Energy Survey (Gal-WiggleZ), used by [68]. The first two data sets are actually extracted from the same big survey, the Sloan Digital Sky Survey (SDSS).

For sufficiently deep galaxy surveys, it is possible to separate galaxies into redshift bins and compute different correlation functions at different redshifts. This technique, called tomography, can be very useful for constraining the scale-dependent growth factor induced by neutrino masses. In that case, the data can be reduced to a set of two-dimension power spectra in different shells, each of them related to in a narrow redshift range. Recently, such a tomographic analysis was used by [69] for constraining neutrinos, using galaxies from the Canada-France-Hawaii-Telescope Legacy Survey (CFHTLS), split in three redshift bins covering the ranges , and (this dataset will be denoted as Gal-CFHTLS).

###### 7.2.2. Cluster Mass Function

Instead of probing directly the matter power spectrum from the spatial distribution of objects, it is possible to constrain integrated quantities of the type , where stands for a given window function. One such quantity is related to the histogram of cluster masses. If the mass of a significant number of galaxy clusters within a given redshift bin is known, this histogram gives an estimate of the so-called cluster mass function, , with being the number of clusters of redshift approximately equal to , and with a mass in the range . This function is related to , the variance of the density in spheres enclosing a mass , itself derived from the convolution of the power spectrum with an appropriate window function. In the next subsection we will refer to bounds derived from cluster abundances probed by X-ray observations from the ROSAT survey, presented by [70] (denoted later as Clus-ROSAT) and by optical observations from the MaxBCG catalogue, presented by [71] (denoted later as Clus-MaxBCG).

###### 7.2.3. Galaxy Weak Lensing

The image of observed galaxies is distorted by gravitational lensing effects, caused by density fluctuations along the line of sight. One of these effects is called cosmic shear. It corresponds to the squeezing of an image in one direction in the sky and its stretching in the orthogonal direction. Because such distortions are coherent over the angular size of the lensing potential wells responsible for lensing, they tend to align slightly the apparent major axis of galaxies in a given patch of the sky. Hence the average cosmic shear in a given direction can be estimated statistically by averaging over major axis orientations. The analysis of a catalogue of images leads to a map of the lensing potential, itself related to the matter power spectrum through the Poisson equation. If the number of source galaxies is sufficient, it is possible to split the catalogue in several redshift bins and to obtain a three-dimensional reconstruction of the gravitational potential in our past-line cone and of at various redshifts. Cosmic shear tomography is particularly useful for measuring neutrino masses, because it can probe the scale-dependent growth factor induced by neutrino masses over an extended range of redshifts. Current cosmic shear surveys already allow to put bounds on neutrino parameters: we will refer later to results from the CFHTLS presented by [72] and denoted as WL-CFHTLS.

###### 7.2.4. Lyman Alpha Forests

The most luminous and distant compact objects that we can observe are quasars. Some of the photons emitted by quasars interact along the line-of-sight. In particular, a fraction of photons are absorbed at the Lyman alpha wavelength by hydrogen atoms located in the Interstellar Galactic Medium (IGM). The absorbed fraction in a given point of the photon trajectory is proportional to the local density of neutral hydrogen. Because photons are continuously redshifted, absorption in a given point is seen by the observer as a depletion of the spectrum at a given frequency. Hence, inside a limited range called the * Lyman alpha forest*, the frequency dependence of quasar spectra is a tracer of the spatial fluctuations of the hydrogen density along the line of sight. Lyman alpha forests in quasar spectra offer an opportunity to reconstruct the hydrogen density fluctuation along several line-of-sights in a given redshift range. After Fourier expanding each spectrum and averaging over many spectra, one gets an estimate of the so-called flux power spectrum that can be related to the total matter power spectrum . Unfortunately, the flux power spectrum does not probe linear scales but mildly nonlinear scales. In order to relate to , it is necessary to perform N-body simulations with a hydrodynamical treatment of baryons, accounting for the complicated thermodynamical evolution of the IGM (which depends on star formation). Also, a limitation of this technique comes from the fact that the emitted quasar spectra already have a nontrivial frequency dependence and that photons are affected by several other effects than Lyman alpha absorption along the line-of-sight. Nevertheless, a careful modeling of all relevant effects allows to obtain interesting constraints. The fact that Lyman alpha forests probe mildly nonlinear scales rather than strongly nonlinear ones is of course crucial for keeping systematic errors under control. Lyman alpha observations typically constrain the matter power spectrum in the wavenumber range /Mpc and in the redshift range . We mention below some neutrino mass bounds inferred from quasar spectra obtained by the SDSS and presented in [73], denoted as Ly--SDSS.

##### 7.3. Large-Scale Structure Bounds

Using LSS observations in combination with CMB data offers an opportunity to observe (or to bound) the step-like suppression of the matter power spectrum in presence of neutrino masses, as explained in Section 6.3 and illustrated in Figure 7. The use of CMB data is crucial in order to constrain parameters like the baryon density, the primordial spectrum amplitude, the tilt, and a combination of and . Without such constraints, there would be too much freedom in the matter power spectrum fitted to LSS data for identifying a smooth step-like suppression.

We summarize in Table 1 the main constraints on available at the time of writing, obtained from combinations of CMB plus homogeneous cosmology data, galaxy power spectrum data, and cluster abundance data. Current data set are far from reaching the sensitivity required to probe the mass splitting of the total mass between different species. The constraints mentioned below have been derived in the case of three degenerate neutrinos with mass , but they roughly apply to the total mass of any scenario. Also, the bounds shown in the first column of Table 1 have been obtained under the assumption of a minimal CDM model with massive neutrinos, featuring seven free parameters. More conservative bounds are sometimes derived for basic extensions of this model, with one or two more parameters. The constraints do not change significantly when assuming, for instance, a primordial spectrum with a running of the tilt or a significant contribution to the CMB of primordial gravitational waves [71]. Parameters known to be slightly degenerate with neutrino masses and leading to weaker bounds are , the equation-of-state parameter of a Dark Energy component (substituting the cosmological constant), and , the effective number of neutrinos discussed in Section 3. The degeneracy with , explained in [74], is illustrated by the last column in Table 1. We did not include in this table current limits from Lyman alpha forest data: this is a delicate matter for the reasons mentioned previously, and a careful modeling of all systematic effects leads to rather weak neutrino mass bounds. The conservative analysis of [73], based on Ly--SDSS, gives a bound eV (95% CL) from Ly--SDSS data alone.

In conclusion, the combination of available data sets consistently indicates that the total neutrino mass is below 0.3 eV at the 95% CL (0.5 eV if we allow for Dark Energy with arbitrary ). This means that the “degenerate scenario” in which all neutrinos share roughly the same mass is almost excluded. The data is about to probe the region in which masses are different from each other and are ranked according to the NH or IH scenario. Other recent summaries of existing bounds have been recently presented in [6, 17, 50, 71, 75].

#### 8. Future Sensitivities on Neutrino Masses from Cosmology

Future CMB observations will have increasing sensitivity to neutrino masses, not only thanks to smaller error bars on the temperature and polarization spectra. They will make it possible to probe the large-scale structure of the universe with a new technique: CMB lensing extraction. The weak lensing of the last scattering surface by nearby galaxy clusters induces specific non-Gaussian patterns in CMB maps that can be extracted using some nonlinear estimators. This method allows to measure the lensing potential up to and to infer constraints on at such high redshifts.

In the very close future, significant improvements on the neutrino mass bounds will be triggered by the Planck CMB satellite. The forecasts presented in [76] predict a neutrino mass sensitivity of eV from Planck alone, using the lensing extraction technique of [77]. This would be twice better than without lensing extraction.

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 (http://cosmology.lbl.gov/BOSS/) (BOSS), the Dark Energy Survey (http://www.darkenergysurvey.org/) (DES), the Large Synoptic Survey Telescope (http://www.lsst.org/) (LSST), or the Euclid satellite (http://sci.esa.int/euclid). Also, in [78] it was pointed out that in the future accurate measurements could be inferred from cluster surveys. Because clusters are more luminous than galaxies, they can be mapped up to higher redshift. Concerning cosmic shear surveys, spectacular improvements are expected from Pan-STARRS (http://pan-starrs.ifa.hawaii.edu/) or the DES, LSST, and the Euclid surveys already mentioned above.

In a near future, the prediction of [79] is that the combination of Planck (with lensing extraction) with BAO scale information from BOSS could lower the error down to eV. In addition, the authors of [80] 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.

With better tomographic data (for either galaxy clustering or cosmic shear), it will become possible to probe the scale dependence of the growth factor induced by neutrino masses (or in other words, the fact that the step-like suppression has an amplitude increasing with time) and to reach spectacular sensitivities. We present below the typical sensitivities expected for a collection of planned surveys (not all approved). These numbers should be taken with care because forecasts are based on an idealization of each experiment, as well as on several assumptions like the underlying cosmological model or even the fiducial value of the neutrino mass itself.

In [81] it was found that the measurement of the galaxy harmonic power spectrum in seven redshift bins by the DES should lead to a sensitivity of eV when combined with Planck data (without lensing extraction). Similar bounds were derived in [82] for another combination of comparable experiments. This shows that at the horizon of 2014 or 2015, a total neutrino mass close to 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 .

The sensitivity of cosmic shear data from a satellite experiment comparable to Euclid was calculated in [83], where it was found that it would shrink to eV in combination with Planck data (without lensing extraction). The forecast of [84] based on galaxy clustering data also from Euclid (completed at small redhsift by similar data from BOSS) gives comparable numbers. Constraints based on the ground-based Large Synoptic Survey Telescope should be slightly weaker [85]. 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, eV, at the 2- level. If the total mass is instead close to eV, it will be detected at the 4- level. However, in that case, available experiments would not have enough sensitivity for making the difference between an inverted and normal hierarchy scenario with the same .

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 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 ) and to probe the matter power spectrum deep inside the matter-dominated regime on linear scales [86]. 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 eV for the combination of Planck with the Square Kilometer Array (SKA) project, or eV with the Fast Fourier Transform Telescope (FFTT), was found in [61]. 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 [87] shows that for a CMB satellite of next generation one could get eV alone, thanks to a very precise reconstruction of the CMB lensing potential, while [78] discusses the potential of cluster surveys. Finally, the authors of [60] 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.

#### 9. 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 it can provide 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.

#### Acknowledgments

S. Pastor was supported by the Spanish Grants FPA2011-22975, Multidark Consolider CSD2009-00064 (MINECO), and PROMETEO/2009/091 (Generalitat Valenciana), as well as by the EC Contract UNILHC PITN-GA-2009-237920. He thanks the Galileo Galilei Institute for Theoretical Physics for the hospitality and the INFN for partial support during the completion of this work.