Isotopes of hydrogen — protium, deuterium, tritium: Visual walkthrough

The parent note told you that a C–H bond breaks 6–7 times faster than a C–D bond, even though deuterium is chemically "the same" as protium. That is the strangest claim in the whole chapter: identical chemistry, but one reacts several times slower. This page builds that number — the kinetic isotope effect — one picture at a time, starting from what a bond even is.

We only need one big idea: a chemical bond behaves like a mass hanging on a spring, and a heavier mass on the same spring wobbles more slowly. Everything follows from that.

Step 1 — A bond is a ball on a spring

WHAT. Picture two atoms joined by a chemical bond. Pull them apart slightly and they snap back; push them together and they push out. That "snap back to the middle" behaviour is exactly what a spring does.

WHY this picture. We need a model simple enough to do maths on. The spring (called a harmonic oscillator) is the simplest object that has a natural resting length and pushes back when disturbed — which is all a stiff bond does for small stretches. We borrow it from the physics of moving particles because it is the one model where we can compute everything exactly.

A first simplification. Carbon is about twelve times heavier than hydrogen, so in this first picture we pin it in place and let only the light H atom slide. This is a deliberate approximation to get the idea across — we will relax it in Step 2 and let carbon move a little too.

PICTURE. In the figure, a heavy carbon atom (violet, on the left, treated for now as a fixed wall) is joined by a coiled spring to a light hydrogen atom (orange). The bold arrow shows the direction the H atom is free to slide: toward and away from carbon.

Step 2 — Let carbon move too: the reduced mass $μ$

WHAT. Step 1 pinned carbon as a wall — but in a real molecule both atoms recoil against the spring, carbon just moves less because it is heavy. To keep the honest "both move" picture while still using the simple one-ball spring formula, physics replaces the two real masses with a single effective mass called the reduced mass $μ$ .

$μ = \frac{m _{1} m _{2}}{m _{1} + m _{2}}$

$m_{1}$ — mass of the first atom (the carbon), doing its job in the top of the fraction and the bottom.
$m_{2}$ — mass of the second atom (the H or D), same.
$μ$ — the single "pretend" mass we hang on the spring instead of the two real ones.

WHY this tool, and how it connects to Step 1. A spring formula needs one mass, not two. The reduced mass is the honest way to collapse a two-body wobble into a one-body wobble. Notice it contains the Step 1 approximation: when $m_{1}$ (carbon) is huge, $μ \to m_{2}$ — the fixed-wall picture is just the limit of the reduced mass when the wall is infinitely heavy. So Step 1 was not wrong, it was the $μ$ formula with carbon taken as immovable; Step 2 puts the small correction back.

PICTURE. The figure shows the two-ball system on the left (both free to recoil) collapsing into an equivalent single-ball system on the right, with $μ$ labelled on the lone ball.

Step 3 — How fast the ball wobbles: the vibration frequency $ν$

WHAT. A ball of mass $μ$ on a spring of stiffness $k$ wobbles back and forth at a definite rate. That rate is the vibrational frequency $ν$ :

$ν = \frac{1}{2 π} \frac{k}{μ}$

$k$ — spring stiffness (top of the root): stiffer spring → faster wobble.
$μ$ — reduced mass (bottom of the root): heavier ball → slower wobble.
$ν$ — wobbles per second.

WHY the square root, and why $μ$ downstairs. This is the one place a tool enters that we must justify. Why does mass sit under a square root? Because a heavier ball has more inertia — it accelerates less for the same spring push — so it takes longer to complete a swing. The square root is exactly how Newton's $F = ma$ balances against a spring's linear pull-back; you can't get any other power. The key consequence: since $k$ is identical for H and D, the only thing changing $ν$ is $μ$ .

PICTURE. Two springs side by side: the light H ball (orange) blurred into a fast wide wobble, the heavy D ball (magenta) a slow narrow wobble. Same spring, different speeds.

$\frac{ν _{D}}{ν _{H}} = \frac{μ _{H}}{μ _{D}} = \frac{1}{2} = \frac{1}{2} \approx 0.707$

The D bond wobbles about 30% slower — the ratio flips (D on top, but $μ_{D}$ downstairs) so the heavier ball comes out slower.

Step 4 — The ball can never fully stop: zero-point energy

WHAT. Classically, at absolute zero temperature the ball should sit perfectly still at the bottom of the spring's energy valley, with zero energy. Quantum mechanics forbids this. The lowest energy the ball is allowed to have is not zero but

$E_{0} = \frac{1}{2} h ν$

