Case Study 2 — From Alpha Emission to Cluster Radioactivity

A Prediction Fulfilled

In 1980, Aurel Sandulescu, Dorin N. Poenaru, and Walter Greiner published a theoretical paper with a bold prediction: certain heavy nuclei should be able to emit fragments much heavier than alpha particles -- nuclear clusters such as ${}^{14}\text{C}$, ${}^{20}\text{O}$, ${}^{24}\text{Ne}$, and ${}^{28}\text{Mg}$ -- via the same quantum tunneling mechanism that governs alpha decay. They modeled the process as extremely asymmetric fission, using the liquid-drop energy surface with shell corrections to identify candidates with favorable Q-values and manageable barriers. Their calculations predicted that ${}^{223}\text{Ra}$ should emit ${}^{14}\text{C}$ with a partial half-life of approximately $10^{15}\,\text{s}$ and a branching ratio of roughly $10^{-9}$ relative to alpha decay.

Four years later, in 1984, H.J. Rose and G.A. Jones at the University of Oxford set out to test this prediction experimentally.

The Experiment

Rose and Jones used a simple but elegant setup. A thin source of ${}^{223}\text{Ra}$ (half-life 11.43 days, produced from ${}^{227}\text{Ac}$) was placed inside a vacuum chamber facing a silicon semiconductor detector. The detector recorded the energy of every charged particle emitted by the source. Alpha particles from ${}^{223}\text{Ra}$ appeared as a prominent peak near 5.7 MeV. The ${}^{14}\text{C}$ cluster, if emitted, would carry away a kinetic energy of approximately 30 MeV -- well separated from the alpha peak but arriving at the detector as a much rarer event.

The challenge was statistics. With a branching ratio predicted to be $\sim 10^{-9}$, the experimenters needed to detect roughly 1 carbon-14 event for every billion alpha particles. Over a counting period of several weeks, they accumulated enough statistics to identify a small but statistically significant cluster of events at the predicted energy.

The result: ${}^{14}\text{C}$ emission from ${}^{223}\text{Ra}$ was confirmed, with a measured branching ratio of $(8.5 \pm 2.5) \times 10^{-10}$ -- in remarkable agreement with the Sandulescu-Poenaru-Greiner prediction.

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

The decay left the daughter ${}^{209}\text{Pb}$, one neutron beyond the doubly magic ${}^{208}\text{Pb}$ -- a crucial detail for understanding why this particular decay channel was observable.

The Physics: Why ${}^{208}\text{Pb}$ Matters

The discovery of cluster radioactivity raised an immediate question: why ${}^{14}\text{C}$ from ${}^{223}\text{Ra}$ and not, say, ${}^{12}\text{C}$ or ${}^{16}\text{O}$? The answer lies in the interplay between the Q-value (determined by nuclear masses) and the preformation/tunneling probabilities.

The Q-value for any particle emission is:

$$Q = [M(\text{parent}) - M(\text{daughter}) - M(\text{cluster})]c^2$$

For cluster emission to have a measurable rate, three conditions must be met simultaneously:

1. $Q > 0$: The decay must be energetically allowed.
2. $Q$ must be large: A larger Q-value gives a thinner barrier and higher tunneling probability.
3. The preformation factor must not be vanishingly small: The cluster must have a reasonable probability of forming inside the parent nucleus.

The doubly magic nucleus ${}^{208}\text{Pb}$ ($Z = 82$, $N = 126$) has an extraordinarily high binding energy per nucleon due to shell closure effects. When the daughter nucleus is at or near ${}^{208}\text{Pb}$, the Q-value receives a bonus of several MeV from the shell-correction energy -- enough to partially compensate for the much smaller tunneling probability of a heavy cluster compared to an alpha particle.

This explains the striking empirical pattern: every observed cluster radioactivity decay produces a daughter within a few nucleons of ${}^{208}\text{Pb}$. The shell effect is not merely a refinement; it is the enabling mechanism.

The Expanding Landscape

