La respuesta es la misma razón por la que un vaso de agua a temperatura ambiente, se evapora. Aunque la mayoría de las partículas va a estar por debajo del punto de ebullición, el equilibrio que uno espera no es del todo en la fase líquida. La ocasional particular energía de la molécula de agua se evapora, así como el ocasional neutral átomo de hidrógeno será golpeado por una particular energía de los fotones y se ionizan.
En lugar de una brusca transición que ocurre cuando la temperatura es igual al potencial de ionización, una mejor aproximación se obtiene con la ecuación de Saha, que explica el hecho de que las distribuciones de partículas que tienen un rango de energías.
Contabilidad para la tirada, degeneraciones correctamente (obligado de hidrógeno tiene un factor de $2$ más de los estados disponibles de aislados de protones), el Saha ecuación se convierte en
$$ \frac{n_\mathrm{e}n_\mathrm{p}}{n_\mathrm{H}} = \frac{(2\pi m_\mathrm{e}kT)^{3/2}}{h^3} \mathrm{e}^{-E_\mathrm{bind}/kT}. $$
Aquí $\mathrm{e}$, $\mathrm{p}$, y $\mathrm{H}$ stand para liberar electrones, protones libres, y el hidrógeno neutro.
Debido a que el universo es neutral, $n_\mathrm{e} = n_\mathrm{p}$. Supongamos que queremos una fracción $f$ de ionización: $n_\mathrm{p}/n_\mathrm{H} = f$. Si el universo tiene una densidad de $\rho$,$n_\mathrm{p} + n_\mathrm{H} = \rho/m_\mathrm{p}$. La combinación de todo esto, el lado izquierdo puede escribirse
$$ \frac{n_\mathrm{e}n_\mathrm{p}}{n_\mathrm{H}} = \frac{f^2}{f+1} \cdot \frac{\rho}{m_\mathrm{p}}. $$
Ahora tenemos sólo dos incógnitas en la ecuación, $T$$\rho$. Sin embargo, la previsible expansión del universo nos dice
$$ T = T_0 (1+z) $$
para $T_0 = 2.725\ \mathrm{K}$, y que
$$ \rho = \Omega_{\mathrm{b},0} \rho_\mathrm{crit} (1+z)^3 $$
para $\Omega_{\mathrm{b},0} = 0.04628$. Aquí la densidad crítica es
$$ \rho_\mathrm{crit} = \frac{3H_0^2}{8\pi G} $$
para $H_0 = 69.32\ \mathrm{km/s/Mpc}$.
Poner todo junto, podemos reescribir la ecuación de Saha
$$ \frac{f^2}{f+1} \cdot \frac{\Omega_{\mathrm{b},0}\rho_\mathrm{crit}}{m_\mathrm{p}} \left(\frac{T}{T_0}\right)^3 = \frac{(2\pi m_\mathrm{e}kT)^{3/2}}{h^3} \mathrm{e}^{-E_\mathrm{bind}/kT}. $$
La solución de los rendimientos $T = 3660\ \mathrm{K}$.
Nota: este análisis se hace algunas aproximaciones, tales como el olvidar de helio.