Calculated fission-fragment mass yields and average total kinetic energies of heavy and superheavy nuclei

Fission-fragment mass and total-kinetic-energy (TKE) distributions following fission of even-even nuclides in the region 74≤Z≤126\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$74 \le Z \le 126$$\end{document} and 92≤N≤230\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$92 \le N \le 230$$\end{document}, comprising 896 nuclides have been calculated using the Brownian shape-motion method. The emphasis is the region of superheavy nuclei. To show compatibility with earlier results the calculations are extended to include earlier studied regions. An island of asymmetric fission is obtained in the superheavy region, 106≤Z≤114\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$106\le Z\le 114$$\end{document} and 162≤N≤176\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$162\le N\le 176$$\end{document}, where the heavy fragment is found to be close to 208\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{208}$$\end{document}Pb and the light fragment adjusts accordingly. Most experimentally observed α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}-decay chains of superheavy nuclei with Z>113\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z > 113 $$\end{document} terminate by spontaneous fission in our predicted region of asymmetric fission. In these cases, the pronounced large asymmetry is accompanied by a low TKE value compatible with measurements.


Introduction
The seven new elements in the range 107 ≤ Z ≤ 113 were all identified through α-decay chains ending in previously observed α decays. However, for the still heavier elements created in 48 Ca-induced fusion-evaporation reactions it is more involved to establish the specific isotope created since most α-decay chains end in spontaneous fission and not in a previously known α decay. For an overview, see Refs. [1,2]. To contribute to the interpretation of these experimental results we calculate fission-fragment mass and kinetic-energy distributions and average total kineticenergies (TKE) for 896 even-even nuclides in the region 74 ≤ Z ≤ 126 and 92 ≤ N ≤ 230. We use the Brownian shape-motion (BSM) method [3] which has been extena e-mail: martin.albertsson@matfys.lth.se sively compared to experimental data [4] and used in studies of the "new region of asymmetry" in the neutron-deficient Pb region [5,6].

