<div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><div><b>Dear Molpro users,</b><br><b>    I got a problem when i was doing qsdpath calculation.</b><br><b>    I have already done the saddle search and calculated the frequences ,and the piont has only one imaginary vibration frequence which means </b><br><b>it is a transation state.</b><br><b>   But when i run the qsd module,I got this error:</b><br>  <br>  Task list generated. Total number of displacements:    156<br>   15 tasks completed, CPU=1h 10m 4s Elapsed=1h 17m 34s<br>   30 tasks completed, CPU=2h 20m 8s Elapsed=2h 34m 52s<br>   45 tasks completed, CPU=3h 30m 11s Elapsed=3h 51m 58s<br>   60 tasks completed, CPU=4h 40m 14s Elapsed=5h 9m 25s<br>   75 tasks completed, CPU=5h 50m 22s Elapsed=6h 26m 39s<br>   90 tasks completed, CPU=7h 0m 25s Elapsed=7h 44m 25s<br>  105 tasks completed, CPU=8h 10m 30s Elapsed=9h 1m 53s<br>  120 tasks completed, CPU=9h 20m 36s Elapsed=10h 19m 30s<br>  135 tasks completed, CPU=10h 30m 43s Elapsed=11h 36m 58s<br>  150 tasks completed, CPU=11h 40m 44s Elapsed=12h 54m 29s<br> Tasks on processor    0 finished.  CPU=  43722.84 sec, Elapsed=  48318.93 sec. All tasks on all processors finished.<br><br> CCSD(T)-F12 hessian saved to record  5301.2<br><br> Combined Powell-Murtagh-Sargent Update of Hessian<br><br> Quadratic Steepest Descent -  Reaction Path Following using exact Hessian<br><br> Hessian eigenvalues:    -0.033385  0.000257  0.005602  0.015645  0.026065  0.068585  0.128212  0.208396  0.265296  0.288874<br>                          0.305995  0.560727<br><br> IDIR= 1 requested but starting point is not a stationary point. Gradient norm=  0.40D-04  Step length=  0.18D-02<br> ? Error<br> ? Not a stationary point<br> ? The problem occurs in QSDPATH1<br><br><b>So why isn't it a stationary point??</b><br><b>Here is my input:</b><br><br> ***,  Reaction Path<br> memory,250,m<br> charge=-1<br> basis=vtz<br> r1=1.95891 ang<br> r2=1.12844 ang<br> r3=1.13476 ang<br> r4=2.0873  ang<br> r5=0.92512 ang<br> a1=97.80701 degree<br> a2=98.16029 degree<br> a3=138.52460 degree<br> a4=173.67715 degree<br> a5=-167.40668 degree<br> a6=102.85072 degree<br> a7=41.71448 degree<br> geometry={<br>           Cl<br>           C,Cl,r1<br>           H1,C,r2,Cl,a1<br>           H2,C,r3,Cl,a2,H1,a6<br>           H3,C,r4,Cl,a3,H1,a7<br>           F,H3,r5,C,a4,Cl,a5}<br> int<br> rhf<br> ccsd(t)-f12<br> optg,root=2,saveact=geo_ts,rewind            ! Find and store the TS<br> {optg,method=qsdpath,dir=1, numhess=5,hesscentral,saveact=geo_path}         ! find IRC in positive direction<br> readvar,geo_ts.act                           ! Reset geometry to TS<br> {optg,method=qsdpath,dir=-1,numhess=5,hesscentral,saveact=geo_path,append} !find IRC in negative direction<br><br> readvar,geo_path.act<br><br> alpha=alpha*pi/180    !convert angle to radian<br><br> table,irc,r1,r2,r3,r4,r5,a1,a2,a3,a4,a5,e_opt   !tabulate results<br><br> {table,irc,e_opt                !plot energy profile as function of irc<br>  plot,file='geo_eopt.plot'}<br><br> {table,irc,r1,r2,r3,r4,r5,a1,a2,a3,a4,a5        !plot distances and angle as function of irc<br>  plot,file='geo_dist.plot'}<br><br>Please give me a hint.<br>Thank you.<br><br><br><br></div></div><br><br><span title="neteasefooter"><p> </p></span></div><br><br><span title="neteasefooter"><p> </p></span>