Chapter 13 — Alpha Decay: Quantum Tunneling Through the Coulomb Barrier

"It has generally been assumed that the charged particles need to have an energy greater than the height of the potential barrier in the field of the nucleus in order to get in or out... We show, however, that the wave mechanical treatment of this problem gives a finite probability for a particle to go through a potential barrier even when its kinetic energy is less than the height of the barrier." — George Gamow, Zeitschrift für Physik 51, 204 (1928)

Chapter Overview

In 1928, nuclear physics faced a paradox. Alpha particles emitted by heavy nuclei carried kinetic energies of 4--9 MeV. Yet the Coulomb barrier between an alpha particle and a daughter nucleus -- the electrostatic energy required to bring them to the nuclear surface -- stood at 25--30 MeV. Classically, a particle inside a potential well cannot escape if its energy is less than the barrier height. The alpha particle should remain trapped inside the nucleus forever.

Three physicists -- George Gamow in Göttingen, and independently Ronald Gurney and Edward Condon at Princeton -- recognized that the newly formulated quantum mechanics offered a resolution. In quantum mechanics, a particle described by a wave function has a nonzero probability of penetrating a potential barrier even when its kinetic energy is less than the barrier height. This phenomenon, quantum tunneling, not only explained alpha decay but also provided one of the first quantitative triumphs of quantum mechanics applied to the atomic nucleus.

This chapter develops the full theory of alpha decay as barrier penetration:

• We set up the physical problem: an alpha particle confined inside a nucleus, facing a Coulomb barrier it cannot classically surmount.
• We apply the WKB approximation (Chapter 5) to calculate the tunneling probability through the Coulomb barrier, deriving the Gamow factor.
• We derive the Geiger-Nuttall law -- the remarkable linear relationship between log half-life and $Z/\sqrt{E_\alpha}$ -- and show that it is a direct consequence of the tunneling calculation.
• We examine alpha decay energetics: Q-values, systematics across the chart of nuclides, and the role of the SEMF (Chapter 4).
• We explore fine structure in alpha decay: transitions to excited states of the daughter nucleus.
• We extend the tunneling framework to proton radioactivity and cluster radioactivity -- exotic decay modes where the emitted particle is not an alpha but a proton or a heavier nuclear cluster.

🏃 Fast Track: If you are familiar with the WKB approximation from Chapter 5, you can skim Section 13.2 and focus on the Gamow factor derivation in Section 13.3. The essential physics is in Sections 13.3--13.5.

🔬 Deep Dive: The full WKB integral evaluation (Section 13.3.3) and the connection to cluster radioactivity (Section 13.8) reward careful study. The cluster decay discussion connects to modern research in nuclear structure.

13.1 The Alpha Decay Puzzle

13.1.1 What Is an Alpha Particle?

An alpha particle is the nucleus of ${}^4\text{He}$: two protons and two neutrons bound together with exceptional stability. Its binding energy is $B({}^4\text{He}) = 28.30\,\text{MeV}$, giving $B/A = 7.07\,\text{MeV/nucleon}$ -- remarkably high for such a light nucleus. This extraordinary stability is the reason nature preferentially emits alpha particles rather than individual nucleons or other light clusters: the energy gain from forming the tightly bound ${}^4\text{He}$ configuration makes alpha emission energetically favorable for many heavy nuclei.

In alpha decay, a parent nucleus ${}^A_Z\text{X}$ emits an alpha particle and transforms into a daughter nucleus:

$${}^A_Z\text{X} \to {}^{A-4}_{Z-2}\text{Y} + {}^4_2\text{He}$$

The parent loses two protons and two neutrons, reducing both its atomic number by 2 and its mass number by 4.

13.1.2 The Energetics: Q-Value

The energy released in alpha decay is the Q-value:

$$\boxed{Q_\alpha = \left[M({}^A_Z\text{X}) - M({}^{A-4}_{Z-2}\text{Y}) - M({}^4_2\text{He})\right]c^2}$$

where all masses are atomic masses (the electron masses cancel). Alpha decay is energetically allowed if and only if $Q_\alpha > 0$.

The Q-value is shared between the kinetic energy of the alpha particle and the recoil kinetic energy of the daughter nucleus. By conservation of momentum (the parent is at rest):

$$p_\alpha = p_\text{daughter} \implies M_\alpha T_\alpha = M_d T_d$$

where $T$ denotes kinetic energy and we use the nonrelativistic relation $p^2 = 2MT$. Combined with energy conservation $Q = T_\alpha + T_d$:

$$\boxed{T_\alpha = Q \cdot \frac{M_d}{M_d + M_\alpha} \approx Q \cdot \frac{A-4}{A}}$$

For heavy nuclei ($A \gg 4$), the alpha particle carries nearly all the released energy. For example, in the decay ${}^{238}\text{U} \to {}^{234}\text{Th} + \alpha$: - $Q = 4.270\,\text{MeV}$ - $T_\alpha = 4.270 \times 234/238 = 4.198\,\text{MeV}$ - $T_\text{Th} = 4.270 \times 4/238 = 0.072\,\text{MeV}$

The daughter recoil carries only 1.7% of the total energy.

13.1.3 The Classical Paradox

Now consider the Coulomb barrier. After the alpha particle separates from the daughter nucleus and moves beyond the range of the nuclear force (at the nuclear surface $r = R$), it faces a repulsive Coulomb potential:

$$V_C(r) = \frac{z_\alpha Z_d e^2}{4\pi\epsilon_0 r} = \frac{2 Z_d e^2}{4\pi\epsilon_0 r} \quad \text{for } r > R$$

where $z_\alpha = 2$ and $Z_d = Z - 2$. The barrier height at the nuclear surface is:

$$V_C(R) = \frac{2 Z_d e^2}{4\pi\epsilon_0 R}$$

Using $R \approx r_0(A_d^{1/3} + 4^{1/3})$ with $r_0 \approx 1.2\,\text{fm}$:

Numerical example -- ${}^{238}\text{U}$: - $Z_d = 90$ (thorium), $A_d = 234$ - $R \approx 1.2(234^{1/3} + 4^{1/3}) \approx 1.2(6.16 + 1.59) = 9.3\,\text{fm}$ - $V_C(R) = \frac{2 \times 90 \times 1.440\,\text{MeV}\cdot\text{fm}}{9.3\,\text{fm}} \approx 27.9\,\text{MeV}$

The alpha particle emerges with $T_\alpha = 4.2\,\text{MeV}$, but the Coulomb barrier at the nuclear surface is approximately $28\,\text{MeV}$ -- nearly seven times higher. Classically, the alpha particle cannot escape. It would need to be at least 28 MeV to overcome the barrier, yet it emerges with only 4.2 MeV.

💡 Intuition: Imagine a ball rolling in a bowl. Classically, if the ball's energy is less than the height of the rim, it can never escape -- it just oscillces back and forth forever. Yet alpha particles do escape. Something beyond classical physics must be at work.

This was the puzzle that confronted physicists in the late 1920s. The resolution came from quantum mechanics.

13.1.4 Historical Context: The State of Play in 1928

By 1928, quantum mechanics was only two years old. Schrödinger had published his wave equation in 1926; Heisenberg, Born, and Jordan had developed matrix mechanics the year before. The wave-mechanical formalism was being applied with spectacular success to atomic physics — energy levels of hydrogen, the helium atom, the structure of the periodic table — but its application to the nucleus was essentially unexplored.

