The interview season for the "Theoretical and Mathematical Physics" master programme at LMU is approaching quickly. We have to come up with new questions and problems that help us judge our applicants.
It turns out to be easy to find questions in quantum mechanics and those easily lead over to mathematics questions. However, we were always short on good stat mech problems. One possibility is to have an easy start with the harmonic oscillator and then couple that to a heat bath and compute the partition function (with geometric series featuring).
But this time, we thought we could vary this a bit and came to a surprising realization: Hydrogen is unstable! This was news to me but google finds a number of pages where this is discussed. Often wrongly, but the good explanation is in a 2001 paper by Miranda.
The idea is the following: Everybody knows that the energy of the -th level of the hydrogen atom has energy proportional to . This level has degeneracy since runs from 1 to and then runs from to . So the partition function is . First, we thought that this might be a function named after some 19th century mathematician but mathematica told us its name is actually since the exponent quickly approaches 1 for every positive .
The conclusion seems to be that there is something wrong with the hydrogen atom. And we have not even started to consider the positive energy scattering states. Obviously, this problem has an IR divergence and it is probably better to embed it in some cavity of finite radius. But still, you would think that then for a large cavity, most of the statistical weight would be in the highly excited states and the probability to be in the ground state would go to zero as the cavity gets larger.
The conclusion would be that a hydrogen atom at any temperature would almost never be in its ground state but always highly excited or even ionized. And all this only because the density of states diverges at 0. This looks like a situations worse than the Hagedorn transition that strings experience due to the exponentially growing density of states.
The solution in the above mentioned paper is quite simple: Rather than these scaling arguments one should put in some numbers! Let us start with the Bohr radius, which is m and the radius grows like . This means in ameter sized cavity we can only fit states up to roughly . However, at room temperature, Boltzmann exponent and . Thus, to balance the Boltzmann suppression of the higher levels compared to the ground state one has to take into account at the order of states and not just the first . Or put differently, one should use and exponentially large cavity. Otherwise the partition function is essentially cut off at and the probablility to find the ground state is very very very close to 1.