<div dir="ltr"><div><div><div>Thanks for your reply.<br></div>   I have checked the log file and the forces information is not written there. <br><br><br></div>Thanks once again.<br><br></div>Nag<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Sep 18, 2015 at 4:30 AM, george <span dir="ltr"><<a href="mailto:george.iitm@gmail.com" target="_blank">george.iitm@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>hi <br><br></div>Do the calculation in single input file for both opt+freq. After the opt is over in freq state read the orbital from optg calculation. In this way its possible to get same energy as optimal energy in sp calculation. Now reg. forces i guess it is written in log file in the form xyz coords.<br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Thu, Sep 17, 2015 at 5:31 PM, samala nagaprasad reddy <span dir="ltr"><<a href="mailto:snpreddy63@gmail.com" target="_blank">snpreddy63@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div><div><div><div>Hi<br></div><div>  I have a couple of questions.<br><br></div>1)
 I am calculating optimization and frequency calculation of benzene. 
After optimization when I calculate the single point energy calculation 
the SP energy is lower than the optimized energy. I am taking the 
optimized Z matrix to run the SP energy calculation. How it is possible?
 What I am missing? Here I am giving my input of optimization and output
 files of both SP and optimization calculations.<br><br></div><div>Optimization input<br></div><div><br> ***,coupling constants<br>  memory,450,m<br> cc2    =    1.394207,<br> ccc3   =    120.000,<br> dih4   =      0.000,<br> hc7    =    1.082115,<br> dih7   =    180.000<br>  basis,avtz<br>  geomtype = zmat<br>  geometry = {angstrom,<br> c;<br> c,   1, cc2;<br> c,   2, cc2,        1, ccc3;<br> c,   3, cc2,        2, ccc3,          1, dih4;<br> c,   4, cc2,        3, ccc3,          2, dih4;<br> c,   5, cc2,        4, ccc3,          3, dih4;<br> h,   2, hc7,        3, ccc3,          4, dih7;<br> h,   3, hc7,        2, ccc3,          1, dih7;<br> h,   4, hc7,        3, ccc3,          2, dih7;<br> h,   5, hc7,        4, ccc3,         3, dih7;<br> h,   6, hc7,        5, ccc3,         4, dih7;<br> h,   1, hc7,        2, ccc3,         3, dih7;<br> }<br> {hf,ENERGY=1.0E-10,orbital=1.0E-8;}<br> {df-hf,ENERGY=1.0E-10,orbital=1.0E-8;}<br> {df-mp2,ENERGY=1.0E-10,orbital=1.0E-8;}<br> optg;<br> frequencies;<br></div><div><br></div><div>SP energy input<br><br>! Cite this work as:<br> ***,coupling constants<br> memory,450,m<br>cc2    =    1.394247,<br>cc3    =    1.394247,<br>ccc3   =    120.000,<br>cc4    =    1.394247,<br>ccc4   =    120.000,<br>dih4   =      0.000,<br>cc5    =    1.394247,<br>ccc5   =    120.000,<br>dih5   =      0.000,<br>cc6    =    1.394247,<br>ccc6   =    120.000,<br>dih6   =      0.000,<br>hc7    =    1.082198,<br>hcc7   =    120.000,<br>dih7   =    180.000,<br>hc8    =    1.082198,<br>hcc8   =    120.000,<br>dih8   =    180.000,<br>hc9    =    1.082198,<br>hcc9   =    120.000,<br>dih9   =    180.000,<br>hc10   =    1.082198,<br>hcc10  =    120.000,<br>dih10  =    180.000,<br>hc11   =    1.082198,<br>hcc11  =    120.000,<br>dih11  =    180.000,<br>hc12   =    1.082198,<br>hcc12  =    120.000,<br>dih12  =    180.000<br>basis,avtz<br>!gprint,orbitals,civectors<br> geomtype = zmat<br> geometry = {angstrom,<br>c;<br>c,   1, cc2;<br>c,   2, cc3,        1, ccc3;<br>c,   3, cc4,        2, ccc4,          1, dih4;<br>c,   4, cc5,        3, ccc5,          2, dih5;<br>c,   5, cc6,        4, ccc6,          3, dih6;<br>h,   3, hc7,        2, hcc7,          1, dih7;<br>h,   4, hc8,        3, hcc8,          2, dih8;<br>h,   5, hc9,        4, hcc9,          3, dih9;<br>h,   6, hc10,       5, hcc10,         4, dih10;<br>h,   1, hc11,       2, hcc11,         3, dih11;<br>h,   2, hc12,       3, hcc12,         4, dih12;<br>}<br>{hf,ENERGY=1.0E-10,orbital=1.0E-8;}<br>{df-hf,ENERGY=1.0E-10,orbital=1.0E-8;}<br>{df-mp2,ENERGY=1.0E-10,orbital=1.0E-8;}<br><br><br></div>2)
 In the frequency (numerical) calculation, it distort the molecule to 
get Hessian. I want the forces at each displacement. 
When I check the output files, it is not wring the information. How can I
 get the information?<br><br></div>I look forward to hear from you<br><br></div>Thank you</div>
<br></div></div><span class="">_______________________________________________<br>
Molpro-user mailing list<br>
<a href="mailto:Molpro-user@molpro.net" target="_blank">Molpro-user@molpro.net</a><br>
<a href="http://www.molpro.net/mailman/listinfo/molpro-user" rel="noreferrer" target="_blank">http://www.molpro.net/mailman/listinfo/molpro-user</a><br></span></blockquote></div><br></div>
</blockquote></div><br></div>