Fusion Reaction of Weakly Bound Nuclei

Abstract

Semiclassical and full quantum mechanical approaches are used to study the effect of channel coupling on the calculations of the total fusion reaction cross section σfus and the fusion barrier distribution Dfus for the systems 6Li + 64Ni, 11B + 159Tb, and 12C + 9Be. The semiclassical approach used in the present work is based on the method of the Alder and Winther for Coulomb excitation. Full quantum coupled-channel calculations are carried out using CCFULL code with all order coupling in comparison with our semiclassical approach. The semiclassical calculations agree remarkably with the full quantum mechanical calculations. The results obtained from our semiclassical calculations are compared with the available experimental data and with full quantum coupled-channel calculations. The comparison with the experimental data shows that the full quantum coupled channels are better than semiclassical approach in the calculations of the total fusion cross section σfus and the fusion barrier distribution Dfus.

1. Introduction

In recent years, big theoretical and experimental efforts had been dedicated to expertise the effect of breakup of weakly bound nuclei on fusion cross sections [1, 2]. This subject attracts special interests for researchers and scholars, because the fusion of very weakly bound nuclei and exotic radioactive nuclei is reactions that have special interests in astrophysics which play a very vital role in formation of superheavy isotopes for future applications [3, 4, 5, 6, 7, 8]. Since the breakup is very important to be considered in the fusion reaction of weakly bound nuclei, the following should be considered: the elastic breakup (EBU) in which neither of the fragments is captured by the target; incomplete fusion reaction (ICF), which happens when one of the fragments, is captured by the target; and complete fusion following BU (CFBU), which happens in all breakup fragments that are captured by the target, is called the sequential complete fusion (SCF). Therefore, the total breakup cross section is the sum of three contributions: EBU, ICF, and CFBU, whereas the sum of complete fusion (including two body fusions and CFBU) and incomplete fusion is called total fusion (TF) [1, 8, 9, 10]. Fusion reactions with high-intensity stable beams which have a significant breakup probability are good references for testing the models of breakup and fusion currently being developed. The light nuclei such as ⁶Li breakup into ⁴He+²H, with separation energy Sα = 1.48 MeV; 11B breakup into ⁴He+⁷Li with separation energy Sα = 8.664 MeV and ¹²C breakup into three α particles induced by neutrons or protons by ¹²C (p, p′) 3α [3, 11, 12]. The breakup channel is described by the continuum discretized coupled-channel (CDCC) method. The continuum that describes the breakup channel is discretized into bins [13, 14]. To study the coupled-channel problem, this requires a profound truncation of the continuum into discrete bin of energy into equally spaced states. The CDCC method is totally based on this concept. Surrey group extended the discretization procedure discussed in [14] for the deuteron case to study the breakup and fusion reactions of systems involving weakly bound nuclei [15, 16]. Recently, Majeed and Abdul-Hussien [17] utilized the semiclassical approach based on the theory of Alder and Winther. They carried out their calculations to investigate the role of the breakup channel on the fusion cross section ${\sigma }_{\mathit{fus}}$ and fusion barrier distribution $D_{\mathit{fus}}$ for ^6,8H halo [17]. Semiclassical coupled-channel calculations in heavy-ion fusion reactions for the systems ⁴⁰Ar + ¹¹⁰Pd and ¹³²Sn + ⁴⁸Ca were carried out by Majeed et al. [18]. They argued that including the channel coupling between the elastic channel and the continuum enhances the fusion reaction cross section ${\sigma }_{\mathit{fus}}$ and the fusion barrier distribution $D_{\mathit{fus}}$ calculations quite well below and above the Coulomb barrier for medium and heavy systems. This study aims to employ a semiclassical approach by adopting Alder and Winther (AW) [19] theory originally used to treat the Coulomb excitation of nuclei. The semiclassical approach has been implemented and coded using FORTRAN programming codenamed (SCF) which is written and developed by Canto et al. [20]. The fusion cross section ${\sigma }_{\mathit{fus}}$ and fusion barrier distribution $D_{\mathit{fus}}$ are calculated here utilizing the semiclassical approach. The results from the present study are compared with the quantum mechanical calculations using the FORTRAN code (CC) [21] and with the experimental data for the three systems ⁶Li + ⁶⁴Ni, ¹¹B + ¹⁵⁹Tb, and ¹²C + ⁹Be.

2. The semiclassical theory

2.1 The single-channel description

