cmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmc
      subroutine LiHF_AP_pes(x,v)
      implicit real*8 (a-h,o-z)
      dimension r(3),x(3)
      r(1)=x(3)
      r(2)=x(2)
      r(3)=x(1)
      call LIFHPAN_LiHF_AP(r,v)
      return
      end

cmmmmmm  Paniagua's LiFH PES fitted by Paniagua
************************************************************************
*     From this point all in atomic units
************************************************************************
      SUBROUTINE LIFHPAN_LiHF_AP(R,ENER)
      IMPLICIT REAL * 8 (A-H,O-Z)
      DIMENSION R(3),DER(3)
      Y=R(1)
      Z=R(2)
      X=R(3)
      CALL LIFH7_LiHF_AP(X,Y,Z,ENER,DER)
      CALL LIF_LiHF_AP(X,E1,D1)
      CALL FH_LiHF_AP(Y,E2,D2)
      CALL LIH_LiHF_AP(Z,E3,D3)
      ENER = ENER + E1 + E2 + E3
C     equilibrium distance  HF = 1.7413
C     energy in a.u. -0.201791344
      ENER=ENER+0.201791344d0
***** DERIVATIVES CALCULATION *****
      DER(1) = DER(1)+D1
      DER(2) = DER(2)+D2
      DER(3) = DER(3)+D3
      RETURN
      END

************************************************************************
*                                                                      *
*                                                                      *
*  ESTA SUBRUTINA CALCULA LAS ENERGIAS DE UNA CURVA DIATOMICA          *
*  CON UN RMS DE        .0715 Kcal/mol                                 *
*  Y UN EMAX  DE        .1764 Kcal/mol                                 *
*                                                                      *
************************************************************************
      SUBROUTINE LIF_LiHF_AP(X,ENER,DER)
      IMPLICIT REAL*8 (A-H,O-Z)
      PARAMETER (CF0= .259060355634939D+02)
      PARAMETER (CF1=-.571402848087689D+00)
      PARAMETER (CF2= .708433808104887D+00)
      PARAMETER (CF3=-.918329923793955D+00)
      PARAMETER (ALFA=.183971702322879D+01)
      PARAMETER (BETA=.585612557426368D+00)
      AUX = DEXP(- ALFA*X)
      BUX = 1.D0/X
      CUX = DEXP(- BETA*X)
***** CALCULO DE LA ENERGIA *****
      ENER = 0.D0
      ENER=ENER+CF0*AUX*BUX
      RHO =X*CUX
        ENER=ENER+CF1*RHO
        ENER=ENER+CF2*RHO**2
        ENER=ENER+CF3*RHO**3
***** CALCULO DE LA DERIVADA *****
        DER =CF1 +2.D0*CF2*RHO
        DER =DER +3.D0*CF3*RHO**2
      DER = DER*(1.D0-BETA*X)*CUX
      DER = DER - CF0*BUX*(ALFA+BUX)*AUX
      RETURN
      END


************************************************************************
*                                                                      *
*                                                                      *
*  ESTA SUBRUTINA CALCULA LAS ENERGIAS DE UNA CURVA DIATOMICA          *
*  CON UN RMS DE        .2535 Kcal/mol                                 *
*  Y UN EMAX  DE        .7348 Kcal/mol                                 *
*                                                                      *
************************************************************************
      SUBROUTINE FH_LiHF_AP(X,ENER,DER)
      IMPLICIT REAL*8 (A-H,O-Z)
      PARAMETER (CF0= .146251212596914D+02)
      PARAMETER (CF1=-.875424365172117D-02)
      PARAMETER (CF2=-.508168931989300D+01)
      PARAMETER (CF3= .750286609001793D+01)
      PARAMETER (ALFA=.308363644148830D+01)
      PARAMETER (BETA=.104065052653433D+01)
      AUX = DEXP(- ALFA*X)
      BUX = 1.D0/X
      CUX = DEXP(- BETA*X)
***** CALCULO DE LA ENERGIA *****
      ENER = 0.D0
      ENER=ENER+CF0*AUX*BUX
      RHO =X*CUX
        ENER=ENER+CF1*RHO
        ENER=ENER+CF2*RHO**2
        ENER=ENER+CF3*RHO**3
***** CALCULO DE LA DERIVADA *****
        DER =CF1 +2.D0*CF2*RHO
        DER =DER +3.D0*CF3*RHO**2
      DER = DER*(1.D0-BETA*X)*CUX
      DER = DER - CF0*BUX*(ALFA+BUX)*AUX
      RETURN
      END


