text:
# Neutrinos from the Big Bang

## 1 Relic neutrinos and the cosmological mass bound

Several years before neutrinos had even been experimentally detected, Alpher, Follin & Hermann  noted that they would have been in thermal equilibrium in the early universe “…through interactions with mesons” at temperatures above $`5MeV`$; below this temperature the neutrinos “…freeze-in and continue to expand and cool adiabatically as would a pure radiation gas”. These authors also observed that the subsequent annihilation of $`e^\pm `$ pairs would heat the photons but not the decoupled neutrinos, so by conservation of entropy $`T_\nu /T`$ would decrease from its high temperature value of unity, down to $`(4/11)^{1/3}`$ at $`Tm_e`$ .

(Thus the present density of massless relic neutrinos (neglecting a possible chemical potential i.e. lepton–antilepton asymmetry) would be

$$\frac{n_\nu }{n_\gamma }=\left(\frac{T_\nu }{T}\right)^3\left(\frac{n_\nu }{n_\gamma }\right)_{T=T_{dec}}=\frac{4}{11}\left(\frac{3}{4}\frac{g_\nu }{2}\right),$$
(1)
where $`g_\nu =2`$ corresponds to left-handed neutrinos and right-handed antineutrinos, and the factor 3/4 reflects Fermi versus Bose statistics. This would also be true for massive neutrinos if the neutrinos are relativistic at decoupling i.e. for $`m_\nu T_{dec}`$. For a present cosmic microwave background (CMB) blackbody temperature $`T_0=2.725\pm 0.002K`$ , the abundance per flavour should then be $`\frac{3}{11}\times \frac{2\zeta (3)}{\pi ^2}T_0^3111.9cm^3`$. As long as they remain relativistic, the neutrinos would retain a Fermi-Dirac distribution with phase-space density

$$f_\nu =\frac{g_\nu }{(2\pi )^3}\left[\mathrm{exp}\left(\frac{p}{T_\nu }\right)+1\right]^1,$$
(2)
since the momentum and temperature would redshift identically.)

Subsequently, Chiu & Morrison  found the rate for $`e^+e^{}\nu _e\overline{\nu _e}`$ in a plasma to be $`\mathrm{\Gamma }_\nu G_F^2T^5`$ for the universal Fermi interaction and Zel’dovich  equated this to the Hubble expansion rate in the radiation-dominated era,

$$H=\sqrt{\frac{8\pi G_N\rho }{3}},with\rho =\frac{\pi ^2}{30}g_{}T^4,$$
(3)
where $`g_{}`$ counts the relativistic degrees of freedom, to obtain the ‘decoupling’ temperature below which the neutrinos expand freely without further interactions as $`T_{dec}(\nu _e)2MeV`$. (Neutral currents were then unknown so $`T_{dec}(\nu _\mu )`$ was estimated from the reaction $`\mu e\overline{\nu _e}\nu _\mu `$ to be $`12MeV`$. Later De Graaf  noted that these would keep $`\nu _\mu `$’s coupled to the plasma down to the same temperature as $`\nu _e`$’s.111In fact $`T_{dec}(\nu _\mu ,\nu _\tau )3.5MeV`$ while $`T_{dec}(\nu _e)2.3MeV`$ because of the additional charged current reaction . Actually decoupling is not an instantaneous process so the neutrinos are slightly heated by the subsequent $`e^+e^{}`$ annihilation, increasing the number density (1) by $`1\%`$ .) However Zel’dovich  and Chiu  concluded that relic neutrinos, although nearly as numerous as the blackbody photons, cannot make an important contribution to the cosmological energy density since they are probably massless.

Interestingly enough, some years earlier Pontecorvo & Smorodinski  had discussed the bounds set on the cosmological energy density of MeV energy neutrinos (created e.g. by large-scale matter-antimatter annihilation) using data from the Reines–Cowan and Davis experiments.222These authors  were also the first to suggest searching for high energy neutrinos by looking for upward going muons in underground experiments, the basis for today’s neutrino telescopes.) Not surprisingly these bounds were rather weak so these authors stated somewhat prophetically that “…it is not possible to exclude a priori the possibility that the neutrino and antineutrino energy density in the Universe is comparable to or larger than the average energy density contained in the proton rest mass”. Zel’dovich & Smorodinski  noted that better bounds can be set by the limit on the cosmological energy density in any form of matter $`\rho _m(\mathrm{\Omega }_m\rho _c`$) following from the observed present expansion rate $`H_0`$ and age $`t_0`$ of the universe.333The critical density is $`\rho _c=3H_0^2/8\pi G_N1.879\times 10^{29}h^2gmcm^3`$ where the Hubble parameter $`hH_0/100kmsec^1Mpc^1`$, so the present Hubble age is $`H_0^1=9.778h^1Gyr`$. Of course they were still discussing massless neutrinos. Weinberg  even speculated whether a degenerate sea of relic neutrinos can saturate the cosmological energy density bound and noted that such a sea may be detectable by searching for (scattering) events beyond the end-point of the Kurie plot in $`\beta `$-decay experiments.

Several years later, Gershteĭn & Zel’dovich  made the connection that if relic neutrinos are massive, then a bound on the mass follows from simply requiring that

$$m_\nu n_\nu <\rho _m.$$
(4)
Using the general relativistic constraint $`\mathrm{\Omega }t_0^2H_0^2<(\pi /2)^2`$, they derived $`\rho _m<2\times 10^{28}gmcm^3`$ (just assuming that $`t_0>5Gyr`$, i.e. that the universe is older than the Earth) and inferred that $`m_{\nu _e},m_{\nu _\mu }<400`$ eV for a present photon temperature of $`3K`$. Their calculation of the relic neutrino abundance was rather approximate — they adopted $`g_\nu =4`$ i.e. assumed massive neutrinos to be Dirac particles with fully populated right-handed (RH) states (even though they commented en passim that according to the $`VA`$ theory such states are non-interacting and would thus not be in equilibrium at $`T_{dec}`$) and moreover they did not allow for the decrease in the neutrino temperature relative to photons due to $`e^+e^{}`$ annihilation. Nevertheless their bound was competitive with the best laboratory bound on $`m_{\nu _e}`$ and $`10^4`$ times better than that on $`m_{\nu _\mu }`$, demonstrating the sensitivity (if not the precision!) of cosmological arguments.

