arxiv: v2 [nucl-th] 9 Feb 2017

Transcription

1 EPJ manuscript No. (wi be inserted by the editor) arxiv: v2 [nuc-th] 9 Feb 2017 Improved WKB approximation for quantum tunneing: Appication to heavy ion fusion A. J. Toubiana 1,2,3, L. F. Canto 4,5, M. S. Hussein 6,7,8 1 Departamento de Engenharia Nucear, Escoa Poitécnica, Universidade Federa do Rio de Janeiro, C.P , , Rio de Janeiro, RJ, Brazi, 2 Écoe CentraeSupec, Pateau du Mouon 3, rue Joiot-Curie GIF-SUR-YVETTE - France, 3 Paris Sacay, Espace Technoogique, Bat. Discovery - RD 128-2e t., Saint-Aubin, France, 4 Instituto de Física, Universidade Federa do Rio de Janeiro, CP 68528, Rio de Janeiro, Brazi, 5 Instituto de Física, Universidade Federa Fuminense, Av. Litoranea s/n, Gragoatá, Niterói, R.J., , Brazi, 6 Instituto de Estudos Avançados, Universidade de São Pauo C. P , São Pauo-SP, Brazi, and Instituto de Física, Universidade de São Pauo, C. P , São Pauo,-SP, Brazi, 7 Departamento de Física Matemática, Instituto de Física, Universidade de São Pauo, C.P , , São Pauo, SP, Brazi, 8 Departamento de Física, Instituto Tecnoógico de Aeronáutica, CTA, São José dos Campos, São Pauo, SP, Brazi 1 Introduction February 10, 2017 Abstract. In this paper we revisit the one-dimensiona tunneing probem. We consider Kembe s approximation for the transmission coefficient. We show how this approximation can be extended to above-barrier energies by performing the anaytica continuation of the radia coordinate to the compex pane. We investigate the vaidity of this approximation by comparing their predictions for the cross section and for the barrier distribution with the corresponding quantum mechanica resuts. We find that the extended Kembe s approximation reproduces the resuts of quantum mechanics with great accuracy. PACS Eq Bc 25.60Gc Quantum tunneing (QT) has been with us for as ong as the Schrödinger equation (SE). It describes a purey quantum phenomenon of a passage through a barrier. One of the first to popuarize the essentias of QT was Gamow in his description of α decay of atomic nucei. For this purpose, Gamow used an approximation of the soution of the SE based on the short wave ength imit. This approximation came to be known as the WKB approximation and it corresponds to using the cassica action as the phase of the wave function. In the deveopment of the tunneing theory, severa authors introduced concepts and recipes in order to get the tunneing probabiity as cose as possibe to the exact resut obtained from the soution of the SE (for a recent review, see Ref. [1]). One such attempt was made by Kembe [2] who introduced a form for the probabiity which works at energies beow the top of the barrier, attains the vaue 1/2 at the top of the barrier as the exact resut dictates, but fais at energies above the barrier. Hi and Wheeer [3], used a paraboic form of the barrier which aows an exact anaytica soution of the Schrödingier equation, which coincides with the resut obtained using the Kembe formua for the same barrier for energies beow the top of the barrier. The HW form does work both beow and above the barrier. However it has the major shortcoming that the paraboic approximation is ony vaid at energies in the vicinity of the top of the Couomb barrier in nuceus-nuceus coisions. In the current paper we discuss the approximation of the Couomb barrier by a paraboa, which is frequenty adopted in the description of fusion of nucei. In this case, the tunneing probabiity is caed the transmission coefficient. The purpose of this paper is to demonstrate that the genera WKB-based Kembe formua for the tunneing probabiity can be extended to the cassicay aowed region encountered at energies above the top of the barrier. With this extension we have a WKB formua for the tunneing probabiity vaid for any we-behaved, compex nucear potentia at a energies. The paper is organized as foows. In section II, we discuss fusion reactions in heavy ions coisions. In Section III we discuss the WKB approximation as used in nucear physics, and introduce the concept of fusion barrier distributions, which has attracted considerabe interest over the ast two decades. In section IV, we deveop the extension of the Kembe formua to energies above the barrier top. This deveopment