The alpha decay puzzle had been recognized since at least 1911, when Rutherford demonstrated that alpha particles scattered from heavy nuclei experienced a Coulomb potential consistent with a point-like nuclear charge. The Geiger-Nuttall law, discovered the same year, established that half-lives and alpha energies were correlated, but no one understood why. For 17 years, the mechanism of alpha emission remained a mystery.

The breakthrough came nearly simultaneously from two groups. George Gamow, working at the Institute for Theoretical Physics in Göttingen under Max Born, attacked the problem in the summer of 1928. He modeled the alpha particle as a quantum-mechanical wave inside a potential well surrounded by the Coulomb barrier and applied the WKB semiclassical approximation to compute the tunneling probability. His paper, "Zur Quantentheorie des Atomkernes," appeared in Zeitschrift für Physik in 1928.

Meanwhile, Ronald Gurney and Edward Condon at Princeton University were working on the same problem, unaware of Gamow's efforts. Their paper, "Wave Mechanics and Radioactive Disintegration," appeared as a letter to Nature in September 1928. Both groups arrived at the same conclusion: alpha decay is quantum tunneling through the Coulomb barrier.

📜 Historical Note: The simultaneous discovery by Gamow and by Gurney-Condon is one of the most striking examples of independent co-discovery in physics. It reflects the fact that by 1928, the mathematical tools (WKB, Schrödinger equation) and the physical problem (alpha energies vs. barrier heights) were both sufficiently well established that the solution was, in some sense, "in the air."

The impact was immediate and profound. Gamow's calculation was the first successful application of quantum mechanics to a nuclear phenomenon. It demonstrated that the nucleus obeyed the same quantum laws as the atom, provided the first quantitative explanation of the Geiger-Nuttall law, and established the WKB approximation as a practical computational tool. Within a year, Gamow and independently Atkinson and Houtermans applied the "inverse" of the tunneling idea — charged particles tunneling into nuclei — to explain thermonuclear reactions in stars, laying the foundation for nuclear astrophysics.

13.1.5 Additional Q-Value Examples

To build facility with alpha decay energetics, let us compute Q-values for several important alpha emitters. All atomic masses are from the AME2020 evaluation.

Example 1: ${}^{210}\text{Po} \to {}^{206}\text{Pb} + \alpha$

Using $M({}^{210}\text{Po}) = 209.982874\,\text{u}$, $M({}^{206}\text{Pb}) = 205.974465\,\text{u}$, $M({}^4\text{He}) = 4.002603\,\text{u}$:

$$Q = (209.982874 - 205.974465 - 4.002603) \times 931.494 = 0.005806 \times 931.494 = 5.407\,\text{MeV}$$

$$T_\alpha = 5.407 \times \frac{206}{210} = 5.304\,\text{MeV}$$

This is the classic polonium alpha decay studied by Marie Curie and many others. The half-life is $138.4\,\text{days}$.

Example 2: ${}^{226}\text{Ra} \to {}^{222}\text{Rn} + \alpha$

Using $M({}^{226}\text{Ra}) = 226.025410\,\text{u}$, $M({}^{222}\text{Rn}) = 222.017578\,\text{u}$:

$$Q = (226.025410 - 222.017578 - 4.002603) \times 931.494 = 0.005229 \times 931.494 = 4.871\,\text{MeV}$$

$$T_\alpha = 4.871 \times \frac{222}{226} = 4.784\,\text{MeV}$$

The half-life of ${}^{226}\text{Ra}$ is 1,600 years — long enough that radium persists in Earth's crust from the uranium decay series, but short enough that it is intensely radioactive (specific activity $3.66 \times 10^{10}\,\text{Bq/g}$). This is the isotope that Marie and Pierre Curie isolated in 1898.

Example 3: ${}^{148}\text{Sm} \to {}^{144}\text{Nd} + \alpha$

Using $M({}^{148}\text{Sm}) = 147.914829\,\text{u}$, $M({}^{144}\text{Nd}) = 143.910093\,\text{u}$:

$$Q = (147.914829 - 143.910093 - 4.002603) \times 931.494 = 0.002133 \times 931.494 = 1.987\,\text{MeV}$$

This is one of the lightest and longest-lived alpha emitters: $t_{1/2} = 6.3 \times 10^{15}\,\text{yr}$ — about 450,000 times the age of the universe. The tiny Q-value translates, through the Gamow factor, to an astronomically long half-life.

🔄 Check Your Understanding: Compute the Q-value for ${}^{232}\text{Th} \to {}^{228}\text{Ra} + \alpha$ using $M({}^{232}\text{Th}) = 232.038055\,\text{u}$ and $M({}^{228}\text{Ra}) = 228.031070\,\text{u}$. You should find $Q = 4.082\,\text{MeV}$, giving $T_\alpha = 4.012\,\text{MeV}$.

13.2 The Quantum Tunneling Mechanism

13.2.1 The One-Body Model

The Gamow model treats alpha decay as a one-body tunneling problem. The alpha particle is assumed to exist as a preformed entity inside the nucleus, bouncing back and forth against the potential walls. Each time it strikes the barrier, there is a small but nonzero probability that it tunnels through.

The potential energy experienced by the alpha particle as a function of its distance $r$ from the center of the daughter nucleus is modeled as:

$$V(r) = \begin{cases} -V_0 & r < R \quad \text{(nuclear interior: square well)} \\ \displaystyle\frac{2Z_d e^2}{4\pi\epsilon_0 r} & r \geq R \quad \text{(Coulomb barrier)} \end{cases}$$

where $V_0 \approx 30$--$50\,\text{MeV}$ is the depth of the nuclear potential well and $R$ is the nuclear radius at which the alpha particle separates from the daughter. The key features of this potential are:

1. Inside the nucleus ($r < R$): The alpha particle is bound in an attractive nuclear potential of depth $\sim V_0$.
2. At the surface ($r = R$): The potential jumps to the Coulomb value $V_C(R) \gg E_\alpha$, where $E_\alpha = T_\alpha$ is the kinetic energy of the emitted alpha.
3. The barrier region ($R < r < b$): The Coulomb potential $V_C(r) > E_\alpha$. This is the classically forbidden region.
4. Beyond the barrier ($r > b$): $V_C(r) < E_\alpha$, and the alpha particle propagates freely.

The classical turning point $b$ is where $V_C(b) = E_\alpha$:

$$b = \frac{2Z_d e^2}{4\pi\epsilon_0 E_\alpha}$$

For ${}^{238}\text{U}$: $b = 2 \times 90 \times 1.440 / 4.2 \approx 61.7\,\text{fm}$.

The barrier extends from $R \approx 9.3\,\text{fm}$ to $b \approx 62\,\text{fm}$ -- a width of over 50 fm. It is tall (28 MeV at the surface) and wide. The tunneling probability is exponentially sensitive to both the height and width.

13.2.2 The Decay Rate

In the one-body model, the decay constant (inverse of the mean lifetime) is:

$$\boxed{\lambda = f \cdot P \cdot S}$$