A better bound of $`m_{\nu _\mu }<130`$ eV was quoted by Marx & Szalay  who numerically integrated the cosmological Friedmann equation from $`\nu _\mu `$ decoupling down to the present epoch, subject to the condition $`t_0>4.5Gyr`$. Independently Cowsik & McClleland  used direct limits on $`\mathrm{\Omega }_m`$ and $`h`$ to obtain an even more restrictive bound of $`m_\nu <8`$ eV, assuming that $`m_\nu =m_{\nu _e}=m_{\nu _\mu }`$; however they too assumed that $`T_\nu =T`$ and that RH states were fully populated. As Shapiro, Teukolsky & Wasserman  first emphasized, even if massive neutrinos are Dirac rather than Majorana, the RH states have no gauge interactions so should have decoupled much earlier than the left-handed ones. Then subsequent entropy generation by massive particle annihilations would have diluted their relic abundance to a negligible level.444Although spin-flip scattering (at a rate $`(m_\nu /T)^2`$) can generate RH states, this can be neglected for $`m_\nu 1MeV`$. Even if RH neutrinos have new (superweak) interactions, their relic abundance can be no more than $`10\%`$ of LH neutrinos so the bound on their mass (5) is relaxed to $`1`$ keV . Even if the (post-inflationary) universe was not hot enough to bring such interactions into equilibrium, a cosmologically interesting abundance can still be generated if the RH states have small mixings with LH states . Thus we arrive at the modern version of the ‘Gershteĭn-Zel’dovich’ bound’ : the conservative limits $`t_0>10Gyr`$ and $`h>0.4`$ imply $`\mathrm{\Omega }_mh^2<1`$ i.e. $`\rho _m<10.54keVcm^3`$ ; combining this with the relic neutrino number density, which is $`1\%`$ larger  than in eq.(1), gives:555If neutrinos are non-relativistic at decoupling, then they drop out of chemical equilibrium with an abundance inversely proportional to their self-annihilation cross-section so $`\mathrm{\Omega }_\nu h^2(m_\nu /2GeV)^2`$ for $`m_\nu m_Z`$ . Thus neutrinos with a mass of $`𝒪(GeV)`$ can also account for the dark matter; however LEP has ruled out such (possible 4th generation) neutrinos upto a mass $`m_Z/2`$ . (Conversely $`\mathrm{\Omega }_\nu h^2>1`$ for the mass range $`100eV2GeV`$, which is thus cosmologically forbidden for any stable neutrino having only electroweak interactions .) For heavier masses upto the highest plausible value of $`𝒪`$(TeV), the relic abundance decreases steadily due to the increasing annihilation cross-section and remains cosmologically uninteresting .

$$\mathrm{\Omega }_\nu h^2=\underset{i}{}\left(\frac{m_{\nu _i}}{93eV}\right)\left(\frac{g_{\nu _i}}{2}\right)<1.$$
(5)
This is a rather conservative bound since galaxy surveys indicate much tighter constraints on the total amount of gravitating (dark) matter in the universe, e.g. the observed ‘redshift space distortion’  suggests $`\mathrm{\Omega }_m0.3`$, averaged over a volume extending several hundred Mpc. Together with the Hubble Key Project determination of $`h=0.72\pm 0.08`$ , this implies that the sum of all neutrino masses cannot exceed about 15 eV.

Although this has historically been the most restrictive constraint on neutrino masses, it is no longer competetive with the direct laboratory bound on the electron neutrino mass from the Mainz and Troitsk tritium $`\beta `$-decay experiments :

$$m_\nu <2.2eV(95\%\mathrm{c}.\mathrm{l}.).$$
(6)
Although the kinematic mass limits on the other neutrino flavours are much weaker (viz. $`m_{\nu _\mu }<190`$ keV, $`m_{\nu _\tau }<18.2`$ MeV ), the bound above now applies in fact to all eigenstates  given the rather small mass-differences indicated by the oscillation interpretation of the Solar ($`\mathrm{\Delta }m^27\times 10^5eV^2`$) and atmospheric ($`\mathrm{\Delta }m^23\times 10^3eV^2`$) neutrino anomalies . This means that the laboratory limit on the sum of neutrino masses (for comparison with the new cosmological bounds to be discussed) is presently 6.6 eV; the sensitivity of planned future experiments is at the eV level . Moreover relic neutrinos contribute at least $`\mathrm{\Omega }_\nu h^20.07/938\times 10^4`$ to the cosmic budget (assuming a mass hierarchy), nearly as much as visible baryons .

## 2 Neutrinos as the galactic ‘missing mass’

The cosmological bound (5) assumes conservatively that neutrinos constitute all of the (dark) matter permitted by the dynamics of the universal Hubble expansion. Further constraints must be satisfied if they are to cluster on a specified scale (e.g. galactic halos or galaxy clusters) and provide the dark matter whose presence is inferred from dynamical measurements. Cowsik & McClleland  were the first to suggest that neutrinos with a mass of a few eV could naturally be the ‘missing mass’ in clusters of galaxies. This follows from the relation $`m_\nu ^81/G_N^3r_{cl}^3M_{cl}`$ (reflecting the Pauli principle) which they obtained by modeling a cluster of mass $`M_{cl}`$ as a square potential well of core radius $`r_{cl}`$ filled with a Fermi-Dirac gas of neutrinos at zero temperature. Subsequently Tremaine & Gunn  noted that this provides a lower bound on the neutrino mass. Although the microscopic phase-space density (2) is conserved for collisionless particles, the ‘coarse-grained’ phase-space density in bound objects can decrease below its maximum value of $`g_\nu /2(2\pi )^3`$ during structure formation. Modeling the bound system as an isothermal sphere with velocity dispersion $`\sigma `$ and core radius $`r_{cl}^2=9\sigma ^2/4\pi G_N\rho (r_{cl})`$ then gives