************************************************************************
*                                                                      *
*                                                                      *
*  ESTA SUBRUTINA CALCULA LAS ENERGIAS DE UNA CURVA DIATOMICA          *
*  CON UN RMS DE        .3207 Kcal/mol                                 *
*  Y UN EMAX  DE        .8803 Kcal/mol                                 *
*                                                                      *
************************************************************************
      SUBROUTINE LIH_LiHF_AP(X,ENER,DER)
      IMPLICIT REAL*8 (A-H,O-Z)
      PARAMETER (CF0= .326698672830820D+02)
      PARAMETER (CF1= .224436053091229D-01)
      PARAMETER (CF2=-.140464674839071D+01)
      PARAMETER (CF3= .215282944314996D+01)
      PARAMETER (ALFA=.286994963677707D+01)
      PARAMETER (BETA=.658031223064788D+00)
      AUX = DEXP(- ALFA*X)
      BUX = 1.D0/X
      CUX = DEXP(- BETA*X)
***** CALCULO DE LA ENERGIA *****
      ENER = 0.D0
      ENER=ENER+CF0*AUX*BUX
      RHO =X*CUX
        ENER=ENER+CF1*RHO
        ENER=ENER+CF2*RHO**2
        ENER=ENER+CF3*RHO**3
***** CALCULO DE LA DERIVADA *****
        DER =CF1 +2.D0*CF2*RHO
        DER =DER +3.D0*CF3*RHO**2
      DER = DER*(1.D0-BETA*X)*CUX
      DER = DER - CF0*BUX*(ALFA+BUX)*AUX
      RETURN
      END