$h$ — Planck's constant, the fixed conversion between frequency and energy.
$ν$ — the wobble rate from Step 3.
$E_{0}$ — the zero-point energy: the leftover wobble that can never be removed.
$\frac{1}{2}$ — the specific floor set by quantum mechanics for a spring.

WHY it can't be zero. If the ball sat perfectly still at a perfectly known spot, we'd know both its position and its (zero) speed exactly — which quantum mechanics does not allow. So the ball must keep a minimum jiggle. This is the crucial physics: the bond always keeps a floor of energy, and that floor depends on $ν$ , which depends on mass. Note this floor is temperature-independent — it survives even at absolute zero, which is exactly why the isotope effect does not vanish as things get cold (we return to this in Step 7).

PICTURE. The energy valley (a U-shaped curve) with two horizontal floor lines drawn inside it: the H floor (orange) sits higher up the valley walls, the D floor (magenta) sits lower. Both are above the very bottom.

Combining Steps 3 and 4:

$\frac{E _{0, H}}{E _{0, D}} = \frac{ν _{H}}{ν _{D}} = \frac{μ _{D}}{μ _{H}} = 2 \approx 1.41$

The H bond keeps $2$ times more leftover energy than the D bond. The heavier ball settles into a deeper, quieter hole.

Step 5 — Why the deeper hole means slower reaction

WHAT. To break a bond, the atoms must climb out of the energy valley up to a summit called the transition state, where the bond is half-broken. The height they must climb is the activation energy $E_{a}$ — the "push" a reaction needs.

WHY this matters for isotopes. Both isotopes climb to almost the same summit, and here is the mechanistic reason. At the summit the C–H bond being broken is stretched nearly to breaking, so its spring stiffness $k$ has collapsed toward zero. By Step 3, $ν \propto k$ , so as $k \to 0$ the vibration frequency along the breaking bond $\to 0$ , and by Step 4 its zero-point energy $\frac{1}{2} h ν \to 0$ for both isotopes. When an energy is near zero for both, the difference between them is near zero too. That is the quantitative justification: the summit zero-point-energy difference is small because the mode that carried the difference has been converted into the reaction path and no longer stores vibrational energy. So the isotopic energy gap that survives to the summit is only the small share held in the bending modes, which is why the reactant gap dominates:

$Δ E_{a} = E_{0, H} - E_{0, D} = \frac{1}{2} h (ν_{H} - ν_{D}) > 0$

$E_{0, H}, E_{0, D}$ — starting floors from Step 4.
$Δ E_{a}$ — the extra activation energy the D bond needs. It is positive because H starts higher, so H needs less climb.

PICTURE. The full reaction landscape: a valley on the left rising to a summit. Two starting floors (orange H higher, magenta D lower) share one summit, with two arrows showing climb heights — the D arrow visibly taller.

Step 6 — Turning the energy gap into a rate ratio

WHAT. Reaction rate depends on activation energy through the Arrhenius equation:

$k = A e^{- E_{a} / R T}$

$k$ — the reaction rate (how fast the bond breaks).
$A$ — a prefactor (how often atoms attempt to cross).
$E_{a}$ — activation energy (the climb from Step 5).
$R$ — the gas constant; $T$ — absolute temperature.
$e^{- E_{a} / R T}$ — the fraction of attempts with enough energy to make the climb.

A qualification about $A$ . To isolate the effect we assume $A_{H} \approx A_{D}$ — that H and D attempt the crossing equally often. This is a simplification, and it is honest to flag where it can break down: isotopic substitution shifts the vibrational and rotational energy levels, which changes the entropy of the transition state and hence the prefactor, and — more dramatically at low temperature — it changes how much the light atom can tunnel through the barrier rather than climb over it (Step 7). For a first, temperature-moderate estimate near room temperature the $A$ 's nearly cancel and the exponential dominates; keep in mind the real ratio can drift from this idealised value.

Divide the two rates ( $A$ cancels under the assumption above):

$\frac{k _{H}}{k _{D}} = \frac{A e ^{- E_{a, H} / R T}}{A e ^{- E_{a, D} / R T}} = e^{(E_{a, D} - E_{a, H}) / R T} = e^{Δ E_{a} / R T}$

PICTURE. A plot of $k_{H} / k_{D}$ against $Δ E_{a}$ , rising steeply — showing how a tiny energy gap (a few kJ/mol) balloons into a rate ratio of 6–7.

Step 7 — The edge cases (hot, heavy, spectator, and cold tunnelling)

Every case matters, so let us walk the limits.