The semiclassical theory is used to estimate the fusion cross section in the one-dimensional potential model which assumes that one can describe the degree of freedom only of the relative motion between the colliding heavy ions [19, 20]. The semiclassical theory deals with the Schrödinger equation assuming independent energy and angular momentum and the potential energy for the radial part of the relative motion through quantum tunneling:

$$\left[-{\mathrm{\hslash }}^2{\nabla }^2/2\mu +V\left(r\right)-E\right]\mathrm{\Psi }\left(\mathrm{r}\right)=0$$E1

where $\mu$ is the reduced mass of the system and $V\left(r\right)$ is the total potential energy of the system. Semiclassical reaction amplitudes can be evaluated as a function of time, assuming the particle trajectory to be determined by classical dynamics, including Coulomb and real nuclear and centrifugal potentials, which can be written in the form

$$V\left(r\right)=V_C\left(r\right)+V_N\left(r\right)+V_l\left(r\right)$$E2

In coupled-channel effects on the elastic channel, the imaginary part should be added to the nuclear potential, represented by complex potential as

$$V_N\left(r\right)=U_N\left(r\right)-\mathit{iW}\left(r\right)$$E3

The method can be extended to describe interference of different $l$ waves due to strong nuclear attraction and absorption caused by the imaginary nuclear potential [6, 20, 22]. When the two nuclei come across the potential barrier into the inner region, the fusion occurs according to the semiclassical theory, and the penetrability probability below barrier can be evaluated using WKB approximation [5, 6, 19, 23, 24]:

$$P_{\mathit{fus}}^{\mathit{WKB}}\left(l,E\right)=\frac{1}{1+\mathit{exp}\left[2{\int }_{r_b^{\left(l\right)}}^{r_a^{\left(l\right)}}{\kappa }_l\left(r\right)\mathit{dr}\right]}$$E4

Then, the latter can be rewritten as follows:

$$P_{\mathit{fus}}^{\mathit{WKB}}\left(l,E\right)=\frac{1}{1+\mathit{exp}\left[\frac{2\pi }{\mathrm{\hslash }{\mathrm{\Omega }}_l}\left(V_b\left(l\right)-E\right)\right]}$$E5

where ${\kappa }_l\left(r\right)$ is the local wave number and $r_b^{\left(l\right)}$ and $r_a^{\left(l\right)}$ are the inner and outer classical turning points at the fusion barrier potential. If one approximates the fusion barrier by a parabolic function, then the penetrability probability above barrier is given by the Hill-Wheeler formula [3, 19]:

$$P_{\mathit{fus}}^{\mathit{WH}}\left(l,E\right)=\frac{1}{1+\mathit{exp}\left[\frac{2\pi }{\mathrm{\hslash }{\mathrm{\Omega }}_l}\left(E-V_b\left(l\right)\right)\right]}$$E6

where $V_b\left(l\right)$ and ${\mathrm{\Omega }}_l$ are the height and the curvature parameter of the fusion barrier for the partial wave, respectively, and $E$ is the bombarding energy. Ignoring the $l$ dependence of $\omega$ and of the barrier position $R_b$ and assuming that the $l$ dependence of $V_b\left(l\right)$ is given only by the difference of the centrifugal potential energy, one can obtain Wong’s formula which is given in Section 5. The fusion cross sections can be estimated by the one-dimensional WKB approximation by the following relations [21, 24]:

$${\sigma }_{\mathit{fus}}\left(E\right)=\frac{\pi }{{\kappa }^2}\sum \left(2l+1\right)P_{\mathit{fus}}^{\mathit{WKB}}\left(l,E\right)$$E7

$$P_{\mathit{fus}}^{\gamma }\left(l,E\right)=\frac{4k}{E}\int \mathit{dr}{\left|u_{\mathit{\gamma l}}\left(k_{\gamma },r\right)\right|}^2{W}_{\mathit{fus}}^{\gamma }\left(r\right)$$E8

where $u_{\mathit{\gamma l}}\left(k_{\gamma },r\right)$ represents the radial wave function for the partial wave $l$ in channel $\gamma$ and $W_{\mathit{fus}}^{\gamma }\left(r\right)$ is the absolute value of the imaginary part of the potential associated to fusion in that channel.