$$m_\nu >120eV\left(\frac{\sigma }{100kmsec^1}\right)^{1/4}\left(\frac{r_{cl}}{kpc}\right)^{1/2}.$$
(7)
This is indeed consistent with the cosmological upper bound (5) down to the scale of galaxies, however there is a conflict for smaller objects, viz. dwarf galaxies which require a minimum mass of $`100`$ eV . In fact the central phase space density of observed dark matter cores in these structures decreases rapidly with increasing core radius, rather than being constant as would be expected for neutrinos . Moreover since neutrinos would cluster more efficiently in larger potential wells, there should be a trend of increasing mass-to-light ratio with scale. This was indeed claimed to be the case initially  but later it was recognised that the actual increase is far less than expected . Thus massive neutrinos are now disfavoured as the constituent of the ‘missing mass’ in galaxies and clusters.

## 3 The rise and fall of ‘hot dark matter’

Nevertheless such cosmological arguments became of particular interest in the 1980’s after the ITEP tritium $`\beta `$-decay experiment claimed a $`30`$ eV mass for the electron neutrino. The attention of cosmologists turned to how the large-scale structure (LSS) of galaxies, clusters and superclusters would have formed if the universe is indeed dominated by such massive neutrinos. The basic picture  is that structure grows through gravitational instability from primordial density perturbations; these perturbations were first detected by the COsmic Background Explorer (COBE) via the temperature fluctuations they induce in the CMB . On small scales ($`<10Mpc`$) structure formation is complicated by non-linear gravitational clustering as well as non-gravitational (gas dynamic) processes but on large scales gravitational dynamics is linear and provides a robust probe of the nature of the dark matter.

Density perturbations in a medium composed of relativistic collisionless particles are subject to a form of Landau damping (viz. phase-mixing through free streaming of particles from high to low density regions) which effectively erases perturbations on scales smaller than the free-streaming length $`41Mpc(m_\nu /30eV)^1`$ . This is essentially the (comoving) distance traversed by a neutrino from the Big Bang until it becomes non-relativistic, and corresponds to the scale of superclusters of galaxies. Thus huge neutrino condensations (generically in the shape of ‘pancakes’), containing a mass $`3\times 10^{15}M_{}(m_\nu /30eV)^2`$, would have begun growing at a redshift $`z_{eq}7\times 10^3(m_\nu /30eV)`$ when the universe became matter-dominated and gravitational instability set in. This is well before the epoch of (re)combination at $`z_{rec}10^3`$ so the baryons were still closely coupled to the photons, while the neutrinos were mildly relativistic ($`v/c0.1`$) hence ‘hot’. After the universe became neutral, baryonic matter would have accreted into these potential wells, forming a thin layer of gas in the central plane of the pancakes. Thus superclusters would be the first objects to condense out of the Hubble flow in a ‘hot dark matter’ (HDM) cosmogony, and smaller structures such as galaxies would form only later through the fragmentation of the pancakes.

The gross features of such a ‘top-down’ model for structure formation are compatible with several observed features of LSS, in particular the distinctive ‘voids’ and ‘filaments’ seen in large galaxy surveys. It was also noted that since primordial density perturbations can begin growing earlier than in an purely baryonic universe, their initial amplitude must have been smaller, consistent with extant limits on the isotropy of the microwave background. Detailed studies  found however that galaxies form too late through the breakup of the pancakes, at a redshift $`z<1`$, counter to observations of galaxies and quasars at $`z>4`$. (Another way of saying this is that galaxies should have formed last in an HDM universe, whereas our Galaxy is in fact dynamically much older than the local group .) There are other difficulties such as too high ‘peculiar’ (i.e. non-Hubble) velocities , excessive X-ray emission from baryons which accrete onto neutrino clusters , and too large voids  (although detailed simulations  showed later that some of these problems had perhaps been exaggerated).

Therefore cosmologists abandoned HDM and turned, with considerably more success, to cold dark matter (CDM) , i.e. particles which were non-relativistic at the epoch of matter-domination . Detailed studies of CDM universes gave excellent agreement with observations of galaxy clustering and a ‘standard CDM’ model for large-scale structure formation was established, viz. a critical density CDM dominated universe with an initially scale-invariant spectrum of density perturbations . Moreover particle physicists provided plausible candidate particles, notably the neutralino in supersymmetric models with conserved $`R`$-parity which naturally has a relic abundance of order the critical density .

## 4 COBE and the advent of ‘mixed dark matter’

Nevertheless neutrinos were resuscitated some years later as a possible (sub-dominant) component of the dark matter when the CDM cosmogony itself ran into problems. To appreciate the background to this it is necessary to recapitulate the essential ingredients of a model for cosmic structure formation. A key assumption made concerns the nature of the primordial density perturbations which grow through gravitational instability in the dark matter. Cosmologists usually assume such fluctuations to have a power spectrum of the scale-free form:

$$P(k)=|\delta _k|^2=Ak^n,$$
(8)
with $`n=1`$ corresponding to the scale-invariant ‘Harrison-Zel’dovich’ spectrum. Here $`\delta _k\frac{\delta \rho (\stackrel{}{x})}{\overline{\rho }}\mathrm{e}^{i\stackrel{}{k}\stackrel{}{x}}d^3x`$ is the Fourier transform of spatial fluctuations in the density field (of wavelength $`\lambda =2\pi /k`$). Moreover the perturbations are assumed to be gaussian (i.e. different phases in the plane-wave expansion are uncorrelated) and to be ‘adiabatic’ (i.e. matter and radiation fluctuate together). Concurrent with the above studies concerning the nature of the dark matter, powerful support for this conjecture was provided by the development of the ‘inflationary universe’ model . Here the perturbations arise from quantum fluctuations of a scalar field $`\varphi `$, the vacuum energy of which drives a period of accelerated expansion in the early universe. The classical density perturbation thus generated has a spectrum determined by the ‘inflaton’ potential $`V(\varphi )`$, with a power-law index which is dependent on $`k`$ :

