<p dir="ltr">Let me make a comment on this, with the caveat to everyone that Peter & I do not completely agree on this (which is OK). A few subtle, but important points.</p>
<p dir="ltr">With DFT we are looking for a fixed point, i.e. the density (and if relevant density matrix or orbital potential) is the same at the start and finish of a single iteration. For simple functionals this is garunteed to exist from the KS theorem.</p>
<p dir="ltr">However, with more complex functionals it is not garunteed to exist in all cases. If Wien2k did a variational minimum one would at least get the lowest energy, but a  fixed point method is not completely the same.</p>
<p dir="ltr">There is a whole world of more complex phenomena which can occur with fixed point methods, and a start is the Wikipedia page on attractors at <a href="https://en.m.wikipedia.org/wiki/Attractor">https://en.m.wikipedia.org/wiki/Attractor</a> . The topic is really rather complex.</p>
<p dir="ltr">The take home message is that you have to check not just :DIS, but also other terms to make sure that you really have found the fixed point and are not stuck at some more complex quasi solution. There is an undocumented parameter :MV which can be checked for orbital dependent parameters; it should be less than 1% for a decent solution in my estimation.</p>
<p dir="ltr">N.B., the next mixer release will have an additional option "PUSH" which takes some steps to try and reduce some of these issues. Unfortunately if the underlying problem is chaotic (in the mathematical sense) it becomes very complicated and I don't have simple solutions.</p>
<p dir="ltr">---<br>
Professor Laurence Marks<br>
Department of Materials Science and Engineering<br>
Northwestern University<br>
<a href="http://www.numis.northwestern.edu">http://www.numis.northwestern.edu</a><br>
Corrosion in 4D <a href="http://MURI4D.numis.northwestern.edu">http://MURI4D.numis.northwestern.edu</a><br>
Co-Editor, Acta Cryst A<br>
"Research is to see what everybody else has seen, and to think what nobody else has thought"<br>
Albert Szent-Gyorgi</p>
<div class="gmail_quote">On Aug 3, 2015 00:52, "Peter Blaha" <<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Then the dmat and the scf-files do not fit together !<br>
<br>
I've noted myself that sometimes the mixer does not properly "mix" the<br>
dmat files, i.e. the new dmat files and the mixed ones disagree and do<br>
not converge to the same numbers.<br>
<br>
To check if this happens in your case:<br>
<br>
compare the dmats of case.scfdmup, case.dmatup/dn_unmixed  (they MUST<br>
agree) and case.dmatup/dn.<br>
<br>
If the unmixed dmats contain values close to 1, but the normal (mixed)<br>
dmats don't, you have such a case. Then do:<br>
<br>
save_lapw calc_not_fully_converged<br>
cp case.dmatup_unmixed case.dmatup   and dn<br>
another runsp -so<br>
<br>
Of course, l=0 means angular momentum l=0 ("s").<br>
<br>
<br>
Am 03.08.2015 um 02:01 schrieb Bin Shao:<br>
> The qtl lines for spin up is<br>
><br>
> :QTL013: 1.0036 2.8012 0.2016 6.8976 0.9330 0.9328 0.9340 0.0384 0.0388<br>
> 0.0388 0.0436 0.0436<br>
>          Q-s-low E-s-low   Q-p-low E-p-low   Q-d-low E-d-low   Q-f-low<br>
> E-f-low<br>
><br>
> There is 7 f electrons. Why the occupation number in dmat files is less<br>
> than 7 (about 4)?<br>
><br>
> The calculation got a magnetic solution with 7 muB and the dmat was from<br>
> a GGA+SOC.<br>
><br>
> "I usually use the magnetic field only together with l=0."<br>
><br>
> what's the meaning of l? the quantum number of angular moment?<br>
><br>
><br>
> --<br>
> Bin Shao<br>
> Postdoc<br>
> Department of Physics, Tsinghua University<br>
> Beijing 100084, P. R. China<br>
> Email: <a href="mailto:binshao1118@gmail.com">binshao1118@gmail.com</a> <mailto:<a href="mailto:binshao1118@gmail.com">binshao1118@gmail.com</a>><br>
><br>
><br>
> _______________________________________________<br>
> Wien mailing list<br>
> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.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/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
><br>
<br>
--<br>
Peter Blaha<br>
Inst.Materials Chemistry<br>
TU Vienna<br>
Getreidemarkt 9<br>
A-1060 Vienna<br>
Austria<br>
+43-1-5880115671<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.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/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>