2 2 A.J. Toubiana et a.: An improved WKB approximation for heavy ion fusion... is made for a typica potentia for nuceus-nuceus coisions. In section V, we use the approximations discussed in the previous section in the cacuations of fusion cross sections and barrier distributions in the case of the ight heavy ion system, 6 Li + 12 C. We chose a ight system to stress the imitation of the paraboic approximation discussed in the present work. For heavy systems, this approximation works better. They do break down but this occurs at ower energies (in comparison to the Couomb barrier). Finay, in Section VI severa concuding remarks are presented. j T W F V B j in r 1 r 2 E j R 2 Fusion reactions in heavy ion coisions In a fusion reaction the projectie merges with the target, forming a compound nuceus (CN). The densities of the coision partners overap so strongy that they ose their initia identities. Thus, the fraction of the incident current that reaches sma projectie-target distances does not emerge in the eastic channe. In this way, any description of the coision based excusivey on the eastic channe vioates the continuity equation. In muti-channe descriptions of the coision, this is not a probem because the current removed from the eastic channe re-emerges (after a time-deay) when the CN decays. However, if one adopts a singe-channe approach, it is necessary to simuate the oss of fux in some way. This coud be achieved through the introduction of a negative imaginary part in the projectie-target interaction potentia, W F (r), or by adopting an ingoing wave boundary condition (IWBC). In a heavy ion coision the rea part of the potentia has the genera form, V (r) = V C (r) + V N (r) + 2 ( + 1), (1) 2µ r2 where µ is the reduced mass of the projectie-target system, and V C and V N are respectivey the Couomb and the nucear potentias. The third term in Eq. (1) is the centrifuga barrier, which appears in the radia equations. For the sake of definiteness, we consider the 6 Li 12 C coision throughout this paper. As it is frequenty done, we approximate the Couomb potentia as, V C (r) = Z PZ T e 2, for r R C, (2) r = Z PZ T e 2 ) (3 r2, for r < R C. (3) 2 R C Above, Z P and Z T are respectivey the atomic numbers of the projectie and the target, and e is the absoute vaue of the eectronic charge. The Couomb radius, R C, is expressed in terms of the mass numbers of ( the projectie ) and the target (A P and A T ) as R C = r 0C A 1/3 P + A 1/3 T, with r 0C 1 fm. R 2 C Fig. 1. (Coor onine) Schematic representation of the rea and the imaginary potentias in a typica heavy ion coision. The cassica turning points are indicated by r 1 and r 2. For V N (r), we adopt the Akyüz-Winther potentia [4]. This potentia is evauated through a doube foding integra of the nucear densities with a M3Y [5] nuceonnuceon interaction [6]. For practica purposes, it is approximated by the Woods-Saxon (WS) function V N (r) = 1 + exp [(r R 0 ) /a 0 ], (4) ( ) where R 0 = r 0 A 1/3 P + A 1/3 T. The WS parameters of the Akyüz-Winther potentia are expressed as functions of the mass numbers of the coision partners. For the 6 Li+ 12 C system they have the vaues: V 0 = MeV, r 0 = 1.16 fm and a 0 = 0.56 fm. In singe-channe descriptions of fusion reactions based on Quantum Mechanics, one frequenty uses a compex potentia (the optica potentia), whose rea part is given by Eq. (1). The sum of the attractive nucear potentia with the repusive Couomb and centrifuga potentias gives rise to a potentia barrier, with height V B. The potentia has aso a negative imaginary part, W F, very intense and with a short range, that accounts for the incident fux ost to the fusion channe. The situation is schematicay represented in Fig. 1, for a coision energy (E) beow the Couomb barrier. The figure shows the incident current, j in, the refected current, j R, and the current transmitted through the barrier, j T, which reaches the strong absorption region. For practica purposes, the radia equation for each anguar momentum is soved numericay. The numerica integration starts at r = 0, where the wave functions vanish and their derivatives are chosen arbitrariy. This choice ony sets the normaization of the wave function, which has no infuence on the cross sections. Fusion absorption can aso be simuated by an ingoing wave boundary condition (IWBC), with a rea potentia. In the IWBC [7,8,9] one makes the assumption that the radia wave function at some distance r = R in ocated in the inner region of the barrier (usuay the minimum of V (r)) behaves as a wave propagating towards the origin. One then uses the WKB approximation to evauate the V 0

3 A.J. Toubiana et a.: An improved WKB approximation for heavy ion fusion... 3 radia wave function and its derivative at r = R in, and starts the numerica integration from this point. The fusion cross section, can be evauated by the partiawave series σ F = π k 2 =0 (2 + 1) P F P CN (5) where P F is the absorption probabiity at the th partiawave, given by P F = 1 S 2, (6) and P CN is the probabiity for the formation of the composite system, the compound nuceus, once the barrier is overcome. In most appications this probabiity is set equa to unity as the density of states of the CN is quite arge. In ight systems, such as 12 C + 12 C, it was demonstrated by Jiang et a. [10] that at ow energies the compound nuceus formation probabiity is significanty smaer than unity due to the ow density of states of the CN, 24 Mg at the excitation energies invoved. However, this phenomenon ies beyond the scope of the present paper. Thus, we wi set P CN = 1 for a vaues of the orbita anguar momentum. Owing to the imaginary potentia, or to an IWBC, the S-matrix is not unitary. Instead, it satisfies the condition: S 2 1. The deviation from unitarity, 1 S 2, measures the absorption (fusion) probabiity. If absorption is simuated by the IWBC, there is no outgoing wave for 0 < r < R in. This means that the fraction of the incident current that reaches this region eads to fusion. Therefore, the fusion probabiity is given by the transmission coefficient through the potentia barrier. Namey, P F T = j T j in. (7) In principe, the IWBC is equivaent to the effects of a deep imaginary potentia acting excusivey in the inner region of the barrier. One frequenty assumes that the detais of the imaginary potentia do not have significant infuence on fusion probabiities, as ong as the potentia is strong and has a short range. In fact, this assumption is questionabe. It may be vaid for heavier systems but it is not accurate in coisions of ight heavy ions, ike 6 Li 12 C. This fact is iustrated in Fig. 2, where the s-wave fusion probabiities obtained with the IWBC and with frequenty adopted vaues of the radius parameter and the diffusivity are shown. Besides studying the energy dependence of the fusion cross section, it was found usefu to go further and study the quantity D F (E) = d2 [Eσ F (E)] de 2, (8) which is referred to as the barrier distribution [11, 12]. Evauating the second derivative numericay, experimenta barrier distributions can be extracted directy from the T 0 (E) T 0 (E) (b) (a) V B V 0 r 0 a IWBC (-dep) E (MeV) Fig. 2. (Coor onine) The transmission coefficient obtained with the IWBC and fusion probabiities evauated with WS imaginary potentias with different r 0 and a 0 parameters (see the egend). The depth of the imaginary potentia is kept fixed at V 0 = 50 MeV. The resuts are shown in inear (pane (a)) and ogarithmic (pane (b)) scaes. data and comparing them with barrier distributions resuting from couped channe cacuations invoving different channes, one can obtain invauabe information about the reaction mechanisms invoved in the coision [13, 14, 15]. 2.1 The Wong formua and barrier distribution It has been a common practice in nucear physics to use approximate forms of the potentia barriers for the purpose of obtaining an anaytica soution of the scattering Schrödinger equation using the IWBC, and thus getting an anaytica form for the transmission coefficient. This has been first done by Hi and Wheeer [3], who approximated V (r) = V B 1 2 µω (r R B) 2, In the above equations, R B is the barrier radius, V B is the barrier height and ω is the curvature parameter, ω = 2 V (R B ). (9) µ

4 4 A.J. Toubiana et a.: An improved WKB approximation for heavy ion fusion... The transmission coefficient through this barrier, known as the Hi-Wheeer transmission coefficient, is given by the exact expression T HW (E) = exp [2 Φ HW (E)], with Φ HW (E) = 2π ω (V B E). This approximation coud be used for each partia wave, eading to anguar momentum dependent barrier radii, heights and curvatures. In this way, one woud have exact transmission coefficient for each and, using them in the partia-wave expansion, one coud obtain the fusion cross section. We discuss the WKB approximation in the cacuation of the transmission coefficients used to determine fusion cross sections (Eqs. (5) and (7)). Let us consider the 6 Li 12 C coision at a sub-barrier energy E, as represented in Fig. 1. Within the WKB approximation, the transmission coefficient through the barrier of the th partia-wave is given by Above, Φ WKB where T WKB (E) = exp [ 2 Φ WKB (E)]. (14) is the integra Φ WKB (E) = r2 r 1 κ (r) dr, (15) 2µ [ V (r) E ] Back in 1973, Wong [16] went one step further. He negected the -dependences of the barrier radii and curvatures, taking for a anguar momenta the vaues for = 0, denoted by R B and ω. In this way, the paraboic approximation for the effective potentia is V (r) = B 1 2 µω (r R B) 2, (10) where B is the barrier of the effective potentia for the th partia-wave, given by 2 B = V B + 2µR 2 ( + 1). (11) B He then treated the anguar momentum quantum number as a continuous variabe, λ = + 1/2, and approximated the sum of partia-waves by an integra over λ. In this way, he got the cosed expression, σ W F (E) = ω R2 B 2E n [ ( )] 2π 1 + exp ω (E V B). (12) Eq. (12), known as the Wong formua, has been quite popuar in ow energy heavy-ion fusion and reactions. More recenty, some improved versions of the Wong formua have been proposed [17, 18]. The fusion barrier distribution associated with the Wong formua can be evauated easiy. Using the cross section of Eq. (12) in Eq. (8), one gets κ (r) =. (16) In the above equations, r 1 and r 2 are the cassica turning points, indicated in Fig. 1. They are the soutions of the equation on r, V (r) = E. (17) Eq. (14) works very we for coision energies we beow the Couomb barrier, but it becomes very poor near the Couomb barrier and above. At E = B, the two cassica turning points coaesce, so that Φ WKB (E) = 0. In this way, one gets T WKB (E) = 1, whereas the quantum mechanica resut is T QM (E) = 1/2. Brink and Smiansky have shown that the situation can be improved if one takes into account mutipe refections under the barrier [21]. In 1935, Kembe [2] showed that the WKB approximation can be improved if one uses a better connection formua. He got the expression T K (E) = exp [2 Φ WKB (E)]. (18) At energies we beow the Couomb barrier, the two turning points are far apart and κ(r) reaches appreciabe vaues within the integration imits. In this way, Φ WKB (E) becomes very arge, so that the unity can be negected in the denominator of Eq. (18). This equation then reduces to Eq. (14). Therefore, the two approximations are equivaent in this energy region. On the other hand, Kembe s approximation remains vaid as the energy approaches the Couomb barrier, so that at E = B one gets Φ WKB (E) = 1/2, which is the correct resut. D W F (E) = π 2 R 2 B/ ω 1 + cosh [2π (E V B ) / ω]. (13) 3 The WKB approximation in Nucear Physics The WKB approximation is a short waveength imit of Quantum Mechanics. Since it is extensivey discussed in text books on Quantum Mechanics and Scattering Theory [19,20,6], its derivation wi not be presented here. We consider ony some appications in scattering theory. Unfortunatey, Eq. (18) gives wrong resuts at abovebarrier energies. In this energy range, there are no cassica turning points, since Eq. (17) has no rea soutions. Thus, Φ WKB (E) = 0. In this way, the transmission coefficient is frozen at the vaue 1/2. Therefore, it has the wrong high energy imit. However, Kembe [2] pointed out that the vaidity of Eq. (18) coud be extended to above-barrier energies by an anaytica continuation of the variabe r to the compex pane (see aso Ref. [22]), athough he did not eaborate on it. This point is considered in detai in the next section.

5 A.J. Toubiana et a.: An improved WKB approximation for heavy ion fusion Kembe transmission coefficient at above-barrier energies Eq. (18) can be extended to above-barrier energies if one soves Eq. (17) in the compex r-pane and evauates the integra of Eqs. (15) and (16) between the compex turning points. Compex turning points have been previousy used by different authors to evauate WKB-phase shifts [23, 24, 25,26]. In this section, we appy this procedure to a typica potentia in heavy ion scattering, and show that it actuay works. First, we introduce the dimensioness coordinate r x = r R B a 0, (19) where R B is the barrier radius and a 0 is the diffusivity parameter of the WS potentia. Next, this coordinate is extended to the compex pane. That is: x z = x + i y. The effective potentia of Eq. (1) then becomes, V (z) = Above, V exp (z + d) + Z PZ T e 2 a 0 z + R B d = R B R 0 a ( + 1/2) 2 2µ (a 0 z + R B ) 2. (20) is a dimensioness constant, associated with the difference between the barrier radius and the radius of the nucear potentia. Note that in Eq. (20) we repaced ( + 1) ( + 1/2) 2. This is a common procedure in semicassica theories. Then, writing the now compex effective potentia as, V (z) V (x + i y) = U (x, y) + i W (x, y), one imposes that its imaginary part vanishes. That is, { } W (x, y) Im V (z) = 0. (21) We shoud remind that the compexity of the potentia is of purey mathematica origin and shoud not be confused with absorption of fux. Using the expicit form of the effective potentia (Eqs. (1), (2) and (4)) in Eq. (21), one gets a compicated equation, that describes ines on the compex z-pane. Soving Eq. (21) by a numerica procedure, one gets the curve Γ represented by a bue dashed ine on pane (a) of Fig. 3. We remark that there are other curves on the compex pane where the potentia is rea. First, there is a curve Γ, at the eft of Γ. That is, for x < d. It ooks ike a distorted refection of Γ around a vertica axis at x = d. In this way, joining this curve with Γ, one gets a cosed contour. The potentia is aso rea on other Fig. 3. (Coor onine) The anaytic continuation of the reaistic potentia of Eq. (1). Pane (a) shows the curves Γ (bue dashed ine) and Γ e (red soid ine), and pane (b) shows the corresponding rea potentias on these ines. For detais, see the text. curves with y > π. The shapes of these curves are expected to be different from Γ, because the Couomb and the centrifuga potentias are not periodic in y. However, the physics invoved in the coision is fuy described by the anaytica continuation on Γ. Therefore, we restrict our discussion to this curve. We find that the curve Γ is very cose to the eipse Γ e given by the equation, ( ) 2 x + d ( y ) 2 + = 1. (22) d π This eipse is represented in Fig. (3) by a red soid ine. Ceary, the eiptica approximation for Γ is indeed quite good. Now we consider the potentia over the contour Γ. Owing to Eq. (21), we can express x as a function of y, x Γ (y), so that the corresponding compex coordinate becomes, z Γ (y) = x Γ (y) + i y. (23) In this way, the potentia depends excusivey on the coordinate y. To simpify the notation of the potentia over Γ, where it is rea, we write U Γ (y) = V (z Γ (y)). (24) Since Eq. (21) is rather compicated, the function x Γ (y) cannot be obtained anayticay. Therefore, the potentia U Γ (y) must be determined numericay.

6 6 A.J. Toubiana et a.: An improved WKB approximation for heavy ion fusion... The situation is much simper within the eiptica approximation for Γ. Using Eq. (22), we get [ ] x Γe (y) = d 1 1 y2 π 2. Inserting the above equation into Eq. (23), and using the resut in Eq. (20), we get an anaytica expression for the potentia V (z Γe (y)). However, this potentia is not exacty rea. Owing to the approximation of Γ by Γ e, V (z Γe (y)) has a sma imaginary part, which shoud be discarded. We then write, { } U Γe (y) = Re V (z Γe (y)). (25) This potentia, evauated for = 0, corresponds to the red soid ine in Fig. (3) (b). Ceary it is very cose to the exact potentia. Inspecting U Γ (y), one immediatey notices that it goes to infinity at the two extremes of the pot, which correspond to y = ±π. These divergences can be traced back to the poes of the WS potentia, ocated at r = R 0 + inπ, where n is any positive or negative integer. In terms of z, the poes are ocated at z = d+i nπ. The infinite vaue of the potentia at y = ±π has two important consequences. The first is that it confines the compex pane to the region π < y < π, which corresponds to d < x < 0. In this way, the curve Γ and curves with y > π are not accessibe to the system. Thus, such curves can be competey ignored. The second consequence is that the minima vaue of x on Γ is x = d, which corresponds to r min = R 0 and this vaue is arger than R C. Therefore, the Couomb potentia for points on Γ are aways given by Eq. (2). This justifies our choice of the Couomb potentia in Eq. (20). From Fig. 3 we concude that eipse is very cose to Γ, and the potentias over the two curves are neary the same. Thus, we henceforth repace the contour Γ by the approximate curve Γ e. This simpifies our cacuations consideraby. Now we evauate Kembe s transmission factor. For this purpose, one has to cacuate the WKB integra between the compex turning points, z ±, Φ WKB (E) = a 0 2µ z+ z dz V (z Γe ) E. (26) We remark that the factor a 0 resuts from the change of variabe of Eq. (19). The integra of Eq. (26) is independent of the integration path on the compex-pane. It depends ony on the integration imits, z and z +. However, for practica purposes, it is convenient to evauate the integra aong the contour Γ e. On this contour, V (x, y) reduces to the rea potentia of Eq. (25), U Γe (y), and the differentia dz can be written as [ ] dxγe (y) dz = + i dy, dy where dx Γe (y)/dy is a rea function of y. For energies above the barrier, U Γe (y) E is negative so that, dz ( ) dxγe (y) U Γe (y) E = + i i E U Γe (y) dy dy = G R (y) dy + i G I (y) dy, (27) where G R (y) and G I (y) are the rea functions, G R (y) = E U Γe (y), (28) G I (y) = dx (y) Γe E UΓe (y). (29) dy The vaues of y corresponding to the integration imits of Eq. (26) are given by the equation, U Γe (y) = E. They are schematicay represented in Fig. (3), for the coision energy E = 2 V B. Since the potentia U Γe (y) is symmetric with respect to y = 0, the turning points have the property, y = y +. (30) Thus, Eq. (30) indicates that the integration imits are symmetric. This property can be formay proved. Since F (z) E V (z) is an anaytic function of z in the vicinity of Γ (or Γ e ), we can write F (z) = F (z ). Thus, if z + is a compex turning point, that is, it satisfies the equation F (z + ) = 0, z = z + wi aso be a turning point. Then, writing z + = x + +i y +, the other turning point, wi be z = x + i y, with x = x + and y = y +. Expressing the integra of Eq. (26) in terms of the variabe y, one gets Φ WKB (E) = a 0 2µ { y+ Re G R (y) dy y + i y+ y G I (y) dy }. (31) Fig. 3(a) indicates that the potentia U Γe (y) is an even function of y, and so is the function G R (y) of Eq. (28). On the other hand, inspecting Fig. 3(b) one concudes that dx Γe (y)/dy is an odd function of y. For this reason, the function G I (y) of Eq. (29) is odd. Since the integration imits in Eq. (31) are symmetric, the second integra in this equation vanishes. Then, the WKB integra reduces to, Φ WKB (E) = a 0 2µ VB y+ y dy Above, we have introduced the notation, ε = E V B. ε U Γe (y) V B. Inserting Eq. (4) into (Eq. (18)), one gets Kembe s transmission coefficient at above-barrier energies.

7 A.J. Toubiana et a.: An improved WKB approximation for heavy ion fusion (b) T 0 (E) QM Kembe HW (a) 0.8 T 0 (E) V B E (MeV) Fig. 4. (Coor onine) Comparison the S-wave transmission coefficients obtained with exact quantum mechanics (bue soid ines), with the Kembe s approximation (stars), and the Hi- Wheeer transmission coefficient of Eq. (9) (red dashed-ines), corresponding to a paraboic approximation to the Couomb barrier. The system is 6 Li + 12 C. The resuts are shown in inear (pane (a)) and ogarithmic (pane (b)) scaes. 5 Study of fusion in the 6 Li 12 C coision In this section we assess the vaidity of the Kembe approximation extended to the compex r-pane, discussed in the previous section. In Fig. 4 we compare the S-wave transmission coefficients obtained with this approximation (bue soid ines) with the corresponding quantum mechanica resuts (stars). The figure shows aso the Hi- Wheeer transmission coefficients, corresponding to the paraboic fit to the Couomb barrier (red dashed ines). Ceary, the Kembe approximation reproduces accuratey the exact resuts, above and beow the Couomb barrier. On the other hand, the Hi-Wheeer transmission coefficients are cose to the exact ones at energies above the barrier, but it overestimates them by orders of magnitude at sub-barrier energies. This probem can be traced back to the paraboic approximation to the Couomb barrier at arge radia distances. Athough the paraboa is cose to the Couomb barrier at r R B, it fas off much faster as r increases. Then, at ow coision energies, the cassica turning point for the paraboic barrier is much arger than the one for the exact potentia. This makes the Hi- Wheeer transmission coefficient unreaisticay arge. Now we discuss the use of the same approximations to predict observabe quantities. First we ook at fusion cross Fig. 5. (Coor onine) Comparison of the exact quantum mechanica fusion cross section for the 6 Li + 12 C systems with resuts of approximations discussed in the previous sections. The resuts are shown in inear (pane (a)) and ogarithmic (pane (b)) scaes. sections. We compare the exact cross section, obtained with fu quantum mechanics (stars), with the cross section using Kembe s transmission coefficients through the exact barriers (bue soid ines). In a Kembe s cacuations, we used the eiptica approximation for Γ. We consider aso the Wong cross section, which is evauated with exact transmission coefficients but through the paraboae fitting the exact -dependent barriers. The resuts are shown in Fig. 5, in inear (pane (a)) and ogarithmic (pane (b)) scaes. One observes that the three cross sections of Fig. 5 show the same trends of the S-wave transmission coefficients of the previous figure. First, one notices that Kembe s approximation reproduces accuratey the quantum mechanica cross section in the whoe energy range. The cross section obtained with Wong s approximations is reasonaby cose to its quantum mechanica counterpart in the vicinity of the Couomb barrier, and it is sighty arger at higher energies (E > 4.5 MeV). The reason is that, as the energy increases, Wong s cross section gets reevant contributions from higher partia-waves, where the radius and the curvature parameters become significanty different from those for = 0, used for a partia-waves. On the other hand, the situation is very bad at sub-barrier energies. The agreement with the quantum mechanica cross section is poor, becoming progressivey worse as the energy decreases. At the owest energies in the pot, E 1

8 8 A.J. Toubiana et a.: An improved WKB approximation for heavy ion fusion... D F (mb / MeV) Li + 12 C Hi-Wheeer Morse Kembe QM r (fm) Fig. 6. (Coor onine) Barrier discributions for the fusion cross sections of the previous figure. MeV, it overestimates the quantum mechanica cross section by more than three orders of magnitude. The probem with Wong s cross section in this energy region resuts from the paraboic approximation to the -dependent barriers. As we expained in the discussion of the previous figure, the transmission coefficients through paraboae for E B are much arger than those through the corresponding exact barriers. Now we discuss the predictions of the same cacuations for the fusion barrier distributions. The resuts are presented in Fig. 6. The notation is the same as in the previous figure. Since barrier distributions are aways cacuated near the barrier, Wong s barrier distribution is not too bad. It is in quaitative agreement with the quantum mechanica fusion barrier distribution. On the other hand, the resuts of Kembe s WKB are extremey good. They reproduce the exact barrier distribution throughout the energy region of the figure. An interesting aspect of the figure is that the maximum of the quantum mechanica barrier, and aso of Kembe s barrier distribution, occurs at E max = 2.95 MeV. Athough is not far from the barrier of the bare potentia (V B = 3.25 MeV), it is about 0.3 MeV ower. 6 Concusions In this paper we considered the more genera, and potentiay more usefu, Kembe approximation for the tunneing probabiity, which is an improved version of the WKB approximation. We considered in detai the extension of Kembe s approximation to energies above the barrier, which can be achieved by judicious contour integration in the compex r-pane. By a judicious choice of a contour in the compex r-pane we have extended the vaidity of the Kembe formua to energies above the barrier. We have accompished this for barriers of a typica optica potentia in heavy ion coisions. We have discussed aso cross sections and barrier distributions, as given by the second derivative of Eσ F with respect to energy, and shown that Kembe s approximation gives an exceedingy good description of the quantum mechanica resuts at near-barrier energies. Thus we have, in this paper, finay obtained a WKB tunneing probabiity which is vaid for an arbitrary potentia and at a energies, without reying on the paraboic or other approximation. We are gratefu to Kouichi Hagino for providing us with a singe-channe code using the ingoing wave boundary condition, and to Rau Donangeo for criticay reading the manuscript. Partia support from the Braziian funding agencies, CNPq, FAPESP, and FAPERJ is aso acknowedged. M. S. H. acknowedges support from the CAPES/ITA Senior Visiting Professor Feowship Program A.J.T. acknowedges financia support from the French committee of the Brafitec program. References 1. K. Hagino, N. Takigawa, Prog. Theor. Phys. 128, 1061 (2012) 2. E.C. Kembe, Phys. Rev. 48, 549 (1935) 3. D.L. Hi, J.A. Wheeer, Phys. Rev. 89, 1102 (1953) 4. R.A. Brogia, A. Winther, Heavy Ion Reactions (Westview Press, 2004) 5. G.F. Bertsch, J. Borysowicz, H. McManus, W.G. Love, Nuc. Phys. A284, 399 (1977) 6. L.F. Canto, M.S. Hussein, Scattering Theory of Moecues, Atoms and Nucei (Word Scientific Pubishing Co. Pte. Ltd., 2013) 7. G.H. Rawitscher, Phys. Rev. 135, B605 (1964) 8. G.H. Rawitscher, Nuc. Phys. 85, 337 (1966) 9. K. Hagino, N. Rowey, A.T. Kruppa, Comp. Phys. Commun. 123, 143 (1999) 10. C.L. Jiang, B.B. Back, H. Esbensen, R.V.F. Janssens, K.E. Rehm, R.J. Charity, Phys. Rev. Lett. 110, (2013) 11. N. Rowey, G.R. Satcher, P.H. Steson, Phys. Lett. B 254, 25 (1991) 12. K. Hagino, N. Rowey, Phys. Rev. C69, (2004) 13. M. Dasgupta, D. Hinde, N. Rowey, A. Stefanini, Ann. Rev. of Nuc. Part. Sci. 48, 401 (1998) 14. L.F. Canto,, P.R.S. Gomes, R. Donangeo, M.S. Hussein, Phys. Rep. 424, 1 (2006) 15. L.F. Canto, P.R.S. Gomes, R. Donangeo, J. Lubian, M.S. Hussein, Phys. Rep. 596, 1 (2015) 16. C.Y. Wong, Phys. Rev. Lett. 31, 766 (1973) 17. C.Y. Wong, Phys. Rev. C 86, (2012) 18. N. Rowey, K. Hagino, Phys. Rev. C 91, (2015) 19. E. Merzbacher, Quantum Mechanics, 3rd edn. (John Wiey & Sons, Inc., 1998) 20. C.J. Joachain, Quantum Coision Theory (North Hoand, 1983) 21. D.M. Brink, U. Smiansky, Nuc. Phys. A 405, 301 (1983) 22. N. Fröman, P. Fröman, JWKB Approximation: Contributions to the Theory, 1st edn. (North-Hoand, Amsterdam, 1965) 23. D.M. Brink, N. Takigawa, Phys. Rev. C47, R2470 (1977) 24. S.Y. Lee, N. Takigawa, C. Marty, Nuc. Phys. A 308, 161 (1978) 25. S.Y. Lee, N. Takigawa, Nuc. Phys. A308, 189 (1978) 26. E. Vigezzi, A. Winther, Ann. Phys. 192, 432 (1989)