Cluster decay

The first information about the atomic nucleus was obtained at the beginning of the 20th century by studying radioactivity. For a long period of time only three kinds of nuclear decay modes (alpha, beta, and gamma) were known. They illustrate three of the fundamental interactions in nature: strong, weak, and electromagnetic. Spontaneous fission became better studied soon after its discovery in 1940 by Konstantin Petrzhak and Georgy Flyorov because of both the military and the peaceful applications of induced fission. This was discovered circa 1939 by Otto Hahn, Lise Meitner, and Fritz Strassmann.

There are many other kinds of radioactivity, e.g. cluster decay, proton emission, various beta-delayed decay modes (p, 2p, 3p, n, 2n, 3n, 4n, d, t, alpha, f), fission isomers, particle accompanied (ternary) fission, etc. The height of the potential barrier, mainly of Coulomb nature, for emission of the charged particles is much higher than the observed kinetic energy of the emitted particles. The spontaneous decay can only be explained by quantum tunneling in a similar way to the first application of the Quantum Mechanics to Nuclei given by G. Gamow for alpha decay.

In 1980 A. Sandulescu, D.N. Poenaru, and W. Greiner described calculations indicating the possibility of a new type of decay of heavy nuclei intermediate between alpha decay and spontaneous fission. The first observation of heavy-ion radioactivity was that of a 30-MeV, carbon-14 emission from radium-223 by H.J. Rose and G.A. Jones in 1984.

— Encyclopædia Britannica, [4]

Usually the theory explains an already experimentally observed phenomenon. Cluster decay is one of the rare examples of phenomena predicted before experimental discovery. Theoretical predictions were made in 1980,[5] four years before experimental discovery.[6]

Four theoretical approaches were used: fragmentation theory by solving a Schrödinger equation with mass asymmetry as a variable to obtain the mass distributions of fragments; penetrability calculations similar to those used in traditional theory of alpha decay, and superasymmetric fission models, numerical (NuSAF) and analytical (ASAF). Superasymmetric fission models are based on the macroscopic-microscopic approach[7] using the asymmetrical two-center shell model[8][9] level energies as input data for the shell and pairing corrections. Either the liquid drop model[10] or the Yukawa-plus-exponential model[11] extended to different charge-to-mass ratios[12] have been used to calculate the macroscopic deformation energy.

Penetrability theory predicted eight decay modes: ¹⁴C, ²⁴Ne, ²⁸Mg, ^32,34Si, ⁴⁶Ar, and ^48,50Ca from the following parent nuclei: ^222,224Ra, ^230,232Th, ^236,238U, ^244,246Pu, ^248,250Cm, ^250,252Cf, ^252,254Fm, and ^252,254No.

The first experimental report was published in 1984, when physicists at Oxford University discovered that ²²³Ra emits one ¹⁴C nucleus among every billion (10⁹) decays by alpha emission.

The quantum tunneling may be calculated either by extending fission theory to a larger mass asymmetry or by heavier emitted particle from alpha decay theory.[13]

Both fission-like and alpha-like approaches are able to express the decay constant ${\displaystyle \lambda =\ln 2/T_{\text{c}}}$ , as a product of three model-dependent quantities

${\displaystyle \lambda =\nu SP_{\text{s}}}$

where ${\displaystyle \nu }$ is the frequency of assaults on the barrier per second, S is the preformation probability of the cluster at the nuclear surface, and P_s is the penetrability of the external barrier. In alpha-like theories S is an overlap integral of the wave function of the three partners (parent, daughter, and emitted cluster). In a fission theory the preformation probability is the penetrability of the internal part of the barrier from the initial turning point R_i to the touching point R_t.[14] Very frequently it is calculated by using the Wentzel-Kramers-Brillouin (WKB) approximation.

A very large number, of the order 10⁵, of parent-emitted cluster combinations were considered in a systematic search for new decay modes. The large amount of computations could be performed in a reasonable time by using the ASAF model developed by Dorin N Poenaru, Walter Greiner, et al. The model was the first to be used to predict measurable quantities in cluster decay. More than 150 cluster decay modes have been predicted before any other kind of half-lives calculations have been reported. Comprehensive tables of half-lives, branching ratios, and kinetic energies have been published, e.g.[15][16] Potential barrier shapes similar to that considered within the ASAF model have been calculated by using the macroscopic-microscopic method.[17]

Previously[18] it was shown that even alpha decay may be considered a particular case of cold fission. The ASAF model may be used to describe in a unified manner cold alpha decay, cluster decay, and cold fission (see figure 6.7, p. 287 of the Ref. [2]).

