subroutine fit3d(r12,r13,r23,e,der) ************************************************************* * PES of the excited triplet state ( 1 ^3A', first of its symmetry) of OH+F system * * S. Gomez-Carrasco, O. Roncero, L. Gonzalez-Sanchez, M. Luz Hernandez, J.M. Alvariño, M. Paniagua and A. Aguado * J. Chem. Phys. 123 (2005), 114310 * * INPUT/OUTPUT * r12 --> OH * r13 --> OF * r23 --> HF * e : energy * der(3): derivatives * * All units arr in atomic units ************************************************************ implicit real * 8 (a-h,o-z) dimension der(3) call diat12(r12,e12,d12) call diat13(r13,e13,d13) call diat23(r23,e23,d23) call triabc(r12,r13,r23,e123,der) e=e12+e13+e23+e123 der(1)=d12+der(1) der(2)=d13+der(2) der(3)=d23+der(3) return end ************************************************************************ subroutine diat12(r,ener,der) ************************************************************************ * This subroutine computes the energies of a diatomic potential * fitted to 99 points * rms = 0.03663237 kcal/mol * emax = 0.15496720 kcal/mol ************************************************************************ implicit real*8 (a-h,o-z) dimension cf( 7) data cf( 1)/0.166099112432E+02/ data cf( 2)/-.234610012185E-01/ data cf( 3)/-.903833157870E+00/ data cf( 4)/-.520290888529E+02/ data cf( 5)/0.360070444062E+03/ data cf( 6)/-.975763500496E+03/ data cf( 7)/0.102060053232E+04/ e0= 0.0000000E+00 der=0.d0 vex1= 0.1038258E+01 vex2= 0.3230380E+01 aux = 1.d0/r bux = dexp(-vex2*r)*aux cux = dexp(-vex1*r) ener=e0+cf(1)*bux dux=1.d0 eux=r*cux do 1 i=2, 7 der=der+(i-1)*cf(i)*dux dux=dux*eux ener=ener+cf(i)*dux 1 continue der=der*(1.d0-vex1*r)*cux der=der-cf(1)*(vex2+aux)*bux return end ************************************************************************ subroutine diat23(r,ener,der) ************************************************************************ * This subroutine computes the energies of a diatomic potential * fitted to 76 points * rms = 0.05060986 kcal/mol * emax = 0.17810149 kcal/mol ************************************************************************ implicit real*8 (a-h,o-z) dimension cf( 7) data cf( 1)/0.501366364227E+01/ data cf( 2)/-.380767904959E-01/ data cf( 3)/-.250639850172E+02/ data cf( 4)/0.173181530914E+03/ data cf( 5)/-.697554822298E+03/ data cf( 6)/0.133214204883E+04/ data cf( 7)/-.105113999439E+04/ e0= 0.0000000E+00 der=0.d0 vex1= 0.1212931E+01 vex2= 0.1543465E+01 aux = 1.d0/r bux = dexp(-vex2*r)*aux cux = dexp(-vex1*r) ener=e0+cf(1)*bux dux=1.d0 eux=r*cux do 1 i=2, 7 der=der+(i-1)*cf(i)*dux dux=dux*eux ener=ener+cf(i)*dux 1 continue der=der*(1.d0-vex1*r)*cux der=der-cf(1)*(vex2+aux)*bux return end ************************************************************************ subroutine diat13(r,ener,der) ************************************************************************ * This subroutine computes the energies of a diatomic potential * fitted to 54 points * rms = 0.14129135 kcal/mol * emax = 0.46773777 kcal/mol ************************************************************************ implicit real*8 (a-h,o-z) dimension cf( 7) data cf( 1)/0.328464606788E+03/ data cf( 2)/-.751930083968E-01/ data cf( 3)/0.389150997186E+01/ data cf( 4)/-.193735989756E+03/ data cf( 5)/0.160255155466E+04/ data cf( 6)/-.523667435431E+04/ data cf( 7)/0.742190036246E+04/ e0= 0.0000000E+00 der=0.d0 vex1= 0.1094583E+01 vex2= 0.4382704E+01 aux = 1.d0/r bux = dexp(-vex2*r)*aux cux = dexp(-vex1*r) ener=e0+cf(1)*bux dux=1.d0 eux=r*cux do 1 i=2, 7 der=der+(i-1)*cf(i)*dux dux=dux*eux ener=ener+cf(i)*dux 1 continue der=der*(1.d0-vex1*r)*cux der=der-cf(1)*(vex2+aux)*bux return end ************************************************************* subroutine triabc(r12,r13,r23,ener,der) ************************************************************* * This subroutine computes the energies of a 3D PES * for the ABC system class fitted to 7372 points * rms = 0.84399196 kcal/mol * emax = 5.91125528 kcal/mol ************************************************************* implicit real*8(a-h,o-z) dimension i1( 192),i2( 192),i3( 192),cf( 192) dimension f12(-1: 8),f13(-1: 8),f23(-1: 8) dimension der(3) data cf( 1)/0.6852261766779042E+01/ data i1( 1)/ 0/,i2( 1)/ 1/,i3( 1)/ 1/ data cf( 2)/-.3824344395727603E+00/ data i1( 2)/ 1/,i2( 2)/ 0/,i3( 2)/ 1/ data cf( 3)/-.8867357748869557E+01/ data i1( 3)/ 1/,i2( 3)/ 1/,i3( 3)/ 0/ data cf( 4)/-.1097073994299347E+03/ data i1( 4)/ 0/,i2( 4)/ 1/,i3( 4)/ 2/ data cf( 5)/-.1797578895680513E+03/ data i1( 5)/ 0/,i2( 5)/ 2/,i3( 5)/ 1/ data cf( 6)/0.2196580964093944E+02/ data i1( 6)/ 1/,i2( 6)/ 0/,i3( 6)/ 2/ data cf( 7)/0.1773539530772989E+03/ data i1( 7)/ 1/,i2( 7)/ 1/,i3( 7)/ 1/ data cf( 8)/0.2421311613837463E+03/ data i1( 8)/ 1/,i2( 8)/ 2/,i3( 8)/ 0/ data cf( 9)/0.1318302483659165E+02/ data i1( 9)/ 2/,i2( 9)/ 0/,i3( 9)/ 1/ data cf( 10)/0.3716109902220446E+03/ data i1( 10)/ 2/,i2( 10)/ 1/,i3( 10)/ 0/ data cf( 11)/0.2178043510455885E+04/ data i1( 11)/ 0/,i2( 11)/ 1/,i3( 11)/ 3/ data cf( 12)/0.5388212547328319E+04/ data i1( 12)/ 0/,i2( 12)/ 2/,i3( 12)/ 2/ data cf( 13)/0.2937569723449801E+04/ data i1( 13)/ 0/,i2( 13)/ 3/,i3( 13)/ 1/ data cf( 14)/-.4861235940698659E+03/ data i1( 14)/ 1/,i2( 14)/ 0/,i3( 14)/ 3/ data cf( 15)/-.2072986704161895E+04/ data i1( 15)/ 1/,i2( 15)/ 1/,i3( 15)/ 2/ data cf( 16)/0.8491152241816341E+02/ data i1( 16)/ 1/,i2( 16)/ 2/,i3( 16)/ 1/ data cf( 17)/-.5434332886328948E+04/ data i1( 17)/ 1/,i2( 17)/ 3/,i3( 17)/ 0/ data cf( 18)/-.9937078480550882E+02/ data i1( 18)/ 2/,i2( 18)/ 0/,i3( 18)/ 2/ data cf( 19)/-.2640233573267568E+04/ data i1( 19)/ 2/,i2( 19)/ 1/,i3( 19)/ 1/ data cf( 20)/-.7792590684080530E+04/ data i1( 20)/ 2/,i2( 20)/ 2/,i3( 20)/ 0/ data cf( 21)/-.1166775840021859E+03/ data i1( 21)/ 3/,i2( 21)/ 0/,i3( 21)/ 1/ data cf( 22)/-.3501743516089569E+04/ data i1( 22)/ 3/,i2( 22)/ 1/,i3( 22)/ 0/ data cf( 23)/-.2786835423743504E+05/ data i1( 23)/ 0/,i2( 23)/ 1/,i3( 23)/ 4/ data cf( 24)/0.2517399810438231E+05/ data i1( 24)/ 0/,i2( 24)/ 2/,i3( 24)/ 3/ data cf( 25)/-.1231577084682017E+06/ data i1( 25)/ 0/,i2( 25)/ 3/,i3( 25)/ 2/ data cf( 26)/-.1114829843557744E+06/ data i1( 26)/ 0/,i2( 26)/ 4/,i3( 26)/ 1/ data cf( 27)/0.3921776442811302E+04/ data i1( 27)/ 1/,i2( 27)/ 0/,i3( 27)/ 4/ data cf( 28)/0.1533882746679526E+05/ data i1( 28)/ 1/,i2( 28)/ 1/,i3( 28)/ 3/ data cf( 29)/0.1015316156984282E+05/ data i1( 29)/ 1/,i2( 29)/ 2/,i3( 29)/ 2/ data cf( 30)/0.5676084258349147E+05/ data i1( 30)/ 1/,i2( 30)/ 3/,i3( 30)/ 1/ data cf( 31)/0.1228963098116260E+06/ data i1( 31)/ 1/,i2( 31)/ 4/,i3( 31)/ 0/ data cf( 32)/0.3193248063796376E+04/ data i1( 32)/ 2/,i2( 32)/ 0/,i3( 32)/ 3/ data cf( 33)/-.5864402685091711E+04/ data i1( 33)/ 2/,i2( 33)/ 1/,i3( 33)/ 2/ data cf( 34)/-.1154055804852693E+05/ data i1( 34)/ 2/,i2( 34)/ 2/,i3( 34)/ 1/ data cf( 35)/-.6892204112535980E+05/ data i1( 35)/ 2/,i2( 35)/ 3/,i3( 35)/ 0/ data cf( 36)/0.2402358203708709E+04/ data i1( 36)/ 3/,i2( 36)/ 0/,i3( 36)/ 2/ data cf( 37)/0.2913743623198661E+05/ data i1( 37)/ 3/,i2( 37)/ 1/,i3( 37)/ 1/ data cf( 38)/0.1434070575655124E+06/ data i1( 38)/ 3/,i2( 38)/ 2/,i3( 38)/ 0/ data cf( 39)/0.9619313230107892E+03/ data i1( 39)/ 4/,i2( 39)/ 0/,i3( 39)/ 1/ data cf( 40)/0.9311449862141881E+04/ data i1( 40)/ 4/,i2( 40)/ 1/,i3( 40)/ 0/ data cf( 41)/0.1565390347221617E+06/ data i1( 41)/ 0/,i2( 41)/ 1/,i3( 41)/ 5/ data cf( 42)/-.2735333345889045E+06/ data i1( 42)/ 0/,i2( 42)/ 2/,i3( 42)/ 4/ data cf( 43)/0.2591229372382639E+06/ data i1( 43)/ 0/,i2( 43)/ 3/,i3( 43)/ 3/ data cf( 44)/0.4980072332775927E+06/ data i1( 44)/ 0/,i2( 44)/ 4/,i3( 44)/ 2/ data cf( 45)/0.2039681031931866E+07/ data i1( 45)/ 0/,i2( 45)/ 5/,i3( 45)/ 1/ data cf( 46)/-.1287239795662128E+05/ data i1( 46)/ 1/,i2( 46)/ 0/,i3( 46)/ 5/ data cf( 47)/-.7489255489946195E+05/ data i1( 47)/ 1/,i2( 47)/ 1/,i3( 47)/ 4/ data cf( 48)/0.9841936357743079E+05/ data i1( 48)/ 1/,i2( 48)/ 2/,i3( 48)/ 3/ data cf( 49)/0.6810031064196413E+06/ data i1( 49)/ 1/,i2( 49)/ 3/,i3( 49)/ 2/ data cf( 50)/-.1431057501361903E+07/ data i1( 50)/ 1/,i2( 50)/ 4/,i3( 50)/ 1/ data cf( 51)/-.9619631501269611E+06/ data i1( 51)/ 1/,i2( 51)/ 5/,i3( 51)/ 0/ data cf( 52)/-.4396272251556070E+05/ data i1( 52)/ 2/,i2( 52)/ 0/,i3( 52)/ 4/ data cf( 53)/0.1121025378394929E+06/ data i1( 53)/ 2/,i2( 53)/ 1/,i3( 53)/ 3/ data cf( 54)/-.6260047103285310E+06/ data i1( 54)/ 2/,i2( 54)/ 2/,i3( 54)/ 2/ data cf( 55)/-.1430707831951109E+06/ data i1( 55)/ 2/,i2( 55)/ 3/,i3( 55)/ 1/ data cf( 56)/0.1148061713146735E+07/ data i1( 56)/ 2/,i2( 56)/ 4/,i3( 56)/ 0/ data cf( 57)/-.3314977006481269E+04/ data i1( 57)/ 3/,i2( 57)/ 0/,i3( 57)/ 3/ data cf( 58)/0.4604280170527119E+05/ data i1( 58)/ 3/,i2( 58)/ 1/,i3( 58)/ 2/ data cf( 59)/0.2312465362896609E+06/ data i1( 59)/ 3/,i2( 59)/ 2/,i3( 59)/ 1/ data cf( 60)/-.3732767577562881E+06/ data i1( 60)/ 3/,i2( 60)/ 3/,i3( 60)/ 0/ data cf( 61)/-.1591108863305906E+05/ data i1( 61)/ 4/,i2( 61)/ 0/,i3( 61)/ 2/ data cf( 62)/-.1453247935426672E+06/ data i1( 62)/ 4/,i2( 62)/ 1/,i3( 62)/ 1/ data cf( 63)/-.7784609019256184E+06/ data i1( 63)/ 4/,i2( 63)/ 2/,i3( 63)/ 0/ data cf( 64)/-.8502192976119250E+04/ data i1( 64)/ 5/,i2( 64)/ 0/,i3( 64)/ 1/ data cf( 65)/0.1477826132323723E+05/ data i1( 65)/ 5/,i2( 65)/ 1/,i3( 65)/ 0/ data cf( 66)/-.4175666815424968E+06/ data i1( 66)/ 0/,i2( 66)/ 1/,i3( 66)/ 6/ data cf( 67)/0.1359260123048554E+06/ data i1( 67)/ 0/,i2( 67)/ 2/,i3( 67)/ 5/ data cf( 68)/0.4034649368368340E+07/ data i1( 68)/ 0/,i2( 68)/ 3/,i3( 68)/ 4/ data cf( 69)/-.8267007728951159E+07/ data i1( 69)/ 0/,i2( 69)/ 4/,i3( 69)/ 3/ data cf( 70)/0.3378356190140494E+07/ data i1( 70)/ 0/,i2( 70)/ 5/,i3( 70)/ 2/ data cf( 71)/-.1209723151473749E+08/ data i1( 71)/ 0/,i2( 71)/ 6/,i3( 71)/ 1/ data cf( 72)/0.1732914245605234E+05/ data i1( 72)/ 1/,i2( 72)/ 0/,i3( 72)/ 6/ data cf( 73)/0.3690726711391444E+06/ data i1( 73)/ 1/,i2( 73)/ 1/,i3( 73)/ 5/ data cf( 74)/-.8465960410072956E+06/ data i1( 74)/ 1/,i2( 74)/ 2/,i3( 74)/ 4/ data cf( 75)/-.1409941682268315E+07/ data i1( 75)/ 1/,i2( 75)/ 3/,i3( 75)/ 3/ data cf( 76)/-.2702380512339350E+07/ data i1( 76)/ 1/,i2( 76)/ 4/,i3( 76)/ 2/ data cf( 77)/0.3612155239662247E+07/ data i1( 77)/ 1/,i2( 77)/ 5/,i3( 77)/ 1/ data cf( 78)/-.7493433507206135E+05/ data i1( 78)/ 1/,i2( 78)/ 6/,i3( 78)/ 0/ data cf( 79)/0.2013249493656697E+06/ data i1( 79)/ 2/,i2( 79)/ 0/,i3( 79)/ 5/ data cf( 80)/-.9334316619316234E+06/ data i1( 80)/ 2/,i2( 80)/ 1/,i3( 80)/ 4/ data cf( 81)/0.2906821554296402E+07/ data i1( 81)/ 2/,i2( 81)/ 2/,i3( 81)/ 3/ data cf( 82)/-.2961849169306458E+07/ data i1( 82)/ 2/,i2( 82)/ 3/,i3( 82)/ 2/ data cf( 83)/0.1173166939082443E+08/ data i1( 83)/ 2/,i2( 83)/ 4/,i3( 83)/ 1/ data cf( 84)/0.7045812902993235E+06/ data i1( 84)/ 2/,i2( 84)/ 5/,i3( 84)/ 0/ data cf( 85)/0.1162771332287755E+05/ data i1( 85)/ 3/,i2( 85)/ 0/,i3( 85)/ 4/ data cf( 86)/-.1198259780925069E+06/ data i1( 86)/ 3/,i2( 86)/ 1/,i3( 86)/ 3/ data cf( 87)/0.1543679755537522E+07/ data i1( 87)/ 3/,i2( 87)/ 2/,i3( 87)/ 2/ data cf( 88)/-.3241497230602823E+07/ data i1( 88)/ 3/,i2( 88)/ 3/,i3( 88)/ 1/ data cf( 89)/-.9245658463421097E+07/ data i1( 89)/ 3/,i2( 89)/ 4/,i3( 89)/ 0/ data cf( 90)/0.6449777074674296E+05/ data i1( 90)/ 4/,i2( 90)/ 0/,i3( 90)/ 3/ data cf( 91)/-.2761217289553591E+06/ data i1( 91)/ 4/,i2( 91)/ 1/,i3( 91)/ 2/ data cf( 92)/0.7867568822470627E+05/ data i1( 92)/ 4/,i2( 92)/ 2/,i3( 92)/ 1/ data cf( 93)/0.5522599900426950E+07/ data i1( 93)/ 4/,i2( 93)/ 3/,i3( 93)/ 0/ data cf( 94)/0.3360240965071740E+05/ data i1( 94)/ 5/,i2( 94)/ 0/,i3( 94)/ 2/ data cf( 95)/0.2690836292860286E+06/ data i1( 95)/ 5/,i2( 95)/ 1/,i3( 95)/ 1/ data cf( 96)/0.1595697967711335E+07/ data i1( 96)/ 5/,i2( 96)/ 2/,i3( 96)/ 0/ data cf( 97)/0.3936004198431874E+05/ data i1( 97)/ 6/,i2( 97)/ 0/,i3( 97)/ 1/ data cf( 98)/-.1077584974629257E+06/ data i1( 98)/ 6/,i2( 98)/ 1/,i3( 98)/ 0/ data cf( 99)/0.5421621450810836E+06/ data i1( 99)/ 0/,i2( 99)/ 1/,i3( 99)/ 7/ data cf(100)/0.1966253087587640E+07/ data i1(100)/ 0/,i2(100)/ 2/,i3(100)/ 6/ data cf(101)/-.1750737162001819E+08/ data i1(101)/ 0/,i2(101)/ 3/,i3(101)/ 5/ data cf(102)/0.2258383487827523E+08/ data i1(102)/ 0/,i2(102)/ 4/,i3(102)/ 4/ data cf(103)/-.1140482688169004E+07/ data i1(103)/ 0/,i2(103)/ 5/,i3(103)/ 3/ data cf(104)/0.1262452110542107E+08/ data i1(104)/ 0/,i2(104)/ 6/,i3(104)/ 2/ data cf(105)/-.2913294913437732E+07/ data i1(105)/ 0/,i2(105)/ 7/,i3(105)/ 1/ data cf(106)/-.1644725979465426E+05/ data i1(106)/ 1/,i2(106)/ 0/,i3(106)/ 7/ data cf(107)/-.1193906354254684E+07/ data i1(107)/ 1/,i2(107)/ 1/,i3(107)/ 6/ data cf(108)/0.1227274805993167E+07/ data i1(108)/ 1/,i2(108)/ 2/,i3(108)/ 5/ data cf(109)/0.1847910307377122E+08/ data i1(109)/ 1/,i2(109)/ 3/,i3(109)/ 4/ data cf(110)/-.3594103415995298E+08/ data i1(110)/ 1/,i2(110)/ 4/,i3(110)/ 3/ data cf(111)/0.8265461648934915E+07/ data i1(111)/ 1/,i2(111)/ 5/,i3(111)/ 2/ data cf(112)/0.2882830357625307E+08/ data i1(112)/ 1/,i2(112)/ 6/,i3(112)/ 1/ data cf(113)/0.2028072332452891E+08/ data i1(113)/ 1/,i2(113)/ 7/,i3(113)/ 0/ data cf(114)/-.3202624632733013E+06/ data i1(114)/ 2/,i2(114)/ 0/,i3(114)/ 6/ data cf(115)/0.3487063093711384E+07/ data i1(115)/ 2/,i2(115)/ 1/,i3(115)/ 5/ data cf(116)/-.1082039145629694E+08/ data i1(116)/ 2/,i2(116)/ 2/,i3(116)/ 4/ data cf(117)/0.4394705604203492E+07/ data i1(117)/ 2/,i2(117)/ 3/,i3(117)/ 3/ data cf(118)/0.4467215583239708E+08/ data i1(118)/ 2/,i2(118)/ 4/,i3(118)/ 2/ data cf(119)/-.4872079607148981E+08/ data i1(119)/ 2/,i2(119)/ 5/,i3(119)/ 1/ data cf(120)/-.2473154242456895E+08/ data i1(120)/ 2/,i2(120)/ 6/,i3(120)/ 0/ data cf(121)/-.2517260385267277E+06/ data i1(121)/ 3/,i2(121)/ 0/,i3(121)/ 5/ data cf(122)/-.1667918780357699E+06/ data i1(122)/ 3/,i2(122)/ 1/,i3(122)/ 4/ data cf(123)/0.3811631935374358E+06/ data i1(123)/ 3/,i2(123)/ 2/,i3(123)/ 3/ data cf(124)/-.9514426509255016E+07/ data i1(124)/ 3/,i2(124)/ 3/,i3(124)/ 2/ data cf(125)/-.3948629192600297E+08/ data i1(125)/ 3/,i2(125)/ 4/,i3(125)/ 1/ data cf(126)/0.3898762315123350E+08/ data i1(126)/ 3/,i2(126)/ 5/,i3(126)/ 0/ data cf(127)/0.6652773053979500E+05/ data i1(127)/ 4/,i2(127)/ 0/,i3(127)/ 4/ data cf(128)/0.7837071533570838E+06/ data i1(128)/ 4/,i2(128)/ 1/,i3(128)/ 3/ data cf(129)/-.2925539210427737E+07/ data i1(129)/ 4/,i2(129)/ 2/,i3(129)/ 2/ data cf(130)/0.2523967223395173E+08/ data i1(130)/ 4/,i2(130)/ 3/,i3(130)/ 1/ data cf(131)/0.5716330294967476E+07/ data i1(131)/ 4/,i2(131)/ 4/,i3(131)/ 0/ data cf(132)/-.1620133521072305E+06/ data i1(132)/ 5/,i2(132)/ 0/,i3(132)/ 3/ data cf(133)/0.4276724131664031E+06/ data i1(133)/ 5/,i2(133)/ 1/,i3(133)/ 2/ data cf(134)/-.4283859458459591E+07/ data i1(134)/ 5/,i2(134)/ 2/,i3(134)/ 1/ data cf(135)/-.1377878984599167E+08/ data i1(135)/ 5/,i2(135)/ 3/,i3(135)/ 0/ data cf(136)/-.5175883725204592E+05/ data i1(136)/ 6/,i2(136)/ 0/,i3(136)/ 2/ data cf(137)/0.1559281282913527E+06/ data i1(137)/ 6/,i2(137)/ 1/,i3(137)/ 1/ data cf(138)/-.1015045042967999E+07/ data i1(138)/ 6/,i2(138)/ 2/,i3(138)/ 0/ data cf(139)/-.7902549609591832E+05/ data i1(139)/ 7/,i2(139)/ 0/,i3(139)/ 1/ data cf(140)/0.1553517712259731E+06/ data i1(140)/ 7/,i2(140)/ 1/,i3(140)/ 0/ data cf(141)/-.3068972503998835E+06/ data i1(141)/ 0/,i2(141)/ 1/,i3(141)/ 8/ data cf(142)/-.2606881069293094E+07/ data i1(142)/ 0/,i2(142)/ 2/,i3(142)/ 7/ data cf(143)/0.1403139252997388E+08/ data i1(143)/ 0/,i2(143)/ 3/,i3(143)/ 6/ data cf(144)/0.1957534613570661E+08/ data i1(144)/ 0/,i2(144)/ 4/,i3(144)/ 5/ data cf(145)/-.1096697384425099E+09/ data i1(145)/ 0/,i2(145)/ 5/,i3(145)/ 4/ data cf(146)/0.1529420881180830E+09/ data i1(146)/ 0/,i2(146)/ 6/,i3(146)/ 3/ data cf(147)/-.1586590592582733E+09/ data i1(147)/ 0/,i2(147)/ 7/,i3(147)/ 2/ data cf(148)/0.1495886413107959E+09/ data i1(148)/ 0/,i2(148)/ 8/,i3(148)/ 1/ data cf(149)/0.2613284402260471E+05/ data i1(149)/ 1/,i2(149)/ 0/,i3(149)/ 8/ data cf(150)/0.1443371219504169E+07/ data i1(150)/ 1/,i2(150)/ 1/,i3(150)/ 7/ data cf(151)/0.6576436556610378E+06/ data i1(151)/ 1/,i2(151)/ 2/,i3(151)/ 6/ data cf(152)/-.4133389418926948E+08/ data i1(152)/ 1/,i2(152)/ 3/,i3(152)/ 5/ data cf(153)/0.5537360178547990E+08/ data i1(153)/ 1/,i2(153)/ 4/,i3(153)/ 4/ data cf(154)/-.1049826333710249E+08/ data i1(154)/ 1/,i2(154)/ 5/,i3(154)/ 3/ data cf(155)/0.4580450896890437E+08/ data i1(155)/ 1/,i2(155)/ 6/,i3(155)/ 2/ data cf(156)/-.2146916937769951E+09/ data i1(156)/ 1/,i2(156)/ 7/,i3(156)/ 1/ data cf(157)/-.3122490540937454E+08/ data i1(157)/ 1/,i2(157)/ 8/,i3(157)/ 0/ data cf(158)/0.9369338956631406E+05/ data i1(158)/ 2/,i2(158)/ 0/,i3(158)/ 7/ data cf(159)/-.4367918387697926E+07/ data i1(159)/ 2/,i2(159)/ 1/,i3(159)/ 6/ data cf(160)/0.1574085481348679E+08/ data i1(160)/ 2/,i2(160)/ 2/,i3(160)/ 5/ data cf(161)/0.9775479807388103E+07/ data i1(161)/ 2/,i2(161)/ 3/,i3(161)/ 4/ data cf(162)/-.9548837026566038E+07/ data i1(162)/ 2/,i2(162)/ 4/,i3(162)/ 3/ data cf(163)/-.7414973089524232E+08/ data i1(163)/ 2/,i2(163)/ 5/,i3(163)/ 2/ data cf(164)/0.1793364640748591E+09/ data i1(164)/ 2/,i2(164)/ 6/,i3(164)/ 1/ data cf(165)/0.6442134323693813E+07/ data i1(165)/ 2/,i2(165)/ 7/,i3(165)/ 0/ data cf(166)/0.5769621005587326E+06/ data i1(166)/ 3/,i2(166)/ 0/,i3(166)/ 6/ data cf(167)/0.1088546719787212E+06/ data i1(167)/ 3/,i2(167)/ 1/,i3(167)/ 5/ data cf(168)/-.6614520948461380E+07/ data i1(168)/ 3/,i2(168)/ 2/,i3(168)/ 4/ data cf(169)/-.1541201787947690E+08/ data i1(169)/ 3/,i2(169)/ 3/,i3(169)/ 3/ data cf(170)/-.1774431140749521E+07/ data i1(170)/ 3/,i2(170)/ 4/,i3(170)/ 2/ data cf(171)/-.6818325615773022E+07/ data i1(171)/ 3/,i2(171)/ 5/,i3(171)/ 1/ data cf(172)/0.1533861808341136E+08/ data i1(172)/ 3/,i2(172)/ 6/,i3(172)/ 0/ data cf(173)/-.4207752521631964E+06/ data i1(173)/ 4/,i2(173)/ 0/,i3(173)/ 5/ data cf(174)/0.1831982531281996E+06/ data i1(174)/ 4/,i2(174)/ 1/,i3(174)/ 4/ data cf(175)/0.8081588808120503E+07/ data i1(175)/ 4/,i2(175)/ 2/,i3(175)/ 3/ data cf(176)/0.1489286882606027E+08/ data i1(176)/ 4/,i2(176)/ 3/,i3(176)/ 2/ data cf(177)/0.4233140403408662E+08/ data i1(177)/ 4/,i2(177)/ 4/,i3(177)/ 1/ data cf(178)/-.6694307558274058E+08/ data i1(178)/ 4/,i2(178)/ 5/,i3(178)/ 0/ data cf(179)/0.3587582749128680E+06/ data i1(179)/ 5/,i2(179)/ 0/,i3(179)/ 4/ data cf(180)/-.1430330263889948E+07/ data i1(180)/ 5/,i2(180)/ 1/,i3(180)/ 3/ data cf(181)/-.1813752994229190E+07/ data i1(181)/ 5/,i2(181)/ 2/,i3(181)/ 2/ data cf(182)/-.3532432170832552E+08/ data i1(182)/ 5/,i2(182)/ 3/,i3(182)/ 1/ data cf(183)/0.1931755255191875E+08/ data i1(183)/ 5/,i2(183)/ 4/,i3(183)/ 0/ data cf(184)/-.6250653617654167E+05/ data i1(184)/ 6/,i2(184)/ 0/,i3(184)/ 3/ data cf(185)/0.1579711027633776E+06/ data i1(185)/ 6/,i2(185)/ 1/,i3(185)/ 2/ data cf(186)/0.7906592362440074E+07/ data i1(186)/ 6/,i2(186)/ 2/,i3(186)/ 1/ data cf(187)/0.8085573915353918E+07/ data i1(187)/ 6/,i2(187)/ 3/,i3(187)/ 0/ data cf(188)/0.8944396015574322E+05/ data i1(188)/ 7/,i2(188)/ 0/,i3(188)/ 2/ data cf(189)/-.7293579279942845E+06/ data i1(189)/ 7/,i2(189)/ 1/,i3(189)/ 1/ data cf(190)/-.1181132367272004E+06/ data i1(190)/ 7/,i2(190)/ 2/,i3(190)/ 0/ data cf(191)/0.5461367919265255E+05/ data i1(191)/ 8/,i2(191)/ 0/,i3(191)/ 1/ data cf(192)/-.5360217261156529E+05/ data i1(192)/ 8/,i2(192)/ 1/,i3(192)/ 0/ vex1=0.9118542142893855E+00 vex2=0.1350010350000000E+01 vex3=0.1003241840000000E+01 f12(0)=1.d0 f13(0)=1.d0 f23(0)=1.d0 bux12=r12*dexp(-vex1*r12) bux13=r13*dexp(-vex2*r13) bux23=r23*dexp(-vex3*r23) do 1 i=1, 8 f12(i)=f12(i-1)*bux12 f13(i)=f13(i-1)*bux13 f23(i)=f23(i-1)*bux23 1 continue ener = 0.d0 der12 = 0.d0 der13 = 0.d0 der23 = 0.d0 do 2 l=1,192 aux=f12(i1(l))*f13(i2(l))*f23(i3(l)) dux12=i1(l)*f12(i1(l)-1)*f13(i2(l))*f23(i3(l)) dux13=i2(l)*f12(i1(l))*f13(i2(l)-1)*f23(i3(l)) dux23=i3(l)*f12(i1(l))*f13(i2(l))*f23(i3(l)-1) ener=ener+cf(l)*aux der12=der12+cf(l)*dux12 der13=der13+cf(l)*dux13 der23=der23+cf(l)*dux23 2 continue der(1)=der12*(1.d0-vex1*r12)*dexp(-vex1*r12) der(2)=der13*(1.d0-vex2*r13)*dexp(-vex2*r13) der(3)=der23*(1.d0-vex3*r23)*dexp(-vex3*r23) return end