The complete fusion cross section in heavy ions evaluated using semiclassical theory is based on the classical trajectory approximation $r$. And the relevant intrinsic degrees of freedom of the projectile, represented by ξ with applying the continuum discretized coupled-channel (CDCC) approximation of Alder and Winther (AW) theory [16], have been proposed [17]. The projectile Hamiltonian is then given by

$$h=h_0\left(\xi \right)+V\left(\xi ,r\right)$$E9

where ${\mathrm{h}}_0\left(\mathrm{\xi }\right)$ is the intrinsic Hamiltonian and $\mathrm{V}\left(\mathrm{\xi },\mathrm{r}\right)={\mathrm{V}}_N\left(\mathrm{\xi },\mathrm{r}\right)+{\mathrm{V}}_C\left(\mathrm{\xi },\mathrm{r}\right)$, $\mathrm{V}\left(\mathrm{\xi },\mathrm{r}\right)$ is the interaction between the projectile and target nuclei. The Rutherford trajectory depends on the collision energy, $E$, and the angular momentum, $l$. In this case the trajectory is the solution of the classical motion equations with the potential $V\left(r\right)=⟨{\mathrm{\Psi }}_0|V\left(r,\xi \right)|{\mathrm{\Psi }}_0⟩$, where ${\mathrm{\Psi }}_0$ is the ground state (g.s.) of the projectile. In this way, the interaction becomes time-dependent in the ξ-space $V_l\left(\xi .t\right)=V\left(r_{l\left(t\right)},\xi \right)$, and the eigenstates of the intrinsic Hamiltonian $|{\mathrm{\Psi }}_{\gamma }⟩$ satisfy the Schrödinger equation [25, 26]:

$$\left.h\right|{\mathrm{\Psi }}_{\gamma }⟩=\left.\epsilon \right|{\mathrm{\Psi }}_{\gamma }⟩$$E10

After expanding the wave function in the basis of intrinsic eigenstates

$$\mathrm{\Psi }\left(\xi ,t\right)=\sum a_{\gamma }\left(l,t\right){\mathrm{\Psi }}_{\gamma }\left(\xi \right)e^{-i{\epsilon }_{\gamma }t/\mathrm{\hslash }}$$E11

and inserting Eq. (11) into the Schrödinger equation for $\mathrm{\Psi }\left(\xi ,t\right)$, the AW equations can be obtained:

$$i\mathrm{\hslash }{\dot{a}}_{\gamma }\left(l,t\right)=\sum _{ϵ}⟨{\mathrm{\Psi }}_{\gamma }|V\left(\xi ,t\right)|{\mathrm{\Psi }}_{\gamma }⟩e^{\raisebox{1ex}{$i\left({\epsilon }_{\gamma }-{\epsilon }_{ϵ}\right)t$}\!\left/ \!\raisebox{-1ex}{$\mathrm{\hslash }$}\right.}\mathit{\gamma ϵ}\left(l,t\right)$$E12

These equations should be solved with initial conditions $a_{\gamma }\left(l,t\to -\infty \right)={\delta }_{\gamma 0}$ which mean that before the collision $\left(t\to -\infty \right)$, the projectile was in its ground state. The final population of the channel $\gamma$ in a collision with angular momentum $l$ is $P_{\mathit{fus}}^{\gamma }\left(l,E\right)={\left|a_{\gamma }\left(l,t\to -\infty \right)\right|}^2$. Eq. (8) gives the general expression for the fusion cross section in multichannel scattering [26].

2.2 The coupled-channel description

The variables employed to describe the collision are the projectile-target separation vector $\overrightarrow{r}$ and the relevant intrinsic degrees of freedom of the projectile $\xi$. For simplicity, we neglect the internal structure of the target. The Hamiltonian is then given by [27]

$$H=H_0\left(\xi \right)+V\left(\overrightarrow{r},\xi \right)$$E13

where $H_0\left(\xi \right)$ is the intrinsic Hamiltonian of the projectile and $V(\overrightarrow{r},\xi )$ represents the projectile-target interaction. Since the main purpose of the present work is to test the semiclassical model in calculations of sub-barrier fusion, the nuclear coupling is neglected. Furthermore, for the present theory-theory comparison, only the Coulomb dipole term is taken into account. The eigenvectors of $H_0\left(\xi \right)$ are given by the equation [27]

$$H_0|{\phi }_{\beta }⟩={\epsilon }_{\beta }|{\phi }_{\beta }⟩$$E14

