<div dir="ltr">Hi Chris,<div><br></div><div> Adding this time to Jürg's answer, rigorously the total energy _is_ only defined at self-consistency, and then the Hohenberg-Kohn theorem holds and you should have a the lowest possible energy. The mixing causes the "estimates" of the total energy (what is printed out at each SCF step) to be a bit more "confused" because the density and orbitals are even less consistent. The problem is usually worst in metallic systems where the density oscillations occur the easiest (smallest "excitation" energies for mixing up the electronic structure).</div>
<div><br></div><div> Another sources of small errors would be from the fractional occupation numbers (again, only in metallic systems) or due to numerical inaccuracies (that is why I suggested checking those parametres like the cut-off energy).</div>
<div><br></div><div> Greetings,</div><div><br></div><div> apsi</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2014/1/14 Chris E <span dir="ltr"><<a href="mailto:christoph...@gmail.com" target="_blank">christoph...@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Marcella and Ari,<div><br></div><div>Thank you for the information. I have already tried tightening the EPS_SCF cutoff but not the EPS_DEFAULT. I will try that as well. </div>
<div>Maybe I should review my DFT books to find this, but its worth asking: if tightened cutoffs produce the same results, is the higher energy calculated earlier in the SCF run a false energy due to mixing? Or does the mixing simply allow the energy to be non-variational? </div>
<div><br></div><div>Thanks again,<br>Chris<br><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><br></div></div></blockquote></div></div><span class="HOEnZb"><font color="#888888">
<p></p>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/groups/opt_out" target="_blank">https://groups.google.com/groups/opt_out</a>.<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-<br> Ari P Seitsonen / <a href="mailto:Ari.P.S...@iki.fi">Ari.P.S...@iki.fi</a> / <a href="http://www.iki.fi/~apsi/">http://www.iki.fi/~apsi/</a><br>
Physikalisch-Chemisches Institut der Universität Zürich<br> Tel: +41 44 63 55 44 97 / Mobile: +41 79 71 90 935
</div>