Dear Nicola<br>Is it related to nstep<br>I chose nstep&nbsp; = 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 &lt;<a href="mailto:marzari@mit.edu">marzari@mit.edu</a>&gt; 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>
&gt; dear Nicola<br>
&gt; this is ekinc for end of minimization I found it small after verlet dynamics<br>
&gt; &nbsp;ekinc &nbsp; &nbsp; &nbsp; &nbsp; : &nbsp; &nbsp; &nbsp; &nbsp;0.10064 &nbsp; &nbsp; &nbsp; 0.10064 (AU)<br>
&gt; &nbsp;ekin &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;: &nbsp; &nbsp; &nbsp;144.72166 &nbsp; &nbsp; 144.72166 (AU)<br>
&gt; &nbsp; &nbsp;epot &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;: &nbsp; &nbsp; -473.07863 &nbsp; &nbsp;-473.07863 (AU)<br>
&gt; &nbsp; &nbsp;totel energy &nbsp;: &nbsp; &nbsp; -280.01088 &nbsp; &nbsp;-280.01088 (AU)<br>
&gt;<br>
&gt; why they are different with<br>
&gt; &nbsp;&gt; scf: &nbsp;total energy &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;-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&#39;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>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;nicola<br>
<font color="#888888"><br>
--<br>
</font><div><div></div><div class="Wj3C7c">---------------------------------------------------------------------<br>
Prof Nicola Marzari &nbsp; Department of Materials Science and Engineering<br>
13-5066 &nbsp; MIT &nbsp; 77 Massachusetts Avenue &nbsp; 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>