where ${\epsilon }_{\beta }$ is energy of internal motion. The AW method is implemented in two steps. First, one employs classical mechanics for the time evolution of the variable $\overrightarrow{r}$. The ensuing trajectory depends on the collision energy, $E$, and the angular momentum,$\mathcal{\hslash l}$. In its original version, an energy symmetrized Rutherford trajectory ${\overrightarrow{r}}_{\ell }\left(t\right)$ was used. In our case, the trajectory is the solution of the classical equations of motion with the potential $V(\overrightarrow{r})=⟨{\phi }_0|V(\overrightarrow{r},\xi ){\phi }_0⟩$, where $|{\phi }_0⟩$ is the ground state of the projectile. In this way, the coupling interaction becomes a time-dependent interaction in the $\xi$-space, $V_{\ell }\left(\xi ,t\right)\equiv V({\overrightarrow{r}}_{\ell }\left(t\right),\xi )$. The second step consists in treating the dynamics in the intrinsic space as a time-dependent quantum mechanics problem. Expanding the wave function in the basis of intrinsic eigenstates [19]

$$\psi \left(\xi ,t\right)=\sum _{\beta }{a}_{\beta }\left(\ell ,t\right){\phi }_{\beta }\left(\xi \right)e^{-i{\epsilon }_{\beta }t/\mathrm{\hslash }}$$E15

and inserting this expansion into the Schrödinger equation for $\psi \left(\xi ,t\right)$, one obtains the AW equations [27]

$$i\mathrm{\hslash }{\dot{a}}_{\beta }\left(\ell ,t\right)=\sum _{\alpha }{a}_{\gamma }\left(\ell ,t\right)⟨{\phi }_{\beta }|V_{\ell }\left(\xi ,t\right)|{\phi }_{\beta }⟩e^{-i\left({\epsilon }_{\beta }-{\epsilon }_{\gamma }\right)t/\mathrm{\hslash }}$$E16

These equations are solved with the initial conditions $a_{\beta }\left(\ell ,t\to -\infty \right)={\delta }_{\beta 0}$, which means that before the collision $\left(t\to -\infty \right)$, the projectile was in its ground state. The final population of channel $\beta$ in a collision with angular momentum $\ell$ is $P_{\ell }^{\left(\beta \right)}={\left|a_{\beta }\left(\ell ,t\to +\infty \right)\right|}^2$ and the angle-integrated cross section is [19]

$${\sigma }_{\beta }=\frac{\pi }{k^2}\sum _{\ell }\left(2\ell +1\right)P_{\ell }^{\left(\beta \right)}$$E17

To extend this method to fusion reactions, we start with the quantum mechanical calculation of the fusion cross section in a coupled-channel problem. For simplicity, we assume that all channels are bound and have zero spin. The fusion cross section is a sum of contributions from each channel. Carrying out partial-wave expansions, we get [28]

$${\sigma }_F=\sum _{\beta }\left[\frac{\pi }{k^2}\sum _{\ell }\left(2\ell +1\right)P_{\ell }^F\left(\beta \right)\right]$$E18

with

$$P_{\ell }^F\left(\beta \right)=\frac{4k}{E}\int W_{\beta }^F\left(r\right){\left|u_{\beta \ell }\left(k_{\beta },r\right)\right|}^2\mathit{dr}$$E19

Above, $u_{\beta \ell }\left(k_{\beta },r\right)$ represents the radial wave function for the $\ell$th partial wave in channel $\beta$, and $W_{\beta }^F$ is the absolute value of the imaginary part of the optical potential associated to fusion in that channel.

To use the AW method to evaluate the fusion cross section, we make the approximation [27]

$$P_{\ell }^F\left(\beta \right)\simeq {\overline{P}}_{\ell }^{\left(\beta \right)}{T}_{\ell }^{\left(\beta \right)}\left(E_{\beta }\right)$$E20

where ${\overline{P}}_{\ell }^{\left(\beta \right)}$ is the probability that the system is in channel $\beta$ at the point of closest approach on the classical trajectory and $T_{\ell }^{\left(\beta \right)}\left(E_{\beta }\right)$ is the probability that a particle with energy $E_{\beta }=E-{\epsilon }_{\beta }$ and reduced mass $\mu =M_P{M}_T/\left(M_P+M_T\right)$, where $M_P,M_T$ are the masses of the projectile and target, respectively, tunnels through the potential barrier in channel $\beta$ [19].

