Semiclassical treatment of fusion and breakup processes of 6 ; 8 He halo nuclei

A semiclassical approach has been used to study the effect of channel coupling on the calculations of the total fusion reaction cross section r fus, and the fusion barrier distribution D fus for the systems 6 He þ 238 U and 8 He þ 197 Au. Since these systems invloves light exotic nuclei, breakup states channel play an important role that should be considered in the calculations. In semiclassical treatment, the relative motion between the projectile and target nuclei is approximated by a classical trajectory while the intrinsic dynamics is handled by time-dependent quantum mechanics. The calculations of the total fusion cross section r fus, and the fusion barrier distribution D fus are compared with the full quantum mechanical calculations using the coupled-channels calculations with all order coupling using the computer code and with the available experimental data.


Introduction
The effects of channel coupling in fusion reactions induced by light weakly bound projectiles have attracted great interest over the last decade [1][2][3]. Some theoretical studies predict strong influence of the breakup channel over the complete fusion (CF) cross section [4][5][6][7][8][9][10]. This has primarily been motivated by the present availability of radioactive ion beams, some of which exhibit unusual features like halo/skin structure and large breakup probabilities. A critical understanding of the fusion mechanism with radioactive ion beams is very significant for the understanding of reactions of astrophysical interest and for the production of new nuclei near the drip lines [11][12][13][14][15][16][17][18][19]. The natural theoretical tool for their description is the coupled-channel method. However, when a weakly bound collision partner is involved, breakup states channel play an important role. Since they are in the continuum, an infinite number of channels should be taken into account. To handle this situation, the continuum can be discretized by the continuum discretized coupled-channel method (CDCC). However, the method is very complicated and requires considerable computer power. For practical purposes, it becomes necessary to approximate the continuum by a finite set of states, as in the CDCC method [20][21][22][23]. This procedure has been extended to the case of fusion reactions in Refs. [8][9][10]. A semiclassical treatment alternative based on the classical trajectory approximation of Alder and Winther (AW) [24] has been proposed by Marta et al. [25].
The aim of the present work focuses on using a semiclassical approach by adopting Alder and Winther (AW) theory originally used to treat the Coulomb excitation of nuclei. The approximation of AW consists of using classical mechanics to describe the projectile-target relative motion while the excitation of the extrinsic degrees of the nuclei are described by quantum mechanics. This semiclassical approximation has been implemented and coded in FORTRAN codename (SCF) [26] which has been used to calculate the total fusion cross section r fus (mb) and the fusion barrier distribution D fus (mb/MeV) for the systems involving light exotic nuclei 6 He þ 238 U and 8 He þ 197 Au. The results of the present study are compared to the full quantum mechanical calculations using the coupled-channel calculations (CC) using the computer code CCFULL [27] and with the available experimental data.

