<div dir="ltr"><br><div class="gmail_extra">Dear Ramon,<br><br></div><div class="gmail_extra">I thought, the problem arising in IRC was because of negative Hessian eigenvalues but, now I realized that's not the case. Earlier I have done transition state optimization with Hessian calculated by forward differences and IRC with central differences so that, it was giving 4 negative Hessian values in the case of IRC.<br>
<br></div><div class="gmail_extra">Now I ran IRC with calculating Hessian using forward differences and  all positive Hessian eigenvalues except one got too. Even now also surprisingly giving the same error which I got earlier. Here, I also have posted the echo of output having error.<br>
<br>***********Partial out put************<br><br> Geometry optimization using default procedure for command MULTI<br><br> Geometry written to block  1 of record 700<br><br> Generated   9 symmetrical displacements<br><br>
 Numerically approximating hessian using forward gradient differences<br><br> Task list generated. Total number of displacements:      9<br><br>    5 tasks completed, CPU=1h 3m 15s Elapsed=1h 10m 13s<br><br> Numerical MULTI000 hessian completed. CPU-time:  6877.21 sec, Elapsed:  7645.49 sec<br>
<br> MULTI hessian saved to record  5300.2<br><br> Quadratic Steepest Descent -  Reaction Path Following using exact Hessian<br><br> Hessian eigenvalues:    -0.012734  0.000786  0.011948  0.018747  0.075415  0.141616  0.303108  1.646533  1.810536<br>
<br> IDIR= 1 requested but starting point is not a stationary point. Gradient norm=  0.98D-04  Step length=  0.36D-01<br> ? Error<br> ? Not a stationary point<br> ? The problem occurs in QSDPATH1<br>***********<br><br></div>
<div class="gmail_extra">What's wrong I am doing here....??<br><br></div><div class="gmail_extra">Thanks,<br><br></div><div class="gmail_extra">Mahi. <br></div><div class="gmail_extra"><br><br></div><div class="gmail_extra">
<br><div class="gmail_quote">On Wed, Aug 21, 2013 at 9:18 PM, Ramon Hernandez <span dir="ltr"><<a href="mailto:ramon@uaem.mx" target="_blank">ramon@uaem.mx</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Dear Mahi,<br>
<br>
I suggest you double-check that indeed your initial CASSCF calculation<br>
corresponds<br>
to a transition state, perhaps run a frequency calculation with tight<br>
parameters. Then<br>
when you run the IRC make sure you indeed converge to the same energy<br>
as the original<br>
CASSCF calculation. If you are sure about both these points then the<br>
problem might really<br>
have to do with a bug in IRC.<br>
<br>
Also if you run again the initial frequency calculation you can save<br>
that information and<br>
then read it directly in the next IRC job so you avoid calculating it twice.<br>
<br>
Cheers<br>
<br>
Ramon<br>
<br>
2013/8/21, G Mahesh <<a href="mailto:gmahesh@students.iiserpune.ac.in">gmahesh@students.iiserpune.ac.in</a>>:<br>
<div class=""><div class="h5">> Dear all,<br>
><br>
> I am familiar with MOLPRO from last two months. Since then I am trying to<br>
> find a transition state for the reaction of reactants and products having<br>
> multi reference character so, I am using CASSCF method for my problem.<br>
> After a long struggle, I have found the transition state with one imaginary<br>
> frequency (150 cm-1) and this normal mode having imaginary frequency seems<br>
> like reaction mode which will connect reactants to products.<br>
><br>
> Now the problem is with IRC, here, I started IRC with the geometry which I<br>
> got from the transition state optimization with QSDPATH method and IDIR=1<br>
> then, after Hessian calculation, it's complaining that, the given structure<br>
> is not stationary with 4 negative Hessian eigenvalues. Does this means, my<br>
> initial structure is not real transition state structure.....??  Should I<br>
> have to try for the transition state structure again ..?? Any ideas to<br>
> solve the problem....?? Pls help me...!<br>
><br>
> Below additional details may be required to understand the problem.<br>
><br>
> We are using the MOLPRO version 2012.<br>
> Our system has 5 atoms. We are using 13 active orbitals and 15 active<br>
> electrons as my cas space and AVTZ as basis.<br>
><br>
> Thanks,<br>
> Mahi<br>
><br>
</div></div></blockquote></div><br></div></div>