We now proceed to study the CF cross sections in reactions induced by weakly bound projectiles. For simplicity, we assume that the g.s. is the only bound state of the projectile and that the breakup process produces only two fragments, $F_1$ and $F_2$. In this way, the labels $\beta =0$ and $\beta \ne 0$ correspond, respectively, to the g.s. and the breakup states represented by two unbound fragments. Neglecting any sequential contribution, the CF can only arise from the elastic channel. In this way, the cross section ${\sigma }_{\mathit{CF}}$ can be obtained from Eq. (20), dropping contributions from $\beta \ne 0$. That is [19],

$${\sigma }_{\mathit{CF}}=\frac{\pi }{k^2}\sum _{\ell }\left(2\ell +1\right)P_{\ell }^{\mathit{\text{Surv}}}{T}_{\ell }^{\left(0\right)}\left(E\right)$$E21

where

$$P_{\ell }^{\mathit{\text{Surv}}}\equiv {\overline{P}}_{\ell }^{\left(0\right)}={\left|a_0\left(\ell ,t_{\mathit{ca}}\right)\right|}^2$$E22

is usually called survival (to breakup) probability [19].

3. Fusion barrier distribution

Nuclear fusion is related to the transmission of the incident wave through a potential barrier, resulting from nuclear attraction plus Coulomb repulsion. However, the meaning of the fusion barrier depends on the description of the collision. Coupled-channel calculations include static barriers, corresponding to frozen densities of the projectile and the target. Its most dramatic consequence is the enhancement of the total fusion reaction cross section ${\mathrm{\sigma }}_{\mathit{fus}}$ at Coulomb barrier energies $V_b$, in some cases by several orders of magnitude. One of the possible ways to describe the effect of coupling channels is a division of the fusion barrier into several, the so-called fusion barrier distribution $D_{\mathit{fus}}$ and given by [20, 29]

$$D_{\mathit{fus}}\left(E\right)=\frac{d^{2}F\left(E\right)}{d{E}^2}$$E23

where $F\left(E\right)$ is related to the total fusion reaction cross section through [29]

$$F\left(E\right)=E{\sigma }_{\mathit{fus}}\left(E\right)$$E24

The experimental determination of the fusion reaction barrier distribution has led to significant progress in the understanding of fusion reaction. This comes about because, as already mentioned, the fusion reaction barrier distribution gives information on the coupling channels appearing in the collision. However, we note from Eq. (24) that, since fusion reaction barrier distribution should be extracted from the values of the total fusion reaction cross section, it is the subject to experimental as well as numerical uncertainties [29, 30, 31]:

$$D_{\mathit{fus}}\left(E\right)\approx \frac{F\left(E+∆E\right)+F\left(E-∆E\right)-2F\left(E\right)}{∆E^2}$$E25

where $∆E$ is the energy interval between measurements of the total fusion reaction cross section. From Eq. (25), one finds that the statistical error associated with the fusion reaction barrier distribution is approximately given by [30]

$${\mathit{\delta D}}_{\mathit{fus}}^{\mathit{\text{stat}}}\left(E\right)\approx \frac{\sqrt{{\left[\mathit{\delta F}\left(E+∆E\right)\right]}^2+{\left[\mathit{\delta F}\left(E-∆E\right)\right]}^2+4{\left[\mathit{\delta F}\left(E\right)\right]}^2}}{{\left(∆E\right)}^2}$$E26

where $\mathit{\delta F}\left(E\right)$ means the uncertainty in the measurement of the product of the energy by the total fusion reaction cross section at a given value of the collision energy. Then the uncertainties can approximately be written as [30]

$${\mathit{\delta D}}_{\mathit{fus}}^{\mathit{\text{stat}}}\left(E\right)\approx \frac{\sqrt{6}\mathit{\delta F}\left(E\right)}{{\left(∆E\right)}^2}$$E27

4. Results and discussion

In this section, the theoretical calculations are obtained for total fusion reaction ${\sigma }_{\mathit{fus}},$ and the fusion barrier distribution $D_{\mathit{fus}}$ using the semiclassical theory adopted the continuum discretized coupled channel (CDCC) to describe the effect of the breakup channel on fusion processes. The semiclassical calculations are carried out using the (SCF) code, while the full quantum mechanical calculations are performed by using the code (CC) for the systems ⁶Li + ⁶⁴Ni, ¹¹B + ¹⁵⁹Tb and ¹²C + ⁹Be. The values of the height $V_c$, radius $R_c$, and curvature $\mathrm{\hslash }\omega$ for the fusion barrier are displayed in Table 1.

