<div dir="auto">Probably overlapping response from Peter...<div dir="auto"><br></div><div dir="auto">Based upon the Check-mixing output, the density is very well converged and did not change with the jump in energy. The most likely way for :DEN alone to change, and nothing else, is that there was a change in how EXC is calculated, for instance using a metaGGA or changing case.in0, without a change in VXC.</div><div dir="auto"><br></div><div dir="auto">Do "grep :POT *scf" to check this.<br><br><div data-smartmail="gmail_signature" dir="auto">_____<br>Professor Laurence Marks<br>"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Gyorgi<br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Dec 7, 2017 7:19 AM, "Pavel Ondračka" <<a href="mailto:pavel.ondracka@email.cz" target="_blank">pavel.ondracka@email.cz</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Peter Blaha píše v Čt 07. 12. 2017 v 13:34 +0100:<br>
> As the charge distance is not going bad, I don't expect it has to do<br>
> with mixing, nor with any scf problem.<br>
><br>
> I'd guess, only E-tot is calculated wrongly in mixer or the run<br>
> script<br>
> has a problem .<br>
> It could be that you have 2 times the same iteration number in the<br>
> scf<br>
> file (:ITE) or that for some reasons once contribution to :ENE is<br>
> not<br>
> listed/listed twice, ...<br>
><br>
> Please check the scf file for :ITE, but also manually for the<br>
> different<br>
> E-tot contributions like (grep :ENE case.scf)<br>
><br>
> :ENE<br>
> :DEN<br>
> :SUM<br>
> and the that the core eigenvalues are listed only once.<br>
><br>
> ALSO: save the calculation, and continue with another run_lapw<br>
> Eventually it should stop after 3 cycles, because everything is<br>
> converged.<br>
><br>
Thanks for the advice,<br>
so the :ITE check did not spot any problems, the :SUM also stays the<br>
same, however the :DEN show the same fluctuations as the total energy,<br>
like this...<br>
<br>
:DEN  : DENSITY INTEGRAL  =        -15670.17810081   (Ry)<br>
:DEN  : DENSITY INTEGRAL  =        -15614.53403108   (Ry)<br>
:DEN  : DENSITY INTEGRAL  =        -15614.55064266   (Ry)<br>
:DEN  : DENSITY INTEGRAL  =        -15614.55963292   (Ry)<br>
:DEN  : DENSITY INTEGRAL  =        -15670.17573346   (Ry)<br>
:DEN  : DENSITY INTEGRAL  =        -15670.17527628   (Ry)<br>
:DEN  : DENSITY INTEGRAL  =        -15670.17574871   (Ry)<br>
<br>
<br>
I'm not sure how to check that "core eigenvalues are listed only once"<br>
but there is nothing suspicious at the first glance.<br>
<br>
Also attached is the output from the Check-mixed as requested by prof.<br>
Marks.<br>
<br>
I saved the calculations as advised and restarted. However I've done<br>
this already several times (I've also tried clean_lapw) without any<br>
effect.<br>
<br>
Best regards<br>
Pavel</blockquote></div></div>