What is the numerical relation between energy density ( $\Omega_\text{s}, \Omega_\alpha$), and the number density of neutrinos (sterile - $n_s$ , active - $n_\alpha$)?

Background info that might help better understand my question:

This is based on an equation derived through the relationship of sterile neutrinos distribution function and the active neutrinos distribution function:

$$f_s = k \sin^2 2\theta f_\alpha \tag{1}$$

integrating in terms of momentum gives:

$$\int \frac{d^3 p}{(2\pi)^3} f_s= \int \frac{d^3 p}{(2\pi)^3} k \sin^2 2\theta f_\alpha \tag{2}$$

where $k$ is a numerical constant, $f_i$ is the distribution function of $i$ neutrinos and $\theta$ is the collision angle of neutrinos.

Based on the definition of number density, $n_i= \int \frac{d^3 p}{(2\pi)^3} f_i$ , we get:

$$\frac{n_s}{n_\alpha}= \gamma \tag{3}$$

where $\gamma = \frac{1}{n_\alpha}\int \frac{d^3 p}{(2\pi)^3} \sin^2 2\theta f_\alpha$

**I'm trying to connect this last equation to $\Omega$ but I need to understand the relation in my question first. **

I've seen in the following article the following statement, based on $(4)$ below:

$$f_s = k f_\alpha \tag{4}$$

$f_s$ has the same functional form as $f_\alpha$ and therefore $\Omega_s/ \Omega_\alpha = (m_s/m_\alpha) (f_s/f_\alpha)$.

I still don't understand where this equation was deduced from.


If I'm not mistaken, the quantities $\Omega$ represent the energy density of the corresponding particle relative to the critical density $\rho_\text{crit} = 3 H_0^2 m_\text{pl}^2$. Thus, to understand the claim, let's write down the number density $$n_i = \int\frac{{\rm d}^3p}{(2\pi)^3}f_i\,.$$ Once whatever process is generating the particles shuts off, the number of species $i$ per co-moving volume is fixed as long as it doesn't decay. When the universe cools below the mass of the species, it becomes non-relativistic, and its energy denisty scales like $m_i n_i$. Therefore, $$\Omega_i = \frac{m_i n_i}{\rho_\text{crit}}\,.$$ Now let's take the ratio of species for which we know $f_a\propto f_s$. Suppose that $\alpha \equiv f_a/f_s$ is the constant of proportionality. Then $$\frac{n_a}{n_s} = \frac{\int\frac{{\rm d}^3 p}{(2\pi)^3}f_a}{\int\frac{{\rm d}^3 p}{(2\pi)^3}f_s} = \frac{\alpha\int\frac{{\rm d}^3 p}{(2\pi)^3}f_s}{\int\frac{{\rm d}^3 p}{(2\pi)^3}f_s} = \alpha\,.$$ Hence $$\frac{\Omega_a}{\Omega_s} = \frac{m_a n_a}{m_s n_s} = \frac{m_a}{m_s}\alpha = \frac{m_a f_a}{m_s f_s}\,.$$


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.