ARS - 2008
Permanent URI for this collectionhttp://repository.kln.ac.lk/handle/123456789/166
Browse
5 results
Search Results
Item On the systematic of anomalous absorption of partial waves by nuclear optical potential(Faculty of Graduate Studies, University of Kelaniya, 2008) Amarasinghe, D.; Munasinghe, J.M.; Piyadasa, R.A.D.An interesting phenomenon relating to the nuclear optical potential was discovered (Kawai M & Iseri Y,(1985)) [1] which is called the anomalous absorption of partial waves by the nuclear optical potential. They found, by extensive computer calculations, that, for a special combinations of the total angular momentum (j) ,angular momentum(/) ,energy (E) and the target nuclei(A), the elastic S-matrix elements corresponding to nucleon elastic scattering become zero. This phenomenon is universal for light ion elastic scattering on composite nuclei. [2]. It is very interesting that this phenomenon occurs for the realistic nuclear optical potential and it exhibits striking systematic in various parameter planes. For example, all nuclei which absorb a partial l waves of a definite node lie along a straight in the plane (Re, A 3 ) as shown in the figure , where Re is the closest approach and A is mass number of the target nucleus. Theoretical description of this systematic has been actually very difficult, though attempts have been made by the Kyushu group in Japan. In this contribution, we explain mathematically the most striking systematic of this phenomenon. Explanation of the systematic Partial wave· u 1 ( k, r) of angular momentum I and incident wave number k satisfies the Schrodinger equation d21 + [ k2 _ l (l : l ) _ 21 {V(r)+iW(r)}] u,(k,r)= 0 dr r 1i , where V (r) is the total real part and W (r) is the total imaginary part of the optical potential. Starting from this equation , one obtains (1) lu1(k,r)I2=2 XI du, 12 -g(ru1(r2Jdr (2) dr dr 0 where g(r)= [k2-�� V(r)-l(l r:1)J. If u1(k, r) is the anomalously absorbed partial wave, the corresponding S-matrix element is zero and hence in the asymptotic region I u1(k,r) I is almost constant. Therefore [1;1' -g(ru1(k,rt ] o (3) for large r. Now, from (1) and (3), it is not difficult to obtain[3] the equation _1 !!I 12 -- g'( r) wh (r) 'Jw ( )J 12d (4) 2 u1 - ( ) + 2 h r,., u, r lu,l dr 2g r g(r u,l 0 166 Proceedings of the Annual Research Symposium 2008- Faculty of Graduate Studies University of Kelaniya which is valid for large r , and has been numerically tested in case of an anomalously absorbed partial waves , where Wh(r) =- 2 W(r) . If W(r) decays much more rapidly n than V(r) in case of a partial wave under consideration lu,l2 =- g'( ( r) ) and by lu,l dr 2g r integrating this equation with respect to r, we obtain I iu1 (k, r)i2 (g(r) 2 = C (5) ,where C is a constant, and the equation (5) is valid for large values of r. In case of anomalous absorption of the partial wave, I u 1 ( k, r ) I is constant in the asymptotic region and therefore g(r) is also constant. We have found that for all partial waves corresponding to a straight line of definite node, g(r) is constant at the respective I [l(l + 1)]2 closest approach. For example, at Re = k , g(r) is constant for all partial waves lying on a straight line in case of anomalous absorption of neutron partial waves by the nuclear optical potential. Therefore, neglecting the spin-orbit potential , we get -I I 21tV0[1+exp[([l(/+l)F -1.17A3)]/arr1 =C0 n2 k where V0 is depth of the real potential and A is the target mass and the optical potential parameter ar = 0.75 and C0 is a constant. Therefore, in case of neutron, we get the linear relation [I (I+ 1)]2 = 1.11 A + C1 k . (6) where C1 is again a constant. This relation has found to be well satisfied in the cases we have tested numerically. The equation (6) well accounts for the anomalous absorption of neutron partial waves by the Nuclear Optical Potential as shown in the figure below. : [:::::r�;;�����i����:::]::::::::::::::::::::::::] :::::::::--::::::::::::: 1 I I I I I I 7 L----------J------------L---------__ J _ ----------L----------- : ----- ------L----------- ::N :: 6 : ----------: ------------: -----------: _ ___________ :L _ -------i: --------..:---: ----------- ;::::" 5 ,1- ------ l l l : : ! - ---,------------r---- -------, ---- ------r--- --------,------------r----------- + : : : I : : : ::-::::: 4 :I- -----------;I ------ I I I I I ------:----- ---"t-------- ----:------------ 1------------:------------ 1....1 I I I I I I I 3 lI- ----------1I ------------I -----------1I ------------rI -----------1I ------------rI ----------- 2 ::- ----------1: ------------: -----------1: ------------: -----------1: ------------: ----------- I I I I I I I ,a , 2 3 4 5 6 All3 Gradient of straight line predicted by ( 6) is 1.1 7 and the actual value is 1.1828 . Very small discrepancy is due to the negligence of the spin-orbit potential.Item Singularities of the elastic S-matrix element,(Faculty of Graduate Studies, University of Kelaniya, 2008) Jayasinghe, W.J.M.L.P.; Piyadasa, R.A.D.It is well known that the standard conventional method of integral equations is not able to explain the analyticity of the elastic S-matrix element for the nuclear optical potential including the Coulomb potential. It has been shown[1],[2] that the cutting down of the potential at a large distance is essential to get rid of the redundant poles of the S-matrix element in case of an attractive exponentially decaying potential. This method has been found [3] to be quite general and it does not change the physics of the problem. Using this method , analiticity and the singularities of the S-matrix element is discussed. Singularities of the elastic S-matrix element Partial wave radial wave equation of angular momentum l corresponding to elastic scattering is given by, [ d2 - 2 /(/+ 1)] 2p [ . ] 2 + k - 2 u1(k,r)=-2 V(r)+ Vc(r)+ zW(r) u1(k,r) M r n (1) where V (r) is the real part of nuclear potential, W (r) is the imaginary part of the optical potential · (r) is the Coulomb potential, and k is the incident wave number. Energy dependence of the optical potential is usually through laboratory energy E1ah and hence it depend on k2 and therefore k2- 2 [V(r) + Vc(r) + i W(r)] is depending on k n through e . It is Well known that . (r) is independent Of k Jn Order tO make U 1 ( k, r) an entire function of k , we impose k independent boundary condition at the origin . Now, we can make use of a well known theorem of Poincare to deduce that the wave function is an entire function of k2 and hence it is an entire function of k as well. We cut off the exponential tails of the optical potential at sufficiently large Rm and use the relation -1 --d u1 =.:u ;<--l -(k,-r)--s--'-r (--k,R--m----') u-'-;(+l--- ( k,-r) u, dr u/-l(k ,r)-sr(k,Rm) uj+l(k,r) (2) to define St(k,Rm),where u,<-l(k,r) and u,<+l(k,r) stand for incoming and outgoing Coulomb wave functions respectively which are given by I (±l(k )-+· [r(/+ l+i 17) ]2 [J[2"+iU+nJ w (-2 'u1 ,r - _l . e 1 1k r ) r(Z+I-z17) +i,,/+2 (3) where Ware the Whittaker functions. In the limit Rm oo St (k,Rm) ,the nuclear part of the S-matrix element , becomes St { k) and the redundant poles removed[1 ],[2].Now, the nuclear S-matrix element , in terms of the Whittaker functions is given by 143 Proceedings of the Annual Research Symposium 2008- Faculty of Graduate Studies University of Kelaniya where w' 1 (2ikr)-(k,r) W I (2ikr) IIJ, 1+-2 in'" 1+-2 , ,r 2 Rm W. 1 (-2ikr)-(-k,r)W 1 (-2ikr) -IIJ, I+-2 -ill ' 1+-2 P1(k,r) = u;(k,r) ,and St (k) has an essential singularity at k = 0, which u1(k,r) , (4) is apparent from the Wister's definition of the Gamma function l(z) smce z= l+ 1 ±i lJ .However, this singularity has no any physical meaning and is an outcome of treating 21Jk as well defined quantity for all k including k = 0 in the corresponding r Schrodinger equation .The infinite number of zeros and poles of S- matrix element due to the Gamma functions associated with S - matrix element have to be interpreted 1 carefully. S;'(k)=O at the zeros of ---- f(l+1+i1J ) and then the total wave function reduces to [ . ] I J( '7 +i(/+I)ZZ"j uj-l(k,r)=-i f(l+1-'7) e l 2 2 W (2ikr) f(l+1+zlJ) i11,1+l2 which is also zero. Even though the corresponding energies of these states are negative since the corresponding wave number is given by ? k= z · z,z2 e- 2 n= 0' 1' 2 , ... 11 (n+l+1) they are not physically meaningful bound states as found in[1],[2] long ago. These states are unphysical since poles are redundant poles. This fact is clearly understood by the fact that all these poles are absent in the physically meaningful total S - matrix element. For large 1k1, Sin' (k) (- ) { e-21k r S(k), where S(k) = [-ik+(k)] • +2k [ik-(-k)]. smce W = e- " for large k. Therefore the S-matrix element has an essential singularity at infinity, which is on the imaginary axis. It is clear that there are no redundant poles in the total S-matrix element is free from redundant poles sinceSJ (k) =SeS t , where Se = f(l + 1 + i7J) f(l+1-i1]) .Item Simple theorem on the integral roots of special class of prime degree polynomial equations(University of Kelaniya, 2008) Piyadasa, R.A.D.Even in case of a simple polynomial x3 + l5xb + 28 = 0 , where (3, b) = 1 , it may be extremely difficult to discard the integral solutions without knowing the number b exactly .In this case, one can make use of the method of Tartaglia and Cardan [Archbold J.W.1961] and its solutions can be written as u + v,Ul:o + vw2 ,uw2 + vw, where u3, v3 are the roots of the equation x2 + 28x -125b3 = 0, and {J) is the cube root I . . [-28 ± .J282 + 500b3 J 3 of umty .Also, u or v can be wntten as 2 and this expression is obviously zero only when b = 0. Therefore if b :f. 0 , it is very difficult to determine that I k [-28±.J282+500b3J3. . = Th h 'lib I' d · h 2 IS an mteger or not . e t eorem w1 e exp rune m t e . following , is Capable of discarding all integral solutions of this equation using only one condition (3, b) = 1 . The theorem in its naive form discards all integral solutions of the polynomial . xP + pbx-cP = 0, where p is a prime and (p,b) = (p,c) = 1 Theorem xP + pbx-cP = 0 has no integral solutions if (p,b) = (p,c) = 1 , where b,c are any integers and p is any prime . Proof Proof of the theorem is based on the following lemma Lemma If (a,p) = (b,p) = 1, and ifs= aP-bP is divisible by p , then p2 dividess. This is true even when s = a P + bP and p is odd. Proof of the Lemma s=aP-a-(bP-b)+a-b and since s is divisible by p and aP-a ,bP-b are divisible by p due to Fermat's little theorem, it follows that a-b is divisible by p. aP -bP= (a-b)[(ap-l -bP-1)+b(aP-2 -bP-2)+ .. ·+bP-3(a-b)+ pbP-1)] (1) From ( 1 ), it follows that s is divisible by p2 • Proof of the lemma for a P +bP is almost the above. It is well known that the equation xP + pbx-cP = 0 (2) has either integral or irrational roots. If this equation has an integral root l, let x = l and (p,l) = 1. Then , lP -cP+ pbl= 0 . From the Lemma, it follows that p2 I CF - cP) .Therefore pI b , and this is a contradiction. Therefore equation has no integral roots which are not divisible by p .If it has an integral solution which is divisible by p , then let x=pf3k,(p, k)=l. Then we have, (p/3 k)P + pbp/3 k- cP = 0 ,and hence pI c ,which is again a contradiction since(p, c) = 1 which completes the proof. As an special case of the theorem , consider the equation x3 + 15xb + 28 = 0 (3) which can be written as x3 +1+15xb+33 =0 (4) and it is clear that this equation has no integral root l = O(mod 3) since 1 is not divisible by 3. If this equation has an integral root k which is not divisible by 3 , then k3 + 1 + 15kb + 32 = 0 from which it follows that 31 b due to the Lemma(in case of negative c ) and is a contradiction . Therefore the equation has no integral roots. In case of p = 2 , it follows from the theorem that the equation x2 - 2 bx- q2 =0 where (2, q) = 1 = (2, b) 'has no integral roots. Again from the theorem it follows that x P - pcx - p f3p a P - bP = 0 , where (p, c) = 1 = ( b, p) and p is a prime, has no integral solutions .. In particular here, p f3p a 3, bP are two components o( F ermat triples. It is easy to deduce that this equation has no integral roots. This theorem may hold for some other useful forms of polynomial equations.Item Analytical proof of Fermat's last theorem for n=4(University of Kelaniya, 2008) Piyadasa, R.A.D.Fermat's last theorem for n = 4 is usually proved [1] using the famous mathematical tool of the method of infinite descent of F ermat. In this contribution, it will be shown that the parametric solution of the polynomial equation d4 = e4 + g4, (e, g) = 1 can be obtained using a simple mathematical technique and thereby the proof of the theorem can be done, without depending on the sophisticated structure of primitive Pythagorean triples of Fermat[1] given by X= 2lm, y =!2-m2 , z = !2 + m2 ' where l >m> 0 and l, m are of opposite parity. The main objective of this contribution is to introduce a new simple mathematical technique which may be very useful in some other problems as well. If the equation (1) has a non-trivial integral solution for (x, y, z) , then one of e, g is even and we can assume that d, e, g are positive. If gis even, (d2-g2)(d2+g2)=e4 and terms in the brackets are eo-prime and hence , one writes d2 +g2 =x4 d2-g2 =y4 (2a) (2b) From these two equations, we get 2d2 = x4 + y4. (2c) Therefore (x2 -d)(x2+ d)=( d-y2 )(d + y2) and it is easy to deduce that terms in the brackets on the left-hand side or on the right-hand side of this equation may have only factor 2 in common since all numbers are odd and x, d, y are eo-prime to one another. In the following, a new simple mathematical technique is used to obtain the parametric solution for x, y ,d, g from this single equation. If x2 -d = d-y2, 2d = x2 + y2 and therefore 4d2 = x4 + y4 + 2x2 y2 , which means d2 = x2y2, and it leads to a contradiction since (d, e) = 1. Similarly we can easily show thatx2-d:t=d+y2. Now, let (d-y2)= a (x2-d) ,to obtain x2-d = ba-1(d-y2), b where(a, b) = l.Then k (x2 +d )= (d+ y2), x2+ d = ab-1( d + y2) .Now, let us form the a following two simultaneous equations, to obtain, x2-d =ba-1(d-y2) (a) x2 + d = ab-1(d + y2) (b) (3) Proceedings of the Annual Research Symposium 2008- Faculty of Graduate Studies University of Kelaniya (4) Since ( d, y) = 1 , b 2 + a 2 = dk , where k has to be determined. Then , one easily obtains 2 2ab+b2-a2 a2+b2 . 2ab = a2 -b2 + y2 k , y = , d = .Now , from(3), It follows that k k 2abx2 = (a2 +b2)d +(a2 -b2)y2 = (a2 +b2)2 +(a2 -b2)(2ab+b 2 -a2) k (5) It is clear from (5) that a and b cannot be of opposite parity since then k2 x2 y2 cannot be either odd or even. Hence a and b are both odd. and therefore k2 = 4 or 41 k2. a2 +b2 Thereforex 2y2 =(4a2b2 -(a2 -b2)2)14 =e2, d = , ab (a2 -b2)=g2 as given 2 below , which is the parametric solution of the equation (1 ), where a, bare parameters. 2 2ab+a2-b2-a2-b2 2a (a-b) 2 2 Now, x -d = = and k is a factor of a +b and if k k it is a factor of a-b , one deduces that k is 2 or a factor of a or b. Since (a, b)= 1, we conclude that k = 2 . Therefore x2 y2 = a2 b2 -( a2 -b2 )2 /4 Since (x2-y2)(x2 + y2) = x4-y4 = 2g2, which follows from(2a),(2b), it IS easy to deduce (6) Therefore a, b, (a2-b2) should be perfect squares. Now, if a=r2, b =s2, then r4 -s4 = t2 for some integers r, s, t. The famous and the only theorem that Fermat has proved is that there are no integers r, s, t satisfying r4-s4 = t2. Hence the Fermat's last theorem for n = 4 can be deduced. It is quite interesting that applying the mathematical technique used in this contribution ,we have shown[2] very easily that the equation r4-s4 = t2 has no non- trivial integral solution for r,s,t ,and then the Fermat's last theorem for n = 4 follows at once[1].Item Simple and analytical proof of Fermat's last theorem for n=3(University of Kelaniya, 2008) Piyadasa, R.A.D.It is well known that Fermat's last Theorem, in general, is extremely difficult to prove although the meaning of the theorem is very simple. It is surprising that the proof of theorem for n = 3 ,the smallest, corresponding number, given by Leonard Euler, which has been recommended for amateurs[!], is not only difficult but also has a gap in the proof. Paulo Rebenboim claims[1] that he has patched up Euler's proof, which is very difficult to understand, however. It was shown[2] that the parametric solution for the x,y,z in the equation z3 = y3 +x3,(x,y) = 1 could be obtained easily with one necessary condition that must be satisfied by the parameters. Fairly simple analytical proof of the Fermat's last theorem for n = 3 was given [2] using this necessary condition. In this contribution very much simpler proof is given ,which is very suitable for amateurs. Fermat's last theorem for n = 3 can be stated as the equation z3 = y3 + x3, (x,y) ·= 1. (1) has no non-zero integral solution for (x,y,z). If we assume a non-zero( xyz * 0) solution for (x,y,z) , then one of (x,y,z) is divisible by 3 . Since we can assume for (1) for negative integers, without loss of generality one can assume that y is divisible by 3 . Then if y = 3fJ ay, the parametric solution of (1), can be expressed as X = 3fJ a()o + 03 y = 3fJ a()o + 33fJ-1a 3 z = 33fJ-1 a3 + 3fJ a()o + 83 and a necessary condition satisfied by the parameters is ()3 -83 -2.3fJ aoB -33[3-1a 3 =0 In this equation ,() is a factor of z, 8 is a factor of x and r = B5 + 32P-I a2 Proof of the Fermat's last theorem for n = 3 (a) (b) (c) (d) Expressing 33/3-I a3 as 33/3-3 a3 + 8.33/3-3 and substituting e = 3fJ-1 g + o m (d),one gets (g- 2a)(o2 +3fl-1go+32fl-3(g2 +2ag+4a2 ))=32fl-3a3 (2) the condition,()= 3/3-l g + 8 is due to a simple lemma used in [2] from which it follows that f3 > 1 .It is easy to deduce that g-2a is divisible by 3 2/3-3 since (3, o) = 1 and fJ > 1 . If a= ± 1 , then ,g = ±2+32P-3>0. Now, (82 + 3/3-l g8 + 32f3-\g2 +. 2ag + 4a2 ) = ±1 is never satisfied since (2) can be expressed as 3/3-1 2 (8 +-- g )2 +32P-3(L+2g +4)= ±1 2 4 (3) (4) If a* ±1, we deduces from (2) that g-2a= 3 213-3s 3 where s* 1 and s is a factor of a . This is because factor of g cannot be a factor of both 8 and a since B = 3/3-l g + 8 and a is a factor of y . Hence, (5) where a= sq, (s, q) = 1. Now, this quadratic equation ing must be satisfied by 2a +3 213-3s 3. If the roots of this quadratic are a1 and a 2 =2a + 32/3-3 s3, then 82 + 4a2 .32/3-3 _ q3 2a.32f3-3 + 3p-t8 ala 2 = 3213_3 , a1 + a 2 =- 3213_3 It is easy to obtain from these two relations that (6) Hence, 2a + 3213-3s 3 is a factor of 82 + 4.a2 .32/3-3- q3 . In other words, 2a +3 2/3-3s 3 is a factor of the expression (2a)3 -(2a)2 .8s33 2/3-3 -882s3 or, -32/3-3s 3 is an integral root of the equation x3 -8s33213-3x2 -882s3= 0 (7) (8) This means that 32/3-3s 3 is a factor of 882. s3 which contradicts (3, 8)= 1 , that is (3,x) = 1, and the proof of theorem is complete. It should be emphasized that the necessary condition needed for our proof can be obtained without obtaining the parametric solution of (1 ), making the proof given here much shorter.