1 Introduction
In recent years, the synthesis of superheavy elements (SHEs) has been developed considerably, both experimentally and theoretically [1-3] and the Periodic Table of elements has been significantly extended. After the first successful synthesis of element 117 in 2010, all the elements in the seventh row of the Periodic Table have been discovered by mankind [3-10], which marks a great advancement in man’s cognition of the micro-world and a milestone in the evolution of the synthesis of SHEs. At present, efforts are being made to synthesize superheavy nuclei (SHN) experimentally. The elements from 107 to 113 were synthesized at GSI and RIKEN in cold fusion reactions [3-6], and the elements from 114 to 118 were produced at Dubna in hot fusion reactions using a neutron-rich projectile nucleus 48Ca [1, 7-9]. Presently, scientists are interested in the production of SHEs with Z ≥ 119, which will be a tremendous breakthrough in the field of nuclear physics. However, in the synthesis of SHEs, the cross sections decrease with the increase in the charge number of the nuclei, and are close to or even far below the order of pb, which makes the process very difficult.
In addition, the production and measurement of the isotopes of already-known SHEs, which have not yet been synthesized, is important, as they can provide information on the trends of the decay properties along the neutron axis approaching the proposed neutron magic number, including the "island of stability". The synthesis of new SHN can greatly extend the map of the nuclides. According to the theory of nuclear structure [11-13], there is a double magic SHN, the nuclei in the neighborhood of which are relatively stable. In order to accurately locate the "island of stability", more experimental data needs to be available, therefore, the synthesis of more SHN is evidently important. At present, the main method used in the synthesis of SHN is the heavy-ion fusion evaporation reaction, hence, theoretical study of the process is essential for appropriate performance of the experiments [14-27]. Although, the ability to predict the production cross sections and incident energies is limited for theoretical models, it is still beneficial to investigate all the possible combinations systematically and predict the favorable ones for producing SHN [28].
In Fig. 1, the blue crosses represent the nuclei that will be investigated in this work. Due to the transition from the cold fusion reactions to hot fusion reactions, the SHN in this region are not synthesized yet. In this study, the effect of mass asymmetry on the fusion probability and cross sections of the reactions 26Mg + 248Cm, 36S + 238U, 48Ca + 226Ra, which lead to the same compound nucleus 274HS*, are investigated first. In view of the long duration of the experiments and production rate, combinations of stable projectiles and targets with half-lives longer than 30 days are considered to be suitable. All possible combinations of projectiles from O to Ti and targets with half-lives longer than 30 days for producing SHN 264Db, 265Db, 267Sg, 268Bh, 268Sg, 269Bh, 271Hs, 271Mt, 272Hs, 272Mt, 273Mt, 274Ds, 275Ds, 275Rg, 276Ds, 276Rg, 277Rg, 278Cn, 279Cn, and 280Cn are examined and some optimal combinations and incident energies for synthesizing these nuclei are suggested.
-201811/1001-8042-29-11-003/alternativeImage/1001-8042-29-11-003-F001.jpg)
The article is organized as follows: In Sec. 2, a description of the DNS model is given; The Sec. 3 details an analysis of the results of calculation and a prediction of the optimal combinations for synthesizing several SHN in the gap region; The Sec. 4 presents the summary and the prospects.
2 Theoretical description
According to the concept of the DNS model, the complete fusion reaction is described as a diffusion process. All the nucleons of the projectile are transferred to the target and form a compound nucleus. This process is accompanied by the dissipation of energy and angular momentum. The evaporation residue (ER) cross section of the superheavy nucleus under the incident energy Ec.m. in the center-of-mass frame can be written as [21, 24, 29-34]:
where, T(Ec.m., J) is the probability that the collision system overcomes the Coulomb barrier and forms the dinuclear system [35]. PCN(Ec.m., J) is the fusion probability. Wsur(Ec.m., J) is the survival probability of the compound nucleus.
The fusion process described as a diffusion process, can be understood by numerically solving the master equations. The time evolution of the distribution probability function P(Z1, N1, t) for fragment 1 with the mass number A1 and excitation energy E1 at time t is described by the following master equations [17, 24]:
The potential energy surface (PES) of the DNS in the fusion process is defined as:
where Z = Z1 + Z2 and N = N1 + N2. EB(Zi,Ni) and EB(Z,N) are the binding energies of the fragment i and the compound nucleus, respectively. VC and VN are the Coulomb potential and nuclear potential, respectively [24].
After the capture process, the DNS is formed. For the formation of the compound nucleus (CN), the projectile should overcome the inner fusion barrier Bfus, defined as the difference between UBG, the highest point called Businaro–Gallone (B.G.) point on the left of U(≤ta), and the potential at the entrance. Hence, summing up the probability on the left side of the B.G. point, the fusion probability of the DNS is given by [17, 28, 36],
The compound nucleus formed in the fusion reaction is unstable because of its high excitation energy and usually cools down by emitting γ-rays, evaporating particles (neutrons, protons, α and other light charged particles) and fission [37-39]. This process is described by a statistical model. Subsequent to the fission of the composite nucleus, the survival probability of the SHN after evaporation of x neutrons (considering only neutron evaporation and fission), can be written as:
where J and
According to Weisskopf’s theory of evaporation, the evaporation width of the particle ν can be written as [41]:
where sν is the spin, mν is the reduced mass relative to the remaining nucleus, and Bν is the binding energy, of the evaporating particle ν.
The width of fission is given by the Bohr-Wheeler formula [43]:
Here
Here, ED=20 MeV is the shell damping energy.
3 Results and discussion
In order to investigate qualitatively the effects of mass asymmetry on the fusion probability, the driving potential, which is the minimum value of PES for combinations with the same mass asymmetry, is shown in Fig. 2, as a function of mass asymmetry for the reactions 26Mg + 248Cm, 36S + 238U, and 48Ca + 226Ra, leading to the same compound nucleus 274Hs*. Although the systems evolve in both the proton and neutron degrees of freedom, the mass asymmetry can still be a good collective degree for estimating the fusion probability. The B.G. point is shown in this figure. The blue arrow indicates the inner fusion barrier Bfus for the reaction 48Ca + 226Ra. The inner fusion barrier plays a main role in the competition between quasifission and complete fusion. According to the results of the calculation, the inner fusion barriers for the reactions 36S + 238U and 48Ca + 226Ra are Bfus, S = 3.6 MeV and Bfus, Ca = 8.6 MeV, respectively. It is observed that the injection point of the reaction 26Mg + 248Cm is located on the left side of the B.G. point, therefore, the fusion probability for 26Mg + 248Cm would be much larger than that of the other two reactions.
-201811/1001-8042-29-11-003/alternativeImage/1001-8042-29-11-003-F002.jpg)
Figure 3 presents the fusion probabilities as a function of the excitation energy of the compound nucleus in the reactions 26Mg + 248Cm, 36S + 238U, and 48Ca + 226Ra. It can be seen from the calculated results that with the same excitation energy, the fusion probability decreases with decreasing mass asymmetry. This is because a near symmetric system can overcome a high inner fusion barrier for forming a compound nucleus by nucleon transfer. The fusion probability for the reaction 26Mg + 248Cm, is close to 1 and much larger than that of the other two reactions. It can also be seen that the fusion probabilities for all the reactions increase with increasing energy.
-201811/1001-8042-29-11-003/alternativeImage/1001-8042-29-11-003-F003.jpg)
Figure 4 shows the excitation functions for the reactions 26Mg + 248Cm, 36S + 238U, and 48Ca + 226Ra. The available experimental data [44-46], which is indicated by points, is also shown. The curves show the results of calculation for different evaporation channels. It can be seen that the calculations duly authenticate the experimental data.
-201811/1001-8042-29-11-003/alternativeImage/1001-8042-29-11-003-F004.jpg)
The excitation functions for the reactions 22Ne + 254Cf, 26Mg + 250Cm, and 48Ca + 228Ra are shown in Fig. 5. All the reactions lead to the same compound nucleus, 276Hs*. It can be seen that the cross sections decrease with decreasing mass asymmetry.
-201811/1001-8042-29-11-003/alternativeImage/1001-8042-29-11-003-F005.jpg)
However, it cannot be presumed that the reaction with larger value of mass asymmetry is better for producing SHN. For systems with larger values of mass asymmetry, the reactions are usually very hot. Therefore, in the de-excitation process, many neutrons evaporate. This actually limits the neutron richness of the synthesized SHN. Moreover, due to the limitation of the target, the projectile cannot be very light. In addition, the mass asymmetry not only affects the fusion probability, but also has an impact on the probability of the capture and survival stages. Therefore, in this study, all the possible combinations are examined systematically, with stable projectiles of Z=8-22 and targets with half-lives longer than 30 days, for synthesizing the unknown SHN, as shown in Fig. 1.
The cross sections of different projectile-target combinations are compared to obtain the optimal ones for synthesizing the unknown SHN, as shown in Figs. 6, 7, 8, 9 and 10. It can be seen that for producing most of these nuclei, the maximal ER cross sections are larger than 10 pb. For example, for producing the 268Sg, the ER cross section is close to 1 nb in the reaction 18O + 254Es → 268Sg + 4n with an incident energy of 82.87 MeV. Moreover, for producing the 264Db, 265Db, 267Sg, 268Bh, 269Bh, 271Hs, 271Mt, 272Hs, 272Mt, 273Mt, 274Ds, 275Rg, 275Ds, 276Ds, 276Rg, 277Rg, the most promising reactions are 18O + 249Bk → 264Db + 3n (12 pb at 79.36 MeV), 19F + 250Cm → 265Db + 4n (5.6 pb at 89.55 MeV), 17O + 254Cf → 267Sg + 4n (82 pb at 82.01 MeV), 18O + 254Es → 268Bh + 4n (160 pb at 84.22 MeV), 18O + 255Es → 269Bh + 4n (770 pb at 84 MeV), 18O + 257Fm → 271Hs + 4n (520 pb at 84.56 MeV), 22Ne + 252Es → 271Mt + 3n (190 pb at 99.25 MeV), 22Ne + 254Cf → 272Hs + 4n (190 pb at 101.87 MeV), 32Si + 243Am → 272Mt + 3n (51 pb at 134.9 MeV), 22Ne + 255Es → 273Mt + 4n (140 pb at 103.52 MeV), 21Ne + 257Fm → 274Ds + 4n (23 pb at 104.13 MeV), 30Si + 248Bk → 275Rg + 3n (20 pb at 137.64 MeV), 22Ne + 257Fm → 275Ds + 4n (20 pb at 104.92 MeV), 22Ne + 257Fm → 276Ds + 3n (9.2 pb at 99.92 MeV), 32Si + 247Bk → 276Rg + 3n (10 pb at 139.94 MeV), and 32Si + 248Bk → 277Rg + 3n (18 pb at 139.31 MeV), respectively. However, it can be seen that for synthesis of 278-280Cn the ER cross sections are approximately 1 pb, therefore, the reaction is still feasible with the available experimental equipment. The maximal evaporation cross sections, optimal incident energies, and corresponding evaporation channels are deduced and presented clearly in Table. 1. It is found that for synthesizing the isotopes 267Sg and 268Sg, the target 254Cf can be used. The 18O is predicted as a projectile, for producing both the isotopes 268Bh and 269Bh. In addition, the combinations of projectile and target are the same for synthesizing the 275Ds and 276Ds, but in different evaporation channels.
*ON | 2*BC | 2*EC | ECN*(Ec.m.) | 2*σER, max(pb) |
---|---|---|---|---|
264Db | 18O + 249Bk | 3n | 36 (79.36) | 1.2 × 101 |
265Db | 19F + 250Cm | 4n | 44 (89.55) | 5.6 |
267Sg | 17O + 254Cf | 4n | 40 (82.01) | 8.2 × 101 |
268Sg | 18O + 254Cf | 4n | 39 (82.87) | 9.2 × 102 |
268Bh | 18O + 254Es | 4n | 39 (84.22) | 1.6 × 102 |
269Bh | 18O + 255Es | 4n | 39 (84.00) | 7.7 × 102 |
271Hs | 18O + 257Fm | 4n | 38 (84.56) | 5.2 × 102 |
271Mt | 22Ne + 252Es | 3n | 34 (99.25) | 1.9 × 102 |
272Hs | 22Ne + 254Cf | 4n | 39 (101.87) | 1.9 × 102 |
272Mt | 32Si + 243Am | 3n | 32 (134.9) | 5.1 × 101 |
273Mt | 22Ne + 255Es | 4n | 39 (103.52) | 1.4 × 102 |
274Ds | 21Ne + 257Fm | 4n | 43 (104.13) | 2.3 × 101 |
275Rg | 30Si + 248Bk | 3n | 33 (137.64) | 2.0 × 101 |
275Ds | 22Ne + 257Fm | 4n | 39 (104.92) | 2.0 × 101 |
276Ds | 22Ne + 257Fm | 3n | 34 (99.92) | 9.2 |
276Rg | 32Si + 247Bk | 3n | 32 (139.94) | 1.0 × 101 |
277Rg | 32Si + 248Bk | 3n | 32 (139.31) | 1.8 × 101 |
278Cn | 30Si + 251Cf | 3n | 35 (140.57) | 1.7 |
279Cn | 30Si + 252Cf | 3n | 36 (140.22) | 1.1 |
280Cn | 30Si + 254Cf | 4n | 43 (144.86) | 3.0 |
-201811/1001-8042-29-11-003/alternativeImage/1001-8042-29-11-003-F006.jpg)
-201811/1001-8042-29-11-003/alternativeImage/1001-8042-29-11-003-F007.jpg)
-201811/1001-8042-29-11-003/alternativeImage/1001-8042-29-11-003-F008.jpg)
-201811/1001-8042-29-11-003/alternativeImage/1001-8042-29-11-003-F009.jpg)
-201811/1001-8042-29-11-003/alternativeImage/1001-8042-29-11-003-F010.jpg)
4 Summary
The production of several SHN was investigated systematically within the dinuclear system model. The mass asymmetry effects were studied and it was observed that the fusion probabilities decreased with decreasing mass asymmetry. Suitable reactions for producing SHN in the gap region, as shown in Fig. 1, were predicted by reviewing the stable beam induced hot fusion reactions, and 192 possible combinations with projectiles from O to Ti and targets with half-lives longer than 30 days, for producing SHN 264Db, 265Db, 267Sg, 268Bh, 268Sg, 269Bh, 271Hs, 271Mt, 272Hs, 272Mt, 273Mt, 274Ds, 275Ds, 275Rg, 276Ds, 276Rg, 277Rg, 278Cn, 279Cn, and 280Cn were examined. The optimal combinations and incident energies for synthesizing these nuclei were predicted. It was found that the production cross sections, for synthesizing most of these SHN, were larger than 10 pb. The predicted cross section was 920 pb for the production of 268Sg; hence, the reaction could be carried out with the available experimental equipment. In future, it would be necessary to carry out systematic calculations on the production of SHN within other theoretical models so that a comparison of the predicted favorable reactions among the different models can provide valuable information regarding the experimental procedure.
Superheavy nuclei from 48Ca-induced reactions
. Nucl. Phys. A 944, 62 (2015). doi: 10.1016/j.nuclphysa.2015.07.003Cross sections for the production of superheavy nuclei
. Nucl. Phys. A 944, 257 (2015). doi: 10.1016/j.nuclphysa.2015.02.010The discovery of the heaviest elements
. Rev. Mod. Phys. 72, 733 (2000). doi: 10.1103/revmodphys.72.733The reaction 48Ca + 238U → 286112* studied at the GSI-SHIP
. Eur. Phys. J. A 32, 251 (2007), ISSN 1434-601X. doi: 10.1140/epja/i2007-10373-xNew Result in the Production and Decay of an Isotope, 278113, of the 113th Element
. Journal of the Physical Society of Japan 81, 103201 (2012). doi: 10.1143/jpsj.81.103201Review of even element super-heavy nuclei and search for element 120
. Eur. Phys. J. A 52 (2016). doi: 10.1140/epja/i2016-16180-4Synthesis of nuclei of the superheavy element 114 in reactions induced by 48Ca
. Nature 400, 242 (1999 a). doi: 10.1038/22281Synthesis of Superheavy Nuclei in the 48Ca + 244Pu Reaction
. Phys. Rev. Lett. 83, 3154 (1999 b). doi: 10.1103/physrevlett.83.3154Synthesis of a New Element with Atomic Number Z=117
. Phys. Rev. Lett. 104, 142502 (2010). doi: 10.1103/PhysRevLett.104.142502Nuclei in the "Island of Stability" of Superheavy Elements
. J. Phys.: Conf. Ser. 337, 012005 (2012). doi: 10.1088/1742-6596/337/1/012005On the stability of superheavy nuclei against fission
. Z. Phys. A 222, 261 (1969). doi: 10.1007/bf01392125On the spontaneous fission of nuclei with Z near 114 and N near 184. and P. Möller
, Nucl. Phys. A 115, 545 (1968). doi: 10.1016/0375-9474(68)90748-3Stability of heavy and superheavy elements
. J. Phys. G: Nucl. Part. Phys. 20, 1681 (1994). doi: 10.1088/0954-3899/20/11/003Theoretical study of the synthesis of superheavy nuclei with Z=119 and 120 in heavy-ion reactions with trans-uranium targets
. Phys. Rev. C 85, 041601 (2012). doi: 10.1103/PhysRevC.85.041601Influence of proton shell closure on the evaporation residue cross sections of superheavy nuclei
. J. Phys. G: Nucl. Part. Phys. 44, 045105 (2017 a). doi: 10.1088/1361-6471/aa53e8Possibilities for synthesis of new isotopes of superheavy nuclei in cold fusion reactions
. Phys. Rev. C 93, 044615 (2016 b). doi: 10.1103/PhysRevC.93.044615Production of new superheavy Z=108-114 nuclei with 238U, 244Pu, and 248,250Cm targets
. Phys. Rev. C 80, 057601 (2009 a). doi: 10.1103/PhysRevC.80.057601Influence of entrance channel on production cross sections of superheavy nuclei
. Phys. Rev. C 96, 024610 (2017 b). doi: 10.1103/PhysRevC.96.024610Competing fusion and quasifission reaction mechanisms in the production of superheavy nuclei
. Phys. Rev. C 82, 044614 (2010). doi: 10.1103/PhysRevC.82.044614Neutron emission in the fusion process and its effect on the formation of superheavy nuclei
. Phys. Rev. C 89, 024604 (2014). doi: 10.1103/PhysRevC.89.024604Q-value effects on the production of superheavy nuclei
. Phys. Rev. C 74, 057602 (2006). doi: 10.1103/PhysRevC.74.057602ISOSPIN DEPENDENCE OF REACTIONS 48Ca + 243-251Bk
. Int. J. Mod. Phys. E 17, 66 (2008). doi: 10.1142/s0218301308011768Systematics of fusion probability in “hot” fusion reactions
. Phys. Rev. C 84, 061601 (2011). doi: 10.1103/PhysRevC.84.061601Orientation effects on evaporation residue cross sections in 48Ca-induced hot fusion reactions
. Phys. Rev. C 90, 014612 (2014 a). doi: 10.1103/PhysRevC.90.014612Formation of superheavy nuclei in cold fusion reactions
. Phys. Rev. C 76, 044606 (2007). doi: 10.1103/PhysRevC.76.044606Production cross sections of superheavy elements Z=119 and 120 in hot fusion reactions
. Phys. Rev. C 89, 024615 (2014 b). doi: 10.1103/physrevc.89.024615Influence of the neutron numbers of projectile and target on the evaporation residue cross sections in hot fusion reactions
. Phys. Rev. C 93, 064610 (2016). doi: 10.1103/physrevc.93.064610Synthesis of neutron-rich superheavy nuclei with radioactive beams within the dinuclear system model
. Phys. Rev. C 97, 064609 (2018). doi: 10.1103/PhysRevC.97.064609Competition between complete fusion and quasi-fission in reactions between massive nuclei
. The fusion barrier. Phys. Lett. B 319, 425 (1993). doi: 10.1016/0370-2693(93)91746-aProblems in description of fusion of heavy nuclei in the two-center shell model approach
. Nucl. Phys. A 646, 29 (1999). doi: 10.1016/s0375-9474(98)00616-2Effect of the entrance channel on the synthesis of superheavy elements
. Eur. Phys. J. A 8, 205 (2000). doi: 10.1007/s10050-000-4509-7Competition between evaporation channels in neutron-deficient nuclei
. Phys. Rev. C 68, 014616 (2003). doi: 10.1103/PhysRevC.68.014616Possibilities of production of neutron-deficient isotopes of U, Np, Pu, Am, Cm, and Cf in complete fusion reactions
. Phys. Rev. C 78, 044605 (2008). doi: 10.1103/PhysRevC.78.044605Possibilities of production of transfermium nuclei in charged-particle evaporation channels
. Phys. Rev. C 94, 044606 (2016). doi: 10.1103/PhysRevC.94.044606Production cross sections of superheavy nuclei based on dinuclear system model
. Nucl. Phys. A 771, 50 (2006). doi: 10.1016/j.nuclphysa.2006.03.002Influence of nuclear basic data on the calculation of production cross sections of superheavy nuclei
. Phys. Rev. C 92, 014601 (2015 b). doi: 10.1103/PhysRevC.92.014601The conditions for the synthesis of heavy nuclei
. Rep. Prog. Phys. 54, 949 (1991). doi: 10.1088/0034-4885/54/7/002Phenomenological statistical analysis of level densities, decay widths and lifetimes of excited nuclei
. Nucl. Phys. A 543, 517 (1992). doi: 10.1016/0375-9474(92)90278-rSurvival probability of superheavy nuclei
. Phys. Rev. C 65, 024308 (2002). doi: 10.1103/PhysRevC.65.024308Nuclear Ground-State Masses and Deformations
. At. Data Nucl. Data Tables 59, 185 (1995). doi: 10.1006/adnd.1995.1002Statistics and Nuclear Reactions
. Phys. Rev. 52, 295 (1937). doi: 10.1103/physrev.52.295Decay of deformed and superdeformed nuclei formed in heavy ion reactions
. Phys. Rev. C 21, 1770 (1980). doi: 10.1103/physrevc.21.1770The Mechanism of Nuclear Fission
. Phys. Rev. 56, 426 (1939). doi: 10.1103/physrev.56.426Experimental study of the 238U(36S,3-5n)269-271Hs reaction leading to the observation of 270Hs
. Phys. Rev. C 81, 061601 (2010). doi: 10.1103/PhysRevC.81.061601Synthesis and study of decay properties of the doubly magic nucleus 270Hs in the 226Ra + 48Ca reaction
. Phys. Rev. C 87, 034605 (2013 b). doi: 10.1103/PhysRevC.87.034605Observation of the 3n Evaporation Channel in the Complete Hot-Fusion Reaction 26Mg +248Cm Leading to the New Superheavy Nuclide 271Hs
. Phys. Rev. Lett. 100, 132503 (2008). doi: 10.1103/PhysRevLett.100.132503