$\text{System}$	$V_c$	$R_c$	$\mathrm{\hslash }\omega$	$\text{Refs}$
^6Li + ^64Ni	$12.41$	$9.1$	$3.9$	[32]
^11B + ^159Tb	$40.34$	$10.89$	$4.42$	[33]
^12C + ^9Be	4.28	7.43	2.61	[34]

Table 1.

The fusion barrier parameters are height $V_c\left(\mathit{MeV}\right)$, radius $R_c\left( ,\mathit{fm}\right)$, and curvature $\mathrm{\hslash }\omega \left(\mathit{MeV}\right).$

4.1 The reaction ⁶Li + ⁶⁴Ni

The calculations of the fusion cross section ${\sigma }_{\mathit{fus}}$ and fusion barrier distribution $D_{\mathit{fus}}$ are presented in Figure 1 panel (a) and panel (b), respectively, for the system ⁶Li + ⁶⁴Ni. The dashed blue and red curves represent the semiclassical and full quantum mechanical calculations without coupling, respectively. The solid blue and red curves are the calculations including the coupling effects for the semiclassical and full quantum mechanical calculations, respectively. Panel (a) shows the comparison between our semiclassical and full quantum mechanical calculations with the respective experimental data (solid circles).

Figure 1.

The comparison of the coupled-channel calculations of semiclassical treatment (red curves) and full quantum mechanical (blue curves) with the experimental data of complete fusion (black-filled circles) [32] for ⁶Li + ⁶⁴Ni system. Panel (a) refers to the total fusion reaction cross section ${\sigma }_{\mathit{fus}}$ (mb), and panel (b) provides the fusion reaction barrier distribution $D_{\mathit{fus}}$ (mb/MeV).

The experimental data for this system are obtained from Ref. [32]. The real and imaginary Akyüz-Winther potential parameters obtained by using chi-square method are the strength $W_0=50\mathit{MeV}$, radius $r_i=1.0\mathit{fm},$ and diffuseness $a_i=0.25\mathit{fm}$, and for the real part, the depth is $V_0=35.0\mathit{MeV}$, radius is $r_0=1.1\mathit{fm}$, and diffuseness is $a_0=0.8\mathit{fm}$. The ${\chi }^2$ values obtained for the total fusion cross section ${\sigma }_{\mathit{fus}}$ are 1.5057 and 1.1286 in the case of no coupling for semiclassical and quantum mechanical calculations, respectively. The ${\chi }^2$ values obtained for the case of coupling effects included are 0.2431 and 0.3115 for semiclassical and quantum mechanical calculations, respectively. The ${\chi }^2$ values show clearly that semiclassical calculations including coupling effects are more consistent with the experimental data than full quantum mechanical including coupling effects. The ${\chi }^2$ values obtained using single-channel calculations for the fusion reaction barrier distribution $D_{\mathit{fus}}$ are 0.1823 and 1.1914 for semiclassical and quantum mechanical calculations, respectively. The ${\chi }^2$ values obtained when coupled channels are included are 0.1827 and 0.1321 for semiclassical and quantum mechanical calculations, respectively; the fusion barrier distribution $D_{\mathit{fus}}$ has been extracted from the experimental data using Wong fit model along with the three-point difference method. The comparison with the experimental data for $D_{\mathit{fus}}$ shows that the quantum mechanical calculations are in better agreement than the semiclassical calculations including the coupling effects.

4.2 The reaction ¹¹B + ¹⁵⁹Tb