where: - $f$ is the assault frequency -- the number of times per second the alpha particle strikes the barrier wall. For an alpha particle with kinetic energy $\sim 30\,\text{MeV}$ bouncing inside a nucleus of radius $\sim 7\,\text{fm}$: $v \approx \sqrt{2E/m_\alpha} \approx 0.12c$ and $f = v/(2R) \approx 10^{21}\,\text{s}^{-1}$. - $P$ is the barrier penetration factor (tunneling probability), calculated by the WKB method. - $S$ is the preformation factor (also called the spectroscopic factor) -- the probability that the alpha particle exists as a preformed cluster inside the nucleus. For even-even nuclei, $S \sim 10^{-1}$--$10^{-2}$. We will return to this in Section 13.6.

The half-life is:

$$t_{1/2} = \frac{\ln 2}{\lambda} = \frac{\ln 2}{f \cdot P \cdot S}$$

The enormous range of alpha decay half-lives (from $\sim 10^{-7}\,\text{s}$ for ${}^{212}\text{Po}$ to $\sim 10^{17}\,\text{yr}$ for ${}^{148}\text{Sm}$) is overwhelmingly determined by $P$, because $P$ varies over many orders of magnitude while $f$ and $S$ are relatively constant.

13.3 The WKB Calculation: Deriving the Gamow Factor

13.3.1 The WKB Tunneling Formula

Recall from Chapter 5 that the WKB approximation gives the transmission probability through a potential barrier $V(r)$ for a particle with energy $E < V$ as:

$$\boxed{P = \exp\left(-\frac{2}{\hbar}\int_R^b \sqrt{2\mu[V(r) - E]}\,dr\right)}$$

where $\mu$ is the reduced mass of the alpha-daughter system:

$$\mu = \frac{m_\alpha M_d}{m_\alpha + M_d} \approx m_\alpha \left(1 - \frac{4}{A}\right) \approx m_\alpha$$

and the integral runs over the classically forbidden region from the inner turning point $R$ (nuclear surface) to the outer turning point $b$ (where $V(b) = E$).

⚠️ Common Pitfall: The factor of 2 in the exponent comes from the WKB connection formulas. It is not a factor of $2\pi$. Some textbooks absorb the factor differently. Be consistent.

13.3.2 Setting Up the Integral

For the Coulomb barrier, $V(r) = kz_\alpha Z_d e^2/r$ where $k = 1/(4\pi\epsilon_0)$. Define the Coulomb parameter:

$$\eta = \frac{z_\alpha Z_d e^2}{4\pi\epsilon_0 \hbar v} = \frac{Z_d e^2}{2\pi\epsilon_0 \hbar v}$$

where $v = \sqrt{2E/\mu}$ is the asymptotic velocity of the alpha particle. The Sommerfeld parameter $\eta$ is a dimensionless measure of the strength of the Coulomb interaction relative to the quantum of action.

The exponent in the WKB formula (called the Gamow factor $G$) is:

$$G = \frac{2}{\hbar}\int_R^b \sqrt{2\mu\left(\frac{kz_\alpha Z_d e^2}{r} - E\right)}\,dr$$

Using $E = kz_\alpha Z_d e^2/b$, we write:

$$G = \frac{2}{\hbar}\sqrt{2\mu E}\int_R^b \sqrt{\frac{b}{r} - 1}\,dr$$

13.3.3 Evaluating the Integral

This integral has a standard form. Substitute $r = b\cos^2\theta$, so $dr = -2b\cos\theta\sin\theta\,d\theta$:

$$\int_R^b \sqrt{\frac{b}{r} - 1}\,dr = \int_{\theta_0}^{\pi/2}\frac{\sin\theta}{\cos\theta}\cdot 2b\cos\theta\sin\theta\,d\theta = 2b\int_{\theta_0}^{\pi/2}\sin^2\theta\,d\theta$$

where $\cos^2\theta_0 = R/b$, i.e., $\theta_0 = \arccos\sqrt{R/b}$. Using $\sin^2\theta = (1-\cos 2\theta)/2$:

$$2b\int_{\theta_0}^{\pi/2}\sin^2\theta\,d\theta = b\left[\theta - \frac{1}{2}\sin 2\theta\right]_{\theta_0}^{\pi/2} = b\left[\frac{\pi}{2} - \theta_0 + \frac{1}{2}\sin 2\theta_0\right]$$

Since $\sin 2\theta_0 = 2\sin\theta_0\cos\theta_0 = 2\sqrt{1 - R/b}\cdot\sqrt{R/b}$:

$$\int_R^b \sqrt{\frac{b}{r} - 1}\,dr = b\left[\frac{\pi}{2} - \arccos\sqrt{\frac{R}{b}} + \sqrt{\frac{R}{b}\left(1 - \frac{R}{b}\right)}\right]$$

Therefore the Gamow factor is:

$$\boxed{G = \frac{2}{\hbar}\sqrt{2\mu E}\cdot b\left[\arccos\sqrt{\frac{R}{b}} - \sqrt{\frac{R}{b}\left(1-\frac{R}{b}\right)}\right]}$$

where we have used $\pi/2 - \arccos x = \arccos\sqrt{1-x^2}$ ... but let us write this more transparently. Define the dimensionless ratio $\rho = R/b$. Then:

$$\boxed{G = 2\pi\eta\left[\frac{1}{\pi}\arccos\sqrt{\rho} - \frac{1}{\pi}\sqrt{\rho(1-\rho)}\right]}$$

where $\eta = z_\alpha Z_d e^2/(4\pi\epsilon_0\hbar v)$ is the Sommerfeld parameter. Equivalently:

$$G = 2\eta\left[\arccos\sqrt{\rho} - \sqrt{\rho(1-\rho)}\right]$$

This is the exact WKB result for the Gamow factor for a pure Coulomb barrier.

13.3.4 The Limiting Cases

Case 1: Thin barrier ($R/b \to 1$, i.e., $E \to V_C(R)$).

When the alpha energy approaches the barrier height, $\rho \to 1$, and $G \to 0$. The tunneling probability $P = e^{-G} \to 1$ — the particle sails over the barrier classically. This is correct.

Case 2: Thick barrier ($R/b \ll 1$, i.e., $E \ll V_C(R)$).

Expanding for small $\rho$:

$$\arccos\sqrt{\rho} \approx \frac{\pi}{2} - \sqrt{\rho}, \qquad \sqrt{\rho(1-\rho)} \approx \sqrt{\rho}$$

Therefore:

$$G \approx 2\eta\left[\frac{\pi}{2} - 2\sqrt{\rho}\right] = \pi\eta - 4\eta\sqrt{R/b}$$

Using $\eta = z_\alpha Z_d e^2/(4\pi\epsilon_0\hbar v)$ and $b = z_\alpha Z_d e^2/(4\pi\epsilon_0 E)$:

$$G \approx \frac{2\pi z_\alpha Z_d e^2}{4\pi\epsilon_0\hbar v} - \frac{4z_\alpha Z_d e^2}{4\pi\epsilon_0\hbar v}\sqrt{\frac{R}{b}}$$

The first term is the Gamow-Sommerfeld factor $2\pi\eta$ for a pure Coulomb barrier with no inner cutoff. The second term is a correction for the finite nuclear radius. This approximation is excellent for alpha decay, where typically $R/b \sim 0.1$--$0.2$.

