The fusion barrier parameters are height , radius , and curvature
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.
- fusion reaction
- breakup channel
- weakly bound nuclei
- fusion barrier
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 6Li breakup into 4He+2H, with separation energy Sα = 1.48 MeV; 11B breakup into 4He+7Li with separation energy Sα = 8.664 MeV and 12C breakup into three α particles induced by neutrons or protons by 12C (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  for the deuteron case to study the breakup and fusion reactions of systems involving weakly bound nuclei [15, 16]. Recently, Majeed and Abdul-Hussien  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 and fusion barrier distribution for 6,8H halo . Semiclassical coupled-channel calculations in heavy-ion fusion reactions for the systems 40Ar + 110Pd and 132Sn + 48Ca were carried out by Majeed et al. . They argued that including the channel coupling between the elastic channel and the continuum enhances the fusion reaction cross section and the fusion barrier distribution 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)  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. . The fusion cross section and fusion barrier distribution 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)  and with the experimental data for the three systems 6Li + 64Ni, 11B + 159Tb, and 12C + 9Be.
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:
where is the reduced mass of the system and 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
In coupled-channel effects on the elastic channel, the imaginary part should be added to the nuclear potential, represented by complex potential as
The method can be extended to describe interference of different 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]:
Then, the latter can be rewritten as follows:
where is the local wave number and and 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]:
where and are the height and the curvature parameter of the fusion barrier for the partial wave, respectively, and is the bombarding energy. Ignoring the dependence of and of the barrier position and assuming that the dependence of 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]:
where represents the radial wave function for the partial wave in channel and 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 . 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 , have been proposed . The projectile Hamiltonian is then given by
where is the intrinsic Hamiltonian and , is the interaction between the projectile and target nuclei. The Rutherford trajectory depends on the collision energy, , and the angular momentum, . In this case the trajectory is the solution of the classical motion equations with the potential , where is the ground state (g.s.) of the projectile. In this way, the interaction becomes time-dependent in the ξ-space , and the eigenstates of the intrinsic Hamiltonian satisfy the Schrödinger equation [25, 26]:
After expanding the wave function in the basis of intrinsic eigenstates
and inserting Eq. (11) into the Schrödinger equation for , the AW equations can be obtained:
These equations should be solved with initial conditions which mean that before the collision , the projectile was in its ground state. The final population of the channel in a collision with angular momentum is . Eq. (8) gives the general expression for the fusion cross section in multichannel scattering .
2.2 The coupled-channel description
The variables employed to describe the collision are the projectile-target separation vector and the relevant intrinsic degrees of freedom of the projectile . For simplicity, we neglect the internal structure of the target. The Hamiltonian is then given by 
where is the intrinsic Hamiltonian of the projectile and 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 are given by the equation 
where 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 . The ensuing trajectory depends on the collision energy, , and the angular momentum,. In its original version, an energy symmetrized Rutherford trajectory was used. In our case, the trajectory is the solution of the classical equations of motion with the potential , where is the ground state of the projectile. In this way, the coupling interaction becomes a time-dependent interaction in the -space, . 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 
and inserting this expansion into the Schrödinger equation for , one obtains the AW equations 
These equations are solved with the initial conditions , which means that before the collision , the projectile was in its ground state. The final population of channel in a collision with angular momentum is and the angle-integrated cross section is 
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 
Above, represents the radial wave function for the th partial wave in channel , and 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 
where is the probability that the system is in channel at the point of closest approach on the classical trajectory and is the probability that a particle with energy and reduced mass , where are the masses of the projectile and target, respectively, tunnels through the potential barrier in channel .
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, and . In this way, the labels and 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 can be obtained from Eq. (20), dropping contributions from . That is ,
is usually called survival (to breakup) probability .
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 at Coulomb barrier energies , 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 and given by [20, 29]
where is related to the total fusion reaction cross section through 
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]:
where 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 
where 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 
4. Results and discussion
In this section, the theoretical calculations are obtained for total fusion reaction and the fusion barrier distribution 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 6Li + 64Ni, 11B + 159Tb and 12C + 9Be. The values of the height , radius , and curvature for the fusion barrier are displayed in Table 1.
4.1 The reaction 6Li + 64Ni
The calculations of the fusion cross section and fusion barrier distribution are presented in Figure 1 panel (a) and panel (b), respectively, for the system 6Li + 64Ni. 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).
The experimental data for this system are obtained from Ref. . The real and imaginary Akyüz-Winther potential parameters obtained by using chi-square method are the strength , radius and diffuseness , and for the real part, the depth is , radius is , and diffuseness is . The values obtained for the total fusion cross section are 1.5057 and 1.1286 in the case of no coupling for semiclassical and quantum mechanical calculations, respectively. The values obtained for the case of coupling effects included are 0.2431 and 0.3115 for semiclassical and quantum mechanical calculations, respectively. The values show clearly that semiclassical calculations including coupling effects are more consistent with the experimental data than full quantum mechanical including coupling effects. The values obtained using single-channel calculations for the fusion reaction barrier distribution are 0.1823 and 1.1914 for semiclassical and quantum mechanical calculations, respectively. The values obtained when coupled channels are included are 0.1827 and 0.1321 for semiclassical and quantum mechanical calculations, respectively; the fusion barrier distribution 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 shows that the quantum mechanical calculations are in better agreement than the semiclassical calculations including the coupling effects.
4.2 The reaction 11B + 159Tb
In similar analysis we compare our theoretical calculations of the fusion cross section and fusion barrier distribution with the corresponding experimental data in panels (a) and (b) of Figure 2, respectively, for the system 11B + 159Tb. The experimental data for this system is obtained from Ref. . The real and imaginary Akyüz-Winther potential parameters are obtained by using chi-square method: , and , , and . The values 0.9473 and 0.2486 are obtained for using semiclassical and quantum mechanical distribution calculations without including the coupling, respectively, while the for using the semiclassical and quantum mechanical distribution calculations including the coupling effects are 0.2681 and 0.1657, respectively. The for the fusion barrier distribution 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 values for and give clear evidence that the quantum mechanical calculations are in better agreement than the semiclassical calculations as compared with experimental data.
4.3 The reaction 12C + 9Be
Figure 3 (panels (a) and (b)) presents the comparison between our theoretical calculations for and using both semiclassical and quantum mechanical calculations with the corresponding experimental data for the system 12C + 9Be. The experimental data for this system are obtained from Ref. . The real and imaginary Akyüz-Winther potential parameters are obtained by using chi-square method: , , , and . The values obtained from the comparison between the results and experimental data for 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 values for 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.
The semiclassical and quantum mechanical calculations for the total fusion reaction and the fusion barrier distribution calculations below and around Coulomb barrier were discussed for the systems 6Li + 64Ni, 11B + 159Tb, and 12C + 9Be. We conclude that the breakup channel is very important to be taken into consideration to describe the total fusion reaction and the fusion barrier distribution 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.
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).
Gomes PRS, Padron I, Crema E, Capurro OA, Fernández Niello JO, et al. Comprehensive study of reaction mechanisms for the 9Be+144Sm system at near- and sub-barrier energies. Physical Review C. 2006; 73:064606
Gomes PRS, Rios JL, Borges JR, Otomar DR. Fusion, breakup and scattering of weakly bound nuclei at near barrier energies. The Open Nuclear & Particle Physics Journal. 2013; 6:10-15
Dasgupta M, Gomes PRS, Hinde DJ, Moraes SB, Anjos RM, Berriman AC, et al. Effect of breakup on the fusion of 6Li, 7Li, and 9Be with heavy nuclei. Physical Review C. 2004; 70:024606
Takigawa N, Sagawa H. Interaction potential and fusion of a halo nucleus. Physics Letters B. 1991; 265:23-28
Hussein MS, Pato MP, Canto LF, Donangelo R. Near-barrier fusion of 11Li with heavy spherical and deformed targets. Physical Review C. 1992; 46:377-379
Dasso CH, Vitturi A. Does the presence of 11Li breakup channels reduce the cross section for fusion processes? Physical Review C. 1994; 50:R12-R14
Canto LF, Donangelo R, Lotti P. Effect of Coulomb dipole polarizability of halo nuclei on their near-barrier fusion with heavy targets. Physical Review C. 1995; 52:R2848-R2850
Martí GV, Gomes PRS, Rodríguez MD, Fernández Niello JO. Fusion, reaction, and breakup cross sections of 9Be on a light mass target. Physical Review C. 2005; 71:027602
Padron I, Gomes PRS, Anjos RM, Lubian J, Muri C, Alves JJ. Fusion of stable weakly bound nuclei with 27Al and 64Zn. Physical Review C. 2002; 66:044608
Wang B, Zhao W, Gomes PRS, Zhao E, Zhou S. Systematic study of breakup effects on complete fusion at energies above the Coulomb barrier. Physical Review C. 2014; 90:034612
Wolfs FLH, White CA, Bryan DC, Freeman CG, Herrick DM, et al. Breakup of 82 MeV 11B. Physical Review C. 1994; 49:2538-2548
Bogatin VI, Novak Z, Ostroumov VI. Breakup of 12C into three alpha particles accompanying inelastic scattering of 80-MeV mesons. Soviet Physics - JETP. 1963; 16:1116-1121
Matsumoto T, Hiyama E, Ogata K, Iseri Y, Kamimura M, Chiba S, et al. Continuum-discretized coupled-channels method for four-body nuclear breakup in 6He+12C scattering. Physical Review C. 2004; 0616061(R):70
Austern N, Iseri Y, Kamimura M, Kawai M, Rawitscher G, Yahiro M. Continuum-discretized coupled-channels calculations for three-body models of deuteron-nucleus reactions. Physics Reports. 1987; 154:125-204
Nunes FM, Thompson IJ. Nuclear interference effects in 8B sub-Coulomb breakup. Physical Review C. 1998; 57:R2818-R2820
Tostevin JA, Nunes FM, Thompson IJ. Calculations of three-body observable in 8B breakup. Physical Review C. 2001; 63:024617
Majeed FA, Abdul-Hussien YA. Semiclassical treatment of fusion and breakup processes of 6,8He halo nuclei. Journal of Theoretical and Applied Physics. 2016; 107:1-6
Majeed FA, Hamodi RSh, Hussian FM. Journal of Computational and Theoretical Nanoscience. 2017; 14:2242-2247
Alder K, Winther A. Electromagnetic Excitations. Amsterdam: North-Holland; 1975
Canto LF, Donangelo R, Marta HD. Physical Review C. 2006; 73:034608 SCF code unpublished
Canto LF, Gomes PRS, Donangelo R, Hussein MS. Fusion and breakup of weakly bound nuclei. Physics Reports. 2006; 424:1-111
Gomes PRS, Rodríguez MD, Martí GV, Padron I, Chamon LC. Effect of the breakup on the fusion and elastic scattering of weakly bound projectiles on 64Zn. Physical Review C. 2005; 71:034608
Canto LF, Gomes PRS, Donangelo R, Hussein MS. Recent developments in fusion and direct reactions with weakly bound nuclei. Physics Reports. 2015; 596:1-86
Cárdenas WHZ, Canto LF, Donangelo R, Hussein MS, Lubian J, Romanelli A. Approximations in fusion and breakup reactions induced by radioactive beams. Nuclear Physics A. 2002; 703:633-648
Hussein MS, Pato MP, Canto LF, Donangelo R. Real part of the polarization potential for 11Li-induced fusion reactions. Physical Review C. 1993; 47:2398-2400
Ibraheem AA, Bonaccors A. Coulomb breakup effects on the optical potentials of weakly bound nuclei. Nuclear Physics A. 2005; 748:414-432
Liang JF, Signorini C. Fusion induced by radioactive ion beams. International Journal of Modern Physics E: Nuclear Physics. 2005; 14:1121
Abramowitz M, Stegun I. Handbook of Mathematical Functions. New York: Dover; 1972
Rowley N, Satchler GR, Stelson PH. On the “distribution of barriers” interpretation of heavy-ion fusion. Physics Letters B. 1991; 254:25
Bertulani CA, Canto LF. Semiclassical calculation of Coulomb break-up of weakly-bound nuclei. Nuclear Physics A. 1992; 539:163
Dasgupta M, Hinde DJ, Rowley N, Stefanini AM. Measuring barrier to fusion. Annual Review of Nuclear and Particle Science. 1998; 48:401
Shaikh MM, Roy S, Rajbanshi S, Pradhan MK, Mukherjee A, et al. Barrier distribution functions for the system 6Li+64Ni and the effect of channel coupling. Physical Review C. 2015; 91:034615
Mukherjee M, Roy S, Pradhan MK, Sarkar MS, Basu P, et al. Influence of projectile breakup threshold on complete fusion. Physics Letters B. 2006; 636:91-95
Cheung HC, High MD, Čujec B. Fusion and neutron transfer cross sections for 9Be+12C at sub coulomb energies. Nuclear Physics A. 1978; 296:333-348