In similar analysis we compare our theoretical calculations of the fusion cross section ${\sigma }_{\mathit{fus}}$ and fusion barrier distribution $D_{\mathit{fus}}$ with the corresponding experimental data in panels (a) and (b) of Figure 2, respectively, for the system ¹¹B + ¹⁵⁹Tb. The experimental data for this system is obtained from Ref. [33]. The real and imaginary Akyüz-Winther potential parameters are obtained by using chi-square method: ${\mathrm{V}}_0=126.1\mathrm{MeV}$, ${\mathrm{r}}_0=1.2\mathrm{fm},\text{and}{\mathrm{a}}_0=0.5\mathrm{fm}$ and ${\mathrm{W}}_0=55.9\mathrm{MeV}$, ${\mathrm{r}}_{\mathrm{i}}=0.986\mathrm{fm}$, and ${\mathrm{a}}_{\mathrm{i}}=0.614\mathrm{fm}$. The ${\mathrm{\chi }}^2$ values 0.9473 and 0.2486 are obtained for ${\mathrm{\sigma }}_{\mathrm{fus}}$ using semiclassical and quantum mechanical distribution calculations without including the coupling, respectively, while the ${\chi }^2$ for ${\sigma }_{\mathit{fus}}$ using the semiclassical and quantum mechanical distribution calculations including the coupling effects are 0.2681 and 0.1657, respectively. The ${\chi }^2$ for the fusion barrier distribution $D_{\mathit{fus}}$ using the semiclassical and quantum calculations are 0.5828 and 1.2329 for no coupling and 4.5969 and 0.0616 including coupling effects, respectively. The ${\chi }^2$ values for ${\sigma }_{\mathit{fus}}$ and $D_{\mathit{fus}}$ give clear evidence that the quantum mechanical calculations are in better agreement than the semiclassical calculations as compared with experimental data.

Figure 2.

The comparison of the coupled-channel calculations of semiclassical treatment (red curves) and full quantum mechanical (blue curves) with the experimental data of complete fusion (black-filled circles) [33] for ¹¹B + ¹⁵⁹Tb system. Panel (a) shows the total fusion reaction cross section ${\sigma }_{\mathit{fus}}$ (mb), and panel (b) gives the fusion reaction barrier distribution $D_{\mathit{fus}}$ (mb/MeV).

4.3 The reaction ¹²C + ⁹Be

Figure 3 (panels (a) and (b)) presents the comparison between our theoretical calculations for ${\sigma }_{\mathit{fus}}$ and $D_{\mathit{fus}}$ using both semiclassical and quantum mechanical calculations with the corresponding experimental data for the system ¹²C + ⁹Be. The experimental data for this system are obtained from Ref. [34]. The real and imaginary Akyüz-Winther potential parameters are obtained by using chi-square method: $V_0=40.3\mathit{MeV}$, $r_0=1.11\mathit{fm},a_0=0.590\mathit{fm}$, $W_0=0\mathit{MeV}$, $r_i=1.1\mathit{fm},$ and $a_i=0.50\mathit{fm}$. The ${\chi }^2$ values obtained from the comparison between the results and experimental data for ${\sigma }_{\mathit{fus}}$ are 1.0633 and 1.1447 without coupling, and 0.4924 and 0.2072 with coupling, for semiclassical and quantum mechanical calculations, respectively. The obtained ${\chi }^2$ values for $D_{\mathit{fus}}$ using semiclassical and quantum mechanical calculations are 1.2383 and 0.6185 without coupling and 0.9875 and 0.1868 with account for coupling, respectively.

Figure 3.

The comparison of the coupled-channel calculations of semiclassical treatment (red curves) and full quantum mechanical (blue curves) with the experimental data of complete fusion (black-filled circles) [34] for ¹²C + ⁹Be system. Panel (a) shows the total fusion reaction cross section ${\sigma }_{\mathit{fus}}$ (mb), and panel (b) gives the fusion reaction barrier distribution $D_{\mathit{fus}}$ (mb/MeV).

5. Conclusion

The semiclassical and quantum mechanical calculations for the total fusion reaction ${\sigma }_{\mathit{fus}}$ and the fusion barrier distribution $D_{\mathit{fus}}$ calculations below and around Coulomb barrier were discussed for the systems ⁶Li + ⁶⁴Ni, ¹¹B + ¹⁵⁹Tb, and ¹²C + ⁹Be. We conclude that the breakup channel is very important to be taken into consideration to describe the total fusion reaction ${\sigma }_{\mathit{fus}}$ and the fusion barrier distribution $D_{\mathit{fus}}$ for reaction of light projectiles. The full quantum mechanical calculations are closer to the experimental data than the semiclassical calculations; however, semiclassical ones can be considered a successful tool for studying fusion reaction of systems involving light projectiles.

Acknowledgments

The author FA Majeed gratefully acknowledges financial assistance from Conselho Nacional de Desenvolvimento Científico e Tecnolόgico (CNPq) (Brazil) and is especially indebted to the World Academy of Sciences for the advancement of science in developing countries (TWAS) (Italy) for a 1-year grant under the scheme (TWAS-CNPq exchange programs for postdoctoral researchers).