In the early Universe, the three flavour active left-handed neutrinos να and their CP conjugated states were kept in thermodynamical equilibrium with e.m. plasma by CC and NC weak interactions with charged leptons and quarks. In this regime the neutrino distributions are just fixed by thermodynamics to the Fermi-Dirac type
By using such distributions we get at the equilibrium
As the temperature decreases, CC and NC weak interactions rate drop down faster than H(T) and hence become not efficient enough in mantaining neutrinos at equilibrium. At this point one says that neutrinos decouple from the e.m. plasma (baryons, leptons, photon). To describe neutrino decoupling in details, one should use the set of Boltzmann equations governing the evolution of their distribution function. Nevertheless in the following we will assume a simplified approach that provides quite precise results.
In order to simplify the Boltzmann equations ruling neutrino distribution functions one can consider the leading processes contributing to the collisional term, namely the scattering over electrons/positrons and pair conversions,
When the temperature is at most few MeV, we can consider the weak reaction rates as due to a 4-fermions interaction term (W and Z not propagating). In this case the corresponding cross section times velocity for charged and neutral current interactions is of the order of
Therefore the weak interaction rate for neutrinos becomes
ne denotes the electron/positron density, which is ≈ T3 as long as electrons are still relativistic. In that epoch the Universe was RD, hence the Hubble parameter was given by
The temperature for which the weak interaction rate, dropping down because of the expansion, reaches the value of H is called the Decoupling Temperature, hereafter denoted by TνD.
From the above definition the decoupling temperature reads
The decoupling temperature is weakly dependent on g* . For temperature of few MeV the relativistic d.o.f. give
One can get a better prediction for the decoupling temperature via the kinetic equation. The Boltzmann equations written in terms of the variables x=m a and y=p a, with m an arbitrary mass scale (usually taken ≈ me) read
During RD epoch H ≈ x-2, in this case it is easy to find that the solutions are function of the combination y/x3, or y T3. In particular one gets a decoupling temperature momentum dependent, namely TνeD= 2.7 y-1/3 and Tν(μ,τ)D= 4.5 y-1/3. By using the average momentum, <y> ≈ 3, one gets TνeD= 1.87 MeV and Tν(μ,τ)D= 3.12 MeV.
With instantaneous decoupling limit we denote the simplifying assumption of a neutrinos decoupling occurring instantaneously at TνD. In reality the decoupling takes place over an extended range of time. In particular more energetic neutrinos will be kept in equilibrium with e.m. plasma longer than low energy neutrinos. This means that even neutrinos to some extent get profit of the entropy release due to e± annihilations, and hence some thermal distortions will be imprinted in the neutrino distribution with respect to a standard Fermi-Dirac function. However these effects change the neutrino energy density at the percent level only, and thus will be neglected in the following.
After decoupling, neutrinos propagate freely, and their distribution only feels the effect of redshift of physical momentum. In the instantaneous decoupling approximation, we can assume the distribution at TνD to be exactly a Fermi-Dirac one
with aD the scale factor at TνD. From that moment on the distribution function satisfying the Liouville equation for ultrarelativistic particles is
which is the initial distribution with a temperature redshifted T = TνD aD/a . Note that this temperature scaling is equivalent to the conservation of the neutrino entropy per comoving volume, namely
The neutrino to photon temperature ratio after neutrino decoupling can be easily obtained by using the conservation of entropy for the e.m. plasma. As long as electron/positron pairs are relativistic, the photon temperature simply scales as a-1 since the number of degrees of freedom gs is constant. When the temperature drops below a value ≈ me, e.m. , the process cannot be efficiently compensated by the inverse pair production. This implies that electron/positron distributions start to be suppressed by the term e-me/T, and practically e± disappear, apart from a tiny relic electron density related by electric neutrality to the baryon number parameter ηB. In practice, the entropy release due to e± annihilations reheat photons, whose temperature is then not decreasing as a-1 for a while. Due to the smallness of weak reactions at this temperature, neutrinos practically get no profit of such reheating, nevertheless even though to a very small extent, such phenomenon leaves a footprint on neutrino distribution functions.
The entropy density of coupled e± and γ has the following expression
The previous expression strongly simplifies for a temperature much larger than me, in fact we get
whereas for T << me we simply have gs = 2. Since the entropy for the e.m. plasma is conserved inside a coomoving volume, namely
At neutrino decoupling T/Tν =1 and gs = 11/2, whereas for T << me we simply have gs = 2, hence the previous expression becomes
Using the known value of the CMB temperature today, T0= 2.725 K, neutrinos at present have a temperature Tν0= 1.945 K, or equivalently 1.7 10-4 eV in natural units.
As we have seen after the e+ - e- annihilations the density number of photons and neutrinos read
with gγ = 2, and
with gν = 6. By using the ratio T/Tν determined in the framework of instantaneous neutrino decoupling, we get nν = 9/11 nγ. This ratio allows to evaluate for the present value of photon number density, namely nγ0 = 410 cm-3, the corresponding neutrino number density, nν0 = 336 cm-3. Absolutely larger than the present aboundance of baryons, namely nBo ≈ 10-7 cm-3. Hence photons and neutrinos are the most abundant species in the present Universe. In particular since neutrinos at the present are non relativistic (their energy is ≈ 10-4 eV << mνi) we can write
Since Ων ≤ ΩDM ≈ 0.2, we get for h ≈ 0.73, ∑i mνi ≤ 10 eV.
If we consider negligible the neutrino chemical potentials, the number density of cosmic neutrinos for each flavour is given by
which leads to an extremely large flux with respect to other astrophysical neutrino sources, including solar neutrinos. By using the ratio between neutrino and photon temperatures, previously obtained, we can compute the total amount of radiation after e± annihilations which results
Such result holds if:
A convenient way to parameterize any deviation from these assumptions consists in defining the effective number of neutrino species Neff
Let us consider the neutrino distribution functions with with x=μ, τ
The Boltzmann Equations
can be casted in the following form
where x = a(t) me and y= a(t) |p|
In addition one must satisfies the covariant energy conservation
where the total energy density reads
whereas the total pressure reads
The change in the electromagnetic plasma equation of state can be evaluated by first considering the corrections induced on the e± and photon masses. They can be obtained perturbatively by computing the loop corrections to the self-energy of these particles. For the electron/positron mass, up to order α ≡ e2/4 π we find the additional finite temperature contribution
The renormalized photon mass in the electromagnetic plasma is instead given, up to order α, by
The previous corrections modify the corresponding dispersion relations as
This in turn affects the total pressure and the energy density of the electromagnetic plasma
Expanding P with respect to δme2 and δmγ2, one obtains the first order correction
and from the previous relation one gets ρint as well.
The presence of the additional contributions Pint and ρint modify the evolution equation for z≡ T a(t), which now reads
The functions K’(ω) and J’(ω) stand for the first derivativeof K(ω) and J(ω) with respect to their argument. Note that neutrinos affect the final value of z through the terms ,which are not vanishing only if neutrinos have a non thermal behaviour.
In case one neglects the finite temperature QED corrections, the functions G1(x/z) and G2(x/z) vanish, and one recovers a very well known expression reported in literature. Notice that the presence of G2(x/z) in the denominator of the r.h.s. of dz/dx makes, at least inprinciple, not correct to simply sum the neutrino contribution with the QED one.
The initial conditions are fixed by the thermodynamical equilibrium which is satisfied for T > 10 MeV
In particular we have
Asymptotically, as x → ∞, the prediction of the instantaneous decoupling approximation gives
zeq → (11/4)1/3
In the Figure we show the departure of z ≡ T/Tν from zeq for different values of the asymmetry parameters ξα
The asymptotic distribution functions are not perfectly thermodynamical distributions (out of equilibrium effects).
The effective number of neutrino species is used to parametrize the energy density stored in relativistic species, ρR, through the relation
where ργ is the energy density of photons, whose value today isknown from the measurement of the Cosmic Microwave Background (CMB) temperature.
The previous exprtession can be also written as
Finally, from the previous definitions it is straightforward to get the following expression
In the Table are reported the values assumed by the several factors when different effects are considered.
We now consider the possibility that, at the stage of neutrino decoupling, there are extra relativistic degrees of freedom, provided by some X field excitations, which are assumed to have a thermal distribution with some temperature TX. Their contribution ρX to the total energy density of the Universe can be parametrized in terms of an additional contribution in the effective number of neutrinos, as previously defined
and we have defined zX= TX a(t). The parameter gX depends on the spin (gX=1 for a real scalar, gX=7/4 for a Weyl spinor, etc.) as well as on the additional internal degrees of freedom of the X particles. Notice that if the X excitations have decoupled between μ± and e± annihilation phases we simply have NX= 4/7 gX . For an earlier decoupling we have instead NX< 4/7 gX.
The presence of ρX, apart from introducing a new contribution to , slightly affects the results previously obtained, namely the relative change of neutrino energy density induced by incomplete neutrino decoupling, as well as the asymptotic photon temperature zfin. In fact since their energy density increases the expansion rate of the Universe, we expect to decrease with growing NX, and the ratio to become closer to unity, since neutrino decouplingprocess starts at earlier times.
If we denote with and the new values for these parameters asfunctions of NX we therefore have
As it is seems more and more clear from neutrino experiments on solar and atmospheric neutrinos, it is unlikely that neutrino mass differences may be greater than ≈ 0.1 eV, unless we enlarge the standard scenario introducing sterile neutrino states. At the same time, it is also quite well established from Tritium decay data that νe mass is bound to be smaller or, at most, of the order of 1 eV. It is therefore clear that in this scenario all neutrino masses are completely negligible as far as their decoupling is concerned. However if their values is as large as ≈ 1 eV, they start to be relevant as the temperature approaches the range relevant for CMB, and the presence of a finite mass modifies of course the neutrino contribution to . It is interesting to consider how the effects of incomplete decoupling and QED thermal effects just studied would now affect the neutrino energy density. This is conveniently parametrized by the (time dependent) quantity
Since for mα ≈ 1 eV active neutrino are fully elativistic for temperatures of the order of MeV, so we can still take at the e± annihilation phase the results previously obtained, it is easy to see that is given by
where we have
6. Dark Matter
11. Cosmic Rays - I
12. Cosmic Rays - II