arxiv: v1 [nucl-th] 25 Nov 2011

Transcription

1 Deveopment of a Cox-Thompson inverse scattering method to charged partices Tamás Pámai 1, Barnabás Apagyi 1 and Werner Scheid 1 Department of Theoretica Physics Budapest University of Technoogy and Economics, H-1111, Budapest, Hungary Institut für Theoretische Physik der Justus-Liebig-Universität, Giessen, Germany arxiv: v1 [nuc-th] 5 Nov 11 Eectronic mai: pamai@phy.bme.hu, apagyi@phy.bme.hu November 8, 11 Abstract Cox-Thompson fixed-energy quantum inverse scattering method is deveoped further to treat ong range Couomb interaction. Depending on the reference potentias chosen, two methods have been formuated which produce inverse potentias with singuar or finite vaue at the origin. Based on the quaity of reproduction of input experimenta phase shifts, it is guessed that the p α interaction possesses an interesting repusive hard core. PACS:.3.Gp,.3.Zz,.6.Cb, 3.65.Nk, 4.1.-i, 5.4.Cm, 5.4.Ep 1 Introduction Inverse quantum scattering methods [1,, 3, 4] represent usefu techniques to assess effective interactions between coiding composite partices in a mode independent way. In nucear physics the modified Newton-Sabatier (mns) method [5, 6] has been used to determine, in genera, compex vaued optica potentias describing various systems, e.g., 1 C 1 C eastic scattering [7, 8] or to guess the spin-orbit potentias arising in the p α and n α coisions [9]. The Newton-Sabatier (NS) method [1, 11] has the property that it requires an infinite set of phase shifts otherwise the first moment of the potentia generated vanishes. This induces an unpeasant osciation of the potentia at arge distances, in addition to the artificia singuarity at the origin which remains to be present in the potentias generated aso by the mns method. In the recent decade another procedure, the Cox-Thompson (CT) inverse quantum scattering method [1] has been investigated [13, 14, 15]. This method has the advantage that it requires a finite set of phase shifts and possesses a non-zero first momentum of the potentia generated. In addition to this, in genera, it aso reproduces test potentias better than the mns method [13]. The CT method has been appied so far to uncharged partice scattering in order to reconstruct short ranged potentias from synthetic phase shifts [13], and aso to construct n α potentia from experimenta data [14]. We sha deveop here the method further with the intention that the CT method can be appied to charged partice scattering when the potentia to be constructed has a ong ranged Couomb tai. The various extensions wi be appied to both synthetic and experimenta phase shifts. The resuts are obtained by using the noninear sover package of MATLAB. Cox-Thompson (CT) method The origina CT method beongs to the fixed-energy procedures when a finite set of N phase shifts {δ } S is converted to the spherica potentia V(r) appearing in the th partia wave radia Schrödinger equation ( m d dr + m (+1) r ) +V(r) m k Ψ (k;r) = (1) 1

