thanks axel and Teo,<br><br>I use the H-H, H-O and O-O parameters for PBE,  taking from the Chem. Phys. 327 (2006) 54. <br><br>
          ATOMS        H O<br>
          RCUT [angstrom]        10<br>
          RMIN [angstrom]         0.2<br>
          RMAX [angstrom]        11.0<br>
          VARIABLES    R<br>
          PARAMETERS   CSIX      RZERO    D<br>
          VALUES      -2.675     3.947    3<br>
          FUNCTION    (CSIX/(R**6))*((1.0 - EXP(-1.0*D*((R/RZERO)**7)))**4)<br><br>          ATOMS        H H<br>          PARAMETERS   CSIX      RZERO    D<br>          VALUES      -2.290     3.538    3<br><br>          ATOMS        O O<br>
          PARAMETERS   CSIX      RZERO    D<br>          VALUES       -9.631    4.233    3<br><br>Best Regards,<br>Limin<br><br><div><span class="gmail_quote">2008/3/10, Teodoro Laino <<a href="mailto:teodor...@gmail.com">teodor...@gmail.com</a>>:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div style="">Dear Limin,<div><br></div><div>Please read carefully the error/warning  messages:</div><span class="q"><div><br></div><div><div><div>On 10 Mar 2008, at 01:36, lmliu wrote:</div><br><blockquote type="cite"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;">*** 01:34:59 WARNING in force_env_methods:mixed_energy_forces err=-300   ***<br>
 *** ASSERTION (cond) failed at line   1151 Error (0.129985E-11) in       ***<br> *** computing numerical derivatives larger then(0.100000E-11) .          ***</span></blockquote></div></div><div><br></div></span><div>First: this is not an error message but it is a warning (WARNING in [..]). </div>
<div>Second: It clearly says that: Error (0.129985E-11) in computing numerical derivatives larger then(0.100000E-11) .</div><div><br></div><div>I guess you can cope with this error on the evaluation of the MIXED forces.. don't you (look at the magnitude of the numbers)?</div>
<div>If you are annoyed from seeing this warning message you can tune the precision of the evaluation of the numerical derivative</div><div>with the keyword:</div><div><br></div><div><a href="http://cp2k.berlios.de/input/InputReference%7E__ROOT__%7EFORCE_EVAL%7EMIXED%7EGENERIC.html#ERROR_LIMIT" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">http://cp2k.berlios.de/input/InputReference~__ROOT__~FORCE_EVAL~MIXED~GENERIC.html#ERROR_LIMIT</a></div>
<div><br></div><div>10^-12 is an extremely safe parameter. I guess 10^-7 will do your work in a very good way as well (please compare these numbers</div><div>with the precision achieved in computing the DFT forces).</div>
<div>Cheers,</div><div>Teo</div><span class="q"><div><br></div><br>
</span></div><br>
</blockquote></div><br>