Commun. Thor. Phys. 66 06 96 00 Vol. 66, No., August, 06 Scattring Stats of l-wav Schrödingr Equation with Modifid Rosn Mors Potntial Wn-Li Chn í,, Yan-Wi Shi á, and Gao-Fng Wi Ôô, Gnral Education Cntr, Xi an Pihua Univrsity, Xi an 70065, China Shaanxi Ky Laboratory of Surfac Enginring and Rmanufacturing, School of Mchanical and Matrial Enginring, Xi an Univrsity, Xi an 70065, China Rcivd Fbruary 9, 06; rvisd manuscript rcivd Jun, 06 Abstract Within a Pkris-typ approximation to th cntrifugal trm, w xamin th approximatly analytical scattring stat solutions of th l-wav Schrödingr quation with th modifid Rosn Mors potntial. Th calculation formula of phas shifts is drivd, and th corrsponding bound stat nrgy lvls ar also obtaind from th pols of th scattring amplitud. PACS numbrs: 03.65.-w, 03.65.Nk Ky words: modifid Rosn Mors potntial, Scattring stats, approximatly analytical solutions Introduction Th Mors and/or original Rosn Mors potntials, [ ] as th important molcular modl potntials, hav bn widly studid du to its validity in dscribing th vibration of actual diatomic molculs spcially for th lowst lctronic stats. [3] Nvrthlss, du to th Mors and/or original Rosn Mors potntials do not considr th innr-shll radii of two atoms, it is thrfor that thy ncountr invitably som problm in rproducing th xprimntal data for som molcul stats such as 6 Π u stat of 7 Li molcul and X 3 Π stat of SiC radical. [4 6] Vry rcntly, in ordr to fit th ffct of innr-shll radii of two atoms for diatomic molculs, [7] Zhang t al. [4] proposd basd on th original Rosn Mors potntial [] a modifid Rosn Mors potntial in th form of V r = D r rij +, r rij + whr paramtrs r ij = r KD /k, D is th dissociation nrgy, r is th quilibrium bond lngth, and k is th quilibrium harmonic vibrational forc constant, and th paramtr K is a dimnsional constant. [8] It is shown that th modifid Rosn Mors potntial is vry suprior to th Mors and original Rosn Mors potntials in fitting xprimntal data for som molcul stats. [4] Morovr, th vibrational nrgy lvls for 6 Π u stat of 7 Li molcul and X 3 Π stat of SiC radical with this potntial modl [5] ar good agrmnt with th xprimntal data dtrmind by th Rydbrg Klin Rss mthod. [9 ] Up to now, th bound stats of this potntial modl hav bn invstigatd within th non-rlativistic and/or rlativistic quantum mchanics. [ 4] Howvr, thr is no any rport on th scattring stats rgarding this potntial modl to our knowldg. Undoubtdly, th scattring stats spcially th phas shifts ar th important issu no mattr in quantum physics or nuclar physics. Thrfor, it is ncssary to study systmatically th scattring stats of this potntial spcially th analytical proprtis of scattring amplitud. Without considring th application of this potntial modl to a spcific molcular stat, w ar going to study th scattring stats of th modifid Rosn Mors potntial in this papr. On th othr hand, w ar going to study th proprtis of scattring stats to show th rlations btwn th nrgy lvls of scattring stats and thos of bound stats at th pols of scattring amplituds. This papr is organizd as follows. In Sc., w driv th scattring stat solutions of th modifid Rosn Mors potntial with cntrifugal trm. In Sc. 3, w discuss th analytical proprtis of scattring amplitud to find th discrt nrgy lvls from th pols of scattring amplitud. Th concluding rmarks ar givn in Sc. 4. Scattring Stats Th quantum bhaviour of a non-rlativistic particl in a tim-indpndnt potntial V r is dscribd by th stationary Schrödingr quation = µ = Ψr, θ, φ + V rψr, θ, φ = EΨr, θ, φ. Hr, th potntial is takn as th modifid Rosn Mors potntial, i.. Eq.. Taking Ψ nlm r, θ, φ = r u nl ry lm θ, φ, and substituting it into Eq., w Supportd by th National Natural Scinc Foundation of China undr Grant No. 4058, and Natural Scinc Basic Rsarch Plan in Shaanxi Provinc of China undr Grant No. 5JK093 E-mail: physwlchn@63.com Corrsponding author, E-mail: fgwi 000@63.com c 06 Chins Physical Socity and IOP Publishing Ltd http://www.iopscinc.iop.org/ctp http://ctp.itp.ac.cn
No. Communications in Thortical Physics 97 obtain th radial Schrödingr quation as follows, u ll + ] r+ [E D r r r ur = 0, 3 whr paramtr λ = rij. For scattring stats, at larg distanc th solutions of Eq. 3 ar asymptotically sinusoidal with phas shifts δ l that carry th scattring information. [5] Howvr, Eq. 3 cannot b solvd analytically du to th cntrifugal trm ll + /r. Thrfor, w hav to us a propr approximation to th cntrifugal trm. On th othr hand, it is noticd that this potntial has a minimum valu at r = r r 0. Thrfor, w can us th Pkris-typ [6 7] approximation to rplac th cntrifugal trm /r. To this nd, w introduc a variabl y = r r 0 /r 0 and xpand th cntrifugal trm in a powr sris around y = 0 r = r0 + = y r0 y + 3y + oy. 4 On th othr hand, w approximatly writ th cntrifugal trm as th following form r c 0 + c r + c r = c 0 + c r0+y + c r0+y. 5 Aftr xpanding th right hand of Eq. 5 in a powr sris around y = 0 and comparing it with Eq. 4, th cofficints c 0, c, and c can b xprssd as th following form Now w introduc a nw variabl [3 c 0 = r r + r r λ + r 3 + r ] r 4, c = r r [ λ3 + r + r 3 + r ] r 4, c = r r 3 [ r r 3 λ3 + r ] r 4. 6 x = r r r, x 7 and insrt Eq. 5 into Eq. 3, thn th radial Schrödingr quation 3 can b transformd into th following form { x x u x + x 3x + x u x + E c 0 ll + x D x [c ll + + r D ] λ + x [c ll + 4 r xd ] } ux = 0. 8 λ W now tak th trial wavfunction in th form of whr ux = x ik/ x η Fx, 9 c ll + + c l + c l η = λ + c 0 l + c 0 l Eλ + r D, 0 λ k = E c 0 l c 0 l D. Insrting Eq. 9 into Eq. 8 yilds th following hyprgomtric quation [ ik ] xxf x + + η + x + η F x k + ik + η η η + + c l + c l + r D λ Fx = 0, its gnral solution is nothing but th hyprgomtric function ã p b p x p Fx = F ã, b; c; x = c p p!, 3 whr x p is th Pochhammr symbol x p = p=0 Γx + p Γx, 4 and paramtrs ã, b and c ar dfind as follows ã = ik σ + η r, b = ik σ + η + r, c = + η, σ = 4c ll + + λ + 8 r D. 5 Now w can writ th analytical radial wav functions of arbitrary l-wav scattring stats for th modifid Rosn Mors potntial as [ ur = N r ] η ] [ r r ik/ r F ã, b; c; r r, 6
98 Communications in Thortical Physics Vol. 66 whr N is th normalizd constant. To dtrmin th normalizd constant N and scattring phas shifts δ l, w us th following rcurrnc formula of hyprgomtric functions [8] F ã, b; c; x = Γ cγ c ã b Γ c ãγ c b F ã, b; ã + b c + ; x + x c ã b Γ cγã + b c F c ã, c b; c ã b + ; x, 7 ΓãΓ b and notic th F a, b; c; 0 = and r / r λ r for r, thn w hav th following asymptotic bhavior for th hyprgomtric function F ã, b; r [ c; r Γ c Γ c ã b r Γ c ãγ c b ik/ ikr Γ c ã b ] Γ c ãγ c b [ = λ ik/ Γ c ã Γ c b Γ c ãγ c b λ ik/ + ikr Γ c ã b ] Γ c ãγ c b λ ik/. 8 By taking Γ c ã b Γ c ãγ c b = Γ c ã b Γ c ãγ c b iθ, λ ik/ = λ ik/ iθ, 9 and insrting it into Eq. 8, w hav F ã, b; r ik/ Γ c ã bλ c; r r λik/ Γ c Γ c ãγ c b ikr [ ikr+θ+θ + ikr+θ+θ ] = NΓ cλ ik/ ikr ik/ Γ c ã bλ Γ c ãγ c b sin kr + π + θ + θ. 0 Aftr insrting Eq. 0 into Eq. 6 and considring r / r and r / r λ r whn r, w can finally xprss th asymptotic form of th radial wav function 6 for r as follows ik/ Γ c ã bλ ur NΓ c r Γ c ãγ c b sin kr + π + θ + θ. Comparing Eq. with th gnral asymptotic solution [5] ur r sinkr /lπ + δ l, w gt th phas shifts and th normalizd constant as th following forms δ l = π l + + θ + θ = π l + + arg Γ c ã b arg Γ c ã argγ c b + argλ ik/, N = ik/ Γ c ãγ c bλ Γ c Γ c ã b. 3 Th substitution of paramtrs of ã, b and c into Eqs. and 3 allows us to xprss th phas shifts and normalizd constant in th following final forms δ l = π l + + θ + θ = π ik l + + argγ + argλ ik/ argγ + ik + η + 4c ll + + λ + 8 r D λ argγ + ik + η 4c ll + + λ + 8 r D, 4 λ N = λik/ Γ Γ + η + ik + η + 4c ll + + λ + 8 r D λ Γ ik Γ + ik + η 4c ll + + λ + 8 r D. 5 λ 3 Analytical Proprtis of Scattring Amplitud Now lt us study th analytical proprtis of th scattring amplitud. It is wll known that th discrt nrgy lvls of th bound stats corrspond to th singl pols of th scattring amplitud. [5] This suggsts that w discuss
No. Communications in Thortical Physics 99 th proprtis of Gamma function Γ + ik + η 4c ll + + λ + 8 r D λ basd on Eq. 4. From th dfinition of th Gamma function Γz = Γz + z = 6 Γz + zz + = Γz + 3 zz + z + =, 7 w know that z = 0,,,... ar th first ordr pols of th Γz. Thrfor, th first ordr pols of Γ + ik + η 4c ll + + λ + 8 r D λ ar situatd at + ik + η 4c ll + + λ + 8 r D = 0,,, 3,... = n r, n r = 0,,... 8 λ from which, w can obtain th corrsponding bound stat nrgy lvls as E nrl = { c 0ll + + D nr + + 4c ll + + 8 r D / λ 4 [ll + c + c λ + r λ } D ] λ n r +, 9 + 4c ll + + 8 r D / λ whr paramtrs c 0, c, and c ar givn by Eqs. 6. By comparing th abov nrgy lvls 9 with th nrgy quation 5 in Rf. [5] and obsrving th rlations of corrsponding paramtrs, it is asy to find that Eq. 9 is in ssnc th sam as Eq. 5 in Rf. [5]. That is to say, th pols of th S-matrix in th complx nrgy plan corrspond to bound stats for ral pols and scattring stats for complx pols in th lowr half of th nrgy plan. Without considring th application of this potntial modl to a spcific molcular stat, w now prsnt som numrical nrgy lvls calculatd from nrgy quation 9 and compar thm with thos obtaind by th MATHEMATICA packag programmd by Lucha and Schöbrl. [9] To this nd, w tak paramtrs of λ= and D =000 in nrgy quation 9 and th MATHEMAT- ICA packag to gt th numrical nrgy lvls as shown in Tabl. It can b found that our prsnt rsults ar in good agrmnt with thos from th MATHEMATICA packag, vn for larg potntial rang paramtr. Tabl Eignvalus E nrl as a function of for p, 3p, 3d, 4p, 4d, 4f, 5p, 5d, 5f, 5g, 6f and 6g stats for λ =, D = 000 and = µ =. Stats r =.6 r =.4 r = 3. Prsnt Lucha t al. Prsnt Lucha t al. Prsnt Lucha t al. p 0.5.88 566 3.304 459 9.63 8.65 36 5.9 573 6.9 8 3 0.0 3.009 30 5 3.08 40 4 3.005 683 3 3.006 46 3 3.75 73 9 3.75 889 0.30 4.65 37 7 4.654 93 8 4.979 69 6 4.979 99 3 5.465 664 8 5.465 87 6 0.45 7.603 754 8 7.605 94 8.637 048 6 8.637 546 5 9.738 866 4 9.739 440 8 0.60. 950.3 097 3.0 586 3.03 69 4.740 3 4.74 453 3p 0.5 6.097 367 6.0 76 0 6.4 369 4 6.50 645 6.460 644 6.46 967 7 0.0 8.49 500 8 8.307 93 5 8.665 98 5 8.669 76 7 9.36 966 8 9.37 796 0.30 3.64 65 3.85 0 4.580 63 4.58 53 6.89 468 6.90 56 0.45.005 88.03 57 5.534 07 5.536 580 8.979 763 8.98 753 0.60 3.536 099 3.54 95 38.646 674 38.65 005 43.9 59 43.99 36 3d 0.5 3.34 955 3.35 74.59 795 4.53 53 5.4 33.4 704 0.0 3.849 03 8 3.857 896 6 3.367 04 3 3.368 67 9 3.375 37 5 3.375 695 0 0.30 5.475 858 6 5.48 467 8 5.335 45 7 5.335 970 7 5.663 70 5 5.663 43 7 0.45 8.4 766 9 8.44 774 8.988 0 8.988 804 9.934 858 3 9.935 447 0 0.60.99 69.9 538 3.37 386 3.37 446 4.935 378 4.936 705 4p 0.5 9.899 76 6 0.39 300 0.4 83 0.39 60 0.695 899 0.699 569 0.0 3.478 678 3.669 705 4.34 374 4.35 98 5.468 79 5.470 390 0.30.653 899.76 770 4.59 4.64 30 6.889 568 6.89 35 0.45 36.357 7 36.378 57 4.374 734 4.38 384 48.50 409 48.58 9 0.60 53.855 356 53.870 803 64.45 706 64.59 533 7.93 7.930 95
00 Communications in Thortical Physics Vol. 66 Tabl continud Stats r =.6 r =.4 r = 3. Prsnt Lucha t al. Prsnt Lucha t al. Prsnt Lucha t al. 4d 0.5 7.88 583 7.67 598 9 6.556 3 4 6.574 7 9 6.676 340 6.679 645 7 0.0 9.70 7 9.355 58 5 9.059 976 9.069 950 3 9.536 043 5 9.537 800 0.30 4.00 698 4.46 096 4.954 835 4.958 777 6.39 66 6.39 9 0.45.876 39.893 766 5.893 638 5.896 676 9.77 88 9.80 0.60 33.370 893 33.380 54 38.998 865 39.004 97 44.07 670 44.4 50 4f 0.30 6.689 7 9 6.69 600 4 5.866 966 9 5.868 06 7 5.959 85 9 5.959 500 5 0.45 9.65 7 8 9.69 030 9.54 403 9 9.55 5 0.8 770 0.9 37 0.60 3. 43 3.4 898 3.894 349 3.895 449 5.8 5 5.9 540 5p 0.30 30.0 67 30.6 30 33.75 098 33.75 34 37.565 954 37.57 70 0.45 50.657 979 50.70 783 59.59 03 59.7 89 67.50 6 67.65 53 0.60 75.079 639 75.09 764 89.59 47 89.545 683 0.74 50 0.775 8 5d 0.30.70 307.846 3 4.55 9 4.56 88 7.096 355 7.099 787 0.45 37.90 359 37.340 554 4.74 78 4.750 737 48.349 59 48.357 3 0.60 54.77 57 54.753 50 64.50 64 64.55 364 73.08 50 73.5 953 5f 0.45 4.73 8 4.99 369 6.43 496 6.436 54 9.473 49 9.476 36 0.60 34.69 06 34.63 883 39.56 9 39.53 487 44.400 349 44.407 78 5g 0.45.04 65.06 693 0.4 976 0.5 793 0.60 54 0.6 3 0.60 4.695 640 4.698 37 4.59 9 4.59 39 5.68 604 5.69 93 6f 0.60 56.030 686 56.069 963 65.034 360 65.048 85 73.400 758 73.48 444 6g 0.60 36.75 899 36.9 73 40.30 58 40.36 4 44.790 58 44.797 346 4 Conclusions In summary, within th framwork of non-rlativistic quantum thory, w hav studid th scattring stats of th modifid Rosn Mors potntial. By using a Pkris-typ approximation to th cntrifugal trm, th normalizd analytical wav functions ar obtaind and corrsponding calculation formula of phas shifts is drivd. Morovr, w hav discussd th analytical proprtis of scattring amplitud and found th discrt nrgy lvls of th bound stats from th pols of th scattring amplitud. Rfrncs [] P.M. Mors, Phys. Rv. 34 99 57. [] N. Rosn and P.M. Mors, Phys. Rv. 4 93 0. [3] G.F. Wi and W.L. Chn, Chin. Phys. B 9 00 090308. [4] G.D. Zhang, J.Y. Liu, L.H. Zhang, W. Zhou, and C.S. Jia, Phys. Rv. A 86 0 0650. [5] H.M. Tang, G.C. Liang, L.H. Zhang, F. Zhao, and C.S. Jia, Can. J. Chm. 9 04 34. [6] C.S. Jia, X.L. Png, and S. H, Bull, Koran Chm. Soc. 35 04 699. [7] A.A. Frost and B. Musulin, J. Am. Chm. Soc. 76 954 045. [8] C.S. Jia and Z.W. Shui, Eur. Phys. J. A 5 05 44. [9] R. Rydbrg, Z. Phys. 80 933 54. [0] O. Klin, Z. Phys. 76 93 6. [] A.L.G. Rs, Proc. Phys. Soc. 59 947 998. [] C.S. Jia, L.H. Zhang, X.T. Hu, H.M. Tang, and G.C. Liang, J. Mol. Spctrosc. 3 05 69. [3] C.S. Jia, J.W. Dai, L.H. Zhang, J.Y. Liu, and G.D. Zhang, Chm. Phys. Ltt. 69 05 54. [4] C.S. Jia, J.W. Dai, L.H. Zhang, J.Y. Liu, and X.L. Png, Phys. Ltt. A 379 05 37. [5] L.D. Landau and E.M. Lifshitz, Quantum Mchanics, Non-Rlativisitc Thory, 3rd d. Prgamon, Nw York 977. [6] C.L. Pkris, Phys. Rv. 45 934 98. [7] G.F. Wi and S.H. Dong, Phys. Ltt. B 686 00 88. [8] I.S. Gradshtyn and I.M. Ryzhik, Tabls of Intgals, Sris, and Products, 5th d. Acadmic Prss, Nw York 994. [9] W. Lucha and F.F. Schöbrl, Int. J. Mod. Phys. C 0 999 607.