<div dir="ltr"><div><div><div><div><div><div>Dear Prof. Marks,<br><br></div>   Thank you very much for your explanation.<br></div>   Unfortunatelly, I am stil using Wien2k 14.2, without the "Check" command.<br></div>   Anyway, I can wait for the charge convergence (for a while). If it takes too long and nothing bad happens, I will "kill" the optimization.<br></div>   Thank you again.<br></div>                All the best,<br></div>                             Luis<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">2017-07-11 13:13 GMT-03:00 Laurence Marks <span dir="ltr"><<a href="mailto:L-marks@northwestern.edu" target="_blank">L-marks@northwestern.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">It is possible for :ENE (Energy, -ec) to seem to converge, but has not<br>
-- so it is always wise IMO to also use the charge converge (:DIS,<br>
-cc). A value of .0004130 is probably OK, I suggest using the "Check"<br>
command which gives more information. (Peter added it to version 16,<br>
presumably 17 as well as a command Check-mixing.)<br>
<br>
>From the variational principle, the density converges more slowly than<br>
the energy. Depending upon what you are trying to do your charge<br>
convergence may be fine. For something like phonons it probably is<br>
not. I will note that very good charge convergence often requires good<br>
k-mesh, RKMAX etc, more than is needed for a routine calculation.<br>
<br>
N.B., with more atoms charge convergence can be harder to achieve --<br>
and less important.<br>
<div class="HOEnZb"><div class="h5"><br>
On Tue, Jul 11, 2017 at 9:40 AM, Luis Ogando <<a href="mailto:lcodacal@gmail.com">lcodacal@gmail.com</a>> wrote:<br>
> Dear Prof. Marks and Wien2k community,<br>
><br>
>    Good afternoon !<br>
>    Just a curiosity about the charge convergence criteria. I am optimizing a<br>
> supercell and the charge convergence limit I am using is 0.001. The same<br>
> numerical value for the energy convergence.<br>
>    The energy is finally converged, but the charge is not despite the fact<br>
> that its change is below the limit in the last 10 iterations.<br>
><br>
> CarbazoleLDAopt.dayfile::<wbr>CHARGE convergence:  1 0.001 -.0002028<br>
> CarbazoleLDAopt.dayfile::<wbr>CHARGE convergence:  0 0.001 .0006007<br>
> CarbazoleLDAopt.dayfile::<wbr>CHARGE convergence:  1 0.001 -.0002590<br>
> CarbazoleLDAopt.dayfile::<wbr>CHARGE convergence:  0 0.001 .0007240<br>
> CarbazoleLDAopt.dayfile::<wbr>CHARGE convergence:  0 0.001 .0008640<br>
> CarbazoleLDAopt.dayfile::<wbr>CHARGE convergence:  0 0.001 .0005344<br>
> CarbazoleLDAopt.dayfile::<wbr>CHARGE convergence:  0 0.001 .0005217<br>
> CarbazoleLDAopt.dayfile::<wbr>CHARGE convergence:  0 0.001 .0002525<br>
> CarbazoleLDAopt.dayfile::<wbr>CHARGE convergence:  0 0.001 .0004560<br>
> CarbazoleLDAopt.dayfile::<wbr>CHARGE convergence:  0 0.001 .0004130<br>
><br>
>    If I am not wrong, the convergence is achieved after 3 consecutive<br>
> iterations with changes below the limit, but in my case this is happening<br>
> for more than 10 iterations  without convergence.<br>
>    Could you please kindly explain what may be happening ?<br>
>    Thank you !<br>
>                     Luis<br>
><br>
<br>
<br>
<br>
</div></div><span class="HOEnZb"><font color="#888888">--<br>
Professor Laurence Marks<br>
"Research is to see what everybody else has seen, and to think what<br>
nobody else has thought", Albert Szent-Gyorgi<br>
<a href="http://www.numis.northwestern.edu" rel="noreferrer" target="_blank">www.numis.northwestern.edu</a> ; Corrosion in 4D: <a href="http://MURI4D.numis.northwestern.edu" rel="noreferrer" target="_blank">MURI4D.numis.northwestern.edu</a><br>
Partner of the CFW 100% program for gender equity, <a href="http://www.cfw.org/100-percent
Co-Editor" rel="noreferrer" target="_blank">www.cfw.org/100-percent<br>
Co-Editor</a>, Acta Cryst A<br>
______________________________<wbr>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.<wbr>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.<wbr>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/<wbr>wien@zeus.theochem.tuwien.ac.<wbr>at/index.html</a><br>
</font></span></blockquote></div><br></div>