Following the Rose-Jones discovery, experimentalists at laboratories around the world searched for additional cluster decay modes. By the mid-2000s, approximately 25 parent-cluster combinations had been observed:

Carbon-14 emission: ${}^{221}\text{Fr}$, ${}^{221,222,223,224,226}\text{Ra}$, ${}^{225}\text{Ac}$

Oxygen and neon emission: ${}^{228}\text{Th} \to {}^{208}\text{Pb} + {}^{20}\text{O}$, ${}^{230}\text{U} \to {}^{208}\text{Pb} + {}^{22}\text{Ne}$, ${}^{231}\text{Pa} \to {}^{207}\text{Tl} + {}^{24}\text{Ne}$, ${}^{234}\text{U} \to {}^{210}\text{Pb} + {}^{24}\text{Ne}$

Magnesium and silicon emission: ${}^{236}\text{Pu} \to {}^{208}\text{Pb} + {}^{28}\text{Mg}$, ${}^{238}\text{Pu} \to {}^{210}\text{Pb} + {}^{28}\text{Mg}$, ${}^{242}\text{Cm} \to {}^{208}\text{Pb} + {}^{34}\text{Si}$

In every case, the daughter is within a few mass units of ${}^{208}\text{Pb}$.

The heaviest observed cluster is ${}^{34}\text{Si}$ from ${}^{242}\text{Cm}$, with a branching ratio of approximately $10^{-16}$ relative to alpha decay. As the emitted cluster gets heavier, two competing effects determine the rate:

• The tunneling probability decreases dramatically (heavier cluster, higher Coulomb barrier charge product $z_c Z_d$).
• The preformation factor decreases even more steeply (assembling 34 nucleons into a coherent cluster is far less probable than assembling 4).

At some point, the cluster emission rate becomes too small to detect. Whether there is a theoretical limit to the cluster mass, or whether this is purely an experimental sensitivity issue, remains an open question.

Theoretical Frameworks

Two complementary theoretical approaches have been applied to cluster radioactivity:

The Preformed Cluster Model

This approach, developed by Buck, Merchant, and Pepin in the 1990s, treats the cluster as a preexisting entity inside the parent nucleus, just as the Gamow model treats the alpha particle. The decay rate is:

$$\lambda = f \cdot P_{\text{WKB}} \cdot S_c$$

where $S_c$ is the spectroscopic factor (preformation probability) for the cluster. For ${}^{14}\text{C}$ emission, $S_c \sim 10^{-8}$ -- eight orders of magnitude smaller than for alpha particles. The WKB tunneling probability $P_{\text{WKB}}$ is calculated exactly as for alpha decay, with the alpha mass and charge replaced by those of the cluster.

This model is conceptually straightforward and computationally tractable. It reproduces experimental branching ratios to within a factor of 10--100 for most observed decays.

The Asymmetric Fission Model

Sandulescu, Poenaru, and Greiner's original approach views cluster radioactivity as an extreme case of asymmetric fission. The parent nucleus is imagined to deform along a fission-like path, but instead of dividing into two roughly equal fragments (symmetric fission), it divides into one very heavy fragment ($\sim 208$ nucleons) and one very light fragment (the cluster). The fission barrier is calculated using the liquid-drop model with shell corrections.

This model has the advantage of providing a unified framework for alpha decay, cluster radioactivity, and conventional fission as points on a continuous spectrum of binary decay processes:

From this perspective, cluster radioactivity fills the gap between alpha decay and fission on the fragmentation spectrum. The continuous variation of barrier and preformation parameters across this spectrum is one of the most elegant features of nuclear decay physics.

The fission model also provides a natural explanation for the role of ${}^{208}\text{Pb}$: the shell correction energy at the doubly magic closure creates a deep valley in the potential energy surface at the highly asymmetric deformation corresponding to a ${}^{208}\text{Pb}$-like fragment. The fission barrier along this asymmetric path is lowered by the shell effect, enhancing the decay rate by many orders of magnitude compared to what the liquid drop model alone would predict. Without this shell effect, all cluster radioactivity half-lives would be unmeasurably long.