💡 Intuition: The Gamow factor is proportional to $Z_d/\sqrt{E_\alpha}$. Since $\eta \propto Z_d/v \propto Z_d/\sqrt{E_\alpha}$, both the barrier height and width increase with $Z_d$ and decrease with $E_\alpha$. This is why the tunneling probability is so extraordinarily sensitive to the Q-value.

13.3.5 Numerical Evaluation

Let us compute the Gamow factor for ${}^{238}\text{U} \to {}^{234}\text{Th} + \alpha$ ($E_\alpha = 4.198\,\text{MeV}$):

Parameters: - $Z_d = 90$, $z_\alpha = 2$ - $\mu = 4 \times 234/(4+234) \times 931.5 = 3,654\,\text{MeV}/c^2 \approx 3.727\,\text{u}$ - $v = \sqrt{2E_\alpha/\mu} = \sqrt{2 \times 4.198/3654} \cdot c = 0.04793c$ - $\eta = 2 \times 90 \times 1.440/(2 \times 0.04793 \times 197.3) = 259.2/(18.91) = 13.71$

Wait — let us be more careful:

$$\eta = \frac{z_\alpha Z_d e^2}{4\pi\epsilon_0 \hbar v} = \frac{2 \times 90 \times 1.440\,\text{MeV}\cdot\text{fm}}{197.3\,\text{MeV}\cdot\text{fm} \times 0.04793} = \frac{259.2}{9.456} = 27.41$$

Turning points: - $R = 1.2(234^{1/3}+4^{1/3}) = 1.2(6.161+1.587) = 9.30\,\text{fm}$ - $b = 2\times90\times1.440/4.198 = 259.2/4.198 = 61.72\,\text{fm}$ - $\rho = R/b = 9.30/61.72 = 0.1507$

Gamow factor:

$$G = 2\times 27.41\left[\arccos\sqrt{0.1507} - \sqrt{0.1507\times 0.8493}\right]$$

$$= 54.82\left[\arccos(0.3882) - \sqrt{0.1280}\right]$$

$$= 54.82\left[1.1727 - 0.3577\right] = 54.82 \times 0.8150 = 44.69$$

So the tunneling probability is:

$$P = e^{-G} = e^{-44.69} \approx 4.3 \times 10^{-20}$$

📊 Putting This in Perspective: The alpha particle strikes the barrier wall approximately $10^{21}$ times per second, yet the probability of tunneling on each attempt is only $4 \times 10^{-20}$. The average number of assaults before escape is $\sim 10^{20}$. The resulting decay rate gives a half-life on the order of $10^9$ years — consistent with the measured value of $t_{1/2}({}^{238}\text{U}) = 4.468 \times 10^9\,\text{yr}$.

13.3.6 Second Example: ${}^{212}\text{Po}$ — A Fast Alpha Emitter

To appreciate the sensitivity of the tunneling calculation, let us repeat the calculation for ${}^{212}\text{Po}$, which has one of the shortest known alpha decay half-lives: $t_{1/2} = 0.299\,\mu\text{s}$.

Parameters: - $Z_d = 82$ (lead), $A_d = 208$, $E_\alpha = 8.784\,\text{MeV}$ - $\mu = 4 \times 208/(4+208) \times 931.5 = 3642\,\text{MeV}/c^2$ - $v = \sqrt{2 \times 8.784/3642}\cdot c = 0.0695c$ - $\eta = 2 \times 82 \times 1.440/(197.3 \times 0.0695) = 236.2/13.71 = 17.23$ - $R = 1.2(208^{1/3} + 4^{1/3}) = 1.2(5.925 + 1.587) = 9.01\,\text{fm}$ - $b = 2 \times 82 \times 1.440/8.784 = 236.2/8.784 = 26.89\,\text{fm}$ - $\rho = 9.01/26.89 = 0.3350$

$$G = 2 \times 17.23 \times [\arccos(0.5789) - \sqrt{0.3350 \times 0.6650}]$$

$$= 34.46 \times [0.9534 - 0.4722] = 34.46 \times 0.4812 = 16.58$$

$$P = e^{-16.58} = 6.3 \times 10^{-8}$$

Compare this to ${}^{238}\text{U}$ where $G = 44.69$ and $P = 4.3 \times 10^{-20}$. The difference $\Delta G = 44.69 - 16.58 = 28.11$ translates to a ratio $P_{212}/P_{238} = e^{28.11} \approx 1.6 \times 10^{12}$. Combined with differences in the assault frequency (factor $\sim 2$) and accounting for the slight difference in preformation factors, this 12 orders of magnitude difference in $P$ is the primary reason why ${}^{212}\text{Po}$ decays in microseconds while ${}^{238}\text{U}$ takes billions of years.

💡 The Lesson: The Gamow factor $G$ changes from 16.6 to 44.7 — a factor of 2.7 — while the tunneling probability changes by 12 orders of magnitude. The exponential amplification is the heart of alpha decay physics.

13.3.7 Comparison of Exact and Approximate Gamow Factors

The following table compares the exact Gamow factor to the thick-barrier approximation for several alpha emitters:

The thick-barrier approximation is excellent (error $< 6\%$) whenever $\rho < 0.2$, which covers most alpha emitters. It fails for very high-energy decays like ${}^{212}\text{Po}$ where the outer turning point is close to the nuclear surface.

13.4 The Geiger-Nuttall Law

13.4.1 The Empirical Discovery

In 1911--1912, Hans Geiger and John Mitchell Nuttall made a remarkable empirical observation. For alpha emitters within a single radioactive decay series, there is a linear relationship between the logarithm of the decay constant (or half-life) and the logarithm of the range of the alpha particle in air (which is a proxy for its kinetic energy). In modern notation:

$$\boxed{\log_{10}\lambda = a + b\,\frac{Z_d}{\sqrt{E_\alpha}}}$$

or equivalently:

$$\log_{10}t_{1/2} = a' + b'\,\frac{Z_d}{\sqrt{E_\alpha}}$$

where $a$, $b$ (or $a'$, $b'$) are constants. The relationship is approximately linear when isotopes of a given element (fixed $Z_d$) are plotted.

This was discovered 17 years before quantum mechanics could explain it. It is one of the oldest empirical laws in nuclear physics.

13.4.2 Derivation from the Gamow Theory

The Geiger-Nuttall law follows directly from the tunneling calculation. In the thick-barrier approximation (Section 13.3.4):

$$G \approx \pi\eta - 4\eta\sqrt{R/b} = \frac{2\pi z_\alpha Z_d e^2}{4\pi\epsilon_0\hbar}\sqrt{\frac{\mu}{2E_\alpha}} - \frac{4}{\hbar}\sqrt{2\mu R \cdot\frac{z_\alpha Z_d e^2}{4\pi\epsilon_0}}$$

Define constants:

$$c_1 = \frac{2\pi z_\alpha e^2}{4\pi\epsilon_0\hbar}\sqrt{\frac{\mu}{2}} = \frac{\pi e^2}{\epsilon_0\hbar}\sqrt{\frac{\mu}{2}} \approx 1.980\,\text{MeV}^{-1/2}\cdot Z_d \text{ (for alpha decay)}$$

$$c_2 = \frac{4}{\hbar}\sqrt{2\mu R\cdot\frac{z_\alpha Z_d e^2}{4\pi\epsilon_0}}$$