$$n(k)=13M^2\left(\frac{V^{}}{V}\right)_{}^2+2M^2\left(\frac{V^{\prime \prime }}{V}\right)_{}$$
(9)
where $`M(8\pi G_N)^{1/2}2.4\times 10^{18}GeV`$ is the Planck mass and $``$ denotes that this is to be evaluated when a mode of wavenumber $`k`$ crosses the Hubble radius $`H^1`$. For a sufficiently ‘flat potential’ (as is necessary to drive enough e-folds of inflation to solve the problems of the standard cosmology), the spectrum indeed has $`n1`$, with corrections $`\mathrm{ln}(k)`$.

Gravitational instability sets in only when the universe becomes matter-dominated and this modifies the spectrum on length scales smaller than the Hubble radius at this epoch, viz. for $`k>k_{eq}^1\mathrm{\hspace{0.17em}80}h^1Mpc`$. Thus the characteristics of the dark matter can be encoded into a ‘transfer function’ $`T(k)`$ which modulates the primordial spectrum; for HDM this is an exponentially falling function while for CDM it is a more gradual power-law. Now the power spectrum inferred from observations may be compared with theoretical models, but another problem arises concerning how we are to normalize the amplitude of the primordial density perturbations,666This is determined by the inflaton potential $`V(\varphi )`$ but there is, as yet, no ‘standard model’ of inflation . particularly since these are in the dark matter and may differ significantly (i.e. be ‘biased’) from the observable fluctuations in the density of visible galaxies. Fortunately, the primordial perturbations have another unique observational signature, viz. they induce temperature fluctuations in the CMB through the ‘Sachs-Wolfe effect’ (gravitational red/blue shifts) on large angular scales $`>2^0`$, corresponding to spatial scales larger than the Hubble radius on the last scattering surface . It was the COBE measurement of these fluctuations a decade ago that initiated the modern era of cosmological structure formation studies.

The quadrupole anisotropy in the CMB measured by COBE  allows a determination of the fluctuation amplitude at the scale, $`H_0^1\mathrm{\hspace{0.17em}3000}h^1Mpc`$, corresponding to the present Hubble radius. With this normalization it became clear that a $`\mathrm{\Omega }_\nu 1`$ HDM universe indeed had too little power on small-scales for adequate galaxy formation.777To save HDM would require new sources of small-scale fluctuations, e.g. relic topological defects , or isocurvature primordial perturbations  — to date however there is no evidence for either. However it also became apparent that the ‘standard CDM model’ when normalized to COBE had too much power on small-scales (see Fig. 1). It was thus a logical step to invoke a suitable mixture of CDM and HDM to try and match the theoretical power spectrum to the data on galaxy clustering and motions .

In fact the possibility that the dark matter may have both a hot and a cold component had been discussed several years earlier by Shafi & Stecker, motivated by theoretical considerations of SUSY GUTs , and the possibilities for solving the problems of the standard CDM model had been noted . In the post-COBE era, a number of detailed studies of mixed dark matter (MDM) universes were performed and a neutrino fraction of about $`20\%`$ was found to give the best match with observations . The implied neutrino mass was $`5eV`$, presumably that of the $`\nu _\tau `$ given the usual hierarchy implied by the see-saw mechanism . The suppression of small-scale power delays the epoch of galaxy formation, so constraints on the HDM component can also be obtained from the abundance of high redshift objects such as QSOs and Ly-$`\alpha `$ systems; this gave an upper bound of 4.7 eV on the neutrino mass . More baroque schemes in which two neutrinos have comparable masses ($`m_{\nu _\mu }m_{\nu _\tau }2.5eV`$) were also constructed  seeking to reconcile the LSND report of neutrino oscillations  with the atmospheric neutrino anomaly which had just been reported by Kamiokande . This was an exciting time for neutrino cosmology as both laboratory data and astronomical observations supported the possibility that a substantial fraction of the cosmological mass density is in the form of massive neutrinos.

However another way to reconcile a CDM universe with the small-scale observations is to relax the underlying assumption that the primordial spectrum is strictly scale-invariant. As shown in Fig.1, a ‘tilted’ spectrum with $`P(k)k^{0.9}`$ also gives a good fit to the data . At first sight this might strike one as simply introducing an additional parameter (although this is arguably no worse than introducing an additional form of dark matter). However one should really ask why the spectrum should be assumed to have a power-law index $`n=1`$ in the first place. As indicated in Eq.(9), $`n`$ in fact varies slowly with $`k`$ and is determined by the slope and curvature of the scalar potential at the epoch when the fluctuation at a specified value of $`k`$ crosses the Hubble radius. The corresponding number of e-folds before the end of inflation is just $`N_{}(k)51+\mathrm{ln}\left(\frac{k^1}{3000h^1Mpc}\right)`$, for typical choices of the inflationary scale, reheat temperature etc. We see that fluctuations on the scales ($`13000Mpc`$) probed by LSS and CMB observations are generated just $`4050`$ e-folds before the end of inflation. It is quite natural to expect the inflaton potential to begin curving significantly as the end of inflation is approached, especially in ‘new inflation’ (small field) models. There are certainly attractive models of inflation in which the spectrum is significantly tilted in this region, in particular an inflationary model based on $`N=1`$ supergravity  naturally gives $`n(k)(N_{}2)/(N_{}+2)0.9`$ at these scales. With such a scale-dependent tilt for the primordial spctrum, the LSS data can be fitted reasonably well, with no need for any HDM component . (Of course in the absence of a ‘standard model’ of inflation, it might be argued that the inflationary spectrum may instead have $`n>1`$ thus allowing a larger HDM component.)