Quantitative Analysis: Why the ${}^{14}\text{C}$ Rate Is What It Is

Let us trace through the physics quantitatively for the ${}^{223}\text{Ra} \to {}^{209}\text{Pb} + {}^{14}\text{C}$ decay.

Kinematics: - $Q = 31.83\,\text{MeV}$ (from AME2020 masses) - $T_{14\text{C}} = 31.83 \times 209/223 = 29.83\,\text{MeV}$ - $T_d = 31.83 \times 14/223 = 2.00\,\text{MeV}$

Barrier parameters: - Daughter charge: $Z_d = 82$ (lead) - Cluster charge: $z_c = 6$ (carbon) - $R = 1.2(209^{1/3} + 14^{1/3}) = 1.2(5.935 + 2.410) = 10.01\,\text{fm}$ - $V_C(R) = 6 \times 82 \times 1.440/10.01 = 70.8\,\text{MeV}$ - $b = 6 \times 82 \times 1.440/29.83 = 23.74\,\text{fm}$ - $\rho = 10.01/23.74 = 0.422$

The barrier height is 70.8 MeV while the kinetic energy is 29.8 MeV — a ratio of 2.4, which is actually smaller than the typical alpha decay ratio of 5--7. However, the reduced mass is much larger: $\mu = 14 \times 209/(14+209) \times 931.5 = 12.22 \times 931.5 = 11,383\,\text{MeV}/c^2$, versus $\mu \approx 3,700\,\text{MeV}/c^2$ for alpha decay. The heavier particle mass dramatically increases the Gamow factor.

Gamow factor: The Sommerfeld parameter is $\eta \approx 42.5$, and:

$$G = 2 \times 42.5 \times [\arccos(\sqrt{0.422}) - \sqrt{0.422 \times 0.578}] \approx 85 \times [0.865 - 0.494] \approx 31.5$$

This gives $P = e^{-31.5} \approx 2 \times 10^{-14}$. For comparison, the alpha decay tunneling probability from ${}^{223}\text{Ra}$ (with $E_\alpha \approx 5.87\,\text{MeV}$) is $P_\alpha \approx 10^{-8}$. The cluster tunneling probability is about $10^6$ times smaller.

However, the preformation factor for ${}^{14}\text{C}$ inside ${}^{223}\text{Ra}$ is estimated at $S_c \sim 10^{-8}$, versus $S_\alpha \sim 10^{-1}$ for the alpha particle. This additional $10^7$ suppression, combined with the $10^6$ reduction in tunneling probability, gives a total suppression of $\sim 10^{13}$. The observed branching ratio is $\sim 10^{-10}$, which is roughly consistent given the uncertainties in the preformation factors and the simplicity of the model.

The quantitative message: the rarity of cluster radioactivity is dominated by the preformation factor (assembling 14 nucleons into a correlated cluster is vastly less probable than assembling 4), with the tunneling suppression playing a secondary but significant role.

Modern Developments

Research on cluster radioactivity continues in several directions:

1. Search for new cluster emitters. Theoretical calculations predict additional cluster decay modes, particularly for neutron-deficient actinides and transactinides. Experimental searches using improved detectors and more intense radioactive sources are ongoing.

2. Connection to superheavy elements. Some superheavy nuclei ($Z \geq 110$) may have competing alpha and cluster decay channels. Observing cluster emission from superheavy elements would provide unique information about nuclear structure far from stability.

3. Nuclear structure probes. The measured branching ratios for cluster emission carry information about the clustering correlations inside the parent nucleus -- how likely are 14 or 24 nucleons to form a correlated cluster? This connects to fundamental questions about nuclear many-body theory.

4. Two-proton radioactivity and beyond. The discovery of two-proton emission from ${}^{45}\text{Fe}$ (Pfützner et al., 2002; Giovinazzo et al., 2002) demonstrated another exotic tunneling mode. The theoretical framework connects directly to cluster radioactivity -- both involve correlated multi-nucleon emission.

