-
The production of superheavy nuclei (SHN) in the laboratory has drawn considerable attention, because it is related to the existence of an island stability in SHN [1-5] and depends on its location. Currently, there is a gap between the SHN synthesized by hot fusion and cold fusion reactions [6-11]. To fully understand the shell effects and other properties, and to develop theoretical models that will be able to predict the properties of SHN located beyond this area, expanding the span of known isotopes of SHN is extremely important.
In recent years, many efforts were devoted to the investigation of the synthesis mechanism of SHN [12–26], and many approaches were proposed to calculate the fusion probability [27–34]. To date, however, none of these approaches are predominant [35]. The dinuclear system (DNS) concept represents one of these approaches [16, 17,36-39]. Based on the DNS model, calculated evaporation residue cross-sections (ERCS) for the cold and hot fusion reactions leading to heavy and SHN are in good agreement with available experimental data [40–60].
Selecting the optimal composition of the colliding nuclei is one of the most important factors for the successful synthesis of SHN. In order to fill the gap, several isotopes of actinide nuclei can be used as targets. The use of beams of 36S, 40Ar, and 44Ca as the projectiles is also interesting. Recently, the first measurement of ERCSs in the complete fusion reaction 36S + 238U and the observation of 270Hs has already been performed in experimental studies [61].
One of the aims of the present study is to investigate several fusion reactions leading to the formation of the same compound nucleus of unknown isotopes of SHN, and those between already obtained SHN in cold and hot fusion. Usually, fusion hindrance is smaller for the mass-asymmetric reaction system, and hence the corresponding cross-section may be enhanced. Once the best projectile is selected, to find the optimal conditions of synthesis, it is necessary to study the dependence of the ERCS on the isospin composition of colliding nuclei. In the present study, the influence of the target neutron number on the capture cross section, fusion probability, and survival probability for the reactions 36S + 236−244Pu is investigated in detail. In addition, we systematically study the ERCSs of the 36S bombarding targets of the actinide isotopic chain. The aim of our study is to predict the ERCSs of unknown neutron-deficient isotopes of SHN.
-
The evaporation residue cross-section (ERCS) in heavy-ion fusion reactions is calculated as the summation over all partial waves
J [16,17],σER(Ec.m.)=πℏ22μEc.m.Jmax∑J=0(2J+1)T(Ec.m.,J)×PCN(Ec.m.,J)Wsur(Ec.m.,J),
(1) where
Ec.m. is the incident energy in the center-of-mass frame.T(Ec.m.,J) is the transmission probability of the two colliding nuclei overcoming the Coulomb potential barrier in the entrance channel to form the DNS. The capture cross-sectionσcap=πℏ22μEc.m.∑J(2J+1)T(Ec.m.,J) is calculated with an empirical coupled-channel approach [23, 62]. ThePCN is the probability that the system evolves from a touching configuration to the compound nucleus in competition with the quasifission process. The last termWsur is the survival probability of the formed compound nucleus, which can be estimated with a statistic method [63]. -
The capture cross-section is:
σcap(Ec.m.)=πℏ22μEc.m.∑J(2J+1)T(Ec.m.,J),
(2) where the transmission probability can be written as
T(Ec.m.,J)=∫f(B)11+exp{−2πℏω(J)[Ec.m.−B−ℏ22μR2BJ(J+1)]}dB,
(3) where
ℏω(J) is the width of the parabolic Coulomb barrier at the positionRB(J) , and an empirical coupled channel method is used via a barrier distribution function, which is taken as an asymmetric Gaussian form [55]. The nucleus-nucleus interaction potential with quadrupole deformation was used, which is addressed in detail in Ref. [52]. -
The fusion dynamics are described as a diffusion process by numerically solving a two-variable master equation (ME) in the corresponding potential energy surfaces [23, 53]. The temporal evolution of the probability distribution function
P(Z1,N1,ε1,t) for fragment 1 withZ1 andN1 with the local excitation energyε1 at timet is described by the following ME:dP(Z1,N1,ε1,t)dt=∑Z′1WZ1,N1;Z′1,N1(t)×[dZ1,N1P(Z′1,N1,ε′1,t)−dZ′1,N1P(Z1,N1,ε1,t)]+∑N′1WZ1,N1;Z1,N′1(t)×[dZ1,N1P(Z1,N′1,ε′1,t)−dZ1,N′1P(Z1,N1,ε1,t)]−[Λqf(Θ(t))+Λfs(Θ(t))]P(Z1,N1,ε1,t),
(4) where
WZ1,N1;Z1,N′1 is the mean transition probability from channel (Z1,N′1 ) to (Z1,N1 ), anddZ1,N1 denotes microscopic dimensions corresponding to the macroscopic state (Z1,N1 ) [18, 27, 64, 65], which is shown later.ε1 denotes the local excitation energy and is likewise shown later. The sum is taken over all possible proton and neutron numbers that fragmentZ′1 ,N′1 may assume, but only one nucleon transfer is considered in the model.The probability
P(Z1,N1,ε1,t) distributed in the bottom of pocket has the possibility to decay out of the DNS, i.e., the evolution of the DNS along the variableR leads to the quasifission of the DNS, with the corresponding quasifission rateΛqfZ,N(Θ) . In order to consider the influence of the DNS decay on the probability distributionP(Z1,N1,ε1,t) , we need to include the effect of the quasifission rate on the ME. The quasifission rateΛqf in Eq. (10) is estimated with the one-dimensional Kramers formula [66, 67].ΛqfZ,N(Θ)=ω2πωBqf(√(Γ2ℏ2)2+(ωBqf)2−Γ2ℏ)×exp(−Bqf(Z,N)Θ(Z,N)).
(5) The quasifission barrier
Bqf measures the depth of the pocket of the nucleus-nucleus interaction potential. The quasifission barrierBqf decreases with increasingZ , and for near symmetric configurations there is no minimum of the nucleus-nucleus potential [46, 48]. In the present study, in the case where the nucleus-nucleus interaction potential has no minimum, we assumed that the height of the quasifission barrierBqf is 0.5 MeV. The temperatureΘ(Z,N) of the DNS is calculated with the expressionΘ(Z,N)=√ε/a whereε depicts the local excitation energy of the DNS. The level density parameter is calculated with the expressiona=A/12 MeV−1 . Here,ω is the frequency of the harmonic oscillator approximating the potential along the internuclear distance around the bottom of the pocket. The frequencyωBqf is the frequency of the inverted harmonic oscillator approximating the interaction potential of two nuclei along the internuclear distance around the top of the quasifission barrier. The quantityΓ denotes the double average width of the contributing single-particle states. In the present study, constant valuesΓ=2.8 MeV,ℏωBqf=2.0 MeV, andℏω=3.0 MeV were employed.Solving Eq. (4) numerically, the temporal evolution of the probability distribution
P(Z1,N1,ε1,t) to find fragment 1 (Z1+N1 ) with excitation energyε1 at timet is obtained. All the components on the left side of the Businaro-Gallone (BG) point contribute to the compound nuclear formation. The fusion probability represents the Z-N configuration at the BG point, beyond which the system falls into the fusion valley in the potential energy surface as a function of the mass-charge asymmetry parameter. Therefore, the fusion probabilityPCN is the summation ofP(Z1,N1,ε1,t) from (Z1=1,N1=1 ) to (ZBG,NBG ) configurations. The compound nucleus formation probability at the minimum of the nucleus-nucleus potential (Bm ), which corresponds to a certain orientation of the colliding nuclei in the entrance channel, and for the angular momentumJ , is given byPCN(Ec.m.,J,Bm)=ZBG∑Z1=1NBG∑N1=1P(Z1,N1,ε1,τint,Bm).
(6) The interaction time
τint (this will be shown later) in the dissipative process of two colliding nuclei is dependent on the incident energyEc.m. ,J andBm , and it is determined by using the deflection function method [68]. Finally, we obtain the fusion probabilityPCN(Ec.m.,J) asPCN(Ec.m.,J)=∫f(Bm)PCN(Ec.m.,J,Bm)dBm,
(7) where the barrier distribution function is provided in an asymmetric Gaussian form [18].
In order to numerically solve Eq. (4), interaction time and local excitation energy are needed as input. The time interval between formation and break of the composite system is defined as the interaction time
τint . As shown in Fig. 1 in Ref. [68], during this process the composite system rotates about its center of mass. On the one hand, for a given valueJi of the incident angular momentum,τint(Ji) is determined by the rotation of the composite system through the angleFigure 1. Mean interaction times are shown as a function of incident angular momentum
J for 36S + 250Cf, 40Ar + 246Cm, and 44Ca + 242Pu reactions with corresponding excitation energyE∗CN=40 MeV.Δϑ(Ji)=π−ϑi−ϑf−Θ(Ji),
(8) where the Coulomb angles
ϑi andϑf are given by Coulomb trajectories in the entrance and exit channels with the corresponding the energiesEi ,Ef and the angular momentaJi ,Jf values, respectively.ϑi(f)=arcsin2bi(f)/R+εi(f)√4+ε2i(f)−arcsin1√(2/εi(f))2+1,
(9) where
εi(f)=α/(Ec.m.bi(f)) ,α=ZPZTe2 , andbi(f)=ℏJi(f)/ √2μEc.m. .The essential ingredient of the model is the determination of the deflection function
Θ(Ji) from the experimental angular distribution. However, this is achieved by introducing the parametrization [69]Θ(Ji)=ΘC(Ji)−βΘCgrJiJgr(δβ)Ji/Jgr.
(10) The first term on the right-hand side is the Coulomb deflection function. The second term describes the deviation from the Coulomb deflection function due to nuclear interaction between the projectile and target. The parameters
δ andβ are determined by a fit of the differential cross-section obtained from experimental data. The initial angular momentum is assumed asJi=J . The details ofδ andβ are given in Ref. [70]. The grazing angular momentumJgr can be expressed:Jgr=0.22Rint[Ared(Ec.m.−V(Rint))]1/2,
(11) where
V(Rint) denotes the interaction barrier at the interaction radiusRint .Ared is the reduced mass.During this process, the composite system rotates about its center of mass. The relation between
Δϑ andτint is given by the integralΔϑ(Ji)=∫τint0dtdϑdt=∫τint0dtℏJ(t)ζrel(t)
(12) with the time-dependent angular momentum
J(t) and relative moment of inertiaζrel(t) . The dissipation of the relative angular momentum<J(t)> is described by<J(t)>=Jst+(Ji−Jst)exp(−t/τJ),
(13) where the limiting value
Jst given by the sticking condition isJst=Jiζ0rel/ζ0tot . The relaxation timeτJ is1.5×10−21 s. For the relative and total moments of inertia, we assume the rigid-body values:ζrel=μR2 (ζ0rel=μR20 ) andζtot=μR2+ 25m1R21+25m2R22 (ζ0tot=μR20+25m1R21+25m2R22) ; wherem1 ,m2 ,μ ,R1 ,R2 are the masses, the reduced mass, and radii of the fragments, respectively.R0 is the radius of the rotating composite system two nuclei form at close contact. The coupled Eqs. (8)–(13) are solved by iteration to obtain the interaction timeτint .For the subsequent three 36S + 250Cf, 40Ar + 246Cm, and 44Ca + 242Pu reactions, the average interaction times are calculated by the deflection function method [68-70]. In Fig. 1, we plot the mean interaction time as a function of the incident angular momentum
J with the corresponding excitation energiesE∗CN=40 MeV. Fig. 1 shows that the interaction time of the composite system is long for partial waves with a small incident angular momentumJ . Moreover, we found that the interaction time decreases with increasingJ . In Fig. 1, we also observe the decrease of interaction time with decreasing mass asymmetry in the entrance channel with the special excitation energy and angular momentumJ . This is because the Coulomb repulsion increases gradually with the decrease in mass asymmetry.The local excitation energy is defined as [18, 71]
ε=Ex−[U(Z1,N1,Z2,N2,β1,β2,J)−U(ZP,NP,ZT,NT,βP,βT,J)],
(14) where the dissipation energy
Ex of the composite system is converted from the relative kinetic energy loss. The dissipation energyEx is related to the minimum of the nucleus-nucleus potential (Bm ) and is determined for each initial relative angular momentumJ by the parametrization method of the classical deflection function.Ex=Ec.m.−Bm−<J(t)>(<J(t)>+1)ℏ22ζrel−<Erad(J,t)>,
(15) <Erad(J,t)>=Eiradexp[−τJintτrad].
(16) τrad denotes the relaxation time of the dissipation of the radial kinetic energy. The quantityEirad denotes the initial radial kinetic energy at the interaction radius. The radial energy at the initial state isEirad(J,0)=Ec.m.−Bm−EJi(Ji+ 1)ℏ2/(2ζrel) . The initial angular momentum is assumed to beJi=J . The value ofτrad is3×10−22 s [69]. Fig. 1 shows that the angular momentum is within the range of our research, and that the interaction timeτJint is much larger thanτrad . Therefore, for the current three reaction systems, the<Erad(J,t)> value atE∗CN=40 MeV is infinitely close to zero.The second term of Eq. (14) is the driving potential energy [16-18] of the system for the nucleon transfer of the DNS, which is:
U(Z1,N1,β1,β2,J)=B(Z1,N1,β1)+B(Z2,N2,β2)−B(Z,N,β)+UC(Z1,Z2,β1,β2)+UN(Z1,N1,Z2,N2,β1,β2,J)=Qgg+UC(Z1,Z2,β1,β2)+UN(Z1,N1,Z2,N2,β1,β2,J),
(17) where
Z=Z1+Z2 andN=N1+N2 , andβi(i=1,2) andβ represent quadrupole deformations of the two fragments and the compound nucleus, respectively. TheB(Z1,N1,β1) ,B(Z2,N2,β2) , andB(Z,N,β) are the binding energies of two deformed nuclei and the compound nucleus [72], respectively. TheQgg (Qgg=B(Z1,N1,β1)+B(Z2,N2,β2)−B(Z,N,β) ) denotes the ground stateQ value. In the present study, the deformation parameters and binding energies are taken from Refs. [73, 74]. Wong's formula [75] is adopted to calculate the Coulomb interaction, and the nuclear potential is calculated with Skyrme-type interaction without considering the momentum and spin dependence [76]. Here, the inner fusion barrier appears on the driving potential energy surface during the evolution of the mass (charge) asymmetry axis. The inner fusion barrier is determined by the difference between the maximum value of the driving potential and its value at the point corresponding to the initial charge asymmetry of the considered reaction. To form a compound nucleus, the inner fusion barrier must be overcome.In Eq. (4),
WZ1,N1;Z1,N′1 ,dZ1,N1 ,Λqf , andΛfs are all dependent on the local excitation energy of the DNS. The transition probability is related to the local excitation energy, and the neutron transition probabilityWZ1,N1,β1,β2;Z1,N′1,β1,β2 can be written as [64, 65]WZ1,N1;Z1,N′1(t)=τmem(Z1,N1,ε1;Z1,N′1,ε′1)ℏ2dZ1,N1dZ1,N′1×∑ii′|<Z1,N′1,ε′1,i′|V(t)|Z1,N1,ε1,i>|2,
(18) where
i denotes all remaining quantum numbers. The memory timeτmem τmem(Z1,N1,ε1;Z1,N′1,ε′1)=ℏ√2π{<V2(t)>Z1,N1,ε1+<V2(t)>Z1,N′1,ε′1}−1/2,
(19) can be interpreted as the coherence time for the transitions between the subsets (
Z1,N1,ε1 ) and (Z1,N′1,ε′1 ) [64, 65], where<V2(t)>Z1,N1,ε1 and<V2(t)>Z1,N′1,ε′1 stand for the average expectation value withZ1,N1,ε1 andZ1,N′1,ε′1 being fixed, respectively. Thus, the memory timeτmem depends on the neutron numberN1 , proton numberZ1 , and the local excitation energyε1 . For the 36S + 250Cf reaction, the memory timeτmem(Z=16,N=20;Z=16, N = 21) is0.75×10−22 s when the excitation energyE∗CN = 40 MeV.The transition probability of Eq.(18) can be written as
WZ1,N1;Z1,N′1(t)=τmem(Z1,N1,ε1;Z1,N′1,ε′1)ℏ2dZ1,N1dZ1,N′1×{[ω11(Z1,N1,ε1;ε′1)+ω22(Z1,N1,ε1;ε′1)]δN′1,N1+ω12(Z1,N1,ε1;ε′1)δN′1,N1−1+ω12(Z1,N1,ε1;ε′1)δN′1,N1+1},
(20) where
ωkk′(Z1,N1,ε1;ε′1)=∑k,k′N′1|<Z1,N′1,ε′1,i′|Vk,k′|Z1,N1,ε1,i>|2=dZ1,N1<Vk,k′V+k,k′>.
(21) The averages in Eqs. (18), (19), and (21) are carried out by using the method of spectral distributions [77, 78]. We obtain
<Vk,k′V+k,k′>=14U2kk′gkgk′Δkk′ΔεkΔεk′×[Δ2kk′+16(Δε2k+Δε2k′)].
(22) which contains some fixed independent parameters
Ukk′(t) andΔkk′(t) . In the present work, the strength parametersUkk′(t) are assumed as [64]Ukk′(t)=g1/31⋅g1/32g1/31+g1/32⋅1g1/3k⋅g1/3k′⋅2γkk′.
(23) In our calculation
Δ11(t) =Δ12(t) =Δ22(t) =Δ21(t) = 2, and the dimensionless strength parametersγ11 =γ12 =γ22 =γ21 = 3 are assumed. Owing to the excitation, a valence spaceΔεk forms symmetrically around the Fermi surface. Only the particles in the states within this valence space are actively involved in the excitation and transfer [64, 65].Δεk=√4εkgk,εk=εAkA,gk=Ak12(k=1,2).
(24) Here
ε deontes the local excitation energy of the DNS. The microscopic dimension is [64, 65]dZ1,N1(m1,m2)=(N1m1)(N2m2).
(25) There are
Nk=gkΔεk valence states andmk=Nk/2 valence nucleons inΔεk . -
The survival probability of the compound nucleus at excitation energies is a probability for the compound system to resist fission decay in the form of emission of light particles and
γ -decay. For the sake of simplicity, the present work as well as other Refs. [24-26] consistently neglect theγ -decay width and other charged particles at high excitation energies of interest in hot fusion reactions, compared with the evaporation of successive emission neutrons. The survival probability of the excited compound nucleus in the de-excitation process, by means of the neutron evaporation in competition with fission, is expressed as follows:Wsur(E∗CN,x,J)=F(E∗CN,x,J)x∏i=1[Γn(E∗i,J)Γn(E∗i,J)+Γf(E∗i,J)]i,
(26) where
F(E∗CN,x,J) is the realization probability of thexn channel at the excitation energyE∗CN(Ec.m.+Q) of the compound nucleus with angular momentumJ ,i the index of evaporation step,Γn andΓf are the partial widths of neutron emission and fission.The partial width for emission of a neutron from a compound nucleus with the excitation energy
E0 is given by the Weisskopf formulaΓn=gmnσinvπ2ℏ2ρ0(E0−δ0)∫E0−Bn−δn0ρn(E0−Bn−δn−ε)εdε,
(27) where
mn andg are the mass and spin degeneracy of the emitted neutron, respectively;σinv is the cross section for the formation of the decaying nucleus in the inverse process;ρ0(E0−δ0) is the level density of the parent nucleus at the thermal excitation energy corrected for its pairing energyδ0 , andρn(E0−Bn−δn−ε) is the corresponding level density of the daughter nucleus after emitting a neutron.Bn andδn are the neutron separation energy and the pairing energy of the daughter nucleus, respectively.The fission width can be expressed in terms of the transition state theory as
ΓBWf=12πρ0(E0−δ0)∫E0−Bf−δf0ρn(E0−Bf−δf−ε)dε,
(28) where
ρn(E0−Bf−δf−ε) is the level density of the fissile nucleus at the saddle configuration. The calculations of the width of the fission channel are performed, taking into account the effects of nuclear viscosity and the fission delay time,Γf=ℏωgsTωsd[√1+(β2ωsd)2−β2ωsd]×ΓBWf,
(29) where the curvatures of the potential at the ground-state (
ωgs ) and saddle point (ωsd ), and the reduced friction parameterβ have been fixed with the default values ofℏωgs=2.0 MeV,ℏωsd=2.4 MeV andℏβ=3.0 MeV, respectively.The back-shift Fermi-gas model at energies of the hot-fusion reaction of interest is used to determine the level density,
ρ(U,J)=(2J+1)exp[2√aU−J(J+1)2σ2]24√2σ3a1/4U5/4,
(30) with
σ2=Θrigidℏ2√Ua ,Θrigid=25muAR2 ,U=E−δ . The back shiftsδ=−Δ (odd-odd), 0 (odd A) andΔ (even-even), respectively, are related to the neutron and proton paring gapΔ=1/2[Δn(Z,N)+Δp(Z,N)] , which is employed from mass differences of the neighboring nuclei [79]. The dependence of the level density parametera on the shell correction and the excitation energy was initially proposed asa(U,Z,N)=˜a(A)[1+Eshf(U)U]
(31) with
˜a(A)=αA+βA2/3 andf(U)=1−exp(−γDU) . It is worth noting that the differences between the corresponding level density parameters are mainly related to different shell corrections, and thus one should use these parameters at the same shell correction energies. In the present study, parametersα=0.1337 ,β=−0.06571 , andγD=0.04884 [79] are determined by fitting to experimental level density data with the help of the microscopic shell correction from FRDM95 [80], which is adopted to calculate the level density using in the evaporation calculations.We calculated the angular momentum dependence of the transmission, fusion, and survival probabilities as shown in Fig. 2 for the reaction 36S + 250Cf at incident energies 169.64 MeV. The values of the three stages decrease significantly with increasing relative angular momentum. Hence, in the following estimation of the ERCSs, we cut off the maximal angular momentum at
Jmax = 30. A similar result is also illustrated in Ref. [81]. -
Very recently, in order to review our calculated abilities on the ERCSs to synthesize superheavy nuclei using DNS model, the hot fusion reactions producing SHN with
Z⩾104 are systematically studied [82]. The results of systematic calculation show that the current theoretical method can describe the ERCS of SHN. In the present work, the calculations for all reactions were performed with the same parameters and assumptions. -
In order to predict the most suitable projectile-target combination among the probable candidates, production cross-sections of new neutron-deficient SHN with charged numbers
Z = 108−114 are analyzed systematically with 36S, 40Ar, and 44Ca projectiles. In the present study, based on the framework of the DNS model, we calculated the ERCSs of the SHN based on the actinide targets 249−252Cf, 245−248Cm, and 241−244Pu with the projectiles 36S, 40Ar, and 44Ca, as shown in Fig. 3. The ERCSs decreases by about one order of magnitude with increasing charge number of projectile fromZ=16 toZ=20 . This is due to the strong decrease in fusion probabilityPCN and the increase in the quasifission with increasing asymmetry in the entrance channel.Figure 3. (color online) Evaporation residue excitation functions in production of isotopes of superheavy nuclei Fl in reactions 36S + 249−252Cf, 40Ar + 245−248Cm, and 44Ca + 241−244Pu.
In Fig. 3(c), the calculated ERCSs are shown for the production of new neutron deficient isotopes of Fl in the fusion reactions of 36S with 251Cf targets and for the 40Ar + 247Cm fusion reaction leading to the same compound nucleus as in the 44Ca + 243Pu reaction. The parameters relevant to the exit channels, such as the neutron separation energies and the fission barrier heights, are nearly identical for both reactions at the same excitation energy (neglecting small differences in angular momentum of the compound nucleus after its formation). However, as shown in Fig. 3(c), the use of an 40Ar beam is less favorable than 36S. This is attributable to a worse fusion probability of the 40Ar + 247Cm fusion reaction, as the dinuclear system becomes more symmetric, and the fusion probability decreases. Our calculations also demonstrated that the use of a 44Ca beam instead of 36S decreases the yield of the same SHN, owing to a worse fusion probability. Calculations were performed for the reactions 36S + 243−246Cm, 40Ar + 239−242Pu, and 44Ca + 235−238U to produce the superheavy nuclei Cn, as shown in Fig. 4. The strong dependence of the calculated ERCSs for the production of SHN on the mass asymmetry in the entrance channel makes the 36S projectile most promising for the further synthesis of SHN.
-
The calculated maximal ERCSs, 3
n and 4n , and the corresponding optimal excitation energies of the compound nuclei in the 3n and 4n evaporation channel are presented in Fig. 5 for the reactions 36S + APu as functions of the mass number A of the target, respectively. In the 3n and 4n emission channel for the 36S + APu reaction, it Fig. 5(b) shows that the maximum ERCSs increases with the increase in neutron number to the maximum value and then decreases with the further increase in the neutron number. The lower part of Fig. 5(c)indicates that the excitation energies of 3 and 4 neutron emission decrease slowly with the increase in neutron number. To analyze the trend of the change above, the whole process of SHN synthesis needs to be investigated in detail. Next, we investigate the influence of the target neutron number on the capture cross-section, fusion probability, and survival probability.Figure 5. Isospin dependence from 36S + APu hot fusion reactions: (a) Q values for fusion reactions 36S + APu; (b) maximal evaporation residue cross sections as functions of target mass number A, for 3n and 4n emission channels; (c) corresponding excitation energies of compound nuclei.
Fig. 6(a) shows that the capture cross-section as a function of the incident energy is quite close for the three above-mentioned reactions owing to a slight difference in Coulomb barriers. Fig. 6(b) shows the capture cross-section
σcap as a function of the excitation energy of the compound nucleus. In the lower excitation energy regionE∗CN<37 MeV, the capture cross-sections for the reactions 36S + 236Pu are larger than those of the reaction systems 36S + 244Pu, because of the large negativeQ values (E∗CN=Ec.m.+Q ) of the former reactions. When the excitation energy increases beyond 37 MeV, the differences caused by the Q values become less significant, and the capture cross-sections tend to be almost consistent.Figure 6. (a) Calculated capture cross-sections as functions of incident energy in center-of-mass frame for reactions 36S + 236, 240, 244Pu. (b) Calculated capture cross-sections are functions of excitation energy of compound nucleus.
Considering only the effects of mass asymmetry on the fusion probability, when the neutron number in the target nucleus increases, the dinuclear system becomes more asymmetric, and the fusion probability increases [30, 83]. Fig. 7(a) shows that our calculated fusion probability is a function of the excitation energy of the compound nucleus for the reaction 36S + 236, 240, 244Pu. However, we are aware that during the calculation for 36S + 236Pu and 36S + 240Pu, the fusion probability
PCN decreases with increasing neutron number in the lower excitation energy regionE∗CN<45 MeV. When excitation energy increases beyond 45 MeV, the differences between fusion probabilities among the three reactions become very small, and the results tend to be consistent. The irregular behavior of the fusion probabilityPCN changes with increasing neutron number of the targets. The fusion probability depends on the details of the driving potential, which is decided by the ground stateQgg value of the nuclei in each DNS and their interactions [54]. Therefore, our results show that in addition to the mass asymmetry, the reactionQgg value plays an important role in the fusion probability.Figure 7. (a) Calculated fusion probabilities are functions of excitation energy of compound nucleus for reactions 36S + 236, 240, 244Pu. (b) Calculated survival probabilities as functions of compound nucleus excitation energy.
The survival probability of the hot compound nucleus is highly sensitive to the value of the neutron separation energy
Bn and fission barrierBf . The fission barriers of the 3n evaporation channel for compound nuclei 272−280Ds basically increase with an increase in the neutron number to the maximum value, and subsequently decrease with the further increase in neutron number [80]. The survival probabilities of 3n evaporation channel are shown in Fig. 7(b) for the compound nuclei 272, 276, 280Ds. The survival probabilities increase with the increase in the neutron number for excitation energies lower than 50 MeV. Therefore, for the 3n evaporation channel, the increase of survival probability of compound nuclei from 272Ds to 275Ds is not cancelled by an decreasing capture cross-section and fusion probability. The calculated maximal ERCSsσ3n is larger in the 36S + 239Pu reaction than in the 36S + 236−238Pu reactions owing to the larger value of survival probability. For the 4n evaporation channel, the change of the target neutron number has little influence on the fusion probability and capture cross-section. The variation trend of 4n ERCS with the neutron number of target nuclei is mainly determined by the survival probability. Therefore, the 4n evaporation channel corresponding to the ERCS basically increases with the increase in neutron number to the maximum value, and then decreases with the further increase in neutron number. -
Figure 8 shows the comparison of the calculated ERCS with the experimental data in the reaction 36S + 238U. The measured ERCSs of the 3
n and 4n channels are denoted by solid squares and open circles [61], respectively. For the 3n channel, calculated results are closer to the experimental data. We also realized that one calculated ERCS of the 4n channel in the reaction 36S + 238U peaks at the excitation energy of approximately 44 MeV, which is smaller than the excitation energy used in the experiment for 238U(36S, 4n )270Hs [61]. However, for the 4n evaporation channel, the peak position is not clear owing to the lack of experimental data. Taking into account the experimental error bars, we assume that the agreement between our calculated ERCS and the experimental value [61] is sufficient for the 36S + 238U reaction.Figure 8. Excitation functions for 36S + 238U reaction. Calculated results and experimental data [61] are denoted by lines and symbols, respectively. Measured ERCSs of 3n and 4n channels are denoted by solid squares and open cycles, respectively. Calculated 3n, 4n, and 5n channels are indicated by solid, dashed, and dotted lines, respectively.
The isospin dependence of the ERCS of some SHN is based on the same assumptions with one set of parameters. The maximum ERCS,
σER (pb), for 3n and 4n emission channels out of 36S bombarding actinide isotopic chains: AU, ANp, AAm, ACm, ABk, and ACf are shown in Fig. 9 as a function of the mass number of the target. From Figure 9 shows that the isotopes of target nucleus with the largest neutron excess are favorable for most cases of hot fusion with the 36S projectile. Except for the 3n emission channel in the 36S + ACm reaction, in all other channels the ERCSs basically increase with increasing neutron numbers, though sometimes not very distinctly.Figure 9. Isospin dependence of maximal evaporation residue cross-sections from hot fusion reactions: (a) 36S + AU, (b) 36S + ANp, (c) 36S + AAm, (d) 36S + ACm, (e) 36S + ABk, and (f) 36S + ACf as functions of target mass number A, for 3n and 4n emission channels.
Next, we investigate the influence of the target neutron number on ERCS. Similarly to the results of the reactions 36S + APu, our calculations show that the capture cross-sections for the neutron-deficient target nucleus are larger than in the neutron-rich target nucleus in the lower excitation energy region for the above projectile target combinations. We are also aware that during the calculation, the fusion probability
PCN for 36S + AU, 36S + ANp, 36S + AAm, 36S + ACm, 36S + ABk, and 36S + ACf , the PCN change with increasing neutron number of the target is not regular. The behavior of the fusion probabilityPCN with increasing neutron number of the targets depends on the details of the driving potential, which is decided by the properties of nuclei in each DNS and their interactions.For 3
n emission, in most cases, the increases inWsur are not cancelled by the decreasingσcap and the irregularly changingPCN with increasing neutron number. Therefore, the 3n evaporation channel corresponding to the ERCS basically increases with the increase in the neutron number. However, for compound nuclei 278−286Cn, the 3n evaporation channel 275−283Cn corresponding to the fission barriers [80] basically decreases to the minimum value with the increase in neutron number, and then increases with a further increase in neutron number. Therefore, the ERCSs of the reactions 36S + ACm decrease to the minimum value with the increase in neutron number, and then increase with the further increase in neutron number.For the 4
n evaporation channel, the change of the target neutron number has little influence on the fusion probability and the capture cross-section. The variation trend of the 4n ERCS with the neutron number of target nuclei is mainly determined by the survival probability. Therefore, the 4n evaporation channel corresponding to the ERCS basically increases with the increase in neutron number. Figure 9 shows that the ERCSs of SHN decrease slowly with the charged numbers of compound nuclei fromZ=108 toZ=114 owing to the increase in survival probabilityWsur , which is not canceled by the decreasingPCN .Currently, there is a gap between the SHN synthesized by cold fusion and those by hot fusion [7]. Considerably large cross-sections (
σ1n⩾ 1 pb) for many reaction channels used to fill the gap. For example, as shown in Figs. 9(b), given the predicted excitation function ofxn ERCSs for the reaction 36S + ANp. For the 36S+235Np reaction, the maximal ERCSs of the 4n channel is 2.74 pb atEc.m.=168.39 MeV (Ec.m.=E∗CN+Q ). It is 17.14 pb atEc.m.=165.11 MeV for the 36S + 237Np reaction. Through the analysis of Fig. 5(c) and Fig. 9, we found that the nucleus 267, 269Mt, 268, 272, 274−276Ds, 273, 275, 276Rg, 274−276, 278−280Cn, 279−281Nh, and 281−283Fl may be produced. Figure 10 shows that at fixed charge asymmetry in the entrance channel, the optimal excitation energyE∗CN of compound nucleus decreases with increasing neutron excess in the target. One can expect large ERCSs in the actinide-based reactions with the 36S beam. Moreover, based on the targets 235−237Np, 241−243Am, and 247−249Bk, the production of odd SHN withZ= 109, 111, and 113, are promising.Figure 10. Isospin dependence of excitation energies of compound nuclei corresponding to maximal evaporation residue cross-sections from hot fusion reactions: (a) 36S + AU, (b) 36S + ANp, (c) 36S + AAm, (d) 36S + ACm, (e) 36S + ABk, and (f) 36S + ACf as functions of target mass number A, for 3n and 4n emission channels.
-
To investigate the conditions for the synthesis of new neutron-deficient SHN, the projectiles 36S, 40Ar, and 44Ca bombarding some actinide isotope chains are systematically studied within the DNS model. Our results demonstrate that the strong dependence of the calculated ERCSs for the production of SHN on mass asymmetry in the entrance channel makes the 36S projectile the most promising for further synthesis of new neutron-deficient SHN. The influence of the target neutron number on ERCSs in hot fusion reactions is also investigated. There is a certain probability to produce new neutron-deficient SHN by using 36S to bombard 232−238U, 236−244Pu, 242−250Cm, and 249−251Cf. Thus, one can expect to produce new neutron-deficient SHN of the Hs, Ds, Cn, and Fl with ERCS ranging from > 1 pb up to 10 pb. Some new nuclides of Mt, Rg, and Nh may be produced by the reaction channels 36S + 235−237Np, 241−243Am, 247−249Bk with ERCS larger than 1 pb. Hopefully, the results will shed light on the experimental synthesis of new nuclides.
Possibilities for synthesis of new neutron-deficient isotopes of superheavy nuclei
- Received Date: 2019-01-28
- Available Online: 2019-05-01
Abstract: This study investigates the optimal projectile/target combination for the production of new neutron-deficient isotopes of superheavy nuclei (SHN). To this end, the dependence of the evaporation residue cross-section (ERCS) used to synthesize SHN on the mass asymmetry and the isospin of colliding nuclei are analyzed within the dinuclear system (DNS) concept. The predicted ERCSs for the production of new neutron-deficient isotopes of SHN were found to be quite large with the 36S projectile, and the cross-section of SHN decreases slowly with the charge of compound nuclei owing to the increase in their survival probability,