Calculations of fusion capture of breakup fragments of the weakly bound projectile 11Be with several targets

Figures(4) / Tables(1)

Get Citation
A. Gómez Camacho and J. Lubian. Calculations of fusion capture of breakup fragments of the weakly bound projectile 11Be with several targets[J]. Chinese Physics C. doi: 10.1088/1674-1137/ae71a4
A. Gómez Camacho and J. Lubian. Calculations of fusion capture of breakup fragments of the weakly bound projectile 11Be with several targets[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ae71a4 shu
Milestone
Received: 2026-02-26
Article Metric

Article Views(36)
PDF Downloads(1)
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:

Calculations of fusion capture of breakup fragments of the weakly bound projectile 11Be with several targets

  • 1. Departamento de Aceleradores y Estudio de Materiales, Instituto Nacional de Investigaciones Nucleares, Ciudad de México 11801, México
  • 2. Instituto de Fisica, Universidade Federal Fluminense, Niteroi 24210-340, Brazil

Abstract: Calculations of the effect of the breakup reactions of the weakly bound projectile 11Be on total, complete, and incomplete fusion with targets 16O, 28Si, 58Ni, 144Sm, and 209Bi are presented. The calculations primarily focus on the relative effects of projectile neutron and 10Be fragments on incomplete fusion. In fact, the energy-dependent contributions to the incomplete fusion of the fragments exhibit highly interesting behaviors when targets of different masses and charges are considered. Incomplete fusion for neutron absorption becomes more important than that for $ ^{10}\mathrm{Be} $ for all targets at energies around and above the barrier. Indeed, neutron capture becomes increasingly important as the target becomes heavier. In addition, it is found that incomplete fusion becomes larger than complete fusion for energies lower than a certain value that depends on the target. Above this value, complete fusion overcomes incomplete fusion. As a final calculation, an energy-dependent systematic comparison of the total incomplete fusion, incomplete fusion of the neutron, and incomplete fusion of 10Be for the different targets is presented. The calculations are performed with two complementary theoretical approaches. The continuum discretized coupled channel model is used to determine the relative projectile-target radial wave functions, which are subsequently used in the angular-momentum-dependent model of fusion probabilities to calculate complete, incomplete, and total fusion cross sections.

    HTML

    I.   INTRODUCTION
    • In recent years, the study of nuclear reactions of weakly bound projectiles with stable targets has been a topic of intense theoretical and experimental research [15]. One of the most important aspects of these studies is to comprehend the effects that the breakup channels of the projectile have on other reaction channels such as elastic, transfer, and fusion. Weakly bound projectiles have the characteristic that, when interacting with a target, they break up into several fragments. Basically, this occurs because of the low binding energy of these nuclei. Therefore, new and different fusion processes can be present for such projectiles. Direct Complete Fusion (DCF) is a process that appears when the entire projectile is captured by the target without a previous breakup. There may be cases where the projectile breaks up but, afterward, all fragments are absorbed by the target. This is known as Sequential Complete Fusion (SCF). Thus, the sum of DCF and SCF constitutes Complete Fusion (CF). Experimentally, DCF and SCF are very difficult to distinguish. Another very important fusion mechanism is Incomplete Fusion (ICF). ICF of a given fragment occurs when, following breakup, this fragment is captured by the target while the others fly away to continuum states. Of course, the total ICF corresponds to the sum of all partial ICFs of the fragments. Hence, CF plus ICF determines the Total Fusion (TF) of the weakly bound projectile to the target. Certainly, it is possible that none of the fragments are captured by the target after breakup, which is known as the Non-capture Breakup (NCBU) process. NCBU may take place in several manners. For instance, elastic breakup occurs when the projectile and target remain in their ground states following breakup. In contrast, inelastic breakup occurs when the entire projectile, projectile fragments, and/or target are in excited states before or after breakup. A full and realistic description of the effects of breakup processes on fusion must consider couplings to all degrees of freedom involved [6]. Furthermore, reactions involving weakly bound nuclei have a higher degree of complexity at sub-barrier energies. This is because, for these projectiles, transfer reactions followed by breakup are essential process that may predominate over direct breakup [710]. From a theoretical perspective, studies that consider all of the above breakup excitations become extremely complex. For this reason, and for the sake of reducing the complexity, most theoretical studies consider only the most important reaction channels that could affect fusion, such as elastic, some inelastic, and direct breakup channels. The continuum discretized coupled channel (CDCC) method is possibly the most widely used model for calculating fusion for weakly bound nuclei [1113]. Using the CDCC model, it is possible to calculate the effects of the continuum breakup states of the projectile on several reaction observables, such as elastic scattering, elastic breakup, and, more recently, fusion. In fact, recent fusion CDCC calculations for weakly bound projectiles with a core-valence structure have been carried out [1415]. In these calculations, fusion of the core-valence fragments is determined using two short-range absorption potentials, $ W_{F_{1}} $ and $ W_{F_{2}} $, which account for the absorption of the fragments $ F_{i} $, $ (i=1,2) $. That is, the ICF$ _{F_{i}} $ of the fragment i corresponds to its capture through the potential $ W_{F_{i}} $. In actual calculations to determine ICF$ _{F_{i}} $, it is supposed that the other fragment is not absorbed, that is, it is transferred to the continuum. In this case, $ W_{F_{i}}{\ne 0} $ but $ W_{F_{j}}=0 $, with $ i,j=1,2 $, and $ i \ne j $. Then, the total ICF of the projectile to the target is $ ICF= ICF_{F_{1}}+ ICF_{F_{2}} $. The CF is determined in a different calculation when the entire projectile is captured without breakup by an absorption potential $ W_{PT} $. Although the calculations [1617] for CF and ICF may be close to experimental data, this approach has a shortcoming, that is, separate calculations are performed for the contributions to the ICF of each fragment and for CF. In this sense, the calculations are not self-consistent. The projected angular-momentum-dependent model of fusion probabilities is a recent model that simultaneously determines the DCF and SCF as well as the ICF of the fragments of valence-core weakly bound projectiles [1820]. In this model, fusion is not obtained from CDCC calculations, but instead, CDCC is used to obtain the radial projectile-target relative wave functions. These wave functions are, in turn, introduced into the model of fusion probabilities, from which fusion cross sections are calculated. In recent studies, this model was used for calculations of the TF, CF, and ICF of the weakly bound 8B projectile in reactions with light to heavy mass targets 16O, 27Al, 28Si, 58Ni, 208Pb, and 209Bi [21]. In this study, it was found that ICF is more important than CF for the targets 16O, 27Al, 28Si, 40Ar, and 58Ni for incident energies around and below the Coulomb barrier. However, this situation is inverted as the energy increases. That is, at a given energy that depends on the target mass, CF overcomes ICF. For the heavier targets 208Pb and 209Bi, CF is the most important component of TF for energies around and above the barrier. In addition, the ICF of the breakup fragments, a proton and 7Be, of the nucleus 8B was determined. In fact, a systematic study was performed on the specific contributions to the ICF of the 7Be and proton fragments for all targets. It was found that the ICF of proton dominates over that of 7Be for the lighter targets. However, as the mass of the target increases, the contrary occurs. For the heavy targets 208Pb and 209Bi, the ICF of 7Be accounts for most of the ICF for low and high energies. Following similar lines, in this study, a systematic investigation is presented for fusion excitation functions of the neutron-halo weakly bound projectile 11Be (cluster structure $ ^{10}\mathrm{Be}+n $, threshold energy $ E\mathrm{_{thre}}=0.5 $ MeV) with targets in a variety of masses. Specifically, CF and ICF cross sections as well as the contributions to ICF from the fragments of 11Be are determined for reactions with the targets 16O, 28Si, 58Ni, 144Sm, and 209Bi. Particular emphasis is placed on determining the relative contributions of CF and ICF to TF as a function of the target mass and incident projectile energy. Similarly, the contributions to ICF from the neutron and 10Be breakup fragments are investigated as a function of the target mass and incident energy. Experimental data for the ICF of the fragments of 11Be are scarce, and thus, a comparison of the calculations is not possible. However, from a theoretical perspective, it is interesting to attempt to understand how the neutron and 10Be fragments may interact with targets of different mass and charge and subsequently be absorbed. The fusion calculations are performed in two steps. (a) The radial wave functions of the relative motion between the incident projectile 11Be and target are obtained from CDCC calculations. (b) These radial wave functions are, in turn, used in the projected angular-momentum model of fusion probabilities [1820], from which the cross sections for the different fusion mechanisms are obtained. In Section II, brief descriptions of the CDCC model, which describes breakup states in the continuum of 11Be, and of the projected angular-momentum model of fusion probabilities are provided. In Section III, details of the construction of the CDCC discrete breakup space of 11Be are presented. In addition, we present the coupled-channel equations from which the radial wave functions of the projectile-target relative motion are obtained. Finally, the fusion cross sections for the different fusion reaction mechanisms are calculated. Section IV is dedicated to a summary of the results and conclusions.

    II.   CDCC AND PROJECTED ANGULAR- MOMENTUM MODEL OF FUSION PROBABILITIES
    • A complete description of the CDCC model for describing continuum breakup states of the projectile is provided in Refs. [1113]. This model is applied to quantify the effect of continuum breakup states on several reaction observables. In the present calculations, the CDCC approach is used to determine the radial angular-momentum wave functions of the projectile-target relative motion. These wave functions are, in turn, used in the angular-momentum-dependent model of fusion probabilities. This latter model is described in Refs. [1821]. A detailed description of the use of the radial wave functions of the projectile-target relative motion obtained by the CDCC model to calculate fusion probabilities is provided in the appendix of Ref. [19].

      In the CDCC model, the Hamiltonian of the projectile-target dynamics is expressed as

      $ \hat{H}({\boldsymbol{R}},{\boldsymbol{r}})=\hat{h}({\boldsymbol{r}})+\hat{T}+\hat{U}_{1}({\boldsymbol{r}}_{1})+\hat{U_{2}}({\boldsymbol{r}}_{2}), $

      (1)

      where $ \hat{h}({\boldsymbol{r}}) $ is the internal Hamiltonian of the projectile, $ \hat{T} $ the kinetic-energy operator of the relative motion, and $ \hat{U}_{1}({\boldsymbol{r}}_{1}) $, $ \hat{U}_{2}({\boldsymbol{r}}_{2}) $ are the interaction potentials between the projectile fragments and target. $ {\boldsymbol{r}}_{1} $ and $ {\boldsymbol{r}}_{2} $ are the radial distances between the centers of mass of the target and fragment $ F_{i},i=1,2 $, whereas $ {\boldsymbol{R}} $ is the distance from the c.m. of the target to the c.m. of the fragments. Then,

      $ {\boldsymbol{R}}={\boldsymbol{r}}_{\boldsymbol{2}}+\frac{A_{1}}{A}{\boldsymbol{r}}, \ \ \ {\boldsymbol{r}}_{\boldsymbol{1}}={\boldsymbol{R}}+\frac{A_{2}}{A}{\boldsymbol{r}}, \ \ \ A=A_{1}+A_{2}, $

      (2)

      where $ A_{1}, A_{2} $ are the fragment masses and A is that of the target. The internal Hamiltonian of the projectile is expressed as

      $ \hat{h}({\boldsymbol{r}})= \frac{{{\boldsymbol{p}}}^2}{2\mu_{\boldsymbol{12}}}\ + \ V_{\boldsymbol{12}}({{\boldsymbol{r}}}), $

      (3)

      where $ V_{\boldsymbol 12}({{\boldsymbol{r}}}) $ includes Coulomb and nuclear potentials. The nuclear interaction usually has a volume Woods-Saxon shape. The ground, inelastic, and scattering breakup states of the projectile $ \Phi_{i} $ satisfy

      $ \hat{h}({\boldsymbol{r}}) \Phi_{\beta} ({\boldsymbol{r}}) =\epsilon_{\beta} \Phi_{\beta} ({\boldsymbol{r}}), \ \ \ {\beta}=0,1,\cdots,N. $

      (4)

      The interaction potentials $ \hat{U}_{1}({\boldsymbol{r}}_{1}) $ and $ \hat{U_{2}}({\boldsymbol{r}}_{2}) $ in Eq. (1) are expressed as

      $ \hat{U}_{i}({\boldsymbol{r}}_{i})=\hat{V_{i}}({\boldsymbol{r}}_{i})-{\mathrm{i}}\hat{W_{i}}({\boldsymbol{r}}_{i}),\ i=1,2, $

      (5)

      where $ \hat{V_{i}}({\boldsymbol{r}}_{i}) $ are the interactions between the projectile fragments and target, and $ {\boldsymbol{r}}_{i} $ is the fragment-target distance. $ \hat{W}_{i}({\boldsymbol{r}}_{i}) $ represents the corresponding absorption potential. The total absorption potential is

      $ \hat{W}({\boldsymbol{r}}_{1},{\boldsymbol{r}}_{2})= \hat{W}_{1}({\boldsymbol{r}}_{1})+ \hat{W}_{2}({\boldsymbol{r}}_{2}). $

      (6)

      The total wave function $ \Psi^{(+)} $ of the projectile-target nuclear system is obtained by expanding it in terms of the eigenstates of the projectile $ \Phi_{\beta}({\boldsymbol{r}}),\ \beta=0\dots N $. That is,

      $ \Psi^{(+)}({\boldsymbol{R}},{\boldsymbol{r}})=\sum\limits_{\beta=0}^{N}\psi_{\beta}\mathbf({\boldsymbol{R}})\otimes \Phi_{\beta}({\boldsymbol{r}}), $

      (7)

      where $ \psi_{\beta}({\boldsymbol{R}}) $ are the radial relative projectile-target wave functions. $ \beta=0 $ refers to the elastic incident channel, whereas $ \beta=1...N $ denotes the inelastic and excited breakup states. The radial relative wave functions $ \psi_{\beta}(R) $ of Eq. (7) are used to calculate the probabilities of the different fusion mechanisms [1820]. These wave functions can be determined from CDCC calculations by solving the coupled-channel equations obtained from Eqs. (1) and (7) for a discretized breakup energy space that includes the bound (B) and continuum (C) states of the projectile. These coupled-channel equations are expressed as

      $\begin{aligned}[b] &\left[ \hat{T}(R)+U_{\beta, \beta}(R)-(E-\varepsilon _{\beta}-\varepsilon_{T})\right] \psi_{\beta}(R)\\&= -\sum_{\beta^{\prime}} U_{\beta,\beta ^{\prime}}(R)\psi_{\beta ^{\prime }}(R), \end{aligned} $

      (8)

      where the coupling matrices are

      $ U_{\beta\beta ^{\prime }}(R)= <\phi_{\beta }\left\vert \hat{U}_{1}({\boldsymbol{r}}_{1})+\hat{U}_{2}({\boldsymbol{r}}_{2})\right\vert \phi_{\beta ^{\prime }}>. $

      (9)

      Here, $ \varepsilon _{\beta } $ is the excitation energy of the projectile in the β-state, $ \varepsilon_{T} $ is the ground state energy of the target (assumed to be $ \varepsilon_{T}=0 $), and $ U_{\beta, \beta} $, $ U_{\beta,\beta ^{\prime}} $ ($ \beta^{\prime}\neq \beta $) are the radial dependent diagonal and non-diagonal coupling matrices of the interaction potentials between the ground, inelastic, and continuum breakup states $ \beta,\beta ^{\prime}=0 \dots N $. The CDCC bin states $ \phi_{\beta }(r) $ are defined by square-integrable wave functions as follows:

      $ \phi_{\beta}(r) = \sqrt{\frac{2}{\pi N_{\beta}}} \int^{k_i}_{k_{i-1}}w_{\beta}(k)f_{\beta}(k,r){\mathrm{d}}k. $

      (10)

      This wave packet is generated by superposing radial scattering wave functions $ f_{\beta}(k,r) $, where $ f_{\beta}(k,r) $ is the radial part of $ \Phi_{\beta}({\boldsymbol{r}} $) with angular momentum l and parity $ \pi=(-1)^l $ in the momentum interval [$ k_{i-1},k_{i} $], or equivalently, the [$ \epsilon_{i-1},\epsilon_i $] energy interval. $ w_{\beta} $ is a weight function and $ N_{\beta} = \displaystyle\int^{k_i}_{k_{i-1}}w_{\beta}(k){\mathrm{d}}k $ is a normalization constant.

      Now, the TF cross-section for a reaction of a weakly bound projectile with a two-particle cluster structure can be written as

      $ \sigma_{TF}=\sigma_{DCF}+\sigma_{SCF}+\sigma_{ICF}, $

      (11)

      where $ \sigma_{DCF} $ is the DCF and $ \sigma_{SCF} $ is the SCF. Then,

      $ \sigma_{CF}=\sigma_{DCF}+\sigma_{SCF} $

      (12)

      defines the total CF. In turn, the ICF $ \sigma_{ICF} $ is expressed as

      $ \sigma_{ICF}=\sigma_{ICF_{1}}+\sigma_{ICF_{2}}, $

      (13)

      where $ \sigma_{ICF_{i}}\ (i=1,2) $ is the ICF of a given fragment. Note that $ \sigma_{DCF} $ is associated with the bound states (B) of the projectile, whereas $ \sigma_{SCF} $ and $ \sigma_{ICF} $ are associated with the continuum (C) states.

      In terms of the bound and continuum states, the total wave function $ \Psi^{(+)} $ of Eq. (7) can be split as follows:

      $ \Psi^{(+)}({\boldsymbol{R}},{\boldsymbol{r}})=\Psi^{(B)}({\boldsymbol{R}},{\boldsymbol{r}})+\Psi^{(C)}({\boldsymbol{R}},{\boldsymbol{r}}), $

      (14)

      where

      $ \Psi^{(B)}(\boldsymbol{R},\boldsymbol{r})=\sum\limits_B^{ }\psi_B\mathbf{(}\boldsymbol{R})\otimes\phi_B(\boldsymbol{r}) $

      (15)

      and

      $ \Psi^{(C)}({\boldsymbol{R}},{\boldsymbol{r}})=\sum\limits_{C}\psi_{C}\mathbf({\boldsymbol{R}})\otimes \phi_{C}({\boldsymbol{r}}), $

      (16)

      in which $ \psi_{B} $, $ \psi_{C} $ are the radial projectile-target relative wave functions for the bound and continuum states, and $ \phi_{B} $, $ \phi_{C} $ are the corresponding bound and breakup states of the projectile. In our calculations, it is assumed that the matrix elements of the imaginary potentials connecting the bound and breakup states of the projectile are negligible [14, 2223]. That is,

      $ U_{\beta\beta ^{\prime }}(R)= <\phi_{\beta }\left\vert \hat{W}_{1}({\boldsymbol{r}}_{1})+\hat{W}_{2}({\boldsymbol{r}}_{2})\right\vert \phi_{\beta ^{\prime }}>=0, $

      (17)

      if β, $ \beta ^{\prime }=B,C $, with $ \beta\ne\beta ^{\prime } $. Then, the following relations hold:

      $ \sigma_{TF}=\sigma_{TF}^{B}+\sigma_{TF}^{C}, $

      (18)

      where

      $ \sigma_{TF}^{B}= \frac{K}{N^{2}E}\sum\limits_{\beta \beta^\prime\epsilon B} <\psi_{\beta} \left\vert W^{(1)}_{\beta \beta^\prime}+ W^{(2)}_{\beta \beta^\prime} \right\vert \psi_{\beta^\prime}>, $

      (19)

      and

      $ \sigma_{TF}^{C}= \frac{K}{N^{2}E}\sum\limits_{\gamma \gamma^\prime\epsilon C} <\psi_{\gamma} \left\vert W^{(1)}_{\gamma \gamma^\prime}+ W^{(2)}_{\gamma \gamma^\prime} \right\vert \psi_{\gamma^\prime}>. $

      (20)

      Therefore, for the bound states,

      $ \sigma_{TF}^{B}=\sigma_{DCF}^{B}. $

      (21)

      This cross-section is constructed by the angular momentum expansion of probabilities to absorb the entire projectile from its bound state, i.e.,

      $ \sigma_{DCF}^{B}=\frac{\pi}{K^2} \sum(2J+1) \ {\cal{P}}^{(DCF)}(J), $

      (22)

      where K is the wave number in the entrance channel and $ {\cal{P}}^{(DCF)}(J) $ is provided in terms of the imaginary potential that depends on the projectile-target coordinates, i.e.,

      $ {\cal{P}}^{(DCF)}(J) = \frac{4K}{E} \int {\mathrm{d}}r\, W_{PT}(R)\, \left| \psi_{0}(R) \right|^2, $

      (23)

      where J is the total angular momentum, $ \psi_{0}(R) $ is the radial ground state wave function, and $ W_{PT}(R) $ is the projectile-target imaginary potential, without considering the cluster structure of the projectile (the physical reasons for this assumption are discussed in Ref. [20]).

      The absorption from continuum states in Eq. (18) is expressed as

      $ \sigma_{TF}^{C}=\sigma_{ICF_{1}}+\sigma_{ICF_{2}}+\sigma_{SCF}. $

      (24)

      In terms of fusion probabilities [20],

      $ \sigma_{SCF} = \frac{\pi}{K^2}\,\sum\limits_J\ (2J+1)\, {\cal{P}}^{SCF}(J), $

      (25)

      $ \sigma_{ICF1} = \frac{\pi}{K^2}\,\sum\limits_J\ (2J+1)\, {\cal{P}}^{ICF1}(J), $

      (26)

      $ \sigma_{ICF2} = \frac{\pi}{K^2}\,\sum\limits_J\ (2J+1)\, {\cal{P}}^{ICF2}(J), $

      (27)

      with

      $ {\cal{P}}^{SCF}(J) = {\cal{P}}^{(1)}(J) \times {\cal{P}}^{(2)}(J), $

      (28)

      $ {\cal{P}}^{ICF1}(J)= {\cal{P}}^{(1)}(J) \times \big[ 1\, -\, {\cal{P}}^{(2)}(J)\big] ,$

      (29)

      $ {\cal{P}}^{ICF2}(J)= {\cal{P}}^{(2)}(J) \times \big[ 1\, -\, {\cal{P}}^{(1)}(J)\big] . $

      (30)

      In these equations, $ {\cal{P}}^{(i)}(J), i=1, 2 $ are the probabilities of the fragment $ F_i $ being absorbed by the part $ \hat{W}_{i} $ of the imaginary potential (6) independently of what occurs with the other fragment. The expressions of these probabilities are provided in the appendices of Ref. [19].

    III.   CALCULATIONS OF CF, ICF, AND TF CROSS SECTIONS

      A.   CDCC calculations

    • In this work, systematic calculations of TF, CF, and ICF of reactions of the weakly bound nucleus 11Be with targets 16O, 28Si, 58Ni, 144Sm, and 209Bi are presented. As a first step, CDCC calculations of the radial projectile-target relative motion wave function $ \psi_{\beta}(R) $, $ \beta=0...N $ are performed, which are, in turn, used to determine the fusion probabilities for a given capture process. The FRESCO code [24] is used to solve the coupled-channel equations (8). As is the norm in CDCC calculations, exhaustive convergence tests are performed to solve the coupled-channel equations. Convergence is checked upon changes in the size of the integration step $ \Delta r $, matching radius $ R\mathrm{_{max}} $, maximum relative angular momentum $ l\mathrm{_{max}} $ between the projectile fragments, maximum potential multipole λ, maximum energy of the bin states $ \epsilon\mathrm{_{max}} $, energy width for the construction of the bin states $ \Delta\epsilon $, and relative angular momentum of the projectile-target $ L_{\mathrm{max}} $. For the calculations of the bound states of 11Be, the parameters of the nuclear potential $ V_{\boldsymbol{12}}({{\boldsymbol{r}}}) $ of Eq. (3), which includes a spin-orbit term, are those of Ref. [25]. These parameters yield the correct binding energy of the ground state $ 1/2^{+} $, $ E_{0}=0.5 $ MeV, first excited state $ 1/2^{-} $, $ E_{x}=0.32 $ MeV, and resonance $ J^{\pi} = 5/2^{+} $, $ E_{\mathrm{res}}=1.28 $ MeV, with $ \Delta E_{\mathrm{res}}=100 $ keV, for which finer discretization is required. Convergence is achieved with the discretized breakup space described in Ref. [26]: $ \Delta r =0.4 $ fm, $ R_{\mathrm{max}}=40 $ fm, bin states constructed up to maximum orbital angular momenta $ l\mathrm{_{max}}=4\hbar $, coupling matrix elements up to multipoles $ \lambda_{\mathrm{max}}=5 $, and maximum energy of bin states up to $ \epsilon\mathrm{_{max}}=8 $ MeV. It should be noted that some variations in the parameters may be necessary. These variations depend on the incident energy of the projectile and target mass. For instance, the maximum breakup energy $ \epsilon\mathrm{_{max}} $ may depend on the incident energy of the projectile. That is, $ \epsilon\mathrm{_{max}} $ should be reduced for energies close to or below the Coulomb barrier with respect to energies above the barrier. High excitation states of the projectile can not be reached for low incident energies. However, a parameter that depends on the target is the relative angular momentum of the projectile-target relative motion $ L\mathrm{_{max}} $. Fusion for heavier targets may require higher $ L\mathrm{_{max}} $ values than that for lighter ones.

    • B.   Absorption and nuclear potentials

    • The fusion cross sections (Eqs. (24)−(27)) are calculated using the CF-ICF computer code [27] that determines fusion probabilities in terms of the angular momentum scattering wave functions obtained from the CDCC calculations. The mathematical expressions of the probabilities for the capture of the fragments are provided in the appendix of Ref. [18]. The absorption potentials of the breakup components of the projectile $ ^{11}\mathrm{Be}\rightarrow^{10}\mathrm{Be}+n $ (Eq. (6)) are assumed to be Woods-Saxon short-range potentials,

      $ W_{i}(r_{i})=\frac{W_{0}}{1+{\mathrm{exp}}[(r_{i}-R_{i})/a_{0}]}, $

      (31)

      where $ i=1,2 $, $ r_{i} $ is the distance between the fragment $ F_{i} $ and center of mass of the target (T), and $ R_{i}=r_{0}(A_{i}^{1/3}+A_{T}^{1/3}) $, where $ A_{i} $ is the mass of fragment $ F_{i} $. The parameters of the potentials $ W_{1}(r) $ and $ W_{2}(r) $ should be selected such that these potentials are well inside the potential barriers of the fragments and target. For the DCF case, the imaginary potential used in this model is not the sum of the fragment-target potentials, but the projectile-target potential $ W_{PT} $. In this case, $ A_i $ is replaced with $ A = A_1+A_2 $. Sequential and incomplete absorption are not zero once either of the fragments, in their trajectories to the target, are inside the absorption potential region. Moreover, the parameters of the fusion absorption potentials $ W_{PT} $, $ W_{1} $, and $ W_{2} $ are fixed to the values $ W_{0}=50 $ MeV, $ r_{0}=1.0 $ fm, and $ a_{0}=0.2 $ fm, independently of the target. These potential parameters guarantee that the potentials $ W_{PT} (r) $, $ W _ {1} (r) $, and $ W_ {2} (r) $ are well inside the Coulomb barriers for all targets. For instance, for the heaviest $ ^ {209} $Bi with the projectile 11Be, the Sao Paulo Potential (SPP) determines a Coulomb barrier positioned at $ R_{B}=11.76 $ fm. This value is larger than the absorption potential radius $ R=8.157 $ fm. This condition is also satisfied for less massive targets. Therefore, the calculations of the fusion cross sections presented later do not depend on the mass of the target. In fact, the independence of fusion cross-sections has been shown in Ref. [28]. The real interactions $ \hat{V_{i}}({\boldsymbol{r}}_{i}) $ of Eq. (5) for the fragments $ F_{i} $, $ i=1,2 $ with the target include nuclear and Coulomb components. With respect to the nuclear components, we use the SPP [2930]. Although the code for calculating the SPP also yields the imaginary part, we consider only the real part. The imaginary part is that described above; that is, a short-range absorption potential. In addition, it should be noted that the calculation of the radial SPP allows the option to determine neutron-target potentials.

    • C.   Fusion cross sections

    • In this systematic study, our objective is to determine the different fusion processes of 11Be with targets ranging from light to heavy masses. Of particular interest is the ICF of the neutron and its role in the total ICF cross section as a function of the incident energy of the projectile. Figures 1(a)−(e) show the calculations of CF, ICF, and TF for 16O, 28Si, 58Ni, 144Sm, and 209Bi, respectively. As can be observed, from low energies, ICF predominates over CF up to a certain value $ E_{t} $ that depends on the target. In fact, $ E_{t} $ is augmented as the mass and charge of the target increase. For the target 16O, calculations below the barrier are challenging owing to numerical divergences. However, Fig. 1(a) shows that $ E_{t} $ should be well below that barrier. For the heavier 209Bi, ICF is larger than CF, so $ E_{t} $ is above the range of energies shown. As the projectile 11Be approaches the target and enters the region of the attractive nuclear potential on the neutron and the net (Coulomb and nuclear) potential on 10Be, projectile dissociation can be produced. The attractive nuclear potential on the neutron and net potential on 10Be become stronger as the mass and charge of the target increase. Consequently, breakup and ICF result in the most important processes for heavier targets. The contrary occurs for light targets. That is, the weaker potentials on the neutron and 10Be can produce dissociation if the interaction with the target is sufficiently slow, which may occur only for low projectile incident energies. Therefore, as observed, particularly in Figs. 1(a)−(b) for the light targets 16O and 28Si, ICF is larger than CF at low energies; however, as the energy increases, CF becomes increasingly important. As shown in Figs.1(c)−(d), as the strengths of the nuclear and Coulomb potentials of the target are augmented, ICF starts to dominate CF for larger energies.

      Figure 1.  (color online) TF, CF, and ICF cross sections for 11Be with targets 16O, 28Si, 58Ni, 144Sm, and 209Bi.

      In addition, it is interesting to investigate the individual contributions to ICF from the neutron and 10Be at low and high energies for the different targets. Figures 2(a)−(e)present the calculations for these contributions to ICF with the targets 16O, 28Si, 58Ni, 144Sm, and 209Bi. The calculations show that the ICF of the neutron accounts for most of the total ICF for all targets for energies above the barrier $ V_{B} $. In fact, the neutron ICF contributes more to the total ICF for more massive targets. We can attempt to understand this behavior as follows. Once the projectile breaks up, and for low energies and $ L=0 $ incident partial waves, the fragment 10Be must penetrate a barrier during the interaction with the target, whereas the neutron is easily captured (as there is no barrier to penetrate). Therefore, it expected that neutron ICF becomes more appreciable at low and high energies. As the energy increases, 10Be surpasses the barrier with the target and can be more easily captured. For $ L\ne 0 $ partial waves, both the neutron and 10Be must pass through the barrier,

      Figure 2.  (color online) Neutron, 10Be, and total ICF cross sections for 11Be with targets 16O, 28Si, 58Ni, 144Sm, and 209Bi. $ V_{B} $ is the barrier energy of the projectile 11Be and the target.

      $ V(r)=-V_{{\mathrm{nucl}}}(r)+V_{{\mathrm{coul}}}(r)+\hbar^{2}L(L+1)/2\mu r^{2}, $

      (32)

      where, of course, $ V_{{\mathrm{coul}}}(r)=0 $ for the neutron, μ is the reduced mass, and r is the fragment-target distance. For low L values, the barrier of 10Be is higher than that of the neutron. Thus, for energies close to but above the barrier, ICF of the neutron is larger. In fact, as seen in Figs. 2(a)−(e), it is even larger for more massive targets. Now, as L takes bigger values, the pocket barrier for the neutron (formed by the nuclear and centrifugal potentials) vanishes at a certain energy. Then, neutron ICF starts to decrease. However, for 10Be, the pocket formed by the potentials of Eq. (32) starts to flatten for even larger values of L, where fusion is still possible for high energies. This is why, in Figs. 2(a)−(e), fusion of 10Be increases for higher energies.

      In Figs. 3(a)−(d), comparative calculations are shown for the (a) total ICF, (b) neutron ICF, (c) 10Be ICF at low energies, and (d) 10Be ICF at energies well above the barrier. The comparison is presented for the different targets as a function of $ E\mathrm{_{c.m.}}-V_B $, with $ V_{B} $ being the Coulomb barrier of the projectile 11Be and the target. It can be observed that Figs. 3(a) and 3(b) show similar behavior, i.e., total ICF is basically determined by neutron capture. That is, ICF of the neutron becomes larger for heavier targets. This is because the attraction on the neutron becomes stronger while repulsion on 10Be becomes stronger. With respect to the ICF of 10Be shown in Figs. 3(c)−(d) for low and high energies, this behavior can be understood in the following manner. For low energies and a low incident angular momentum L, 10Be must pass through the barrier of Eq. (32) to fuse with the target. This barrier becomes higher as the target becomes heavier. Table 1 presents the barrier parameters for the various targets as determined by the SPP code [2930]. Therefore, the ICF of 10Be for the heavier targets is lower than that for the lighter ones, as observed in Fig. 3(c). For higher energies above the barrier, the deeper nuclear potential $ V\mathrm{_{nucl}} $ of massive targets causes the barrier pocket produced by the potentials of Eq. (32) to still be present for a large angular momentum L. Therefore, fusion can occur for heavy targets at higher energies. This is not the case for light targets with weaker nuclear potentials, for which the pocket vanishes for lower values of L. The results of these calculations are shown in Fig. 3(d). As a final calculation, Figs. 4(a)−(b) show the energy-dependence contributions to the total CF from the direct and sequential components for 11Be with the light 28Si and heavy 209Bi targets. For both targets, the basic component of the CF is DCF. CF predominantly occurs through the ground state of the projectile; that is, direct fusion. It is very interesting to observe that the SCF of 209Bi becomes increasingly larger than that for 28Si. This could be related to the stronger potentials that produce breakup of the projectile and then result in capture of the fragments.

      Figure 3.  (color online) Comparison of (a) total ICF, (b) neutron ICF, (c) 10Be ICF for low energies, and (d) 10Be ICF for energies well above the barrier. $ V_{B} $ is the barrier energy of the projectile 11Be and the target.

      $ V_{B} $/MeV $ R_{B} $/fm $ {\hbar}\omega_{B} $/MeV
      16O 5.11 8.28 1.50
      28Si 8.58 8.64 1.93
      58Ni 15.89 9.40 2.65
      144Sm 30.79 10.84 3.62
      209Bi 38.55 11.64 4.58

      Table 1.  Coulomb barrier parameters $ V_{B} $, $ R_{B} $, and $ {\hbar}\omega_{B} $ for 10Be with targets 16O, 28Si, 58Ni, 144Sm, and 209Bi.

      Figure 4.  (color online) Total CF, DCF, and SCF for the weakly bound 11Be with targets (a) 28Si and (b) 209Bi.

    IV.   SUMMARY
    • TF, CF, and ICF cross sections have been systematically studied for the neutron-halo weakly bound projectile 11Be in reaction with targets 16O, 28Si, 58Ni, 144Sm, and 209Bi. The calculations were performed with the CDCC model to determine the radial relative motion wave functions of the projectile-target system. These wave functions were then introduced into the angular momentum-dependent model of fusion probabilities to calculate the excitation functions of the different fusion mechanisms. The results of the calculations led to the following conclusions.

      (a) The energy dependence of total ICF becomes increasingly important with respect to CF as the target charge and mass increase.

      (b) The contributions to the ICF of the neutron and 10Be fragments were determined for the different systems. For all targets, the total ICF is mainly due to neutron absorption.

      (c) A comparison of the neutron ICF for the different targets showed that, for incident energies above the barrier, the neutron ICF becomes increasingly relevant for more massive targets. In contrast, for the 10Be fragment, for energies below and around the barrier, fusion becomes smaller for the heavy targets. However, this situation is inverted for high energies, where the ICF of 10Be increases as the target mass increases.

      (d) CF is basically composed of the DCF component at all incident energies. However, for heavier targets, the SCF gradually increases. For light targets, SCF remains very small with respect to DCF.

Reference (30)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return