Possibilities for synthesis of new neutron-deficient isotopes of superheavy nuclei

Figures(10)

Get Citation
Xiao-Jun Bao. Possibilities for synthesis of new neutron-deficient isotopes of superheavy nuclei[J]. Chinese Physics C, 2019, 43(5): 054105. doi: 10.1088/1674-1137/43/5/054105
Xiao-Jun Bao. Possibilities for synthesis of new neutron-deficient isotopes of superheavy nuclei[J]. Chinese Physics C, 2019, 43(5): 054105.  doi: 10.1088/1674-1137/43/5/054105 shu
Milestone
Received: 2019-01-28
Article Metric

Article Views(2298)
PDF Downloads(38)
Cited by(0)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Possibilities for synthesis of new neutron-deficient isotopes of superheavy nuclei

  • Department of Physics, Collaborative Innovation Center for Quantum Effects, and Key Laboratory of Low Dimensional Quantum Structures and Quantum Control of Ministry of Education, Hunan Normal University, Changsha 410081, China

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, Wsur. Wsur is not canceled by the decreasing probability, PCN, that the system will evolve from a touching configuration to the compound nucleus in competition with the quasifission process.

    HTML

    1.   Introduction
    • 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 [1226], and many approaches were proposed to calculate the fusion probability [2734]. 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 [4060].

      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.

    2.   Theoretical framework
    • 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.JmaxJ=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]. The PCN is the probability that the system evolves from a touching configuration to the compound nucleus in competition with the quasifission process. The last term Wsur is the survival probability of the formed compound nucleus, which can be estimated with a statistic method [63].

    • 2.1.   Capture cross section

    • 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.B22μR2BJ(J+1)]}dB,

      (3)

      where ω(J) is the width of the parabolic Coulomb barrier at the position RB(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].

    • 2.2.   Fusion probability

    • 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 with Z1 and N1 with the local excitation energy ε1 at time t is described by the following ME:

      dP(Z1,N1,ε1,t)dt=Z1WZ1,N1;Z1,N1(t)×[dZ1,N1P(Z1,N1,ε1,t)dZ1,N1P(Z1,N1,ε1,t)]+N1WZ1,N1;Z1,N1(t)×[dZ1,N1P(Z1,N1,ε1,t)dZ1,N1P(Z1,N1,ε1,t)][Λqf(Θ(t))+Λfs(Θ(t))]P(Z1,N1,ε1,t),

      (4)

      where WZ1,N1;Z1,N1 is the mean transition probability from channel (Z1,N1) to (Z1,N1), and dZ1,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 fragment Z1, N1 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 variable R 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 distribution P(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((Γ22)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 barrier Bqf decreases with increasing Z, 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 barrier Bqf 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 expression a=A/12 MeV1. 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 time t 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 probability PCN is the summation of P(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 momentum J, is given by

      PCN(Ec.m.,J,Bm)=ZBGZ1=1NBGN1=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 energy Ec.m., J and Bm, and it is determined by using the deflection function method [68]. Finally, we obtain the fusion probability PCN(Ec.m.,J) as

      PCN(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 value Ji of the incident angular momentum, τint(Ji) is determined by the rotation of the composite system through the angle

      Figure 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 energy ECN=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 energies Ei, Ef and the angular momenta Ji, 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, and bi(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 as Ji=J. The details of δ and β are given in Ref. [70]. The grazing angular momentum Jgr can be expressed:

      Jgr=0.22Rint[Ared(Ec.m.V(Rint))]1/2,

      (11)

      where V(Rint) denotes the interaction barrier at the interaction radius Rint. 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=τint0dtJ(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+(JiJst)exp(t/τJ),

      (13)

      where the limiting value Jst given by the sticking condition is Jst=Jiζ0rel/ζ0tot. The relaxation time τJ is 1.5×1021 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); where m1, 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 energies ECN=40 MeV. Fig. 1 shows that the interaction time of the composite system is long for partial waves with a small incident angular momentum J. Moreover, we found that the interaction time decreases with increasing J. 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 momentum J. 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 energy Ex is related to the minimum of the nucleus-nucleus potential (Bm) and is determined for each initial relative angular momentum J 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 quantity Eirad denotes the initial radial kinetic energy at the interaction radius. The radial energy at the initial state is Eirad(J,0)=Ec.m.BmEJi(Ji+ 1)2/(2ζrel). The initial angular momentum is assumed to be Ji=J. The value of τrad is 3×1022 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 at ECN=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 and N=N1+N2, and βi(i=1,2) and β represent quadrupole deformations of the two fragments and the compound nucleus, respectively. The B(Z1,N1,β1), B(Z2,N2,β2), and B(Z,N,β) are the binding energies of two deformed nuclei and the compound nucleus [72], respectively. The Qgg (Qgg=B(Z1,N1,β1)+B(Z2,N2,β2)B(Z,N,β)) denotes the ground state Q 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,N1, 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 probability WZ1,N1,β1,β2;Z1,N1,β1,β2 can be written as [64, 65]

      WZ1,N1;Z1,N1(t)=τmem(Z1,N1,ε1;Z1,N1,ε1)2dZ1,N1dZ1,N1×ii|<Z1,N1,ε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,N1,ε1)=2π{<V2(t)>Z1,N1,ε1+<V2(t)>Z1,N1,ε1}1/2,

      (19)

      can be interpreted as the coherence time for the transitions between the subsets (Z1,N1,ε1) and (Z1,N1,ε1) [64, 65], where <V2(t)>Z1,N1,ε1 and <V2(t)>Z1,N1,ε1 stand for the average expectation value with Z1,N1,ε1 and Z1,N1,ε1 being fixed, respectively. Thus, the memory time τmem depends on the neutron number N1, proton number Z1 , and the local excitation energy ε1. For the 36S + 250Cf reaction, the memory time τmem(Z=16,N=20;Z=16,N = 21) is 0.75×1022 s when the excitation energy ECN = 40 MeV.

      The transition probability of Eq.(18) can be written as

      WZ1,N1;Z1,N1(t)=τmem(Z1,N1,ε1;Z1,N1,ε1)2dZ1,N1dZ1,N1×{[ω11(Z1,N1,ε1;ε1)+ω22(Z1,N1,ε1;ε1)]δN1,N1+ω12(Z1,N1,ε1;ε1)δN1,N11+ω12(Z1,N1,ε1;ε1)δN1,N1+1},

      (20)

      where

      ωkk(Z1,N1,ε1;ε1)=k,kN1|<Z1,N1,ε1,i|Vk,k|Z1,N1,ε1,i>|2=dZ1,N1<Vk,kV+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,kV+k,k>=14U2kkgkgkΔ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 parameters Ukk(t) are assumed as [64]

      Ukk(t)=g1/31g1/32g1/31+g1/321g1/3kg1/3k2γ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 and mk=Nk/2 valence nucleons in Δεk.

    • 2.3.   Survival probability

    • 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(ECN,x,J)=F(ECN,x,J)xi=1[Γn(Ei,J)Γn(Ei,J)+Γf(Ei,J)]i,

      (26)

      where F(ECN,x,J) is the realization probability of the xn channel at the excitation energy ECN(Ec.m.+Q) of the compound nucleus with angular momentum J, 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π22ρ0(E0δ0)E0Bnδn0ρn(E0Bnδnε)εdε,

      (27)

      where mn and g 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(E0Bnδ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)E0Bfδf0ρn(E0Bfδfε)dε,

      (28)

      where ρn(E0Bfδ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[2aUJ(J+1)2σ2]242σ3a1/4U5/4,

      (30)

      with σ2=Θrigid2Ua, Θ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 parameter a on the shell correction and the excitation energy was initially proposed as

      a(U,Z,N)=˜a(A)[1+Eshf(U)U]

      (31)

      with ˜a(A)=αA+βA2/3 and f(U)=1exp(γ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].

      Figure 2.  Calculated transmission, fusion, and survival probabilities as functions of relative angular momenta in the reaction 36S + 250Cf at excitation energies of compound nucleus of 40 MeV.

    3.   Numerical results and discussions
    • 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 Z104 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.

    • 3.1.   Predictions of probable projectile-target combinations

    • 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 from Z=16 to Z=20. This is due to the strong decrease in fusion probability PCN 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.

      Figure 4.  (color online) Evaporation residue excitation functions in production of isotopes of superheavy nuclei Cn in reactions 36S + 243−246Cm, 40Ar + 239−242Pu, and 44Ca + 235−238U.

    • 3.2.   Influence of the target neutron number on ERCSs

    • The calculated maximal ERCSs, 3nand 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 region ECN<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 negative Q values (ECN=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 region ECN<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 probability PCN 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 state Qgg value of the nuclei in each DNS and their interactions [54]. Therefore, our results show that in addition to the mass asymmetry, the reaction Qgg 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 barrier Bf. 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.

    • 3.3.   Production cross-sections of neutron-deficient SHN

    • Figure 8 shows the comparison of the calculated ERCS with the experimental data in the reaction 36S + 238U. The measured ERCSs of the 3n 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 probability PCN 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 3n emission, in most cases, the increases in Wsur are not cancelled by the decreasing σcap and the irregularly changing PCN 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 4n 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 from Z=108 to Z=114 owing to the increase in survival probability Wsur, which is not canceled by the decreasing PCN.

      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 of xn ERCSs for the reaction 36S + ANp. For the 36S+235Np reaction, the maximal ERCSs of the 4n channel is 2.74 pb at Ec.m.=168.39 MeV (Ec.m.=ECN+Q). It is 17.14 pb at Ec.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 energy ECN 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 with Z=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.

    4.   Conclusions
    • 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.

Reference (83)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return