One can obtain with good approximation one universal curve (UNIV) for any kind of cluster decay mode with a mass number Ae, including alpha decay

${\displaystyle \log T=-\log P_{s}-22.169+0.598(A_{e}-1)}$

In a logarithmic scale the equation log T = f(log P_s) represents a single straight line which can be conveniently used to estimate the half-life. A single universal curve for alpha decay and cluster decay modes results by expressing log T + log S = f(log P_s).[19] The experimental data on cluster decay in three groups of even-even, even-odd, and odd-even parent nuclei are reproduced with comparable accuracy by both types of universal curves, fission-like UNIV and UDL[20] derived using alpha-like R-matrix theory.

In order to find the released energy

${\displaystyle Q=[M-(M_{d}+M_{e})]c^{2}}$

one can use the compilation of measured masses[21] M, M_d, and M_e of the parent, daughter, and emitted nuclei, c is the light velocity. The mass excess is transformed into energy according to the Einstein's formula E = mc².

The main experimental difficulty in observing cluster decay comes from the need to identify a few rare events against a background of alpha particles. The quantities experimentally determined are the partial half life, T_c, and the kinetic energy of the emitted cluster E_k. There is also a need to identify the emitted particle.

Detection of radiations is based on their interactions with matter, leading mainly to ionizations. Using a semiconductor telescope and conventional electronics to identify the ¹⁴C ions, the Rose and Jones's experiment was running for about six months in order to get 11 useful events.

With modern magnetic spectrometers (SOLENO and Enge-split pole), at Orsay and Argonne National Laboratory (see ch. 7 in Ref. [2] pp. 188–204), a very strong source could be used, so that results were obtained in a run of few hours.

Solid state nuclear track detectors (SSNTD) insensitive to alpha particles and magnetic spectrometers in which alpha particles are deflected by a strong magnetic field have been used to overcome this difficulty. SSNTD are cheap and handy but they need chemical etching and microscope scanning.

A key role in experiments on cluster decay modes performed in Berkeley, Orsay, Dubna, and Milano was played by P. Buford Price, Eid Hourany, Michel Hussonnois, Svetlana Tretyakova, A. A. Ogloblin, Roberto Bonetti, and their coworkers.

The main region of 20 emitters experimentally observed until 2010 is above Z = 86: ²²¹Fr, ^221-224,226Ra, ^223,225Ac, ^228,230Th, ²³¹Pa, ^230,232-236U, ^236,238Pu, and ²⁴²Cm. Only upper limits could be detected in the following cases: ¹²C decay of ¹¹⁴Ba, ¹⁵N decay of ²²³Ac, ¹⁸O decay of ²²⁶Th, ^24,26Ne decays of ²³²Th and of ²³⁶U, ²⁸Mg decays of ^232,233,235U, ³⁰Mg decay of ²³⁷Np, and ³⁴Si decay of ²⁴⁰Pu and of ²⁴¹Am.

Some of the cluster emitters are members of the three natural radioactive families. Others should be produced by nuclear reactions. Up to now no odd-odd emitter has been observed.

From many decay modes with half-lives and branching ratios relative to alpha decay predicted with the analytical superasymmetric fission (ASAF) model, the following 11 have been experimentally confirmed: ¹⁴C, ²⁰O, ²³F, ^22,24-26Ne, ^28,30Mg, and ^32,34Si. The experimental data are in good agreement with predicted values. A strong shell effect can be seen: as a rule the shortest value of the half-life is obtained when the daughter nucleus has a magic number of neutrons (N_d = 126) and/or protons (Z_d = 82).

The known cluster emissions as of 2010 are as follows:[22][23][24]

The fine structure in ¹⁴C radioactivity of ²²³Ra was discussed for the first time by M. Greiner and W. Scheid in 1986.[25] The superconducting spectrometer SOLENO of IPN Orsay has been used since 1984 to identify ¹⁴C clusters emitted from ^{222–224,226}Ra nuclei. Moreover, it was used to discover[26][27] the fine structure observing transitions to excited states of the daughter. A transition with an excited state of ¹⁴C predicted in Ref. [25] was not yet observed.

Surprisingly, the experimentalists had seen a transition to the first excited state of the daughter stronger than that to the ground state. The transition is favoured if the uncoupled nucleon is left in the same state in both parent and daughter nuclei. Otherwise the difference in nuclear structure leads to a large hindrance.

The interpretation[28] was confirmed: the main spherical component of the deformed parent wave function has an i_11/2 character, i.e. the main component is spherical.

• National Nuclear Data Center