Dear Hutter,<div><br></div><div>Thank you for your reply.</div><div><br></div><div>Au parameters are correct as reported in paper. I started the calculation with the Ag-O distance of 2.03Angstrom both in CP2K and MOPAC (simply same initial geometry). The uptimized bond length of Ag-O from CP2K is 1.90 A where as from MOPAC is 2.03A which is consistent with PM6 results in "AM1* parameters for palladium and silver, J Mol Model (2011) 17:2585–2600)".</div>
<div><br></div><div>I did not one shot energy calculation with Ag-O distance of 2.03A in CP2K but the heat of formation is 135.7645 kcal/mole but from MOPAC is 98.21kcal/mole.</div><div><br></div><div>I did the test calculation of AgH+, Ag(NH3)4+, AgF from both in CP2K and MOPAC but results are fairly similar( error ~ or less 1 kcal/mole).</div>
<div><br></div><div>In all test calculations, the initial geometries for CP2K and MOPAC are same.</div><div><br></div><div>Thank you.</div><div><br></div><div>Sincerely,</div><div>Bharat Sharma</div><div><br><div class="gmail_quote">
On Fri, Aug 17, 2012 at 7:36 AM,  <span dir="ltr"><<a href="mailto:hut...@pci.uzh.ch" target="_blank">hut...@pci.uzh.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi<br>
<br>
I quickly checked, the Au parameters should be correct.<br>
Can you make some additional tests? For example is the optimized<br>
geometry the same? What is the difference in energy for a given<br>
geometry?<br>
Are the results correct for other molecules containing these<br>
atomic types?<br>
<br>
regards<br>
<br>
Juerg Hutter<br>
<br>
--------------------------------------------------------------<br>
Juerg Hutter                         Phone : <a href="tel:%2B%2B41%2044%20635%204491" value="+41446354491">++41 44 635 4491</a><br>
Physical Chemistry Institute   FAX   : <a href="tel:%2B%2B41%2044%20635%206838" value="+41446356838">++41 44 635 6838</a><br>
University of Zurich               E-mail:  <a href="mailto:hut...@pci.uzh.ch">hut...@pci.uzh.ch</a><br>
Winterthurerstrasse 190<br>
CH-8057 Zurich, Switzerland<br>
---------------------------------------------------------------<br>
<br>
-----<a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a> wrote: -----<br>
To: <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a><br>
From: bharat<br>
Sent by: <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a><br>
Date: 08/16/2012 04:21PM<br>
Subject: [CP2K:3974] CP2K PM6 result does not consistent with MOPAC PM6<br>
<div class="HOEnZb"><div class="h5"><br>
Dear Experts,<br>
<br>
I did the semiempirical calculation with PM6 for the small system of AgO. CP2K gives the Heat of formation of 134.9850424864 kcal/mol whereas MOPAC gives the 98.21606 kcal/mol. The experimental Heat of formation is about 53.0 kcal/mol.<br>

<br>
I do not know why CP2K can not reproduce the same result as MOPAC.<br>
<br>
Can anybody help me to figure out, what wrong is going on in CP2K calculation.<br>
<br>
My simple input is following.<br>
<br>
Thank you in advance.<br>
<br>
Sincerely,<br>
Bharat Sharma<br>
<br>
&FORCE_EVAL<br>
  METHOD Quickstep<br>
  &DFT<br>
   CHARGE 0<br>
   LSD T<br>
    &QS<br>
      METHOD PM6<br>
      &SE<br>
       ANALYTICAL_GRADIENTS F<br>
      &END<br>
    &END QS<br>
    &SCF<br>
      SCF_GUESS ATOMIC<br>
    &END SCF<br>
  &END DFT<br>
  &SUBSYS<br>
    &CELL<br>
      ABC 20.0 20.0 20.0<br>
      PERIODIC NONE<br>
    &END CELL<br>
    &COORD<br>
 O                 -2.48815585    0.31201747    0.00696372<br>
 Ag                -4.60495278    0.41447825    0.06239400<br>
    &END COORD<br>
  &END SUBSYS<br>
&END FORCE_EVAL<br>
&MOTION<br>
  &GEO_OPT<br>
   MAX_ITER 5000<br>
  &END GEO_OPT<br>
&END MOTION<br>
&GLOBAL<br>
  PROJECT AgO<br>
  RUN_TYPE GEO_OPT<br>
&END GLOBAL<br>
  --<br>
 You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
 To view this discussion on the web visit <a href="https://groups.google.com/d/msg/cp2k/-/3fh3gtrWSUoJ" target="_blank">https://groups.google.com/d/msg/cp2k/-/3fh3gtrWSUoJ</a>.<br>
  To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
 To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
  For more options, visit this group at <a href="http://groups.google.com/group/cp2k?hl=en" target="_blank">http://groups.google.com/group/cp2k?hl=en</a>.<br>
<br>
<br>
--<br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
For more options, visit this group at <a href="http://groups.google.com/group/cp2k?hl=en" target="_blank">http://groups.google.com/group/cp2k?hl=en</a>.<br>
<br>
</div></div></blockquote></div><br></div>