As Fig.1 shows, yet another way of suppressing small-scale power in the CDM cosmogony is to decrease the matter content of the universe, since this postpones the epoch of matter-radiation equality and thus shifts the peak of power spectrum to larger scales. Furthermore the spatial geometry can be maintained flat if there is a cosmological constant with $`\mathrm{\Omega }_\mathrm{\Lambda }=1\mathrm{\Omega }_m0.7`$. Evidence for just such a cosmology ($`\mathrm{\Lambda }`$CDM) has come subsequently from observations of the Hubble diagram of Type Ia supernovae which suggest that the expansion is in fact accelerating , coupled with the observation that $`\mathrm{\Omega }_m`$ does not exceed $`0.3`$ even on the largest scales probed . Assuming such a $`\mathrm{\Lambda }`$CDM cosmology and a scale-invariant power spectrum gives a limit of $`f_\nu \mathrm{\Omega }_\nu /\mathrm{\Omega }_m<0.13`$ ($`95\%\mathrm{c}.\mathrm{l}.`$) using the power spectrum of galaxy clustering determined from the Two degree Field galaxy redshift survey (2dFGRS) as shown in Fig. 2. This corresponds to an upper bound of $`m_\nu <1.8eV(95\%\mathrm{c}.\mathrm{l}.)`$ on the sum of neutrino masses, adopting ‘concordance’ values of $`\mathrm{\Omega }_m`$ and $`h`$ . If the spectral index is allowed to vary in the range $`n=1\pm 0.1`$, then this bound is relaxed to 2.1 eV.

## 5 CMB anisotropy and limits on neutrino dark matter

It is clear that there are uncertainties in the determination of the HDM component from LSS data alone. Fortunately it is possible to reduce these substantially by examining an independent probe of the primordial power spectrum, viz. temperature fluctuations in the CMB.

In general a skymap of the CMB temperature can be decomposed into spherical harmonics

$$T(\theta ,\varphi )=\underset{l=0}{\overset{\mathrm{}}{}}\underset{m=l}{\overset{l}{}}a_l^mY_l^m(\theta ,\varphi ),$$
(10)
where the $`l^{\mathrm{th}}`$ multipole corresponds to an angle $`\theta ^{}200/l`$ and probes spatial scales around $`k^16000h^1l^1Mpc`$. In inflationary theories, the fluctuations are gaussian so the co-efficients $`a_l^m`$ are independent stochastic variables with zero mean and variance $`C_l=|a_l^m|^2`$; each $`C_l`$ has a $`\chi ^2`$ distribution with $`(2l+1)`$ degrees of freedom . For an assumed set of cosmological parameters and given the primordial density perturbation spectrum, the $`C_l`$’s can be determined by solution of the Einstein-Boltzmann equations which describe how the different components (photons, ions, electrons, neutrinos …) evolve . Thus theoretical estimates of the power at each multipole can be compared with observations. The low multipoles (large spatial scales) are sensitive to the primordial spectrum alone,888In principle, primordial gravitational waves can also make a contribution here but this is expected to be negligible in ‘small-field’ inflationary models . but the measurements in the region are particularly uncertain, both because of uncertainties in the foreground substraction and also because there are fewer independent measurements for low multipoles (‘cosmic variance’). For example, by measuring the first $``$20 multipoles COBE could only fix $`n=1.2\pm 0.3`$  so could not discriminate between a scale-invariant and a mildly tilted spectrum. However subsequent ground-based experiments with angular resolution far superior to COBE’s have now measured the power at higher multipoles . The dominant features in the power spectrum here are the ‘acoustic peaks’, the most prominent being at $`l200`$, arising from oscillations of the coupled plasma-photon fluids at last scattering . The position of the first peak is a measure of the horizon length at the epoch of (re)combination of the primordial plasma and thus provides a measure of space curvature — observations by the BOOMERanG  and MAXIMA  experiments were the first to show that this is in fact close to zero. Taken together with the earlier recognition that CDM alone cannot make up the critical density, this has led to the widespread adoption of the $`\mathrm{\Lambda }`$CDM model in which $`\mathrm{\Omega }_m0.3`$, $`\mathrm{\Omega }_\mathrm{\Lambda }0.7`$ .

As seen in Fig. 3, the expectations for CMB anisotropy in a MDM universe do not differ significantly from a CDM universe having the same initial perturbation spectrum. Thus it is clear that by combining CMB and LSS data, it would be possible to determine whether the suppression of small-scale power is intrinsic to the primordial spectrum of inflationary fluctuations, or is induced by a HDM component. However there are additional uncertainties in the input values of the other cosmological parameters, in particular the values of the Hubble parameter and of the baryon density which is usually inferred from considerations of primordial nucleosynthesis. In analysing the observational data, various ‘priors’ are often assumed for these quantities on the basis of other observations. For example a detailed likelihood analysis  yielded the neutrino density fraction $`f_\nu `$ to be at most 0.08 from CMB observations alone, decreasing to 0.06 if LSS data from the PSCz survey is added, and further to 0.04 if the value of $`h=0.72\pm 0.08`$  is adopted. Another analysis  used the more precise 2dFGRS data , together with additional constraints on the matter density from the SNIa data ($`\mathrm{\Omega }_m=0.28\pm 0.14`$ assuming a flat universe) , and on the baryon density from primordial nucleosynthesis ($`\mathrm{\Omega }_Bh^2=0.02\pm 0.002(95\%\mathrm{c}.\mathrm{l}.)`$ ); using the analytic result for the suppression of the power spectrum, $`\mathrm{\Delta }P/P8f_\nu `$ , a bound of $`m_\nu <2.53eV`$ was obtained.