Theoretical framework
Since this work is mainly concerned with the reactions induced by weakly bound projectiles, two variables have been employed, namely r and n for the projectile-target separation vector and the relevant intrinsic degrees of freedom of the projectile, respectively. The internal structure for the target nucleus have been neglected for simplicity. The Hamiltonian then reads where H 0 ðnÞ is the intrinsic Hamiltonian of the projectile and Vðr; nÞ represents the projectile-target interaction. The eigenvectors of H 0 ðnÞ are given by the equation The AW theory [24] is implemented in the following sequence: first, classical mechanics has been used for the time evolution of the variable, r. The trajectory depends on the collision energy, E, the angular momentum, h'. The Rutherford trajectory r ' ðtÞ was used in the original AW method. In our case, the solution of the trajectory is the equations of classical motion with the potential, VðrÞ ¼ hw 0 jVðr; nÞjw 0 i where jw 0 i is a wave function of the ground state of the projectile. The coupling interaction becomes a time-dependent interaction at using the n-space, Vðr ' ðtÞ; nÞ. Second, the dynamics were treated in the intrinsic space as a time-dependent quantum mechanics problem. Expanding the wave function in the basis of intrinsic eigen states, and inserting this expansion into the time-dependent Schrödinger equation for Wðn; tÞ, one obtains the AW coupled equations These equations are solved with the initial conditions A m ð'; t À! À1Þ ¼ d m0 , which means that before the collision ðt À! À1Þ the projectile was in its ground state. The final population of channel m in a collision with angular momentum ' is P ðmÞ ' ¼ jA m ð'; t À! þ1Þj 2 and the angle integrated cross section is: For using this method in the fusion processes, we start with the quantum mechanical calculation of the fusion reaction cross section in a coupled-channel problem. It is assumed that all the channels are bound and have spin zero for simplicity. The fusion reaction cross section is a sum of contributions from each channel. Carrying out partial-wave expansions we get Above, U m' ðk m ; rÞ represents the radial wave function for the ' th -partial-wave in channel m and W F m ðrÞ is the absolute value of the imaginary part of the optical potential associated to fusion reaction in that channel.
To evaluate the fusion cross section using AW method, we had adopted used the approximation: ' is the probability that the system is in channel m at the point of closest approach on the classical trajectory, and T ðmÞ ' ðE m Þ is the probability that a particle with energy E m ¼ E À m and reduced mass M ¼ M P M T =M P þ M T where M P ; M T are, respectively, the masses of the projectile and target, tunnels through the potential barrier in channel m.
We now proceed to study the CF cross sections in reactions induced by weakly bound projectiles. For simplicity, we assume that the ground state (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 m ¼ 0 and m 6 ¼ 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 r CF can be obtained from Eq. (9), dropping contributions from m 6 ¼ 0. That is, where P F ' ðmÞ P Fusion barrier distribution The fusion barrier distribution D fus is extracted from the measured data of r fus , by taking the second derivative, with respect to center of mass energy (E c.m. ), of the product (E c.m. r fus ) [28][29][30] as This quantity is calculated numerically using the threepoint difference formula, where with equal increments of the energy center of mass and the statistical error ðdD fus Þ related with the second derivative at energy E is approximately given by, where the ðdr fus Þ are the errors in the fusion interaction cross sections, and ðdD fus Þ is error of the fusion barrier distribution.

Results and discussion
In this section, the numerical results for the total fusion cross section r fus , and the fusion barrier distribution D fus calculated using the semiclassical approach for the systems 6 He þ 238 U and 8 Heþ 197 Au are explained. Our calculated results of r fus and D fus are compared with the corresponding experimental data and with the full quantum mechanical calculations using the CCFULL code. The Aküz-Winther potential parameters used in the present calculations are listed in Table 1.
In the case of 6 He þ 238 U system, our calculated results represented by dashed blue and black curves for the semiclassical and coupled-channel calculations without including coupling, respectively. The solid blue and black curves are the calculations including the coupling effects for the semiclassical and coupled-channel calculations, for the total fusion reaction cross section r fus , respectively, compared with the respective measured data (solid circles) in Fig.1 panel (a). The experimental data for this system are obtained from Ref. [17]. The coupled-channel calculations using CCFULL performed by considering rotational deformations for both the projectile and target nuclei with deformations parameters b 2 , and b 4 adopted from Moller et al. [31] as listed in Table 2.
The position of the experimental Coulomb barrier V b indicated by an arrow at the E c.m. axis. In the case of nocoupling both semiclassical and quantum mechanical calculations underestimate the experiential data below the Coulomb barrier. The semiclassical calculations are markedly enhanced over the coupled-channel calculations below the Coulomb barrier when we had considered coupling effects. It is clearly seen that, the semiclassical calculations are in better agreement with the measured data compared to the coupled-channel calculations. Panel (b) of Fig. 1 shows the fusion barrier distribution D fus calculations; we could not construct the second derivative accurately using the three-point difference method from the measured data due to few experimental data points for this system. In the corresponding results of the barrier distribution in the lower panel of Fig. 1 we are unable to make clear judgment, if our semiclassical calculations are more satisfactorily fitting the measured data than the coupledchannel calculations.
In a similar analysis, we compare our calculated results of r fus and D fus with the corresponding data in panels (a) and (b) of Fig. 2, respectively, for the 8 He þ 197 Au system. The experimental data of this system are obtained from Ref. [32]. In this case also, we had found that the inclusion of the coupling enhances the semiclassical calculations markedly above and below the Coulomb barrier compared with the measured data, while the CC calculations have very slight improvement below the Coulomb barrier, but still far from the measured data. This system also has few experimental data, that the construction of the second derivative is inaccurate; therefore, we cannot make clear judgment if our semiclassical or coupled-channel calculations are able to fit the measured data for the fusion barrier distribution calculations. Table 1 The parameters of Aküz-Winther potential along

Conclusions
The effect of coupled channel between the elastic channel and the continuum is found to be very essential in the semiclassical calculations which leads to improvement in the total fusion reaction cross section r fus and the fusion barrier distribution D fus around and below the Coulomb barrier and brings the theoretical results closer to the (a) (b) Fig. 1 The comparison of the semiclassical and full quantum mechanical calculations using SCF code (blue curves) and CCFull code (black curves) with the experimental data (black filled circles) [17] for 6 He þ 238 U system. a for the total fusion cross section r fus (mb), and b for the fusion barrier distribution D fus (mb/ MeV), and the arrow represents the position of the Coulomb barrier V b  experimental data. The inclusion of the coupling effects by considering both target and projectile have rotational deformations that enhances the full quantum mechanical calculations around and below the Coulomb barrier, but still have shortfall in reproducing data especially below the Coulomb barrier. The semiclassical approach used in the present work proves to be more adequate to describe light exotic nuclei than the full quantum mechanical ones in comparison with the experiment, because the semiclassical approach has been successfully extended to consider the continuum as a discrete set of channels as in Ref. [8]; therefore, the relative motion between the breakup fragments is more accurately described. The analytical continuation of the time variable is used to extend the trajectories to classically forbidden regions which enhances the calculations below the Coulomb barrier [25]. We could not construct the second derivative to calculate the fusion barrier distribution form the measured data accurately, therefore we could not make clear judgment of whether our semiclassical or coupled-channel calculations agreed with the experimental fusion barrier distribution. This work can be extended to study more systems involving halo nuclei, besides one can test its success or failure in the calculations of fusion reactions involving medium and heavy systems.
Acknowledgments The first author F.A.M. would like to acknowledge the financial support from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (Brazil) and is very grateful 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).
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://crea tivecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.