dear Nicola<br>this is ekinc for end of minimization I found it small after verlet dynamics<br>&nbsp;ekinc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.10064&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.10064 (AU)<br>&nbsp;ekin&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 144.72166&nbsp;&nbsp;&nbsp;&nbsp; 144.72166 (AU) <br>&nbsp;&nbsp; epot&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp; -473.07863&nbsp;&nbsp;&nbsp; -473.07863 (AU) <br>
&nbsp;&nbsp; totel energy&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp; -280.01088&nbsp;&nbsp;&nbsp; -280.01088 (AU)<br>&nbsp;<br>why they are different with<br>
&gt; scf: &nbsp;total energy &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;-243.33284101 Ry<br>Thanks<br><br><div class="gmail_quote">On Thu, Jul 10, 2008 at 11:50 AM, 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;">
<br>
<br>
Shouldn&#39;t ekinc be zero ? I forgot the name used, but in a properly<br>
executed CP calculation, in which you do damped dynamics on the<br>
electrons to reach the ground state at fixed ions, the kinetic energy<br>
of the electrons (as opposed to the quantum kinetic energy of your<br>
system, ekin) should be zero.<br>
<br>
I suggest to anyone learning CP calculations (a very worthy<br>
endeavour, and one that would be extremely useful in the long<br>
run) to study carefully the review paper on the topic by Marx and<br>
Hutter,<br>
<br>
<a href="http://www.theochem.ruhr-uni-bochum.de/research/marx/marx.pdf" target="_blank">http://www.theochem.ruhr-uni-bochum.de/research/marx/marx.pdf</a><br>
<br>
CP in Q-E now can also do conjugate gradients minimization on<br>
the electrons (more efficient than damped dynamics) and, for<br>
metals, ensemble-dft (i.e. born-oppenheimer dynamics).<br>
<br>
For cg/edft, there was one bug (memmory misallocation)<br>
for parallel execution, just found by by Brandon Wood. It&#39;ll<br>
be fixed in the cvs in the next few days.<br>
<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;nicola<br>
<div class="Ih2E3d"><br>
Mansoureh Pashangpour wrote:<br>
&gt; Dear Paolo<br>
&gt; I know one of them is in Ry and another one in Ha,but when I compare<br>
&gt; them they are different<br>
&gt; scf: &nbsp;total energy &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;-243.33284101 Ry<br>
&gt;<br>
&gt; CP: &nbsp;ekinc &nbsp; &nbsp; &nbsp; &nbsp; : &nbsp; &nbsp; &nbsp; 10.67409 &nbsp; &nbsp; &nbsp;10.67409 (AU)<br>
&gt; &nbsp; &nbsp;ekin &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;: &nbsp; &nbsp; &nbsp;153.58930 &nbsp; &nbsp; 153.58930 (AU)<br>
&gt; &nbsp; &nbsp;epot &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;: &nbsp; &nbsp; -458.58466 &nbsp; &nbsp;-458.58466 (AU)<br>
&gt; &nbsp; &nbsp;totel energy &nbsp;: &nbsp; &nbsp; -258.18270 &nbsp; &nbsp;-258.18270 (AU)<br>
&gt; &nbsp; &nbsp;temperature &nbsp; : &nbsp; &nbsp; &nbsp; &nbsp;0.00000 &nbsp; &nbsp; &nbsp; 0.00000 (K )<br>
&gt; &nbsp; &nbsp;enthalpy &nbsp; &nbsp; &nbsp;: &nbsp; &nbsp; -258.18270 &nbsp; &nbsp;-258.18270 (AU)<br>
&gt; &nbsp; &nbsp;econs &nbsp; &nbsp; &nbsp; &nbsp; : &nbsp; &nbsp; -258.18270 &nbsp; &nbsp;-258.18270 (AU)<br>
&gt; &nbsp; &nbsp;pressure &nbsp; &nbsp; &nbsp;: &nbsp; &nbsp; &nbsp; -3.47982 &nbsp; &nbsp; &nbsp;-3.47982 (Gpa)<br>
&gt; &nbsp; &nbsp;volume &nbsp; &nbsp; &nbsp; &nbsp;: &nbsp; &nbsp; &nbsp;488.11442 &nbsp; &nbsp; 488.11442 (AU)<br>
&gt; Thanks<br>
<br>
</div><font color="#888888">--<br>
---------------------------------------------------------------------<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>
</font><div><div></div><div class="Wj3C7c">_______________________________________________<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>