It had been estimated  that the CMB data expected from the Wilkinson Microwave Anisotropy Probe (WMAP) , combined with LSS data from the Sloane Digital Sky Survey (SDSS) , will provide sensitivity to neutrino mass at the eV level. Indeed the recently announced first results from WMAP, combined with 2dFGRS, have already set the bound $`\mathrm{\Omega }_\nu h^2<0.0076`$ corresponding to $`m_\nu <0.7eV(95\%\mathrm{c}.\mathrm{l}.)`$, as shown in Fig. 4 . This severely restricts , but does not altogether rule out , a fourth (sterile) neutrino with a mass of $`𝒪(1)`$ eV as suggested by the LSND experiment . Combined with the new KamLAND data, this also restricts the accessible range for the observation of neutrinoless $`\beta \beta `$ decay , which has been claimed to have been seen already implying an effective Majorana mass of $`|m_\nu |=0.39_{0.34}^{+0.45}eV(95\%\mathrm{c}.\mathrm{l}.)`$ .

## 6 The nucleosynthesis limit on $`N_\nu `$

Hoyle & Taylor  as well as Peebles  had emphasized many years ago that new types of neutrinos (beyond the $`\nu _e`$ and $`\nu _\mu `$ then known) would boost the relativistic energy density hence the expansion rate (3) during Big Bang nucleosynthesis (BBN), thus increasing the yield of $`{}_{}{}^{4}\mathrm{He}`$. Shvartsman  noted that new superweakly interacting particles would have a similar effect. Subsequently this argument was refined quantitatively by Steigman, Schramm & Gunn . In the pre-LEP era when the laboratory bound on the number of neutrino species was not very restrictive , the BBN constraint was used to argue that at most one new family was allowed , albeit with considerable uncertainties . Although LEP now finds $`N_\nu =2.994\pm 0.012`$ , the cosmological bound is still important since it is sensitive to any new light particle , not just $`SU(2)_L`$ doublet neutrinos, so is a particularly valuable probe of new physics, e.g. neutrinos coupled to new gauge bosons expected in string models .

BBN limits on neutrinos come mainly from the observational bounds on the primordial $`{}_{}{}^{4}\mathrm{He}`$ abundance, termed $`Y_\mathrm{p}`$ by astronomers. This is proportional to the neutron fraction which ‘freezes out’ at $`n/p\mathrm{exp}[(m_nm_p)/T_{\mathrm{fr}}]`$ when the weak interaction rate ($`G_F^2T^5`$) falls behind the Hubble expansion rate (3), at $`T_{\mathrm{fr}}(g_{}G_N/G_F^4)^{1/6}1`$ MeV. The presence of additional neutrino flavors (or of any other relativistic species) at this time increases $`g_{}`$, hence the expansion rate $`H`$, leading to a larger value of $`T_{\mathrm{fr}}`$, $`n/p`$, and thus $`Y_\mathrm{p}`$. In the Standard Model, the number of relativistic particle species at 1 MeV is

$$g_{}=5.5+\frac{7}{4}N_\nu ,$$
(11)
where 5.5 accounts for photons and $`e^\pm `$, and $`N_\nu `$ is the number of (massless) neutrino flavors. (The energy density of new light fermions $`i`$ is equivalent to an effective number $`\mathrm{\Delta }N_\nu =_i(g_i/2)(T_i/T_\nu )^4`$ of additional doublet neutrinos, where $`T_i/T_\nu `$ follows from considerations of their (earlier) decoupling.) In Fig. 5, we show the expected abundance of $`{}_{}{}^{4}\mathrm{He}`$ taking $`N_\nu =3`$, as a function of the density of baryons normalized to the blackbody photon density:

$$\eta \frac{n_B}{n_\gamma }2.728\times 10^8\mathrm{\Omega }_Bh^2,$$
(12)
with $`\eta _{10}\eta /10^{10}`$. The computed $`{}_{}{}^{4}\mathrm{He}`$ abundance scales linearly with $`\mathrm{\Delta }N_\nu `$ but it also increases logarithmically with $`\eta `$ :

$$\mathrm{\Delta }Y0.012\mathrm{\Delta }N_\nu +0.01\mathrm{ln}\left(\frac{\eta _{10}}{5}\right).$$
(13)
Thus to obtain a bound on $`N_\nu `$ requires an upper limit on $`Y_\mathrm{p}`$ and a lower limit on $`\eta `$. The latter is poorly determined from direct observations of luminous matter  so must be derived from the abundances of the other synthesized light elements, $`D`$, $`{}_{}{}^{3}\mathrm{He}`$ and $`{}_{}{}^{7}\mathrm{Li}`$, which are power-law functions of $`\eta `$ as seen in Fig. 5. The complication is that these abundances are substantially altered in a non-trivial manner during the chemical evolution of the galaxy, unlike $`Y_\mathrm{p}(^4\mathrm{He})`$ which just increases by a few percent due to stellar production. (This can be tagged via the correlated abundance of oxygen and nitrogen which are made only in stars.) Even so, chemical evolution arguments have been used to limit the primordial abundances of $`D`$ and $`{}_{}{}^{3}\mathrm{He}`$ and thus derive increasingly severe bounds on $`N_\nu `$ , culminating in one below 3 ! However a more conservative view  is that there is no ‘crisis’ with BBN if we recognize that chemical evolution arguments are unreliable and consider only direct measurements of light element abundances.

Fig. 5 shows the current status of such measurements, as reviewed in more detail elsewhere . We observe $`{}_{}{}^{4}\mathrm{He}`$ in clouds of ionized hydrogen (H II regions), the most metal-poor of which are in dwarf blue compact galaxies (BCGs). There is now a large body of data on $`{}_{}{}^{4}\mathrm{He}`$ and C, N, O in these systems which confirm that the small stellar contribution to helium is positively correlated with ‘metal’ production; extrapolating to zero metallicity gives the primordial $`{}_{}{}^{4}\mathrm{He}`$ abundance :

