Dear Nicola<br>Is it related to nstep<br>I chose nstep  = 5000,<br>should it be more or less ?<br><br><div class="gmail_quote">On Thu, Jul 10, 2008 at 12:21 PM, Nicola Marzari <<a href="mailto:marzari@mit.edu">marzari@mit.edu</a>> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="Ih2E3d">Mansoureh Pashangpour wrote:<br>
> dear Nicola<br>
> this is ekinc for end of minimization I found it small after verlet dynamics<br>
>  ekinc         :        0.10064       0.10064 (AU)<br>
>  ekin          :      144.72166     144.72166 (AU)<br>
>    epot          :     -473.07863    -473.07863 (AU)<br>
>    totel energy  :     -280.01088    -280.01088 (AU)<br>
><br>
> why they are different with<br>
>  > scf:  total energy              =  -243.33284101 Ry<br>
<br>
<br>
</div>It is actually (very) large, and it needs to be well below 10-5 au.<br>
<br>
If after that it's still different, try something easy - a hydrogen<br>
molecule, and go from there.<br>
<br>
In any case, the fact that the cp calculation has a much lower energy<br>
than pwscf means that probably something incorrect is going on. As<br>
mentioned, try very simple things, and move up from there.<br>
<br>
<br>
                                nicola<br>
<font color="#888888"><br>
--<br>
</font><div><div></div><div class="Wj3C7c">---------------------------------------------------------------------<br>
Prof Nicola Marzari   Department of Materials Science and Engineering<br>
13-5066   MIT   77 Massachusetts Avenue   Cambridge MA 02139-4307 USA<br>
tel 617.4522758 fax 2586534 <a href="mailto:marzari@mit.edu">marzari@mit.edu</a> <a href="http://quasiamore.mit.edu" target="_blank">http://quasiamore.mit.edu</a><br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br>