The Big Picture: A Spectrum of Tunneling

Cluster radioactivity revealed that nature's repertoire of tunneling-mediated nuclear decay modes is far richer than the traditional categories (alpha, beta, gamma, fission) suggest. The continuum of possibilities -- from single-proton emission through alpha decay, cluster emission, and on to asymmetric and symmetric fission -- is unified by the same quantum-mechanical principle: a system confined by a potential barrier has a nonzero probability of tunneling through it, with the rate exponentially sensitive to the barrier height, width, and the mass of the tunneling object.

Experimental Challenges

Detecting cluster radioactivity poses severe experimental challenges that push the boundaries of nuclear instrumentation:

1. Extremely rare events. With branching ratios of $10^{-9}$ to $10^{-16}$ relative to alpha decay, the experimenter must detect one cluster event among billions or trillions of alpha particles. This requires either very long counting times, very strong sources, or both.

2. Background discrimination. The ${}^{14}\text{C}$ signal at $\sim 30\,\text{MeV}$ is well separated in energy from the alpha particles at $\sim 5$--$6\,\text{MeV}$, but pile-up events (two alpha particles arriving simultaneously in the detector) can mimic a single high-energy event. Careful pulse-shape analysis is required to distinguish genuine cluster events from pile-up.

3. Particle identification. Knowing that a 30 MeV particle hit the detector is not enough — one must confirm that it is ${}^{14}\text{C}$ and not some other particle. This is achieved using $\Delta E$-$E$ telescope detectors (two detectors in series, measuring both $dE/dx$ and total energy, which together identify the charge and mass of the particle) or magnetic spectrographs.

4. Source preparation. The radioactive source must be thin enough that the emitted cluster escapes without losing significant energy, yet contain enough atoms to produce a detectable rate. For the rarest decays (branching ratios $\lesssim 10^{-14}$), milligram-scale sources of transuranium elements are needed — requiring access to reactor or accelerator production facilities.

Rose and Jones's original 1984 experiment used a relatively simple Si surface-barrier detector and accumulated data over several weeks. Modern experiments use multi-detector arrays and sophisticated digital signal processing to push the sensitivity to branching ratios below $10^{-15}$.

The Big Picture: A Spectrum of Tunneling

Cluster radioactivity revealed that nature's repertoire of tunneling-mediated nuclear decay modes is far richer than the traditional categories (alpha, beta, gamma, fission) suggest. The continuum of possibilities -- from single-proton emission through alpha decay, cluster emission, and on to asymmetric and symmetric fission -- is unified by the same quantum-mechanical principle: a system confined by a potential barrier has a nonzero probability of tunneling through it, with the rate exponentially sensitive to the barrier height, width, and the mass of the tunneling object.

This unification is one of the most elegant features of quantum mechanics applied to nuclear physics. The same Gamow factor that Gamow derived for alpha decay in 1928 — generalized to different emitted particles — explains the full spectrum of charged-particle radioactivity nearly a century later.

Discussion Questions

1. The prediction of cluster radioactivity preceded its discovery by four years. How does this compare to other predictions in nuclear/particle physics (e.g., the neutrino, the omega-minus, the Higgs boson)? What role did the tunneling framework play in enabling the prediction?

2. Why are all observed cluster radioactivity daughters near ${}^{208}\text{Pb}$? Could cluster radioactivity exist near other magic numbers (e.g., daughters near ${}^{132}\text{Sn}$, which is doubly magic at $Z = 50$, $N = 82$)? What would the Q-values and barriers look like?

3. The heaviest observed cluster emission is ${}^{34}\text{Si}$ from ${}^{242}\text{Cm}$. Is there a fundamental upper limit on the cluster mass, or is the limit purely experimental? How would you estimate the maximum cluster mass that could be observed with foreseeable improvements in detector sensitivity?

4. The fission model views cluster radioactivity as extremely asymmetric fission. At what mass asymmetry does the conceptual distinction between "cluster emission" and "fission" break down? Is this distinction physically meaningful, or is it merely conventional?