$$Y_\mathrm{p}=0.238\pm 0.002\pm 0.005.$$
(14)
Here and subsequently, the first error is statistical, and the second is an estimate of the systematic uncertainty. The latter clearly dominates, and is based on the scatter in different analyses of the physical properties of the H II regions . Other extrapolations to zero metallicity give $`Y_\mathrm{p}=0.244\pm 0.002`$ , and $`Y_\mathrm{p}=0.235\pm 0.003`$ , while the average in the 5 most metal-poor objects is $`Y_\mathrm{p}=0.238\pm 0.003`$ . The value in Eq.(14), shown in Fig. 5, is consistent with all these determinations.

The systems best suited for Li observations are hot, metal-poor stars belonging to the halo population (Pop II) of our Galaxy. Observations have long shown that Li does not vary significantly in such stars having metallicities less than 1/30 of Solar — the ‘Spite plateau’. Recent precision data suggest a small but significant correlation between Li and Fe  which can be understood as the result of Li production from cosmic rays . Extrapolating to zero metallicity one arrives at a primordial value

$$\mathrm{Li}/\mathrm{H}|_\mathrm{p}=(1.23\pm 0.06{}_{0.32}{}^{+0.68}{}_{}{}^{+0.56})\times 10^{10}.$$
(15)
The last error is our estimate of the maximum upward correction necessary to allow for possible destruction of Li in Pop II stars, due to e.g. mixing of the outer layers with the hotter interior . Such processes can be constrained by the absence of significant scatter in the Li-Fe correlation plot , and through observations of the even more fragile isotope 6Li .

In recent years, high-resolution spectra have revealed the presence of D in quasar absorption systems (QAS) at high-redshift, via its isotope-shifted Lyman-$`\alpha `$ absorption. It is believed that there are no astrophysical sources of deuterium, so any measurement of D/H provides a lower limit to the primordial abundance and thus an upper limit on $`\eta `$; for example, the local interstellar value of $`D/H=(1.5\pm 0.1)\times 10^5`$  requires $`\eta _{10}9`$. Early reports of $`D/H>10^4`$ towards 2 quasars (Q0014+813  and PG1718+4807 ) have been undermined by later analyses . Three high quality observations yield $`D/H=(3.2\pm 0.3)\times 10^5`$ (PKS1937-1009), $`(4.0\pm 0.7)\times 10^5`$ (Q1009+2956), and $`(2.5\pm 0.2)\times 10^5`$ (HS0105+1619); their average value

$$D/H=(3.0\pm 0.4)\times 10^5$$
(16)
has been widely promoted as the primordial abundance . Recently, the same group  have provided another measurement: $`D/H=(2.42_{0.25}^{+0.35})\times 10^5`$ (HS 243+3057). However the observed dispersion in the measurements suggests either that systematic uncertainties have been underestimated, or that there is intrinsic dispersion in the D abundance in QAS. Other values have been reported in different (damped Lyman-$`\alpha `$) systems which have a higher column density of neutral H, viz. $`D/H=(2.24\pm 0.67)\times 10^5`$ (Q0347-3819)  and $`D/H=(1.65\pm 0.35)\times 10^5`$ (Q2206-199) . Moreover, allowing for a more complex velocity structure than assumed in these analyses raises the inferred abundance by upto $`50\%`$ . Even the ISM value of D/H now shows unexpected scatter of a factor of 2 . All this may indicate significant processing of the $`D`$ abundance even at high redshift. Given these uncertainties, it is prudent to bound the primordial abundance with an upper limit set by the non-detection of D absorption in a high-redshift system (Q0130-4021) , and the lower limit set by the local interstellar value , both at $`2\sigma `$:

$$1.3\times 10^5<D/H|_\mathrm{p}<9.7\times 10^5.$$
(17)
For $`{}_{}{}^{3}\mathrm{He}`$, the only observations available are in the Solar system and (high-metallicity) H II regions in our Galaxy . This makes inference of the primordial abundance difficult, a problem compounded by the fact that stellar nucleosynthesis models for $`{}_{}{}^{3}\mathrm{He}`$ are in conflict with observations . Such conflicts can perhaps be resolved if a large fraction of low mass stars destroy $`{}_{}{}^{3}\mathrm{He}`$ by internal mixing driven by stellar rotation, consistent with the observed 12C/13C ratios . The observed abundance ‘plateau’ in H II regions then implies a limit on the primordial value of $`{}_{}{}^{3}\mathrm{He}/H<(1.9\pm 0.6)\times 10^5`$ , which is consistent with the other abundance constraints we discuss.

The overlap in the $`\eta `$ ranges spanned by the larger boxes in Fig. 5 indicates overall concordance between the various abundance determinations. Accounting for theoretical uncertainties as well as the statistical and systematic errors in observations, there is acceptable agreement among the abundances when

$$2.6\eta _{10}6.2.$$
(18)
However the agreement is far less satisfactory if we use only the quoted statistical errors in the observations. As seen in Fig. 5, $`{}_{}{}^{4}\mathrm{He}`$ and $`{}_{}{}^{7}\mathrm{Li}`$ are consistent with each other but favor a value of $`\eta `$ which is discrepant by $`2\sigma `$ from the value $`\eta _{10}=5.9\pm 0.4`$ indicated by the D abundance (16) alone. It is important to note that it is the latter which essentially provides the widely-quoted ‘precision’ determination of the baryon density: $`\mathrm{\Omega }_Bh^2=0.02\pm 0.002(95\%\mathrm{c}.\mathrm{l}.)`$ . Additional studies are required to clarify if this discrepancy is real.

