Last week someone (“Cosm”) posted an interesting comment/question on my post “What is the redshift of the Cosmic Microwave Background (CMB)?” . The question asked how we can say that the temperature of the Universe at recombination is about 3000K when the energy of electrons with such a temperature would be
and the ionisation potential of hydrogen is 13.6 eV. This would imply, naively, that recombination would occur at a higher temperature, roughly . But, it does not occur at a temperature of just over 100,000 K, but at about 3000 K. Why is this?
I decided that it was a sufficiently interesting point that I would go through the detail of why recombination (or decoupling as I prefer to call it) happened when the Universe was about 3,000K; even though electrons cease to be able to thermally ionise hydrogen at a much higher temperature of about 100,000K, so one might think that decoupling would happen earlier. The clue is that it has to do with something called equilibrium theory.
The temperature of decoupling based on equilibrium theory
One can get a reasonable estimate of the temperature (and redshift) when recombination (decoupling) took place by using equilibrium theory. This was developed in the 1920s and is based on the Saha ionisation equation, named after Indian astrophysicist Meghnad Saha.
We are going to assume that the process
is the dominant reaction for creating neutral hydrogen. That is, electrons and bare protons combining. When a free electron combines with a proton to create neutral hydrogen it will also create a photon, and that is what the represents (it is not a gamma-ray photon, it is more likely to be a UV or visible-light photon). But, it is a two way process, as obviously electrons can also be ionised to go back to free electrons and bare protons.
We are going to denote the number density of free electrons (the number per unit volume) as , the number density of free protons (ionised hydrogen) as , and the number density of neutral hydrogen as .
We can calculate the relative abundance of the free electrons to protons and neutral hydrogen via the Saha equation
Where is the mass of an electron, is Boltzmann’s constant, is the temperature, is the reduced Planck constant (which is defined as where is Planck’s constant), and is the ionisation potential of hydrogen (13.6 eV). We know from charge neutrality that . We are going to define as the fraction of free electrons, so by definition
Because we can re-write Equation (1) as
and from (2) we can write
We can also write
Dividing Eq. (4) by Eq. (5) we have
But
so we can re-write Eq. (6) as
Most of the quantities on the right hand side are physical constants, and the rest are known functions of redshift . For example, the temperature history of the CMB at any redshift is just given by
where is the CMB’s current temperature. The total number density of hydrogen (neutral and ionised), is also a function of redshift, and is given by
where is the currently observed density. So, we will re-write Eq. (7) as
Equation (10) cannot be solved analytically, only numerically.
We will assume that we need a fraction of free electrons of for decoupling to have occurred (in other words, of the electrons have bound to protons to form neutral hydrogen). We will try different temperatures to see what fraction they give. Note: when we assume a particular temperature, this will fix the redshift, because of Eq. (8).
Let us first try a temperature of . From Eq. (8) this gives . Plugging these values of and into Eq. (10) gives
Solving this quadratic equation gives
This is obviously much less than , by the time the Universe has cooled to 3000K there are very few free electrons, decoupling is essentially complete.
What if we use ? Doing the same thing we find which then gives
Solving this quadratic gives . This means that of the electrons are free when the temperature is 4000K, meaning of the hydrogen atoms have become neutral.
Thirdly, let us try 3800K. This gives . Plugging this into Eq. (9) gives
Solving this gives
At a temperature of T=3800K just under of the electrons are free, so some of the hydrogen atoms are neutral.
So, we can surmise from this that the Universe had decoupled enough to be about transparent to radiation by the time the temperature was a little below . This is why the temperature of decoupling (recombination) is usually given as lying between a temperature of .
Why is this just an approximation?
As I mentioned above, this is just an approximation, but not a bad one. It’s an approximation because it is a simplification of what is actually going on. In 1968, Jim Peebles (who had done the work in the 1965 Dicke etal. paper) and, independently, Yakov Zel’dovich (he of the Sunyaev-Zel’dovich effect) in the USSR, worked out a more complete theory where the hydrogen has three energy levels, rather than what we have done here where we assume the free electrons go straight into the ground-state (i.e. only two energy levels). Their third level was the n=2 state, the energy level just above the ground state. This is an important energy level for hydrogen, as it is transitions down to the n=2 level which give rise to visible-light photons. Using this more complicated 3-level model gives a Universe which is 90% neutral at , which would correspond to a temperature of T=3000K, which is why this temperature is most often quoted. You can read more about their 3-level model here.
Summary
Using equilibrium theory, which is an oversimplification, gives the following fractions of neutral hydrogen for three different temperatures
- At T=3000K the Universe would have been more than 99% neutral
- At T=4000K the Universe would have been about 40% neutral
- At T=3800K the Universe would have been about 70% neutral
Using a more correct 3-level model developed by Peebles and, independently, by Zel’dovich, gives that the Universe would have been about 90% neutral by the time the temperature had dropped to T=3000K. It is this temperature which is usually quoted when we talk about the temperature of the Universe when recombination (decoupling) occurred.
Excellent article and information! Thank you!
Thanks, this is explanation is nice and clear.
Good. That was my hope! 😉
Your present day hydrogen number density is too high by an order of magnitude. The baryon density is 4.6% of the critical density. Did you mistakenly use the matter density instead? I make it 0.19 H atoms per cubic metre.
Yeah, you’re off by a factor of 10 on the present day density of H. 0.14 to 0.13 atoms per cubic m is a more accurate value.