2 with boundary conditions Ψ (k;) r +1 and Ψ (k;r ) sin(kr π +δ ). We sha aso refer to another form of the Schrödinger equation ( ) d x dx +1 q(x) ψ (x) = (+1)ψ (x), () where dimensioness quantities have been introduced for the distance x = kr, for the potentia q(x) = E 1 V(xk 1 ) with E = k /(m), and the dimensioness wave function ψ (x) is reated to the wave function Ψ (k;r) through the asymptotic normaization constant Ψ (k;r) = C ψ (kr). The conversion from phase shifts {δ } to potentia V(r) can be carried out in two steps [13, 14]. First, we sove either of the foowing two systems of highy non-inear equations with e iδ = 1+iK+ 1 ik K ± = L T, S K + +K or tan(δ ) = +i(k + K S, (3) ), [M sin ] L [Mcos 1 ] L )π/ e±i(, S. (4) and { } { ( ) } Msin 1 sin ( L) π = M cos L L(L+1) (+1) cos ( ) ( L) π, S, L T (5) for the set T of the shifted anguar momenta {L} where the reations S T = and T = S = N hod. Next, one cacuates the potentia as V(r) = Eq(kr), (6) where with q(x) = ( ) d K(x,x) xdx x K(x,y) = L T A L (x)u L (y), (8) where the coefficient functions A L (x) are cacuated by soving the system of inear equations W[u L (x),v (x)] A L (x) (+1) L(L+1) = v (x), S (9) L T with πx πx u n (x) = J n+ 1(x), v n(x) = Y n+ 1 (x), n S T (1) being the reguar and irreguar soutions of the dimensioness free radia Schrödinger equation for the nth partia wave (Riccati-Besse functions [16]). Note that the set T of shifted anguar momenta L contains rea or (in case of compex phase shifts) even compex numbers. Using the behaviour of these functions at the origin and infinity, one easiy derives the foowing properties of the Cox-Thompson potentia. At the origin we get [13], in genera, a constant q() = Q (1 Q) L T with the matrix G L = (L ) 1 and the constant quantity Q = L T At the infinity the CT potentia in genera fas off ike an inverse power of two S (7) G 1 L 1. (11) S G 1 L. L+ 3 q(x ) = 1 x (α T cos(x) β T sin(x))+ x 3(α T sin(x)+β T cos(x)+γ T ) (1) where the coefficients depending on the set T are defined by the foowing reations α T = 1 ( a L cosl π b LsinL π ), (13) L T β T = 1 ( b L cosl π +a LsinL π ), (14) L T γ T = 1 ( b L cosl π a LsinL π ), (15) L T

3 with the coefficients a L,b L defined via the asymptotic expansion functions A a L (x) A L(x ) which admits a periodic expansion [15] as because of the form of the asymptotic version of equation (9) A a L(x) = a L cos(x)+b L sin(x) (16) A a cos(( L) π L(x) ) (+1) L(L+1) = cos(x π ), S. (17) L T 3 Extensions of CT method to ong ranged potentia In nucear physics one frequenty encounters charged partice scattering. The appearance of the ong ranged Couomb potentia requires an appropriate extension of the CT method which wi be carried out in various ways in this section. A the different extensions rey upon the two-potentia formaism [17]. Consider a spherica potentia V(r) which can be written as a sum of a short-ranged interior (or nucear) part ˆV and a ong-ranged exterior (asymptotic or reference) part V () (r), i.e., Then the phase shift itsef is aso spit into two parts V(r) = ˆV(r)+V () (r). (18) δ = ˆδ +δ () (19) where δ () is the phase shift due to V () and ˆδ corresponds to the phase shifts caused by ˆV (in the presence of V () ). If, for instance, ˆV(r), ˆδ =. If the short ranged part of the potentia is zero beyond a finite distance r a, ˆV(r) =, r r a () then, accordingy, the radia wave function can be written in this region as ( ) Ψ (k;r r a ) = C F () (kr)+tanˆδ G () (kr). (1) In the atter equation the functions F () and G () mean, respectivey, the reguar and irreguar soutions of the dimensioness radia Schrödinger equation () with potentia q(x) = E 1 V () (xk 1 ). 3.1 Phase transformation method (PCT) We may design a modified potentia Ṽ(r) which is identica with the origina one (up to a constant energy shift) within the interior region < r < r a and zero outside Ṽ(r) = { V(r) V(ra ), r r a,, r r a. () Thus if equation () appies the ˆV interior potentia can be deduced from Ṽ. The radia scattering wave function of this probem at the shifted energy E B = E V(r a ) kb /m is identica with the origina one in the inside region Ψ (k B ;r) = Ψ (k;r), r r a (3) but differs in the outside region Ψ (k B ;r) = C ) (u (k B r) tan δ v (k B r), r r a. (4) The equaity of the ogarithmic derivatives at r = r a ( d 1 dr og r Ψ (k B ;r)) r=r a = d dr og ( 1 r Ψ (k;r)) r=r a (5) 3

4 with Ψ (k;r) from (1) gives the transformed phase shifts δ (k B ) which can be used to (re)construct the short ranged potentia Ṽ(r) (). To get the whoe origina potentia one simpy adds the reference potentia { Ṽ () V(ra ), r r (r) = a, V () (6) (r), r r a, if it is known. The phase shift transformation method enabes the use of the CT method as formuated in the preceding section. It requires that the modified potentia Ṽ(r) be constructed from the set { δ } by the actua inverse method as accuratey as possibe. But no exact fixed energy inverse scattering method has yet been formuated which is abe to produce a short range potentia that is exacty zero beyond a finite radius, if a finite set of phase shifts is empoyed. Therefore, one aways introduces a sma error into the potentia (re)construction when appying the phase shift transformation procedure. The phase shift transformation method has been introduced by May, Münchow and Scheid [6] and is widey appied in case of the modified Newton Sabatier method. It is appied aso to the CT inverse procedure to treat Couomb scattering [18] and we sha ca this combined method as the PCT method. 3. Generaized CT scheme Based on the two potentia formaism one can derive a generaized CT (gct) scheme which empoys the given nucear phase shifts {ˆδ } for constructing the short-ranged (nucear) potentia ˆV(r) but avoids use of a matching radius r a. The derivation starts with the ansatz for the input symmetrica kerne of the Ge fand-levitan-type integra equation g(x,x ) = S γ F () (x < )G () (x > ), (7) where x < (x > ) denotes the esser (greater) of x,x. Then, by proceeding through the usua steps [13, 14], one arrives at a system of noninear equations identica in structure to (3) e iˆδ = 1+iK+ 1 ik or tan(ˆδ ) = K + +K +i(k + K ), S, (8) with and { Msin M cos } L K ± = = L T, S [M sin ] L [M 1 1 L(L+1) (+1) [ cos] L e ±i ( sin cos ( )π/+δ () δ () ( L) π +δ() L δ() ( ( L) π +δ() L δ() ], S, (9) ) ), S, L T (3) for the set T of the shifted anguar momenta L where the reations S T = and T = S hod. From the set T one cacuates the nucear potentia as ˆV(r) = Eq(kr), (31) where with q(x) = ( ) d K(x,x) xdx x K(x,y) = L T (3) A L (x)f () L (y). (33) The coefficient function A L (x) is cacuated by soving the system of inear equations L T A L (x) W[F() L (x),g() (x)] (+1) L(L+1) = G() (x), S, (34) From this gct scheme severa deveopments are possibe as the actua form of the reference potentia has not yet been specified. 4

5 3.3 Couomb reference potentia method (CCT) If one sets the reference potentia to be the bare Couomb potentia V () (r) = E η kr 1 Z 1 Z e, (35) 4πε r with η being the Sommerfed parameter then one arrives at the Couomb CT (CCT) method. In this case the reguar and irreguar reference functions are the reguar and irreguar Couomb functions, the reference phase shift becomes the Couomb phase σ = 1 i n [ Γ(+1+iη) Γ(+1 iη) ], S. (36) Here it shoud be mentioned that in the courseofappication ofcct method it may become necessary to know the reguar and irreguar Couomb functions and Couomb phases for compex orders. For exampe, in case of non-eastic scattering the phases are compex and therefore, as noted before, the L numbers are aso compex vaued. The Couomb functions are we-defined for compex orders and simiary to the rea order case they can be given as power series (for detais see e.g. [19, ] and Appendix). It is interesting that contrary to the fact that the asymptotic form of the reference functions contains the we known ogarithmic term ηnkr in the argument, it does not appear in the gct formuas because of canceation. Using the known power series of the Couomb functions it can be shown that the CCT method gives a potentia which is proportionato the Couomb potentia near the origin. For one term, i.e. T = S = 1, we get [ ] L(1+) η V(r ) = E +O(1). (37) (1+L) kr For arge r we get V(r ) = E η kr E (kr) L T S ( ) cos ] 1 Lcos(Θ L (kr)+θ (kr))+o (kr) 3, (38) [M 1 with Θ L (x) = x ηnx L π +σ L. We see that the CCT method generates an inverse potentia that gives a Couomb-ike singuarity at the origin and produces a damped osciation around the Couomb tai at arge distances. It is free of the matching parameter r a and requires just the nucear phase shifts ˆδ which are derived by the usua phase shift anaysis procedures. 3.4 Modified Couomb reference potentia method (MCT) In order to obtain an inverse potentia that is finite at the origin, instead of being singuar there, we can modify the Couomb reference potentia accordingy. This reference potentia is constant in the interior domain and purey Couombic outside. This modified Couomb potentia is the same as that empoyed by the phase transformation method for Couomb asymptotics and reads as { E V () η (r) = kr a, r r a, E η kr, r r. (39) a To this reference potentia there beong the foowing reguar and irreguar radia wave function ( ) F () u 1 η kr (kr) = a kr, r < r a, α () F F (kr)+β () F G (kr), r > r a, ( ) G () v 1 η kr (kr) = a kr, r < r a, α () G F (kr)+β () G G (kr), r > r a, and reference phase shift δ () = σ +arctan 5 (4) (41) ( ) βf (). (4) α () F

6 The coefficients α F (),β F (),α () G,β () G can be cacuated from the equaity of the inner and outer wave functions and their derivatives at the matching radius r a. For exampe, the two coefficients necessary for cacuating the reference phase shift are α F () = u ( 1 η xa xa ) G (x a) ( ) 1 η xa u 1 η xa xa G (xa) F (x a) G (x, β F (xa) () F a) G (xa) = u ( 1 η xa xa ) F (x a) with x a = kr a and prime denotes derivation with respect to argument. The tota phase shifts can be written in two different ways ( ) 1 η xa u 1 η xa xa F (xa) G (x a) F (x, (43) G (xa) a) F (xa) δ = ˆδ +σ = δ MCT +δ () (44) where ˆδ means the nucear phase shifts given as data and δ MCT is to be used to perform the CT inverse cacuation outined above. The potentia obtained by the MCT method has a finite vaue at the origin. Because this method empoys a simiar reference potentia as the PCT method, the resuts provided by the two methods shoud aso be very simiar, athough quite different functions are used in the cacuations. The advantage of the MCT over the PCT ies in that it does not invove the sma error in the phase shift reproduction inherent to the PCT. 4 Appications to (re)construct effective potentias 4.1 Synthetic phase shifts To iustrate the genera appicabiity of the ong-range CT inversion procedures we sha first use them to reconstruct mode potentias. We mode the α α scattering with two sighty different potentias: the first is finite at the origin and the second is singuar (describing a possibe non-ocaity). The first mode potentia is given by V(r) = U(r)+V C (r), (45) where the nucear interaction is described by a Woods-Saxon form ( ) 1, U(r) = U 1+d r R a (46) and the Couomb interaction is represented by the potentia of a homogeneous charged sphere as { ( ) Z1Z e V C (r) = R C 3 r, r R R C, C Z 1Z e r, r > R C. For the various parameters we choose the foowing vaues: Z 1 =, Z =, A 1 = 4, A = 4, U = MeV, R = R C =.636 fm, a =.5 fm. We have cacuated the nucear phase shifts ˆδ at energies E = 5 and 35 MeV (see tabes 1 and ). These phase shifts are then used as input data for the various CT cacuations. The resuts for the L-vaues obtained using the noninear sover of MATLAB are shown in tabes 1 and. The corresponding potentias are dispayed in figure 1. One can observe that as expected the PCT and MCT procedures give amost the same resuts and the CCT inverse potentia is divergent at the origin. By increasing the scattering energy more phase shifts become avaiabe for the inversion (see tabes 1 and ) and the potentia reproduction becomes better. The second mode is obtained by adding the singuar potentia term to the previous mode, i.e., (47) V sing (r) = d r r (48) V(r) = U(r)+V C (r)+v sing (r). (49) The reconstruction of phases is isted in tabe 3 and the potentias are shown in figure. We see that the CCT potentia foows nicey the mode potentia in the singuar domain near the origin whie the PCT and MCT methods are unabe to reproduce the singuarity athough their phase shift reconstruction is good. 6

7 Tabe1: Mode data (ˆδ orig phaseshifts), inversionresuts(l shifted anguarmomentaand ˆδ differences between the mode phase shifts and the ones given by the various methods ) of the α α scattering at E c.m. = 5 MeV. The matching parameter used in the PCT and MCT procedures was set to r a = 1 fm. L CCT L MCT PCT orig L ˆδ ˆδ CCT ˆδ MCT ˆδ PCT Note that for the sake of comparison the phases given by the inverse potentias were cacuated by cutting-off the non-physica osciations beyond the matching radius r a used in the PCT procedure. Without the cut-off the MCT and CCT potentias reproduce the phase shifts within an error of the numerica precision. Tabe : Mode data and inversion resuts of the α α scattering at E c.m. = 35 MeV. The matching parameter used in the PCT and MCT procedures was set to r a = 1 fm. L CCT L MCT PCT orig L ˆδ ˆδ CCT ˆδ MCT ˆδ PCT Experimenta phase shifts Our goa is to assess the effective centra potentia governing the p α scattering events. Comprehensive data of phase shift anaysis of the p α scattering has been presented by Ai, Ahmad and Ferdous in [1]. Because of the spin-orbit couping both spin-up δ + and spin-down δ phase shifts contribute to the scattering ampitude at each partia wave. In case of weak spin-orbit couping the combined phase shifts ˆδ = 1 +1 [(+1)δ+ +δ ] (5) are characteristic of the underying centra potentia [], and we sha use them as input for the CT procedures. We have inverted the phase shift data of Ai et a. Since the inverse potentias exhibit simiar characteristics we present ony some of the inversion resuts here. In figure 3 the inverse potentias yieded by three CT methods at E ab = MeV proton energy are depicted. The resuts at this energy are representative of the potentias recovered beow the E ab =.94 MeV, α+p d+ 3 He ineastic threshod. As we see the MCT and PCT potentias are amost identica and strongy resembe a Woods-Saxon form. The range and strength of a the three potentias are simiar 7

8 V(r) (MeV) CCT PCT MCT Mode V(r) (MeV) CCT PCT MCT Mode (a) E c.m. = 5 MeV (b) E c.m. = 35 MeV Figure 1: Mode with potentia (45) and inverse potentias yieded by the CCT, MCT and PCT methods (abeed accordingy) at E c.m. = 5 and 35 MeV center of mass energies. The mode potentia is nonsinguar at the origin. but the CCT potentia is different in shape: it possesses a repusive core. Note that a repusion at sma distances can be theoreticay accounted for as the manifestation of the Couombic non-ocaity [3]. It can aso be seen that this repusion core stabiizes the inverse potentia in the sense that the ampitude of the asymptotic osciations is diminished compared to the non-repusive MCT/PCT resuts. This is aso the reason why the phase shift reproduction with a given precision of the CCT potentia is better than that of the MCT/PCT potentias (see tabe??). In figure 4 numerous other CCT and PCT potentias are shown at various energy vaues incuding those above ineastic threshod. (The MCT potentias are not shown because they coincide within the width of ine with the PCT resuts.) Apart from the Couombic singuarity at the origin, the potentias have a simiar range of 3 4 fm and strength of 5 7 MeV (PCT) and 5 16 MeV (CCT). The imaginary part is much ess compared to the rea part. The reproduction of phase shifts (not shown) gets better at higher energy in case of the PCT potentias because of the fixed cut-off radius which in principe does not appy to CCT (and MCT) method. Without use of this radius the CCT (and MCT) potentias give back the input phase shifts exacty. 5 Concusions, summary We have deveoped the Cox-Thompson(CT) fixed energy inverse quantum scattering method into various directions in order to make it appropriate for treating ong range interactions. By expicit cacuation we have shown that the modified (MCT) and phase transformed (PCT) methods yied practicay the same potentias athough quite other functions are invoved in the cacuation. These potentias are generay finite at the origin therefore they can be used for cases with no singuarity at sma distances. The Couomb (CCT) method produces a characteristic Couomb singuarity at the origin. This method can be appied with success when nonocaity or a repusion hard core pays a roe in the interaction. We have appied the methods first to mode cases with and without singuarity. The mode potentias describing α α scattering are we reproduced by the new methods according to their characteristic properties concerning the finiteness or infiniteness at the origin. Then, the experimenta p α phases of Ai et a [1] have been inverted with the resut that the CCT potentias showing up a repusive hard core reproduce better the input phase shifts than the PCT potentias which exhibit a Woods-Saxon shape. Whether or not the repusive core is present between the proton and apha partice when being scattered by each other shoud be carified aso using independent source of information(see e.g. [1, 4, 5, 6, 7]). Acknowedgements The authors thank the DFG for supports through the contract No. 436 UNG 113/1/-1. 8

9 Tabe 3: Mode data and inversion resuts for the divergent mode potentia at E c.m. = 5 MeV. The matching parameter used in the PCT and MCT procedures was set to r a = 1 fm. L CCT L MCT PCT orig L ˆδ ˆδ CCT ˆδ MCT ˆδ PCT Tabe 4: Inversion resuts of the experimenta p α data at E ab = MeV. The matching parameter used in the PCT and MCT procedures was set to r a = 7 fm from which distance aso the MCT and CCT potentias have been repaced by the pure Couombic tai. L CCT L MCT PCT orig L ˆδ ˆδ CCT ˆδ MCT ˆδ PCT A Couomb functions of compex order In certain appications of both the CCT and MCT formuations of the Cox-Thompson inverse scattering method for ong-ranged potentias it is necessary to evauate the reguar and irreguar Couomb wave functions for compex orders. Based on [19, ] the evauation is accompished by using the power series given beow. The reguar Couomb wave function for L C compex order is given by F L (x) = L d πη [Γ(L+1+iη)Γ(L+1 iη)] 1 Γ(L+) and the C L j constants are defined by the recursion x L+1 j=l+1 C L j (η)x j L 1, (51) C L L+1 = 1, (5) C L L+ C L j = η L+1, (53) = ηcl j 1 CL j (j+l)(j L 1) j > L+. (54) This production of the F L (x) is a simpe anaytic continuation of the formuae in [16]. However the irreguar Couomb wave function cannot be given by such a simpe generaization. We utiize therefore the fact that the Couomb wave functions and the Whittaker functions are reated to each other ineary for fixed L and η. By means of the Whittaker functions the reguar Couomb wave 9

10 6 3 V(r) (MeV) CCT PCT MCT Mode Figure : Mode with potentia (49) and inverse potentias yieded by the CCT, MCT and PCT methods (abeed accordingy) at E c.m. = 5 MeV center of mass energy. The mode potentia is singuar at the origin. function is aso given. with F L (x) = A(L,η)M ( iη,l+ 1,ix), (55) G L (x) = ia(l,η)m ( iη,l+ 1,ix) +B(L,η)W ( iη,l+ 1,ix) (56) πη A(L,η) = d i(l+1)π d [Γ(L+1+iη)Γ(L+1 iη)] 1 Γ(L+) (57) πη B(L,η) = d+ilπ d + [Γ(L+1+iη)Γ(L+1 iη)] 1 Γ(L+1+iη). (58) TheWhittaker functions M(a,b,y)and W(a,b,y)arewe-knownforcompex arguments(see [16] chapter 13). The Gamma function was cacuated by a Lanczos series approximation. The η parameter appearing in the above formuae is the Sommerfed parameter given with the quantities discussed earier as References η = k E e 8πε Z 1 Z = e 8πε m E Z 1Z. (59) [1] Chadan K and Sabatier P C 1977 Inverse Probems in Quantum Scattering Theory (New York: Springer Verag) [] von Geramb H (ed) 1994 Quantum Inversion Theory and Appications (Lecture Notes in Physics vo 47) (Berin: Springer) [3] Apagyi B, Endrédi G and Lévay P (eds) 1997 Inverse and Agebraic Quantum Scattering Theory (Lecture Notes in Physics vo 488) (Berin: Springer) [4] Apagyi B and Horváth M(eds) 8 Proceedings of the Internationa Conference on Inverse Quantum Scattering Theory (Specia Issue, Modern Physics Letters B vo, issue 3) (Word Scientific) 1

11 5 4 E ab = MeV 3 V(r) (MeV) CCT MCT PCT Figure 3: Inverse p α potentias V(r) obtained from input phase shifts (given in tabe??) as a function of the radia distance r at energy E ab = MeV (E c.m. = MeV, k =.731 fm 1 ). Curves obtained by the PCT, MCT and CCT method are abeed accordingy. ˆδ orig [5] Münchow M and Scheid W 198 Phys. Rev. Lett [6] May K-E, Münchow M and Scheid W 1984 Phys. Lett. B [7] Apagyi B, Ostrowski A, Scheid W and Voit H 199 J. Phys. G: Nuc. Part. Phys [8] Apagyi B, Schmidt A, Scheid W and Voit H 1994 Phys. Rev. C [9] Aexander N, Amos K, Apagyi B, and Lun D R 1996 Phys. Rev. C [1] Newton R G 196 J. Math. Phys [11] Sabatier P C 1966 J. Math. Phys [1] Cox J R and Thompson K W 197 J. Math. Phys [13] Apagyi B, Harman Z and Scheid W 3 J. Phys. A: Math. Theor [14] Mechert O, Scheid W and Apagyi B 6 J. Phys. G: Nuc. Part. Phys [15] Pámai T, Horváth M and Apagyi B 8 J. Phys. A: Math. Theor [16] Abramowitz M and Stegun I A 197 Handbook of Mathematica Functions (New York: Dover Pubications) [17] Joachain Ch 1975 Quantum Coision Theory (Amsterdam: North-Hoand) chapter 17 [18] Mechert O 5 Das Inverse Streuprobem bei fester Energie mit der Methode von Cox und Thompson, Thesis (University Giessen) [19] Thompson I J, Barnett A R (1986) J. Comput. Phys [] Dzieco A, Yngve S (1999) J. Math. Phys [1] Ai S, Ahmad A A Z and Ferdous N (1985) Rev. Mod. Phys

12 - - Re V(r) (MeV) E ab = MeV E ab = MeV E ab = 3.48 MeV E ab = 3.43 MeV E ab = 39.8 MeV E ab = 48.8 MeV Re V(r) (MeV) E ab = MeV E ab = MeV E ab = 3.48 MeV E ab = 3.43 MeV E ab = 39.8 MeV E ab = 48.8 MeV (a) Rea part of the CCT potentia. (c) Rea part of the PCT potentia Im V(r) (MeV) 5-5 Im V(r) (MeV) E ab = 3.48 MeV E ab = 3.43 MeV E ab = 39.8 MeV E ab = 48.8 MeV E ab = 3.48 MeV E ab = 3.43 MeV E ab = 39.8 MeV E ab = 48.8 MeV (b) Imaginary part of the CCT potentia. (d) Imaginary part of the PCT potentia. Figure 4: Compex-vaued inverse p α potentias yieded by the CCT (a,b) and PCT (c,d) methods at various E ab proton energies beow and above the ineastic threshod using the experimenta phase shifts of [1]. [] Leeb H, Huber H and Fiededey H (1995) Phys. Lett. B [3] Bacheet G B, Hamann D R, Schüter M (198) Phys. Rev. B [4] Gamme J L and Thaer R M (1958) Phys. Rev [5] Satcher G R, Owen L W, Ewyn A J, Morgan G L and Water R L (1968) Nuc. Phys. A 11 1 [6] Thompson G E, Epstein M B and Sawada T (197) Nuc. Phys. A [7] Thompson D R, Reichstein I, McCure W and Tang Y C (1969) Phys. Rev