Hot limit $T \to \infty$ . As $T$ grows, $Δ E_{a} / R T \to 0$ , so $k_{H} / k_{D} \to e^{0} = 1$ . At very high temperature the isotope effect disappears — thermal energy dwarfs the tiny zero-point gap. Picture the exponent's fraction shrinking to nothing.

Cold limit $T \to 0$ (and where the simple picture breaks). Push the formula the other way: as $T \to 0$ , the exponent $Δ E_{a} / R T \to + \infty$ , so the predicted $k_{H} / k_{D} \to \infty$ . Read literally this says the ratio grows without bound as things get cold — and the reason it does not simply vanish is exactly the zero-point floor of Step 4, which is temperature-independent and therefore still tilts H above D even at absolute zero. But at very low $T$ the Arrhenius picture is no longer the whole story: the light H atom can tunnel — pass through the energy barrier instead of over it, a purely quantum shortcut. Because tunnelling is far easier for the lighter, faster-wobbling H than for D, it makes $k_{H} / k_{D}$ climb to very large values (often far above 7, sometimes into the hundreds) at low temperature, rather than diverging to infinity as the bare exponential claims. So the real cold-limit behaviour is a large but finite, tunnelling-dominated ratio — a genuinely quantum effect the spring-plus-Arrhenius model only hints at.

Heavy-element limit. For, say, a C–C bond where we swap $^{12}$ C for $^{13}$ C, the reduced mass barely changes ( $μ$ ratio $\approx 1.04$ , not $2$ ), so $ν_{H} - ν_{D}$ is tiny and $k_{H} / k_{D} \approx 1.02$ . This is why hydrogen is unique (as the parent note stressed): only for H does one neutron nearly double $μ$ .

Bond-not-breaking case. If the C–H bond is not stretched in the rate-determining step (it's just a spectator), the zero-point energies barely change between start and summit, so $Δ E_{a} \approx 0$ and $k_{H} / k_{D} \approx 1$ . A measured effect of ~7 is therefore a fingerprint that the C–H bond breaks in the slow step — this is exactly why chemists measure it.

Recall Check yourself

As temperature rises, does the kinetic isotope effect grow or shrink? ::: It shrinks toward 1 — the exponent $Δ E_{a} / R T \to 0$ . Why doesn't the effect vanish as temperature falls to absolute zero? ::: The zero-point energy floor is temperature-independent, so the H-above-D tilt survives; at very low T quantum tunnelling makes the ratio large but finite. Why is there almost no isotope effect swapping $^{12}$ C for $^{13}$ C? ::: One neutron barely changes carbon's reduced mass, so $ν$ and the zero-point energies hardly shift. If you measure $k_{H} / k_{D} \approx 1$ for a C–H containing reactant, what does it tell you? ::: The C–H bond is not breaking in the rate-determining step.

The one-picture summary

The single figure stacks the entire chain: heavier D ball → slower wobble $ν$ → lower zero-point floor $E_{0}$ → deeper starting hole → taller climb $Δ E_{a}$ → smaller rate through the exponential → $k_{H} / k_{D} \approx 7$ .

Recall Feynman retelling — say it like a story

Imagine a bond as a little ball bouncing on a spring against a heavy wall of carbon. Quantum rules say the ball can never fully stop bouncing — it always keeps a bit of jiggle even at the coldest temperature. Now swap the ordinary hydrogen ball for a deuterium ball that weighs twice as much. Same spring, heavier ball, so it bounces slower — and a slower bounce means less leftover jiggle. Less jiggle means the deuterium ball is sitting lower, resting deeper in its little energy valley. To break the bond, the ball has to climb out of the valley up to a mountain pass. Since deuterium started deeper, it has a taller climb. And here's the magic: reaction speed depends on that climb exponentially, so even a small extra climb makes deuterium react about seven times slower. Down at very cold temperatures a second quantum trick appears — the light hydrogen can sneak through the mountain instead of over it (tunnelling), which lets it pull even further ahead of the sluggish deuterium. That's the whole kinetic isotope effect — heavier ball, quieter jiggle, deeper hole, taller climb, much slower reaction, plus a cold-weather tunnelling bonus. And it only shows up dramatically for hydrogen, because only hydrogen doubles its weight when you add a single neutron.

Connections: this walkthrough leans on isotopes and mass number, borrows the ball-and-spring energy language cousin to kinetic theory, and its punchline drives reaction rate arguments. The stability side of hydrogen isotopes lives in Nuclear Stability and Binding Energy and Beta Decay; applications appear in Heavy Water and Nuclear Reactors, NMR Spectroscopy and Nuclear Fusion.