Case Study 1: Alpha Decay — Gamow's Triumph

Overview

In 1928, a 24-year-old Soviet physicist named George Gamow solved one of the deepest puzzles in nuclear physics: why certain heavy nuclei spontaneously emit alpha particles with well-defined energies, and why the half-lives of alpha emitters vary by more than 24 orders of magnitude — from microseconds to billions of years — while the alpha energies change by less than a factor of three. The answer, quantum tunneling through the Coulomb barrier, was one of the first applications of the newly minted quantum mechanics to nuclear physics. It remains one of the most spectacular quantitative successes of quantum theory.

This case study traces the history, physics, and quantitative predictions of Gamow's model in detail, connecting it to the WKB formalism developed in this chapter.

Part 1: The Puzzle Before Gamow

The Experimental Landscape (1900–1928)

By the 1920s, a large body of experimental data on alpha decay had been accumulated, principally by Ernest Rutherford and his collaborators at the Cavendish Laboratory in Cambridge.

Key experimental facts:

1. Discrete energies. Each radioactive species emits alpha particles of well-defined kinetic energies. $^{226}$Ra, for example, emits mainly 4.78 MeV alphas. This specificity indicated a definite energy difference between parent and daughter nuclear states.

2. Enormous range of half-lives. The shortest-lived alpha emitters ($^{212}$Po: $t_{1/2} = 0.30\,\mu$s) are separated from the longest-lived ($^{232}$Th: $t_{1/2} = 1.4 \times 10^{10}$ years) by a factor of $\sim 10^{24}$.

3. The Geiger-Nuttall law (1911). Hans Geiger and John Nuttall discovered an empirical linear relationship between $\log_{10} t_{1/2}$ and $1/\sqrt{E_\alpha}$ for members of the same radioactive series. A plot of $\log t_{1/2}$ versus $1/\sqrt{E_\alpha}$ yields a straight line. The physical origin of this relationship was completely mysterious.

4. The energy paradox. Rutherford had shown (1911) that the nucleus is a tiny, positively charged core surrounded by electrons. The alpha particle, carrying charge $+2e$, should face a Coulomb barrier of height:

$$V_{\text{barrier}} = \frac{2kZe^2}{R} \sim 25\text{–}35\,\text{MeV}$$

at the nuclear surface $R \sim 7$ fm. But the observed alpha energies were 4–9 MeV — far below the barrier. How could an alpha particle with 5 MeV of energy escape from a 30 MeV barrier?

The Classical Impossibility

Classical mechanics absolutely forbids a particle from traversing a region where its kinetic energy would be negative. If you throw a ball at a wall with energy $E$ and the wall's "potential energy" is $V > E$, the ball bounces back — always and without exception.

For alpha decay, the barrier is the Coulomb repulsion between the alpha particle ($+2e$) and the daughter nucleus ($+Ze$). The barrier height exceeds the alpha energy by factors of 3 to 8. Classically, the alpha particle should never escape.

Some physicists speculated that the alpha particle gained additional energy from unknown nuclear forces during its emission, somehow "boosting" it over the barrier. Others proposed that the Coulomb law broke down at nuclear distances. Neither explanation was satisfactory, and neither could explain the Geiger-Nuttall law.

Part 2: Gamow's Quantum Solution

The Key Insight: Tunneling

In 1928, Gamow — then visiting the University of Göttingen — recognized that the newly developed quantum mechanics provided a natural explanation. The alpha particle does not go over the Coulomb barrier. It goes through it.

Gamow had studied the quantum mechanics of potential barriers and knew that a particle encountering a barrier higher than its energy has a nonzero probability of appearing on the other side. This probability is exponentially small for thick, high barriers — but it is not zero. He applied this insight to the nuclear Coulomb barrier.

The Model

Gamow's model is remarkably simple:

Inside the nucleus ($r < R$): The alpha particle is bound in a deep potential well created by the strong nuclear force. The well depth is $-V_0 \sim -30$ to $-50$ MeV. The alpha particle bounces back and forth inside this well, hitting the walls with a frequency $f \sim v/(2R)$.

At the nuclear surface ($r = R$): The attractive nuclear force gives way to the repulsive Coulomb force. The potential jumps from $\sim -V_0$ to the Coulomb barrier height $V(R) = 2kZe^2/R$.

Outside the nucleus ($r > R$): The potential is purely Coulombic: $V(r) = 2kZe^2/r$. This falls off as $1/r$ and eventually drops below $E_\alpha$ at the outer turning point $r_2 = 2kZe^2/E_\alpha$.

The forbidden region: For $R < r < r_2$, we have $V(r) > E_\alpha$, and the alpha particle is classically forbidden. This is the region through which it must tunnel.

The Calculation

The WKB transmission probability is:

$$T = \exp\!\left(-\frac{2}{\hbar}\int_R^{r_2}\sqrt{2m_\alpha\left(\frac{2kZe^2}{r} - E_\alpha\right)}\,dr\right) = e^{-2\gamma}$$