Then:

$$G = c_1 \frac{Z_d}{\sqrt{E_\alpha}} - c_2$$

Since $\lambda = fS \cdot e^{-G}$:

$$\ln\lambda = \ln(fS) - G = \ln(fS) - c_1\frac{Z_d}{\sqrt{E_\alpha}} + c_2$$

Converting to $\log_{10}$:

$$\boxed{\log_{10}\lambda = A + B\,\frac{Z_d}{\sqrt{E_\alpha}}}$$

where $A = [\ln(fS) + c_2]/\ln 10$ and $B = -c_1/\ln 10$. For a fixed element ($Z_d$ constant), this gives:

$$\log_{10}t_{1/2} = a' - b'\sqrt{E_\alpha}$$

which is the Viola-Seaborg form of the Geiger-Nuttall law, often written as:

$$\log_{10}t_{1/2}(\text{s}) = \frac{aZ_d + b}{\sqrt{E_\alpha(\text{MeV})}} + cZ_d + d$$

with fitted constants. The physically transparent form is the linear relationship between $\log_{10}\lambda$ and $Z_d/\sqrt{E_\alpha}$.

13.4.3 The Extraordinary Sensitivity

The most striking feature of the Geiger-Nuttall law is the extreme sensitivity of the half-life to the alpha energy. Consider two even-even alpha emitters:

The alpha energies range from about 2 to 9 MeV — a factor of only $\sim 4.5$ — while the half-lives span from sub-microsecond to $10^{15}$ years — over 40 orders of magnitude. This enormous range is a direct consequence of the exponential dependence of the tunneling probability on the Gamow factor.

💡 Intuition: A small increase in $E_\alpha$ reduces $Z_d/\sqrt{E_\alpha}$ and thereby reduces the Gamow factor $G$. Since $P \propto e^{-G}$, even a modest decrease in $G$ produces an exponential increase in the tunneling probability and a corresponding decrease in the half-life. The barrier penetration factor is the amplifier that converts small energy differences into enormous half-life differences.

Numerical example: For ${}^{232}\text{Th}$ ($E_\alpha = 4.012\,\text{MeV}$, $Z_d = 88$) versus ${}^{212}\text{Po}$ ($E_\alpha = 8.784\,\text{MeV}$, $Z_d = 82$):

The Gamow factors are approximately $G_{232} \approx 48.6$ and $G_{212} \approx 18.4$, giving $P_{232}/P_{212} \sim e^{-(48.6-18.4)} = e^{-30.2} \approx 8\times 10^{-14}$. This factor of $10^{13}$ in the tunneling probability, combined with differences in the assault frequency and preformation factor, accounts for the $\sim 10^{23}$ ratio in half-lives.

13.5 Alpha Decay Systematics

13.5.1 Where Does Alpha Decay Occur?

Alpha decay is energetically possible ($Q_\alpha > 0$) for essentially all nuclides with $A \gtrsim 150$ and becomes the dominant decay mode for the heaviest elements ($Z \gtrsim 84$). However, for many nuclides between $A \approx 150$ and $A \approx 210$, the Q-values are positive but small, and the resulting half-lives are immensely long — often exceeding the age of the universe.

The lightest alpha emitter observable in the laboratory is ${}^{108}\text{Te}$ ($Z = 52$), and the lightest naturally occurring alpha emitter is ${}^{144}\text{Nd}$ ($Z = 60$), with $Q_\alpha = 1.905\,\text{MeV}$ and $t_{1/2} = 2.29 \times 10^{15}\,\text{yr}$.

An important exception occurs for very light nuclei: ${}^{8}\text{Be}$ ($Z=4$) decays into two alpha particles with a half-life of only $\sim 6.7 \times 10^{-17}\,\text{s}$. This is not tunneling through a Coulomb barrier in the usual sense, but rather the breakup of a barely unbound resonance state.

13.5.2 The Role of the SEMF

The Q-value for alpha decay can be estimated from the semi-empirical mass formula (Chapter 4):

$$Q_\alpha = B(A-4, Z-2) + B({}^4\text{He}) - B(A,Z)$$

Using the SEMF, the Q-value is positive when the binding energy gained by forming the alpha particle exceeds the binding energy lost by the parent. Let us trace how each term in the SEMF contributes to $Q_\alpha$:

Volume term ($a_V A$): The daughter has 4 fewer nucleons, so the volume energy decreases by $4a_V \approx 63\,\text{MeV}$. This opposes alpha emission because the parent has more volume binding.

Surface term ($-a_S A^{2/3}$): The daughter has a smaller surface area. For heavy nuclei, the surface energy decrease is about $2a_S A^{-1/3} \times 4 \approx 12\,\text{MeV}$. This partially compensates the volume loss.

Coulomb term ($-a_C Z^2/A^{1/3}$): This is the key term. The Coulomb energy of the parent scales as $Z^2/A^{1/3}$, while the daughter has $(Z-2)^2/(A-4)^{1/3}$. The difference is approximately $\Delta E_C \approx a_C[4Z - 4]/A^{1/3} \approx 4a_C Z/A^{1/3}$. For ${}^{238}\text{U}$: $\Delta E_C \approx 4 \times 0.711 \times 92/238^{1/3} \approx 42\,\text{MeV}$. This large positive contribution favors alpha emission and dominates for $Z \gtrsim 52$.

Asymmetry term ($-a_A(A-2Z)^2/A$): This term favors alpha emission from neutron-rich nuclei (which are further from the $N = Z$ line) and opposes it from near-stability nuclei.

Pairing term: Enhances Q-values when both parent and daughter are even-even.

The combined effect: for $A \gtrsim 150$, the Coulomb energy release from reducing $Z$ by 2, plus the 28.3 MeV binding energy gained from forming the very tightly bound ${}^4\text{He}$, exceeds the volume energy lost by shrinking the nucleus. This is why $Q_\alpha > 0$ for essentially all nuclei above samarium.

🔗 Connection to Chapter 4: The SEMF Q-value prediction for ${}^{238}\text{U}$ gives $Q_\alpha^{\text{SEMF}} \approx 4.3\,\text{MeV}$, in good agreement with the measured $Q = 4.270\,\text{MeV}$. The SEMF correctly predicts the trend of $Q_\alpha$ with $A$ and $Z$, but misses local shell effects — particularly the enhanced stability of daughters near ${}^{208}\text{Pb}$, which produces anomalously large Q-values and correspondingly short half-lives.

13.5.3 Even-Even Advantage

Alpha decay half-lives exhibit a clear pattern with respect to the nucleon number parity of the parent:

Even-even nuclei have systematically shorter half-lives (faster decay) than neighboring odd-$A$ or odd-odd nuclei with similar Q-values. This is primarily due to the preformation factor: in even-even nuclei, all nucleons are paired, and the four nucleons that will form the alpha particle are more easily assembled into a correlated ${}^4\text{He}$ cluster. Unpaired nucleons disrupt this clustering.

13.5.4 The Natural Decay Series

Three of the four natural radioactive decay series terminate at stable lead isotopes through a sequence of alpha and beta decays:

• Uranium series: ${}^{238}\text{U} \to \cdots \to {}^{206}\text{Pb}$ (8 alpha, 6 beta decays)
• Actinium series: ${}^{235}\text{U} \to \cdots \to {}^{207}\text{Pb}$ (7 alpha, 4 beta decays)
• Thorium series: ${}^{232}\text{Th} \to \cdots \to {}^{208}\text{Pb}$ (6 alpha, 4 beta decays)

Each alpha decay reduces $A$ by 4 and $Z$ by 2. Since $A \mod 4$ is preserved under alpha decay and changes by 0 under beta decay, the three series correspond to $A = 4n+2$ (uranium), $A = 4n+3$ (actinium), and $A = 4n$ (thorium). The fourth series ($A = 4n+1$, the neptunium series starting from ${}^{237}\text{Np}$) is extinct in nature because ${}^{237}\text{Np}$ has a half-life of only $2.14 \times 10^6\,\text{yr}$, much shorter than the age of the Earth.

🔗 Connection to Chapter 12: The Bateman equations governing the activity of each member of a decay series were developed in Chapter 12. The secular equilibrium condition ($\lambda_1 \ll \lambda_2, \lambda_3, \ldots$) is well satisfied in the uranium and thorium series, where the parent half-life vastly exceeds that of any daughter.

All three natural series terminate at stable lead isotopes because ${}^{208}\text{Pb}$, ${}^{207}\text{Pb}$, and ${}^{206}\text{Pb}$ have negative Q-values for both alpha and beta decay. The doubly magic ${}^{208}\text{Pb}$ ($Z = 82$, $N = 126$) is the heaviest stable nucleus — the endpoint of the thorium series and a cornerstone of nuclear structure physics.

13.5.5 Alpha Decay Energies as Precision Nuclear Data

Alpha decay energies are among the most precisely measured quantities in nuclear physics. Modern semiconductor detectors and magnetic spectrographs can determine alpha energies to better than 1 keV precision. This precision has important consequences:

1. Mass determinations. The relation $Q = (M_p - M_d - M_\alpha)c^2$ means that if two of the three masses are known, the third can be determined from the measured Q-value. Alpha decay Q-values have been used to determine the masses of many transuranium isotopes that are difficult to measure by other means.

2. Nuclear level schemes. Fine structure alpha energies (Section 13.7) provide the excitation energies of daughter states to keV precision, complementing gamma-ray spectroscopy.

3. Cosmochronology. The alpha decay half-lives of ${}^{238}\text{U}$, ${}^{235}\text{U}$, and ${}^{232}\text{Th}$ are used in the U-Pb and Th-Pb geochronological dating methods. The precision of radiometric ages depends directly on the accuracy of these half-life measurements. The currently adopted value $t_{1/2}({}^{238}\text{U}) = (4.4683 \pm 0.0048) \times 10^9\,\text{yr}$ anchors the geological timescale.

4. Superheavy element discovery. As noted in Section 13.7.4, the identification of new superheavy elements relies on matching measured alpha decay energies and half-lives to theoretical predictions and known daughter nuclei.

13.6 The Preformation Factor and Beyond the One-Body Model

13.6.1 What Is the Preformation Factor?

The one-body Gamow model assumes the alpha particle exists as a preformed entity bouncing inside the nucleus. In reality, the alpha particle must assemble from two protons and two neutrons within the nuclear medium before it can tunnel out. The probability of finding a pre-existing alpha cluster at the nuclear surface is the preformation factor (or spectroscopic factor) $S$.

For even-even nuclei in the actinide region, microscopic calculations and comparisons with experiment give $S \sim 10^{-1}$--$10^{-2}$. For odd-$A$ nuclei, $S$ is typically an order of magnitude smaller because the unpaired nucleon disrupts the alpha clustering. This explains much of the even-odd staggering in half-lives noted in Section 13.5.3.

13.6.2 Calculating Half-Lives: The Complete Formula

Combining all factors, the alpha decay half-life in the one-body model is:

$$t_{1/2} = \frac{\ln 2}{f \cdot P \cdot S}$$

where:

• Assault frequency: $f = v/(2R)$, with $v = \sqrt{2E_\text{inside}/m_\alpha}$ the velocity of the alpha particle inside the well. Taking $E_\text{inside} \approx Q + V_0 \approx 35\,\text{MeV}$ and $R \approx 8\,\text{fm}$, $f \approx 3 \times 10^{21}\,\text{s}^{-1}$.

• Penetration factor: $P = e^{-G}$ from the WKB calculation.

• Preformation factor: $S \approx 0.1$ for even-even nuclei, $\sim 0.01$ for odd-$A$.

Numerical check for ${}^{238}\text{U}$:

$$t_{1/2} = \frac{0.693}{3\times10^{21}\times 4.3\times10^{-20}\times 0.1} \approx \frac{0.693}{1.3\times10^{1}} \approx 0.053\,\text{s}$$

This is wildly wrong — the measured half-life is $4.47\times10^{9}\,\text{yr} = 1.41\times10^{17}\,\text{s}$. The discrepancy is approximately 18 orders of magnitude.

The resolution is that the simple estimate above is very sensitive to the precise values of $G$, $f$, and $S$. A change in $G$ from 44.7 to $\sim 85$ would bring the prediction in line with experiment. More refined calculations, using better nuclear radii, angular momentum corrections, and properly computed preformation factors, achieve agreement with experiment to within a factor of 2--5 for even-even nuclei — an impressive achievement given that the half-lives themselves span 20+ orders of magnitude.

⚠️ Common Pitfall: Do not be discouraged by large absolute errors in alpha decay calculations. The correct metric is: does the model reproduce the trend across many isotopes and the logarithm of the half-life? A calculation that predicts $\log_{10}t_{1/2}$ to within $\pm 1$ is capturing the essential physics, because the underlying quantity (the Gamow factor) changes continuously and predictably with $Z$ and $E_\alpha$.

13.6.3 Refined Models

Several improvements to the basic Gamow model have been developed:

1. Nuclear potential shape: Replacing the square well with a more realistic Woods-Saxon potential modifies the inner turning point and the assault frequency.

2. Centrifugal barrier: For alpha decay to excited states with angular momentum transfer $\ell > 0$, the effective potential includes a centrifugal term $\ell(\ell+1)\hbar^2/(2\mu r^2)$, which increases the barrier and reduces the tunneling probability.

3. Deformation effects: For deformed nuclei, the barrier height depends on the angle between the alpha emission direction and the nuclear symmetry axis. Prolate deformations generally enhance tunneling along the symmetry axis.

4. Microscopic preformation: Calculations using the nuclear shell model or cluster models to determine $S$ from first principles, rather than treating it as a fitted parameter.

5. Density-dependent cluster formation: The R-matrix approach and the Fliessbach cluster model compute the formation amplitude at the nuclear surface self-consistently.

13.7 Fine Structure in Alpha Decay

13.7.1 The Phenomenon

When alpha emitters are studied with high-resolution spectrometers, the alpha particle energy spectrum often shows not one but several discrete peaks. This fine structure arises because the parent nucleus can decay to the ground state of the daughter or to its excited states.

Consider ${}^{241}\text{Am}$ ($Z = 95$), a well-known alpha emitter used in household smoke detectors:

13.7.2 Why Do Intensities Decrease for Excited States?

Transitions to excited states of the daughter have lower alpha energies because the excitation energy of the daughter must come from the available Q-value:

$$E_\alpha^* = \left(Q - E_x^*\right)\frac{M_d}{M_d + M_\alpha}$$

where $E_x^*$ is the excitation energy of the daughter state. The lower alpha energy means a thicker Coulomb barrier and a smaller tunneling probability. Since the tunneling probability is exponentially sensitive to $E_\alpha$, even a few hundred keV reduction in $E_\alpha$ can reduce the branching ratio by an order of magnitude.

Additionally, if the spin and parity of the daughter excited state require the alpha particle to carry away orbital angular momentum $\ell > 0$, the centrifugal barrier further suppresses the transition. The selection rules are:

$$\vec{I}_\text{parent} = \vec{I}_\text{daughter} + \vec{\ell} + \vec{s}_\alpha$$

Since the alpha particle has spin zero ($s_\alpha = 0$, it is a boson with $J^\pi = 0^+$), the angular momentum transfer equals $\ell$, and:

$$|\Delta I| = |I_p - I_d| \leq \ell \leq I_p + I_d$$

$$\pi_p = \pi_d \cdot (-1)^\ell$$

For even-even parent nuclei ($I_p^\pi = 0^+$), the ground-state-to-ground-state transition has $\ell = 0$ (favored), while transitions to $2^+$, $4^+$, ... excited states require $\ell = 2, 4, \ldots$ (hindered by the centrifugal barrier).

13.7.3 Hindrance Factors

The hindrance factor $\text{HF}$ quantifies how much slower an alpha transition is compared to the prediction of the simple tunneling model (which accounts for the Coulomb and centrifugal barriers but assumes a constant preformation factor). Hindrance factors are defined relative to the ground-state transition:

$$\text{HF} = \frac{t_{1/2}(\text{observed})}{t_{1/2}(\text{predicted from tunneling})}$$

For favored transitions (no change in nuclear structure), $\text{HF} \approx 1$--$4$. For transitions requiring a change in the nuclear configuration (e.g., breaking a pair), $\text{HF}$ can range from $10$ to $10^3$ or more. Hindrance factors provide valuable nuclear structure information: they reveal which daughter states have large overlaps with the parent configuration and which require substantial rearrangement.

13.7.4 Alpha Spectroscopy as a Nuclear Structure Tool

Alpha spectroscopy — measuring the energies and intensities of alpha groups with high precision — has been one of the most productive techniques in nuclear structure physics. The key advantages are:

1. High energy resolution: Semiconductor detectors (silicon surface-barrier or passivated implanted planar silicon, PIPS) achieve energy resolutions of $\sim 15$--$25\,\text{keV}$ FWHM for alpha particles, and magnetic spectrographs can reach $\sim 1$--$3\,\text{keV}$.

2. Absolute energy calibration: Alpha energies from well-known calibration sources (${}^{241}\text{Am}$: 5485.56 keV, ${}^{239}\text{Pu}$: 5156.59 keV, ${}^{244}\text{Cm}$: 5804.77 keV) are known to sub-keV precision.

3. Spin-parity assignments: The pattern of fine-structure intensities, combined with the selection rules of Section 13.7.2, constrains the spin and parity of daughter states.

4. Superheavy element identification: The discovery and identification of superheavy elements ($Z \geq 104$) relies heavily on measuring the alpha decay energies and half-lives of the produced isotopes and their daughters, forming characteristic "alpha decay chains" that link new isotopes to known ones.

🔗 Connection to Chapter 11: The identification of superheavy elements discussed in Chapter 11 depends critically on alpha decay chains. When a new superheavy nucleus is synthesized, it typically undergoes a sequence of alpha decays, each reducing $Z$ by 2 and $A$ by 4, until reaching a known region of the chart of nuclides. The measured $E_\alpha$ and $t_{1/2}$ at each step must be mutually consistent — providing a fingerprint that confirms the discovery.

13.8 Proton Radioactivity

13.8.1 The Concept

Proton radioactivity is the emission of a single proton from a nuclear ground state (or isomeric state). Like alpha decay, it proceeds by quantum tunneling through the Coulomb barrier, but the emitted particle is a single proton ($Z = 1$) rather than an alpha particle ($Z = 2$). The lower charge of the proton means a lower Coulomb barrier, but proton emission is only energetically possible for very proton-rich nuclei near or beyond the proton drip line.

13.8.2 Discovery and Systematics

Proton radioactivity from a nuclear ground state was first observed in 1981 for ${}^{151}\text{Lu}$ (lutetium-151) by Hofmann et al. at GSI Darmstadt. The measured proton energy was $1.233\,\text{MeV}$ and the half-life $85\,\text{ms}$.

Since then, proton radioactivity has been observed for approximately 30 nuclides, spanning the region from $Z = 53$ (${}^{109}\text{I}$) to $Z = 83$ (${}^{185}\text{Bi}$). Key examples:

13.8.3 Theoretical Description

The tunneling calculation is essentially the same as for alpha decay, with $z_\alpha = 1$ and $m_\alpha \to m_p$. The Gamow factor is smaller (lower $Z$ of the emitted particle), and the centrifugal barrier plays a more prominent role because the proton often carries significant angular momentum.

The decay rate is:

$$\lambda = \frac{v}{2R} \cdot P_\ell \cdot S_p$$

where $P_\ell$ includes both Coulomb and centrifugal barrier penetration and $S_p$ is the spectroscopic factor for the proton state. The spectroscopic factor carries direct information about the single-particle wave function of the emitted proton in the parent nucleus, making proton radioactivity a powerful probe of nuclear structure far from stability.

💡 Intuition: Proton radioactivity is to single-particle shell structure what alpha decay is to clustering: measuring the rate tells you about the nuclear wave function at the surface.

13.8.4 Nuclear Structure Information from Proton Emission

Proton radioactivity is remarkably clean as a nuclear structure probe. Unlike alpha decay, where the preformation factor introduces substantial uncertainty, the emitted proton in proton radioactivity was already a proton inside the parent nucleus. The spectroscopic factor $S_p$ is directly related to the occupation probability of the single-particle orbital from which the proton is emitted. Comparing measured half-lives to tunneling calculations with known $\ell$ values yields spectroscopic factors that can be compared to shell model predictions.

For example, the ground-state proton emission of ${}^{151}\text{Lu}$ involves a proton in the $h_{11/2}$ orbital ($\ell = 5$). The measured half-life of 85 ms, combined with the WKB tunneling calculation, gives a spectroscopic factor $S_p \approx 0.5$ — consistent with shell model predictions for a nearly pure single-particle state. In contrast, ${}^{185}\text{Bi}$ emits a proton from an $s_{1/2}$ orbital ($\ell = 0$), and the absence of a centrifugal barrier makes the decay much faster despite the higher Coulomb barrier.

💡 Intuition: Proton radioactivity is to single-particle shell structure what alpha decay is to clustering: measuring the rate tells you about the nuclear wave function at the surface. The directness of this connection — without the complications of preformation — makes proton emission a uniquely valuable probe.

13.8.5 Two-Proton Radioactivity

A more exotic variant is two-proton radioactivity, in which two protons are emitted simultaneously. This occurs in nuclei where single-proton emission is energetically forbidden (the nucleus is still inside the proton drip line for one-proton emission) but emission of a proton pair is energetically allowed because of the pairing energy. The two protons must tunnel through the barrier as a correlated pair — they are emitted simultaneously, not sequentially.

