In general, the primordial plasma does not satisfies the equilibrium conditions, hence we cannot exploit the powerful methods of equilibrium thermodynamics. This means that we cannot assume that a gas of particles be distributed according to Bose-Einstein or Fermi-Dirac functions. If this occurs, one neds to obtain the time evolution of the statistical system in terms of the kinetic theory, namely solving the Boltzmann equation.
The approach to equilibrium for a system made of photons, neutrinos, charged leptons, etc. is driven by interactions processes such as scatterings (crucial to redistribute particle momenta and to achieve the so-called kinetic equilibrium) and interactions where the number of a particular particle specie can change. Since the Universe is expanding the rate of microscopic interactions which tries to thermalize the system must compel with the expansion rate, measured by the Hubble parameter H, which diluting the plasma tries to freeze the species.
If we consider for simplicity the case of two body processes of the form a +b → c+d, the interaction rate for the species a is then given by
where σ stands for the typical cross section, v is the a-b relative velocity and nb the b number density. Hence the equilibrium is maintained if
When for a specie Γa ≤ H we say that it is decoupled. Typically, at high temperature density and mean energy are large enough to ensure that the equilibrium condition be fulfilled. As the expansion proceeds, if the interaction rate decreases faster than H for some value of redshift z=zD
This defines the so-called decoupling redshift. After a specie has decoupled it is called a relic. If the decoupling takes place when the specie is still relativistic it is called a hot relic, like neutrinos. On the other hand, if it decouples when it is already non-relativistic it is defined a cold-relics), like for example cold dark matter (CDM).
In general a specie is characterized by a distribution function in the phase-space f(xμ,Pμ). In absence of particle interactions, the value of f cannot change as we follow a point in phase space along a geodesic motion f(xμ(λ),Pμ(λ)), namely
which defines the Liouville operator as
In FRW the universe is homogeneous and isotropic, hence the distribution function depends on time and either energy E or the modulus of linear momentum P (via the mass shell condition). For this reason we get the following equation for f(t,P)
Changing the variables, namely using the physical momentum p=P a we have
Finally, if we introduce y ≡ (Σi Pi2)1/2, namely the modulus of the comoving momentum, we get
Thus any function depending on the comoving momentum only is a solution of the collisionless Boltzmann equation. Such equation is modified in presence of particle scatterings, annihilations etc. .Using physical momentum one can write a set of equations, one for each specie,
by fi we denote the distributions of all species interacting with the one we are interested in.
The right hand side of the previous equation, namely C, defines the so-called collisional integral, which depends on fj and fi due to the hypothesis of molecular chaos, namely that particle momenta are uncorrelated with their positions.
Let us consider the case of two-particle into two-particle interactions, a + b ↔ c + d. For the species “a” we have the collisional integral for two body process given by
Note that the factors (1 ± f) are due to the Pauli blocking effect for fermions (minus sign), and stimulated emission for bosons (positive sign). If the system is very diluted and the particle chemical potentials are very small, these factors reduce to unity, as in the case of particles obeying classical statistics. It is possible to prove that by denoting with Q a quantity conserved by interactions (like energy or momentum, electric charge, etc.), then
If we have two particle species, namely a and b which can only interact via the scattering process a + b → a + b both na and nb are conserved. In this case taking a general conserved charge qa (or analogously qb) one still gets that the collisional integral for the a particle vanishes and qa na is covariantly conserved.
Among the possible quantity one can choose Q=E, in this case summing over all species which are mutually interacting, we get
that expresses the covariant conservation of the total energy-stress tensor (namely summed over all interacting species).
At this point one may wonder if the equilibrium distributions satisfy the Boltzmann equation for an expanding Universe. Effectively one can prove that also in this case the equilibrium distribution functions,namely the ones which implies C=0 take the standard form
whenever all species are at the same temperature (kinetic equilibrium), and the chemical potentials satisfy the condition (chemical equilibrium)
To yield a vanishing C is a necessary but not sufficient condition for a distribution function to be an equilibrium distribution. In fact one must prove that it is also in the kernel of the Liouville operator. Let us write a generic equilibrium function in the form
with ξ=μ/T and β=1/T.
One can easily prove that the previous expression is in the kernel of Liouville operator if ot satisfies the equation
For massless particles, namely E=p, one can find an exact solution with
This means that if some specie decouples when relativistic (hot relic) it will keep a Fermi-Dirac or Bose-Einstein distribution with the decreasing of the temperature also after decoupling. In fact, untill it is in thermal contact with other species it keeps the equilibrium distribution, which will also be a solution of the collisionless Boltzmann equation after decoupling, being a function of the comoving momentum only.
For non-relativistic particles (m ≠ 0) the Liouville relation does not admits trivial solutions. In particular we can state that in FRW there are no exact equilibrium solution to the Boltzmann equation for massive particles . Nevertheless the distribution can approach closely the equilibrium distribution if interaction rates are sufficiently larger than H.
For values of temperature much larger than the particle mass, T >>m, one can neglect m and if for a specie strongly coupled its distribution is a Fermi-Dirac or Bose-Einstein with T ≈ a-1. When T decreases down to T ≈ m, the equilibrium distribution is no more a solution of the Boltzmann equation, but one can look for an approximate equilibrium form
con |δ(p,t)|<<1. For large interactions rate Γ >> H, and δ(p,t) ≈ H/Γ <<1 .
As a relevant example, let us consider e- and e+ which interact quite efficiently with γ via Compton scattering and pair processes e+ e- ↔ γ γ. When the temperature is much smaller than the electron mass, the Compton cross section is given by the Thomson formula, namely σT = 8 π e4/(3 me2) ≈ 0.67 10-24 cm2 . Photons have a distribution function, satisfying the Boltzmann equation, and whose temperature scales as T ≈ a-1=1+z. Hence when T = me , z = me/T0 ≈ 1010 . From this value of redshift and till the equality point, zeq ≈ 103 the universe is radiation dominated.
Neglecting the photon reheating due to e± annihilations, we have that the number of photons reads
where ζ(3) = 1.202. From this we get
Substituting the numerical values we get
Hence e± are distributed with an almost exact equilibrium Fermi-Dirac function. Since also pair processes have a large rate compared with H, chemical equilibrium holds, and reminding that photons have zero chemical potential, we can write
Hence we can write fot T ≤ me
In many relevant cases one can safely assume that kinetic equilibrium holds due to fast scattering processes. In these cases it is convenient to integrate the Boltzmann to get the evolution of the number densities. Let us consider the example of two body – two body processes. In this case we have
We have obtained the previous result considering the particular case where all particles involved are either non–relativistic or, in case they are relativistic, with small (or vanishing as for photons) chemical potential. Fot the species considered the number density is a simple function of the chemical potential
In this case the energy conservation can be cast in the following form
which represents the product of the cross section σ times the relative velocity for the process a+b → c+d, averaged over the particle distributions. Note that on the r.h.s. of the Eq. for d/da (na a3) we have the ratio between Γa = <σ|v|>nb and the Hubble parameter.If this ratio is large the speci “a” is in chemical equilibrium. The chemical equilibrium condition can be cast in the following form
which is called the Saha equation.
This equation allows to get the order of magnitude of the time (or redshift) at which some important events take place, like the recombination stage. Such event occurs when electrons and protons recombine to form neutral hydrogen, through the two body process
or to compute the instant at which Deuterium starts forming (the onset of primordial nucleosynthesis) due to proton neutron fusion
In equilibrium conditions one can easily compute energy and pressure for a given species “i“. Summing over particles and antiparticles we have
gi is the number of internal degrees of freedom as helicity, colour, etc. In the previous relations we have assumed opposite chemical potentials for particles and antiparticles. This is true if these particles are in chemical equilibrium with photons, by direct interactions or via processes like
For relativistic particles with ξi = 0 one has
If ξi ≠ 0 and we are always in the relativistic case one gets
In this limit, the particle–antiparticle asymmetry of a fermionic specie takes the exact form
whereas one can only provide an expansion for the number density in terms of powers of ξi , namely up to the first order
For relativistic bosons, since the chemical potential cannot be greater than the particle mass (otherwise when ξi approaches the value of the particle mass a Bose condensate develops) we have
where the last result only applies in case μi << mi and no Bose condensate forms.
For non-relativistic particles and μi << mi the distribution is given by the Maxwell-Boltzmann function
so we have
For particle with positive chemical potential the density of their antiparticles is smaller by a factor
For a homogeneous system, such as fluids in a FRW Universe, the entropy density, s=S/V is obtained from the gran potential Ω = -P V as
Integrating by parts one then obtains
The total entropy density due to several relativistic species can be written in terms of the contribution of one polarization state for the photons, multiplied by an effective number of “entropy relativistic degrees of freedom” gs
If we have vanishing chemical potentials we get
where we have assumed that different species can have different temperatures Ti with respect to the photon temperature T. The quantity gs is in general a function of T. When the temperature T reaches the mass threshold mi , the i-th species becomes non-relativistic and its contribution becomes negligible. In this situation gs decreases, and all species which at this stage are still in thermal equilibrium with “i” particles get reheated due to the entropy released by the rapid “i”- anti”i” annihilations. In this phase from the entropy conservation we get
6. Dark Matter
11. Cosmic Rays - I
12. Cosmic Rays - II