The integral can be evaluated by the substitution $r = r_2\cos^2\theta$:

$$\int_R^{r_2}\sqrt{\frac{2kZe^2}{r} - E_\alpha}\,dr = \sqrt{E_\alpha}\, r_2\left[\arccos\sqrt{\frac{R}{r_2}} - \sqrt{\frac{R}{r_2}\left(1 - \frac{R}{r_2}\right)}\right]$$

In the limit $R/r_2 \ll 1$ (which is typically well-satisfied, since $R \sim 7$ fm and $r_2 \sim 30$–$60$ fm for the heaviest nuclei):

$$2\gamma \approx \frac{2\pi k Z e^2}{\hbar v_\alpha} - 4\sqrt{\frac{2m_\alpha k Z e^2 R}{\hbar^2}}$$

where $v_\alpha = \sqrt{2E_\alpha/m_\alpha}$ is the alpha particle velocity at infinity.

The Decay Rate

The alpha particle hits the barrier $f = v_\alpha'/(2R)$ times per second, where $v_\alpha' = \sqrt{2(E_\alpha + V_0)/m_\alpha}$ is the velocity inside the nucleus. Each time, it has probability $T = e^{-2\gamma}$ of escaping. The decay rate is:

$$\lambda = f \cdot T = \frac{v_\alpha'}{2R}\, e^{-2\gamma}$$

The pre-exponential factor $f \sim 10^{21}$ Hz varies slowly from nucleus to nucleus. The exponential factor $e^{-2\gamma}$ varies enormously. Essentially all of the variation in half-lives comes from the exponential — from the Gamow factor.

Part 3: Quantitative Predictions — A Triumph Across 24 Orders of Magnitude

Systematic Comparison

The following table compares Gamow model predictions with experimental data for a selection of alpha emitters spanning the full range of known half-lives:

The WKB Gamow model tracks the experimental half-lives across 24 orders of magnitude — from $\sim 10^{-7}$ seconds to $\sim 10^{18}$ seconds — with typical errors of less than 0.5 in $\log_{10}(t_{1/2})$, i.e., within a factor of 3.

The Geiger-Nuttall Law Explained

The first term in the Gamow exponent, $2\gamma \approx 2\pi kZe^2/(\hbar v_\alpha)$, is proportional to $Z/\sqrt{E_\alpha}$. For a series of isotopes with similar $Z$ (e.g., the uranium series), this gives:

$$\log_{10} t_{1/2} \approx A + \frac{B}{\sqrt{E_\alpha}}$$

which is precisely the Geiger-Nuttall law. Gamow's theory not only explained the empirical law but derived the constants $A$ and $B$ from fundamental physics.

Why Small Energy Changes Produce Enormous Lifetime Changes

The Gamow exponent for $^{238}$U is $2\gamma \approx 88$, and for $^{212}$Po it is $2\gamma \approx 33$. The difference is $\Delta(2\gamma) \approx 55$, which translates to a factor of:

$$\frac{t_{1/2}(^{238}\text{U})}{t_{1/2}(^{212}\text{Po})} \sim e^{55} \sim 10^{24}$$

A change in alpha energy from 4.27 MeV to 8.95 MeV — barely a factor of 2 — produces a factor of $10^{24}$ in half-life. This extreme sensitivity is the signature of quantum tunneling: the transmission probability depends exponentially on the barrier integral, and even modest changes in the integrand (driven by energy changes) produce dramatic changes in the result.

Part 4: Beyond the Simple Model — Refinements and Limitations

The Pre-Formation Factor

Gamow's model assumes the alpha particle exists as a pre-formed entity inside the nucleus, bouncing against the barrier walls. In reality, the four nucleons that will become the alpha particle spend most of their time distributed among the other nucleons. The probability that they happen to cluster together at the nuclear surface is the spectroscopic factor or preformation probability $P_\alpha$.

For even-even nuclei (even $Z$ and even $N$), $P_\alpha \sim 0.1$–$1.0$. For odd-$A$ nuclei, $P_\alpha$ can be much smaller, reflecting the additional difficulty of assembling the alpha cluster when one nucleon does not fit neatly into the configuration.

Including $P_\alpha$ modifies the decay rate:

$$\lambda = P_\alpha \cdot f \cdot e^{-2\gamma}$$

This correction is typically an order-unity factor that improves agreement with experiment.

Angular Momentum Barrier

If the alpha particle carries orbital angular momentum $l$ (relative to the daughter nucleus), the effective barrier includes the centrifugal term:

$$V_{\text{eff}}(r) = \frac{2kZe^2}{r} + \frac{l(l+1)\hbar^2}{2m_\alpha r^2}$$

The centrifugal barrier raises the effective potential, increasing $\gamma$ and suppressing the decay rate. Alpha decays with $l > 0$ (which change the spin/parity of the nucleus) are systematically slower — an effect known as hindrance.

Nuclear Deformation

Some nuclei (especially in the rare earth and actinide regions) are deformed — their charge distribution is elongated rather than spherical. For these nuclei, the barrier height depends on the direction of emission relative to the nuclear symmetry axis. The effective barrier is lower along the long axis, creating a preferred emission direction.

The WKB calculation must be averaged over orientations:

$$\lambda = \frac{1}{4\pi}\int \lambda(\theta, \phi)\, d\Omega$$

This orientation averaging can modify the predicted half-life by factors of 2–5 for strongly deformed nuclei.

Part 5: Legacy and Modern Relevance

Nuclear Astrophysics

Gamow soon recognized that his tunneling calculation could be run "in reverse" — instead of particles tunneling out of a nucleus, he considered particles tunneling into one. This is nuclear fusion.

In the cores of stars, protons at temperatures of $\sim 10^7$ K (thermal energies $\sim 1$ keV) must fuse despite a Coulomb barrier of $\sim 1$ MeV. The tunneling probability for a proton approaching another proton is:

$$T \propto \exp\!\left(-\frac{\pi k e^2}{\hbar}\sqrt{\frac{2m_p}{E}}\right) = \exp\!\left(-\frac{b}{\sqrt{E}}\right)$$

This is the same Gamow factor, now applied to fusion. The interplay between this tunneling factor and the Maxwell-Boltzmann thermal distribution creates the Gamow peak — the narrow energy window where fusion reactions actually occur. Gamow's tunneling theory is the foundation of nuclear astrophysics.

Proton Radioactivity

In the 1980s, proton radioactivity was discovered — the emission of a single proton from proton-rich nuclei. The Gamow model applies directly, with $Z_\alpha \to 1$ and $m_\alpha \to m_p$. The lower charge and mass of the proton (compared to the alpha) mean that proton emission is typically observed only from nuclei very far from stability, where the proton separation energy is small.

Superheavy Elements

The search for superheavy elements (atomic numbers $Z > 103$) relies heavily on the Gamow model. Theoretical predictions of alpha-decay half-lives, computed using WKB tunneling through the Coulomb barrier (with shell-model corrections), guide experimentalists in designing synthesis reactions and estimating how long they have to observe the new nuclei before they decay.

The "island of stability" predicted near $Z = 114$, $N = 184$ is a region where nuclear shell effects are expected to produce extra binding energy, reducing the alpha energy and enormously increasing the half-life. Elements up to $Z = 118$ (oganesson) have been synthesized, and their measured alpha-decay properties are consistent with WKB predictions.

Discussion Questions

1. Gamow's model treats the alpha particle as classical inside the nucleus (bouncing back and forth) but quantum mechanical at the barrier (tunneling). Is this hybrid approach self-consistent? Under what conditions would a fully quantum mechanical treatment of the alpha particle inside the nucleus give different results?

2. The WKB prediction for $^{238}$U is accurate to within a factor of 2 — across a number that spans 40 decimal places. In what sense is this "good agreement," and in what sense is it "off by a factor of 2"? How does this compare with the precision of, say, the hydrogen spectrum calculations in Chapter 18?

3. Gamow published his theory in 1928, just two years after the Schrödinger equation. What does this rapid application tell us about the culture of theoretical physics in the 1920s? Could a comparable theoretical breakthrough happen as quickly today?

4. The Geiger-Nuttall law was known empirically for 17 years (1911–1928) before Gamow explained it. What prevented an earlier explanation? What had to happen in physics between 1911 and 1928 to make Gamow's calculation possible?

5. Alpha decay conserves energy — the total energy (kinetic + rest mass) is the same before and after. Does this mean the alpha particle "borrows" energy from the vacuum to cross the barrier? How should we think about energy conservation during the tunneling process?

Quantitative Exercises

E1. Using the data in the table in Part 3, make a Geiger-Nuttall plot: $\log_{10}(t_{1/2}/\text{s})$ vs. $1/\sqrt{E_\alpha}$ (MeV)$^{-1/2}$. Verify that the relationship is approximately linear for isotopes with similar daughter $Z$. Determine the slope and intercept from a best-fit line.

E2. Calculate the Gamow exponent and estimated half-life for $^{144}$Nd ($Z_{\text{daughter}} = 58$, $E_\alpha = 1.90$ MeV, $R = 6.7$ fm), one of the rarest natural alpha emitters. The experimental half-life is $2.3 \times 10^{15}$ years. How does your estimate compare?

E3. The reaction $p + p \to d + e^+ + \nu_e$ (proton-proton fusion) occurs in the Sun at a core temperature of $1.5 \times 10^7$ K. Calculate the Gamow peak energy $E_0 = (bk_BT/2)^{2/3}$ and compare with $k_BT$. Why does fusion occur at energies much higher than the thermal average?

E4. If the fine structure constant $\alpha = e^2/(4\pi\epsilon_0\hbar c)$ were 5% larger, how would the half-life of $^{238}$U change? (Hint: the Gamow exponent is proportional to $\alpha$.) What would this mean for the existence of natural uranium on Earth?