First confirmed experimentally for ${}^{45}\text{Fe}$ by Pfützner et al. (2002) at GSI and Giovinazzo et al. (2002) at GANIL, two-proton radioactivity has since been observed in several additional nuclides: ${}^{19}\text{Mg}$ (the lightest known two-proton emitter), ${}^{48}\text{Ni}$, ${}^{54}\text{Zn}$, and ${}^{67}\text{Kr}$. The angular correlations between the two emitted protons carry information about the proton-proton interaction inside the nucleus and about the geometry of the decay — whether the two protons are emitted "back-to-back" or in a ${}^2\text{He}$-like correlated pair ("diproton" emission).

This remains an active area of research, with experiments at radioactive beam facilities (FRIB, RIKEN, GANIL) producing ever-more exotic proton-rich nuclei to study.

13.9 Cluster Radioactivity

13.9.1 Discovery

In 1984, H.J. Rose and G.A. Jones at the University of Oxford observed the emission of ${}^{14}\text{C}$ nuclei from ${}^{223}\text{Ra}$:

$${}^{223}\text{Ra} \to {}^{209}\text{Pb} + {}^{14}\text{C}$$

with a branching ratio of approximately $8.5 \times 10^{-10}$ relative to alpha decay. The daughter ${}^{209}\text{Pb}$ is one neutron beyond the doubly magic ${}^{208}\text{Pb}$, and this proximity to shell closure plays a crucial role in making the decay observable.

13.9.2 Systematics

Since 1984, cluster radioactivity has been observed for approximately 25 parent-cluster combinations. The emitted clusters range from ${}^{14}\text{C}$ to ${}^{34}\text{Si}$. A striking pattern emerges: the daughter nucleus is always close to ${}^{208}\text{Pb}$ (or its immediate neighbors). This reflects the enormous binding energy gain from forming a nucleus near the doubly magic shell closure at $Z = 82$, $N = 126$.

Selected examples:

13.9.3 Theoretical Description

Cluster radioactivity can be treated within the same tunneling framework as alpha decay. The decay rate is:

$$\lambda = f \cdot P \cdot S_c$$

where $S_c$ is the preformation probability for the cluster inside the parent nucleus and $P$ is the WKB penetration factor. The key difference from alpha decay is that the preformation factor $S_c$ is much smaller — typically $10^{-5}$--$10^{-15}$ depending on the cluster size — because assembling 14 or more nucleons into a coherent cluster is far less probable than assembling just 4 into an alpha particle.

Two theoretical frameworks have been particularly successful:

1. Preformed cluster model (Buck, Merchant, and Pepin): Treats the cluster as preformed with a spectroscopic factor $S_c$ determined empirically or from microscopic calculations. The tunneling is then a standard WKB problem.

2. Fission-like model (Sandulescu, Poenaru, and Greiner): Treats cluster emission as an extremely asymmetric fission process, using the liquid-drop energy surface plus shell corrections to compute the deformation barrier. This model predicted cluster radioactivity (Sandulescu, Poenaru, and Greiner, 1980) four years before its experimental observation.

💡 Intuition: Cluster radioactivity interpolates between alpha decay (emission of ${}^4\text{He}$) and asymmetric fission (emission of a fragment with $A \sim 100$). As the emitted cluster gets heavier, the tunneling probability decreases (heavier particle, larger Coulomb barrier), but the preformation factor is the dominant suppression mechanism.

13.9.4 The Role of Shell Effects

The invariable proximity of the daughter to ${}^{208}\text{Pb}$ is not a coincidence. The doubly magic shell closure at $Z = 82$, $N = 126$ provides a large Q-value enhancement (the daughter is very tightly bound), which partially compensates for the reduced tunneling probability of the heavy cluster. Without this shell effect, cluster radioactivity would be unobservably rare.

🔗 Connection to Chapters 4 and 6: The role of ${}^{208}\text{Pb}$ in cluster radioactivity beautifully illustrates the interplay between bulk nuclear properties (the liquid drop, Chapter 4) and quantum shell effects (Chapter 6). The liquid drop model alone would not predict the observed clustering around ${}^{208}\text{Pb}$; it requires the shell correction energy.

13.10 Summary of Key Results

The Physics in Five Equations

1. Alpha decay Q-value:

$$Q_\alpha = \left[M({}^A_Z\text{X}) - M({}^{A-4}_{Z-2}\text{Y}) - M({}^4_2\text{He})\right]c^2$$

2. WKB tunneling probability:

$$P = \exp(-G), \quad G = \frac{2}{\hbar}\int_R^b\sqrt{2\mu[V(r)-E]}\,dr$$

3. Gamow factor (Coulomb barrier):

$$G = 2\eta\left[\arccos\sqrt{\rho} - \sqrt{\rho(1-\rho)}\right], \quad \rho = R/b, \quad \eta = \frac{z_\alpha Z_d e^2}{4\pi\epsilon_0\hbar v}$$

4. Thick-barrier approximation:

$$G \approx \frac{2\pi z_\alpha Z_d e^2}{4\pi\epsilon_0\hbar v} - \frac{4}{\hbar}\sqrt{2\mu R\cdot\frac{z_\alpha Z_d e^2}{4\pi\epsilon_0}}$$

5. Geiger-Nuttall law:

$$\log_{10}\lambda = A + B\,\frac{Z_d}{\sqrt{E_\alpha}}$$

The Threshold Concept

The alpha particle does not go "over" the Coulomb barrier. It goes "through" it. Quantum tunneling converts the classically impossible into the quantum-mechanically improbable — and that improbability, governed by the Gamow factor, explains why half-lives for alpha decay range from nanoseconds to times vastly exceeding the age of the universe, while alpha energies vary over a surprisingly narrow range.

13.11 Project Checkpoint: Alpha Tunneling Calculator

This chapter's contribution to the Nuclear Data Analysis Toolkit is alpha_tunneling.py, which:

1. Calculates the WKB tunneling probability through the Coulomb barrier for user-specified parent nuclei.
2. Computes the Gamow factor and estimated half-life for a set of even-even alpha emitters.
3. Reproduces the Geiger-Nuttall law by plotting $\log_{10}\lambda$ vs. $Z_d/\sqrt{E_\alpha}$.
4. Compares calculated half-lives to measured values.

See the code/ directory for the complete, runnable script and code/project-checkpoint.md for documentation.

💻 Computational Note: The script requires numpy and matplotlib. Run with: python alpha_tunneling.py. It produces three plots: (1) the Coulomb barrier potential with tunneling region shaded, (2) calculated vs. measured half-lives, and (3) the Geiger-Nuttall plot.

What's Next

In Chapter 14, we turn to beta decay — the process by which a neutron converts into a proton (or vice versa) inside the nucleus, emitting an electron (or positron) and a neutrino. Unlike alpha decay, which is governed by the strong nuclear and electromagnetic forces, beta decay is a manifestation of the weak interaction — the force responsible for changing quark flavors. Fermi's theory of beta decay (1934), the neutrino hypothesis (Pauli, 1930), and the discovery of parity violation in the Wu experiment (1957) constitute some of the greatest chapters in the history of physics.