Dear Experts,<div><br></div><div>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.</div><div><br></div><div>I do not know why CP2K can not reproduce the same result as MOPAC.</div><div><br></div><div>Can anybody help me to figure out, what wrong is going on in CP2K calculation.</div><div><br></div><div>My simple input is following.</div><div><br></div><div>Thank you in advance.</div><div><br></div><div>Sincerely,</div><div>Bharat Sharma</div><div><br></div><div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  &DFT</div><div>   CHARGE 0</div><div>   LSD T</div><div>    &QS</div><div>      METHOD PM6</div><div>      &SE</div><div>       ANALYTICAL_GRADIENTS F</div><div>      &END</div><div>    &END QS</div><div>    &SCF</div><div>      SCF_GUESS ATOMIC</div><div>    &END SCF</div><div>  &END DFT</div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC 20.0 20.0 20.0</div><div>      PERIODIC NONE</div><div>    &END CELL</div><div>    &COORD</div><div> O                 -2.48815585    0.31201747    0.00696372</div><div> Ag                -4.60495278    0.41447825    0.06239400</div><div>    &END COORD</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div>&MOTION</div><div>  &GEO_OPT</div><div>   MAX_ITER 5000</div><div>  &END GEO_OPT</div><div>&END MOTION</div><div>&GLOBAL</div><div>  PROJECT AgO</div><div>  RUN_TYPE GEO_OPT</div><div>&END GLOBAL</div></div>