The distribution function $f(\mathbf{x}, \mathbf{p}, t)$ gives the number of particles in the universe with position in $(\mathbf{x}, \mathbf{x}+\delta \mathbf{x})$ and momentum in $(\mathbf{p}, \mathbf{p}+\delta \mathbf{p})$ at time $t$. It satisfies the boundary condition that $f \rightarrow 0$ as $|\mathbf{x}| \rightarrow \infty$ and as $|\mathbf{p}| \rightarrow \infty$. Its evolution obeys the Boltzmann equation

$$\frac{\partial f}{\partial t}+\frac{\partial f}{\partial \mathbf{p}} \cdot \frac{d \mathbf{p}}{d t}+\frac{\partial f}{\partial \mathbf{x}} \cdot \frac{d \mathbf{x}}{d t}=\left[\frac{d f}{d t}\right]_{\mathrm{col}}$$

where the collision term $\left[\frac{d f}{d t}\right]_{c o l}$ describes any particle production and annihilation that occurs.

The universe expands isotropically and homogeneously with expansion scale factor $a(t)$, so the momenta evolve isotropically with magnitude $p \propto a^{-1}$. Show that the Boltzmann equation simplifies to

$$\frac{\partial f}{\partial t}-\frac{\dot{a}}{a} \mathbf{p} \cdot \frac{\partial f}{\partial \mathbf{p}}=\left[\frac{d f}{d t}\right]_{\mathrm{col}}$$

The number densities $n$ of particles and $\bar{n}$ of antiparticles are defined in terms of their distribution functions $f$ and $\bar{f}$, and momenta $p$ and $\bar{p}$, by

$$n=\int_{0}^{\infty} f 4 \pi p^{2} d p \quad \text { and } \quad \bar{n}=\int_{0}^{\infty} \bar{f} 4 \pi \bar{p}^{2} d \bar{p}$$

and the collision term may be assumed to be of the form

$$\left[\frac{d f}{d t}\right]_{\mathrm{col}}=-\langle\sigma v\rangle \int_{0}^{\infty} \bar{f} f 4 \pi \bar{p}^{2} d \bar{p}+R$$

where $\langle\sigma v\rangle$ determines the annihilation cross-section of particles by antiparticles and $R$ is the production rate of particles.

By integrating equation $(*)$ with respect to the momentum $\mathbf{p}$ and assuming that $\langle\sigma v\rangle$ is a constant, show that

$$\frac{d n}{d t}+3 \frac{\dot{a}}{a} n=-\langle\sigma v\rangle n \bar{n}+Q$$

where $Q=\int_{0}^{\infty} R 4 \pi p^{2} d p$. Assuming the same production rate $R$ for antiparticles, write down the corresponding equation satisfied by $\bar{n}$ and show that

$$(n-\bar{n}) a^{3}=\mathrm{constant}$$