Thermodynamics – Why Do Lighter Isotopes Evaporate Faster Than Heavier Isotopes?

climate-sciencenuclear-physicspotential energystatistical mechanicsthermodynamics

I have just read that we are able to estimate the temperature of the earth thousands of years ago by measuring the ratio of certain isotopes present in ice cores that froze thousands of years ago.
My understanding of this idea is that $\mathrm{H_2O}$ molecules containing lighter isotopes (i.e. with $\mathrm{^{16}O}$) evaporate quicker than those molecules that contain heavier isotopes (i.e. with $\mathrm{^{18}O}$) and that this difference in evaporation rate exhibits a temperature-dependent relationship.

My problem is that the equipartition theorem says that each molecule should contain an energy of $\frac{1}{2}kT$ per degree of freedom regardless of the molecular mass. If we assume that all the liquid $\mathrm{H_2O}$ molecules roughly obey a Maxwell-Boltzmann distribution, then the heavier isotopes will be moving slower than the lighter isotopes (on average at a given temperature) but their greater mass will mean that all isotopes will still possess the same average kinetic energy (in accordance with the equipartition theorem). So if all the water molecules possess the same kinetic energy on average (regardless of the isotopes they contain), why do those molecules with lighter mass preferentially evaporate? I understand that they are moving quicker than their heavier counterparts but the determination of whether a particle escapes a potential well and evaporates is its kinetic energy not its velocity.

In effect, I am asking why the speed of the molecules determines their proclivity to evaporate (to escape the potential well they inhabit due to intermolecular attraction) when it should be their kinetic energy that determines whether they escape their intermolecular potential well or not. If the evaporation rate is determined by kinetic energy (which I think it should be?) as opposed to velocity, then there shouldn't be any difference in the evaporation rates of different water molecules containing different isotopes because they all contain the same kinetic energy on average.

Best Answer

You're right, the difference in evaporation rate has little to do with the speed of the molecules. The correct explanation is that hydrogen bonding is slightly stronger in heavy water than in normal water. As a result, heavy water molecules are slightly more attracted to one another than molecules of light water, and it is harder for a molecule of heavy water to evaporate.

Why is hydrogen bonding stronger in heavy water? To get some intuition, let's consider the quantum harmonic oscillator. If you solve the Schrödinger equation with a potential energy function of $V(x) = \frac12 kx^2$, you get solutions with energy $\hbar \omega (n + \frac12)$ where $n$ is an integer and $\omega = \sqrt{\frac{k}{m_r}}$. Here $m_r$ is the reduced mass $\frac{m_1m_2}{m_1+m_2}$. Thus, a quantum harmonic oscillator with higher masses at the endpoints (heavy water molecules) will have lower $\omega$ and a lower ground state energy $\hbar \omega/2$. This means that more energy is needed to excite the bonding electrons from the ground state to the dissociation energy, so the bond should be slightly stronger.

I should mention that this explanation is rather oversimplified in a lot of important ways. In fact, surprisingly many important questions about hydrogen bonding in water are still not well understood. It is known that the hydrogen bond in a heavy water dimer is stronger than the one in a light water dimer (as we'd expect based on the argument above). There is some theoretical evidence that this conclusion should generalize to similar small systems. However, it's possible that the relevant physics may be different (and unfortunately, less elegant) for larger amounts of water. One intriguing experimental study found that the hydrogen bond in light water is actually shorter than the hydrogen bond in heavy water, which would imply a stronger bond. The study also reported that the number of hydrogen bonds per water molecule drops from $3.76 \pm 0.1$ in heavy water to $3.62 \pm 0.1$ in light water because the geometric packing of the molecules is slightly different. Perhaps this effect may be more important than the strength of the individual bonds.

Related Question