Calculational details
The first step is to calculate the potential-energy surfaces as functions of five shape parameters. In previous publications [3,4,[6][7][8][9] potential-energy surfaces calculated in Ref. [10] have been used. The current calculations are also done as specified in Ref. [10] but with one difference: In the Strutinsky shell-correction procedure we use a larger smoothing range where B S is the ratio of the surface area of the current shape to that of a spherical shape. This is particularly important for nuclei in the vicinity of fermium (Fm) as discussed in great detail in Refs. [11,12]. In the calculation, most of the CPU time is used to calculate the single-particle levels. Once the levels are determined, the time needed to calculate the shell and pairing corrections and macroscopic contributions to the potential energy is almost negligible. Therefore we use the same set of levels to calculate the shell corrections for several neighbouring nuclei, which still leads to satisfactory accuracy. We use levels calculated for four "center" nuclei, namely 214 Rn, 288 Pu, 258 Fm, and 270 Hs. This strategy is routinely employed when calculating the potential energy for neighbouring nuclei [3,10,[13][14][15][16] and is sufficiently accurate for our overview here. As pointed out in Sect. 3.1, we obtain very similar results as Ref. [6], for the region where our two studies overlap, even though potentials for somewhat differ-ent nuclei are used to calculate the single-particle levels and the corresponding shell corrections.
The critical neck radius, where fragment separation is assumed to occur, is c 0 = 1.5 fm. To ensure that the vast majority of the random walks reach this neck radius before reaching the boundary of the employed shape lattice we have extended the grid with eight additional points in the elongation direction, corresponding to more than one million additional shapes. The largest value of the dimensionless elongation parameter q 2 (defined in Eq. (11) in Ref. [9]) is extended from roughly 16 to 20.

Fission-fragment mass distributions
We calculate the fission-fragment mass distributions with the BSM method on the multi-dimensional potential-energy surfaces as specified in Refs. [3,4]. The walks are started at the second minimum, when such a minimum exists, otherwise at the ground-state minimum. For the cases where the ground-state shape is spherical, the walk is started at the least deformed symmetric shape included in the grid, corresponding to β 2 ≈ 0.12. We follow the prescription of the earlier study [6] and use the effective level density of Ref. [4]. A more refined calculation would employ shapedependent microscopic level densities as pioneered in Ref. [9] but such a huge effort is not computationally feasible yet. Also, we limit the study to even-even nuclei because yields vary insignificantly between neighbouring nuclides except in a few isolated cases. Each distribution is based on 10,000 walks. Excitation energies are chosen just sufficiently above the barrier to obtain reasonable computing times. As in Ref. [6] the bias potential is 60 MeV. The fission-fragment mass yield Y (A) is defined as the percentage of events resulting in fragment-mass number A. The yield is normalized to 200% because each fission event results in two fragments.

Total kinetic energies
The total available energy in the fission process is given by the initial excitation energy, E exc , with respect to the ground state of the fissioning nucleus and the Q value, (2) where M(Z , N ) is the ground-state mass of the parent nucleus, and M(Z L , N L ) and M(Z H , N H ) are the groundstate masses of the light and heavy fragments, respectively. The proton and neutron numbers, Z and N , are determined by requiring the same Z /N ratio as for the fissioning nucleus. In the present study only fragments with even Z and N are considered.
The available energy is divided between the total kinetic energy and the total excitation energy of the fragments, i.e., The total excitation energy, E TXE , shared between the two fragments, is composed by the two parts, intrinsic excitation energy and deformation energy of the light (L) and heavy (H) fragment, The intrinsic excitation energy, E * sc , is the energy difference between the total energy, E tot , and the potential energy at the scission configuration U (χ sc ), where χ denotes the five shape parameters. The deformation energy is released in the fragments after scission when the accelerated fission fragments relax their respective shapes to ground-state deformations. The released energy is thereby added to the internal excitation energy. The deformation energy of each fragment is calculated as the energy difference between the fragment mass at scission [17] and its ground-state mass, where i = L or H. The masses in Eqs. (2) and (6), are calculated in the same macroscopic-microscopic model that was used to obtain the potential-energy surfaces [18]. The fragment shapes at scission are taken as the spheroidal shapes characterized by the ε 2 values in the 5D shape parametrization, while the groundstate shapes include ε 2 , ε 4 , ε 6 .
The Q * value is calculated from Eq. (2) and E TXE from Eq. (4). The TKE value is then obtained from Eq. (3). The TKE yield Y (E TKE ) is the percentage of fisson events per MeV kinetic energy, while the average TKE is the kinetic energy of the relative motion of two fission fragments averaged over all fission events.

Fission-fragment mass distributions
The calculated regions of symmetric and asymmetric fission in the full region of study are shown in Fig. 1. The results agree very well with Fig. 3 in Ref. [6] in which the lower part of the region in Fig. 1 is shown. Figure 1 also shows the upper part of the asymmetric-fissioning actinide region  Fig. 1 but studies of 70 nuclides from Z = 85 to Z = 94 were presented in Ref. [19]. It was suggested there that the transition between symmetric fission in the lighter actinide region and asymmetric for heavier actinides is A ≈ 226. It is further stated that this is somewhat surprising since one would expect both protons and neutrons to affect what regions fission symmetrically or asymmetrically. However, in Ref. [19] fission mass distributions across the line A ≈ 226 are obtained for only a few proton numbers, namely Z = 89, 90 and 91. Here, and in Ref. [6], which covers a larger, contiguous region of nuclides than the experimental work, the results show that both protons and neutrons affect asymmetry. Particularly interesting is that above Z ≈ 88 (N ≈ 132) the calculated transition line is clearly not a constant mass number A, but approximately a constant N − Z for a range of about eight proton numbers. This prediction has yet to be tested experimentally.
The calculated mass difference between the heavy and light fragments is presented in Fig. 2. Soon after the discovery Only nuclides with asymmetric fission and with a symmetric-yield to peak-yield ratio less than 0.2 are included (red squares in Fig. 1) of fission it was observed that the mass of the heavy fragment remains relatively constant at A ≈ 140, as the mass number of the fissioning system evolves from A ≈ 230 towards heavier systems. Consequently the light mass increases so the heavy/light fragment mass difference decreases as the fissioning system becomes heavier. This is particularly the case towards the neutron-deficient region where experimental data exists. Striking in Fig. 2 is the abrupt transition to a small region (around Z ≈ 110, N ≈ 166) of very large differences between the heavy and light fragment masses as well as a very low symmetric to asymmetric yield ratio seen in Fig. 1. This is illustrated in a complementary way in Fig. 3 where heavy and light fragment mass pairs are plotted as coloured symbols. Figure 3 shows that asymmetric fission of actinides with mass number from A ≈ 220 to A ≈ 246 corresponds to divisions into a 132 Sn-like heavy fragment and the corresponding partner. Similarly, there is a highly asymmetric region just below and at Z ≈ 114 corresponding to division into a 208 Pb-like heavy fragment and the corresponding partner. In the superheavy region, a 208 Pb-like fission fragment is compatible with recent results based on density functional theory [20,21], whereas recent calculations using a prescission point model predicts divisions into a 132 Sn-like light fragment and the corresponding partner [22]. It was observed in Ref. [23] that the potential-energy surfaces for some superheavy nuclei, for example 272 Ds (see Fig. 7 in Ref. [23]), exhibit the usual fission valley but also a fusion valley higher in energy than the fission valley and separated from the fission valley by a pronounced ridge. Its existence might be an additional reason that some "cold" fusion reactions have led to evaporation residue creation; the ridge between fusion valley and fission valley hinders the fusion trajectory to deflect into the fission valley before a compound nucleus is created [23]. It might be tempting to predict symmetric fission because the fission valley for this nucleus corresponds to symmetric shapes. On the other hand small neck radii are reached at much lower Q 2 values in the fusion valley than in the fission valley. In this and many other cases it is obviously hard to predict fission mass divisions from our limited ability to understand the details of the calculated multi-dimensional potential-energy surface. How- Fig. 6 Calculated symmetric-yield to peak-yield ratios versus N and Z for nuclides in the region of observed superheavy elements, for eveneven nuclides. We show some representative observed decay chains with a small solid black dot indicating the start of the decay and the larger solid white dot indicating the termination by spontaneous fission. Pairs of thin parallel lines indicate magic neutron and proton numbers Fig. 7 Calculated symmetric-yield to peak-yield ratios versus N and Z for even-even nuclides in the Fm region ever, the BSM method constitutes a well-defined approach to calculate fission mass yields. For many nuclei in this region the result is highly asymmetric fission, paradoxically corresponding to fission in the fusion valley.
Some specific, calculated, yield distributions are presented in Figs. 4 and 5 in which the highly variable and rapidly changing character of the distributions is readily apparent. The calculated mass yields in Fig. 4 show several similarities with the yields obtained using a scission-point model [24,25], though the exact transition points from asymmetric to symmetric fission in this region differ slightly. Nuclides around the superheavy region exhibit modes leading to both sym-  Fig. 5). Due to the subtle competition between these two modes, there can be abrupt changes in yields between neighbouring nuclides (e.g. (Z , N ) = (108, 156)).
There are some observations of fission kinetic energies in the superheavy region, normally in fission events that terminate α-decay chains, but it is hard to draw conclusions about the asymmetry of fission due to the very few events observed. Some observations of such fission events are reported in Ref. [26]. Sometimes a high TKE value is associated with symmetric fission, in particular for some Fm isotopes where the scission configuration is very compact. However, if the fission configuration is "liquid-drop-like" symmetric fission would correspond to elongated scission shapes and a low TKE. Therefore, a TKE value by itself is not sufficient to establish if fission is symmetric or asymmetric. We show in Fig. 6 an enlarged portion of Fig. 1 and two examples of decay chains discussed in Ref. [26] and their termination by fission at the decay-chain endpoints. Figure 7 shows an enlarged display of the symmetric-yield to peak-yield ratio in Fig. 1 to enhance details in the vicinity of 258 100 Fm 158 . The transition from asymmetric to symmetric fission at 258 100 Fm 158 is very visible in this figure as well as in Fig. 4. In this region there are also drastic variations of TKE distributions, so we now discuss the results obtained for these distributions.

Fission-fragment TKE distributions
The calculated average TKE values (see Sect. 2.2 for details about the method used) for the entire region of study are shown in Fig. 8. Obviously the TKE increases for heavier nuclides and towards the neutron drip line. More easily interpretable is the difference between the Viola TKE systematics [27] E (Viola) and the actually calculated average TKE. This is illustrated in Fig. 9. The substantial abrupt local variations seen in the region near 258 Fm are also seen experimentally [28][29][30][31][32]. In their main features, the measurements are well reproduced in the calculations. In the region of very asymmetric fission below Z = 114 (N ≈ 162-174) the TKE is lower than the systematics as would be expected. Also in the heavy neutronrich region (N ≈ 210-220), where we obtain asymmetric fission yields, the average TKE is lower than given by the Viola systematics. Figure 10 shows in detail the TKE distributions for the three isotopes 256 Fm, 258 Fm, and 260 Fm. The transition from "normal" TKE to high TKE is dramatic, with 258 Fm exhibiting a clear "bimodal" structure, as has been seen experimentally [31]. A bimodal structure in the Fm region has also been obtained within the scission-point model [24].

Summary
The BSM method [3], which in systematic calculations was previously applied to nuclides from Z = 74 to Z = 94 [6], has here been used to perform systematic yield calculations in the region 74 ≤ Z ≤ 126 and 92 ≤ N ≤ 230 for 896 even-even nuclides. Where there is overlap with the previous calculations there is good agreement. Results above Z = 94 show for the first time predictions based on the BSM method for the heavier actinide region and the superheavy region. In the vicinity of Z = 114 there is a new, smaller region of asymmetry corresponding to 208 Pb-like heavy fragments. Neutron-rich nuclei near the neutron drip line are also predicted to fission asymmetrically.