Cluster decay

Cluster decay, also named heavy particle radioactivity, heavy ion radioactivity or heavy cluster decay,^[1] is a rare type of nuclear decay in which an atomic nucleus emits a small "cluster" of neutrons and protons, more than in an alpha particle, but less than a typical binary fission fragment. Ternary fission into three fragments also produces products in the cluster size.

Description

The loss of protons from the parent nucleus changes it to the nucleus of a different element, the daughter, with a mass number A_d = A − A_e and atomic number Z_d = Z − Z_e, where A_e = N_e + Z_e.^[2] For example:

²²³
₈₈Ra → ¹⁴
₆C + ²⁰⁹
₈₂Pb

According to "Ronen's golden rule" of cluster decay, the emitted nucleus tends to be one with a high binding energy per nucleon, and especially one with a magic number of nucleons.^[3]

This type of rare decay mode was observed in radioisotopes that decay predominantly by alpha emission, and it occurs only in a small percentage of the decays for all such isotopes.^[4]

The branching ratio with respect to alpha decay is rather small (see the Table below):

${\displaystyle B=T_{a}/T_{c}}$

where T_a and T_c are the half-lives of the parent nucleus relative to alpha decay and cluster radioactivity, respectively.

Cluster decay, like alpha decay, is a quantum tunneling process: in order to be emitted, the cluster must penetrate a potential barrier. This is a different process than the more random nuclear disintegration that precedes light fragment emission in ternary fission, which may be a result of a nuclear reaction, but can also be a type of spontaneous radioactive decay in certain nuclides, demonstrating that input energy is not necessarily needed for fission, which remains a fundamentally different process mechanistically.

In the absence of any energy loss for fragment deformation and excitation, as in cold fission phenomena or in alpha decay, the total kinetic energy is equal to the Q-value and is divided between the particles in inverse proportion with their masses, as required by conservation of linear momentum

${\displaystyle E_{k}=QA_{d}/A}$

where A_d is the mass number of the daughter, A_d = A − A_e.

Cluster decay exists in an intermediate position between alpha decay (in which a nucleus spits out a ⁴He nucleus), and spontaneous fission, in which a heavy nucleus splits into two (or more) large fragments and an assorted number of neutrons. Spontaneous fission ends up with a probabilistic distribution of daughter products, which sets it apart from cluster decay. In cluster decay for a given radioisotope, the emitted particle is a light nucleus and the decay method always emits this same particle. For heavier emitted clusters, there is otherwise practically no qualitative difference between cluster decay and spontaneous cold fission.

History

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, ^[5]

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,^[6] four years before experimental discovery.^[7]

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^[8] using the asymmetrical two-center shell model^[9][10] level energies as input data for the shell and pairing corrections. Either the liquid drop model^[11] or the Yukawa-plus-exponential model^[12] extended to different charge-to-mass ratios^[13] 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.^[14]

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.

Theory

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.^[15]

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.^[16] 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.^[17][18] Potential barrier shapes similar to that considered within the ASAF model have been calculated by using the macroscopic-microscopic method.^[19]

Previously^[20] 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).^[21] 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^[22] 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^[23] 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².

Experiments

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 in the range 87 ≤ Z ≤ 96: ²²¹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:^[24][25][26]

Isotope	Emitted particle	Branching ratio	log T(s)	Q (MeV)
221Fr	14C	8.14×10−13	14.52	31.290
221Ra	14C	1.15×10−12	13.39	32.394
222Ra	14C	3.7×10−10	11.01	33.049
223Ra	14C	8.9×10−10	15.04	31.829
224Ra	14C	4.3×10−11	15.86	30.535
225Ac	14C	4.5×10−12	17.28	30.476
226Ra	14C	3.2×10−11	21.19	28.196
228Th	20O	1.13×10−13	20.72	44.723
230Th	24Ne	5.6×10−13	24.61	57.758
231Pa	23F	9.97×10−15	26.02	51.844
	24Ne	1.34×10−11	22.88	60.408
230U	22Ne	4.8×10−14	19.57	61.388
232U	24Ne	9.16×10−12	20.40	62.309
	28Mg	< 1.18×10−13	> 22.26	74.318
233U	24Ne	7.2×10−13	24.84	60.484
	25Ne			60.776
	28Mg	<1.3×10−15	> 27.59	74.224
234U	28Mg	1.38×10−13	25.14	74.108
	24Ne	9.9×10−14	25.88	58.825
	26Ne			59.465
235U	24Ne	8.06×10−12	27.42	57.361
	25Ne			57.756
	28Mg	< 1.8×10−12	> 28.09	72.162
	29Mg			72.535
236U	24Ne	< 9.2×10−12	> 25.90	55.944
	26Ne			56.753
	28Mg	2×10−13	27.58	70.560
	30Mg			72.299
236Pu	28Mg	2.7×10−14	21.52	79.668
237Np	30Mg	< 1.8×10−14	> 27.57	74.814
238Pu	32Si	1.38×10−16	25.27	91.188
	28Mg	5.62×10−17	25.70	75.910
	30Mg			76.822
240Pu	34Si	< 6×10−15	> 25.52	91.026
241Am	34Si	< 7.4×10−16	> 25.26	93.923
242Cm	34Si	1×10−16	23.15	96.508

Fine structure

The fine structure in ¹⁴C radioactivity of ²²³Ra was discussed for the first time by M. Greiner and W. Scheid in 1986.^[27] 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^[28][29] the fine structure observing transitions to excited states of the daughter. A transition with an excited state of ¹⁴C predicted in Ref. ^[27] 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^[30] 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.