The recent WMAP data on CMB fluctuations implies a baryon abundance of $`\eta _{10}=6.5_{0.3}^{+0.4}`$ if the primordial perturbations are assumed to have a power-law form ; the value decreases by $`9\%`$ if the spectral index is allowed to vary with scale, as is now indicated by the observations. Although broadly consistent with the determinations based on primordial abundances, it is clear that this exacerbates the tension with the $`{}_{}{}^{4}\mathrm{He}`$ and $`{}_{}{}^{7}\mathrm{Li}`$ measurements mentioned above. It has been proposed that the CMB determination of $`\eta _{10}`$ can be used as an input into BBN calculations . However it would be advisable to await further LSS data from SDSS  to pin down more precisely the primordial fluctuation spectrum, to which such determinations are sensitive. For example, allowing for a ‘step’ in the primordial spectrum at $`k0.05h`$ Mpc-1, as was indicated by data from the APM galaxy survey , can decrease the baryon density inferred from the WMAP+2DFGRS data by upto $`40\%`$ .

Given this situation, it is still neccessary to be conservative in evaluating bounds on $`N_\nu `$ from BBN. An analysis based on simple $`\chi ^2`$ statistics and taking into account the correlated uncertainties of the elemental yields, gives

$$2<N_\nu <4(95\%\mathrm{c}.\mathrm{l}.).$$
(19)
Tighter bounds can of course be obtained under less conservative assumptions, e.g. adopting the D abundance (16) requires $`N_\nu <3.2`$ . If true, such bounds rule out the possibility of ‘sterile’ neutrinos since these would be brought into equilibrium through mixing with the active neutrinos in all currently viable schemes . The WMAP determination of the baryon density seems to require $`N_\nu `$ less than 3  but is consistent with the Standard Model taking systematic uncertainties in the elemental abundances into account . Moreover as emphasized above, the WMAP determination of $`\eta `$ will have to ve lowered if the primordial spectrum is not scale-free and this will considerably ease the present tension. (The WMAP+2dFGRS data by themselves restrict $`1.8<N_\nu <5.7(95\%\mathrm{c}.\mathrm{l}.)`$ assuming a flat universe .)

The above limits hold for the standard BBN model and one can ask to what extent they can be modified if plausible changes are made to the model. For example, there may be an initial excess of electron neutrinos over antineutrinos, parametrised by a ‘chemical potential’ $`\xi _e`$ in the distribution function: $`f_\nu [\mathrm{exp}(p/T\xi _e)+1]^1`$. Then $`np`$ equilibrium is shifted in favour of less neutrons, thus reducing the $`{}_{}{}^{4}\mathrm{He}`$ abundance by a factor $`\mathrm{exp}(\xi _e)`$ . However the accompanying increase in the relativistic energy density speeds up the expansion rate and increases the $`n/p`$ ratio at freeze-out, leading to more $`{}_{}{}^{4}\mathrm{He}`$, although this effect is smaller. For neutrinos of other flavours which do not participate in nuclear reactions, only the latter effect was presumed to operate, allowing the possibility of balancing a small chemical potential in $`\nu _e`$ by a much larger chemical potential in $`\nu _{\mu ,\tau }`$, and thus substantially enlarging the concordance range of $`\eta _{10}`$ . However the recent recognition from Solar and atmospheric neutrino experiments that the different flavours are maximally mixed, no longer permits such a hierarchy of chemical potentials , thus ruling out this possible loophole. Consequently the relic neutrino abundance cannot be significantly different from the value (1). A small chemical potential in electron neutrinos of $`\xi _e<0.08`$ is however still possible and could be another explanation for why the $`{}_{}{}^{4}\mathrm{He}`$ abundance seems to be lower than the value that would be expected on the basis of the measurements of the D abundance or the CMB anisotropy.

Another possible change to standard BBN is to allow inhomogeneities in the baryon distribution, created e.g. during the QCD (de)confinement transition. If the characteristic inhomogeneity scale exceeds the neutron diffusion scale during BBN, then increasing the average value of $`\eta `$ increases the synthesised abundances such that the observational limits essentially rule out such inhomogeneities. However fluctuations in $`\eta `$ on smaller scales will result in neutrons escaping from the high density regions leading to spatial variations in the $`n/p`$ ratio which might allow the upper limit to $`\eta `$ to be raised substantially . Recent calculations show that D and $`{}_{}{}^{4}\mathrm{He}`$ can indeed be matched even when $`\eta `$ is raised by a factor of $`2`$ by suitably tuning the amplitude and scale of the fluctuations, but this results in unacceptable overproduction of $`{}_{}{}^{7}\mathrm{Li}`$ . A variant on the above possibility is to allow for regions of antimatter which annihilate during or even after BBN; however the $`{}_{}{}^{7}\mathrm{Li}`$ abundance again restricts the possibility of raising the limit on $`\eta `$ substantially . Finally the synthesised abundances can be altered if a relic massive particle decays during or after BBN generating electromagnetic and hadronic showers in the radiation-dominated plasma. Interestingly enough the processed yields of D, $`{}_{}{}^{4}\mathrm{He}`$ and $`{}_{}{}^{7}\mathrm{Li}`$ can then be made to match the observations even for a universe closed by baryons , however the production of 6Li is excessive and argues against this possibility.

In summary, standard BBN appears to be reasonably robust (although non-standard possibilities cannot be definitively ruled out) and consistent with $`N_\nu =3`$, leaving little room for new physics. Systematic uncertainties in the elemental abundance determinations need to be substantially reduced in order to make further progress.

## 7 Conclusions

The study of neutrinos in cosmology has had a long history. However as both theory and observations have improved, the fascinating possibilities that had been raised in early work have gradually been eliminated. In particular neutrino (hot) dark matter is no longer required by our present understanding of large-scale structure. Cosmological bounds on neutrino masses continue to be more restrictive than laboratory limits and all known species are now required to be sufficiently light that the rich phenomenology of unstable neutrinos is now largely of historical interest. Moreover there is no motivation from considerations of primordial nucleosynthesis for invoking new sterile neutrinos. Nevertheless such neutrinos can exist (e.g. in the ‘bulk’ ) and there may yet be surprises from the Big Bang in store for us.