*****************************************************************
*   SUBROUTINE LIFH7:   LiFH 3D potential
*    RMS =      0.65       Kcal/mol
*
*****************************************************************

      SUBROUTINE LIFH7_LiHF_AP(X,Y,Z,ENER,DER)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION IZ( 98), JZ( 98), KZ( 98), CF( 98)
      DIMENSION XXX(-1:6), YYY(-1:6), ZZZ(-1:6)
      DIMENSION DER(3)
      DATA CF(  1)/0.478684594760D+01/,IZ(  1)/1/,JZ(  1)/1/,KZ(  1)/0/
      DATA CF(  2)/0.729385703517D+00/,IZ(  2)/1/,JZ(  2)/0/,KZ(  2)/1/
      DATA CF(  3)/-.156396911454D+02/,IZ(  3)/0/,JZ(  3)/1/,KZ(  3)/1/
      DATA CF(  4)/0.132639116265D+02/,IZ(  4)/1/,JZ(  4)/1/,KZ(  4)/1/
      DATA CF(  5)/-.269964012146D+02/,IZ(  5)/2/,JZ(  5)/1/,KZ(  5)/0/
      DATA CF(  6)/-.210930935654D+01/,IZ(  6)/2/,JZ(  6)/0/,KZ(  6)/1/
      DATA CF(  7)/0.106587035834D+03/,IZ(  7)/0/,JZ(  7)/2/,KZ(  7)/1/
      DATA CF(  8)/0.174814146447D+03/,IZ(  8)/0/,JZ(  8)/1/,KZ(  8)/2/
      DATA CF(  9)/0.237158025327D+02/,IZ(  9)/1/,JZ(  9)/2/,KZ(  9)/0/
      DATA CF( 10)/-.423941830065D+00/,IZ( 10)/1/,JZ( 10)/0/,KZ( 10)/2/
      DATA CF( 11)/-.233169881047D+02/,IZ( 11)/2/,JZ( 11)/1/,KZ( 11)/1/
      DATA CF( 12)/-.600217207207D+03/,IZ( 12)/1/,JZ( 12)/2/,KZ( 12)/1/
      DATA CF( 13)/-.755939394816D+03/,IZ( 13)/1/,JZ( 13)/1/,KZ( 13)/2/
      DATA CF( 14)/0.306634456126D+03/,IZ( 14)/2/,JZ( 14)/2/,KZ( 14)/0/
      DATA CF( 15)/-.734950008434D+01/,IZ( 15)/2/,JZ( 15)/0/,KZ( 15)/2/
      DATA CF( 16)/-.363912831022D+03/,IZ( 16)/0/,JZ( 16)/2/,KZ( 16)/2/
      DATA CF( 17)/0.583491107309D+02/,IZ( 17)/3/,JZ( 17)/1/,KZ( 17)/0/
      DATA CF( 18)/0.317300063172D+02/,IZ( 18)/3/,JZ( 18)/0/,KZ( 18)/1/
      DATA CF( 19)/0.134685929111D+04/,IZ( 19)/0/,JZ( 19)/3/,KZ( 19)/1/
      DATA CF( 20)/-.527904207203D+01/,IZ( 20)/0/,JZ( 20)/1/,KZ( 20)/3/
      DATA CF( 21)/-.950064232384D+03/,IZ( 21)/1/,JZ( 21)/3/,KZ( 21)/0/
      DATA CF( 22)/0.139715468262D+02/,IZ( 22)/1/,JZ( 22)/0/,KZ( 22)/3/
      DATA CF( 23)/0.379581931507D+03/,IZ( 23)/2/,JZ( 23)/2/,KZ( 23)/1/
      DATA CF( 24)/0.132766643990D+04/,IZ( 24)/2/,JZ( 24)/1/,KZ( 24)/2/
      DATA CF( 25)/0.948956875368D+02/,IZ( 25)/1/,JZ( 25)/2/,KZ( 25)/2/
      DATA CF( 26)/0.348535973557D+03/,IZ( 26)/3/,JZ( 26)/1/,KZ( 26)/1/
      DATA CF( 27)/0.211080767780D+04/,IZ( 27)/1/,JZ( 27)/3/,KZ( 27)/1/
      DATA CF( 28)/-.510103186668D+03/,IZ( 28)/1/,JZ( 28)/1/,KZ( 28)/3/
      DATA CF( 29)/-.631369820301D+03/,IZ( 29)/3/,JZ( 29)/2/,KZ( 29)/0/
      DATA CF( 30)/0.661457533096D+02/,IZ( 30)/3/,JZ( 30)/0/,KZ( 30)/2/
      DATA CF( 31)/0.281482564549D+03/,IZ( 31)/0/,JZ( 31)/3/,KZ( 31)/2/
      DATA CF( 32)/0.197806092197D+04/,IZ( 32)/0/,JZ( 32)/2/,KZ( 32)/3/
      DATA CF( 33)/-.105743974604D+03/,IZ( 33)/2/,JZ( 33)/3/,KZ( 33)/0/
      DATA CF( 34)/-.705711919175D+02/,IZ( 34)/2/,JZ( 34)/0/,KZ( 34)/3/
      DATA CF( 35)/-.724837111795D+02/,IZ( 35)/4/,JZ( 35)/1/,KZ( 35)/0/
      DATA CF( 36)/-.112253320655D+03/,IZ( 36)/4/,JZ( 36)/0/,KZ( 36)/1/
      DATA CF( 37)/-.118407656996D+05/,IZ( 37)/0/,JZ( 37)/4/,KZ( 37)/1/
      DATA CF( 38)/0.596460369684D+02/,IZ( 38)/0/,JZ( 38)/1/,KZ( 38)/4/
      DATA CF( 39)/0.497387661368D+04/,IZ( 39)/1/,JZ( 39)/4/,KZ( 39)/0/
      DATA CF( 40)/-.750727331646D+02/,IZ( 40)/1/,JZ( 40)/0/,KZ( 40)/4/
      DATA CF( 41)/0.391466559266D+04/,IZ( 41)/2/,JZ( 41)/2/,KZ( 41)/2/
      DATA CF( 42)/-.120979039019D+03/,IZ( 42)/3/,JZ( 42)/2/,KZ( 42)/1/
      DATA CF( 43)/-.112278614225D+04/,IZ( 43)/3/,JZ( 43)/1/,KZ( 43)/2/
      DATA CF( 44)/-.346070725529D+04/,IZ( 44)/1/,JZ( 44)/3/,KZ( 44)/2/
      DATA CF( 45)/-.763479011563D+04/,IZ( 45)/1/,JZ( 45)/2/,KZ( 45)/3/
      DATA CF( 46)/-.522738333805D+04/,IZ( 46)/2/,JZ( 46)/3/,KZ( 46)/1/
      DATA CF( 47)/0.594817464194D+03/,IZ( 47)/2/,JZ( 47)/1/,KZ( 47)/3/
      DATA CF( 48)/0.224748366931D+04/,IZ( 48)/3/,JZ( 48)/3/,KZ( 48)/0/
      DATA CF( 49)/0.213956398677D+02/,IZ( 49)/3/,JZ( 49)/0/,KZ( 49)/3/
      DATA CF( 50)/0.219941926066D+04/,IZ( 50)/0/,JZ( 50)/3/,KZ( 50)/3/
      DATA CF( 51)/-.658409969143D+03/,IZ( 51)/4/,JZ( 51)/1/,KZ( 51)/1/
      DATA CF( 52)/0.110136079082D+05/,IZ( 52)/1/,JZ( 52)/4/,KZ( 52)/1/
      DATA CF( 53)/0.128421948481D+04/,IZ( 53)/1/,JZ( 53)/1/,KZ( 53)/4/
      DATA CF( 54)/0.400598029012D+03/,IZ( 54)/4/,JZ( 54)/2/,KZ( 54)/0/
      DATA CF( 55)/-.985792793715D+02/,IZ( 55)/4/,JZ( 55)/0/,KZ( 55)/2/
      DATA CF( 56)/0.266636830207D+03/,IZ( 56)/0/,JZ( 56)/4/,KZ( 56)/2/
      DATA CF( 57)/0.980902655287D+02/,IZ( 57)/0/,JZ( 57)/2/,KZ( 57)/4/
      DATA CF( 58)/-.661441078012D+04/,IZ( 58)/2/,JZ( 58)/4/,KZ( 58)/0/
      DATA CF( 59)/0.224980053981D+03/,IZ( 59)/2/,JZ( 59)/0/,KZ( 59)/4/
      DATA CF( 60)/0.531197782165D+02/,IZ( 60)/5/,JZ( 60)/1/,KZ( 60)/0/
      DATA CF( 61)/0.144862253668D+03/,IZ( 61)/5/,JZ( 61)/0/,KZ( 61)/1/
      DATA CF( 62)/0.207519072761D+05/,IZ( 62)/0/,JZ( 62)/5/,KZ( 62)/1/
      DATA CF( 63)/-.342480274672D+03/,IZ( 63)/0/,JZ( 63)/1/,KZ( 63)/5/
      DATA CF( 64)/-.522800362269D+04/,IZ( 64)/1/,JZ( 64)/5/,KZ( 64)/0/
      DATA CF( 65)/0.822327905488D+02/,IZ( 65)/1/,JZ( 65)/0/,KZ( 65)/5/
      DATA CF( 66)/-.490934281115D+04/,IZ( 66)/3/,JZ( 66)/2/,KZ( 66)/2/
      DATA CF( 67)/0.962472020345D+04/,IZ( 67)/2/,JZ( 67)/3/,KZ( 67)/2/
      DATA CF( 68)/0.653133561903D+04/,IZ( 68)/2/,JZ( 68)/2/,KZ( 68)/3/
      DATA CF( 69)/-.193077093203D+04/,IZ( 69)/3/,JZ( 69)/3/,KZ( 69)/1/
      DATA CF( 70)/0.117964603759D+04/,IZ( 70)/3/,JZ( 70)/1/,KZ( 70)/3/
      DATA CF( 71)/-.686465507495D+04/,IZ( 71)/1/,JZ( 71)/3/,KZ( 71)/3/
      DATA CF( 72)/0.712775202220D+03/,IZ( 72)/4/,JZ( 72)/2/,KZ( 72)/1/
      DATA CF( 73)/0.197622992810D+03/,IZ( 73)/4/,JZ( 73)/1/,KZ( 73)/2/
      DATA CF( 74)/-.146287814680D+05/,IZ( 74)/1/,JZ( 74)/4/,KZ( 74)/2/
      DATA CF( 75)/0.982200900423D+03/,IZ( 75)/1/,JZ( 75)/2/,KZ( 75)/4/
      DATA CF( 76)/0.139974105843D+05/,IZ( 76)/2/,JZ( 76)/4/,KZ( 76)/1/
      DATA CF( 77)/-.502881797236D+04/,IZ( 77)/2/,JZ( 77)/1/,KZ( 77)/4/
      DATA CF( 78)/-.862358982023D+02/,IZ( 78)/4/,JZ( 78)/3/,KZ( 78)/0/
      DATA CF( 79)/-.629084324618D+02/,IZ( 79)/4/,JZ( 79)/0/,KZ( 79)/3/
      DATA CF( 80)/0.439889079902D+04/,IZ( 80)/0/,JZ( 80)/4/,KZ( 80)/3/
      DATA CF( 81)/0.178631055845D+04/,IZ( 81)/0/,JZ( 81)/3/,KZ( 81)/4/
      DATA CF( 82)/-.390915475835D+04/,IZ( 82)/3/,JZ( 82)/4/,KZ( 82)/0/
      DATA CF( 83)/0.113426507404D+03/,IZ( 83)/3/,JZ( 83)/0/,KZ( 83)/4/
      DATA CF( 84)/0.347194278814D+03/,IZ( 84)/5/,JZ( 84)/1/,KZ( 84)/1/
      DATA CF( 85)/-.397270220449D+05/,IZ( 85)/1/,JZ( 85)/5/,KZ( 85)/1/
      DATA CF( 86)/0.525784164789D+04/,IZ( 86)/1/,JZ( 86)/1/,KZ( 86)/5/
      DATA CF( 87)/-.193317625223D+03/,IZ( 87)/5/,JZ( 87)/2/,KZ( 87)/0/
      DATA CF( 88)/0.570614540262D+02/,IZ( 88)/5/,JZ( 88)/0/,KZ( 88)/2/
      DATA CF( 89)/0.120635405962D+05/,IZ( 89)/0/,JZ( 89)/5/,KZ( 89)/2/
      DATA CF( 90)/-.160938477427D+04/,IZ( 90)/0/,JZ( 90)/2/,KZ( 90)/5/
      DATA CF( 91)/0.168870872884D+05/,IZ( 91)/2/,JZ( 91)/5/,KZ( 91)/0/
      DATA CF( 92)/-.384974842833D+03/,IZ( 92)/2/,JZ( 92)/0/,KZ( 92)/5/
      DATA CF( 93)/-.132778419344D+02/,IZ( 93)/6/,JZ( 93)/1/,KZ( 93)/0/
      DATA CF( 94)/-.645978626900D+02/,IZ( 94)/6/,JZ( 94)/0/,KZ( 94)/1/
      DATA CF( 95)/0.608615169042D+04/,IZ( 95)/0/,JZ( 95)/6/,KZ( 95)/1/
      DATA CF( 96)/-.261023065095D+04/,IZ( 96)/0/,JZ( 96)/1/,KZ( 96)/6/
      DATA CF( 97)/-.840648972386D+04/,IZ( 97)/1/,JZ( 97)/6/,KZ( 97)/0/
      DATA CF( 98)/0.821495977419D+02/,IZ( 98)/1/,JZ( 98)/0/,KZ( 98)/6/
      DATA XXX(-1)/ 0.D0/,YYY(-1)/ 0.D0/,ZZZ(-1)/ 0.D0/
      DATA XXX(0)/ 1.D0/,YYY(0)/ 1.D0/,ZZZ(0)/ 1.D0/
      VEX1 =  0.450010967388186D+00
      VEX2 =  0.121996183914580D+01
      VEX3 =  0.826972051663734D+00
      X1 = DEXP(-VEX1*X)
      Y1 = DEXP(-VEX2*Y)
      Z1 = DEXP(-VEX3*Z)
      XX = X*X1
      YY = Y*Y1
      ZZ = Z*Z1

      DO 1 L=1,6
         XXX(L) = XX*XXX(L-1)
         YYY(L) = YY*YYY(L-1)
         ZZZ(L) = ZZ*ZZZ(L-1)
    1 CONTINUE
      ENER = 0.D0
      DER1 = 0.D0
      DER2 = 0.D0
      DER3 = 0.D0
      DO 2 L=1,98
         ENER=ENER+CF(L)*XXX(IZ(L))*YYY(JZ(L))*ZZZ(KZ(L))
         DER1=DER1+CF(L)*IZ(L)*XXX(IZ(L)-1)*YYY(JZ(L))*ZZZ(KZ(L))
         DER2=DER2+CF(L)*JZ(L)*XXX(IZ(L))*YYY(JZ(L)-1)*ZZZ(KZ(L))
         DER3=DER3+CF(L)*KZ(L)*XXX(IZ(L))*YYY(JZ(L))*ZZZ(KZ(L)-1)
    2 CONTINUE
      DER(1) = DER1 * (1.D0-VEX1*X)*X1
      DER(2) = DER2 * (1.D0-VEX2*Y)*Y1
      DER(3) = DER3 * (1.D0-VEX3*Z)*Z1
      RETURN
      END
