<p dir="ltr">First, MSR1a is a bit better in version 13, so updating would help. (13 & 14 are about the same, the next one should be better, fingers crossed.)</p>
<p dir="ltr">There are four ways this might happen (maybe)<br>
a) You have old *bro* files around. Try deleting them.<br>
b) Something is wrong with your inorb so your U is not working right.<br>
c) Your case.scf file has some strange characters in it from an MPI problem or OS crash. Do mv case.scf case.scf_old<br>
d) Something else.</p>
<p dir="ltr">The mixer thinks there is something inconsistent between the previous history files and the current cycle, so is reverting to a failsafe Pratt. If none of the above work please send case.scf to my private email.</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">www.numis.northwestern.edu</a><br>
<a href="http://MURI4D.numis.northwestern.edu">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 Oct 15, 2014 6:01 PM, "Kortan, Victoria R" <<a href="mailto:victoria-povilus@uiowa.edu">victoria-povilus@uiowa.edu</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">Dear Wien Users,<br>
<br>
I am running WIEN2k_12.1 and attempting to calculate the total energy of an Er atom at an interstitial position in a 2x2x2 Si supercell with relaxed atomic positions.  This is an LDA+U(for the f orbitals of the Er) spin polarized calculation.  I have run (as
 suggested in the mixer ReadMe (VER5.3)):<br>
"runsp_lapw -orb -cc 0.1 -ec 0.25 -fc 20 -p"<br>
with a 1x1x1 k-mesh and RKMAX = 7 (just to get a feel for the calculation) and it looks like the calculation is headed in the correct direction (FRMS decreasing, DIS decreasing).<br>
<br>
I would like to use the MSR1a mixing method, but am having trouble making that happen.  I have changed MSR1 to MSR1a in case.inm, but when I do "grep :MIX case.scf" I get lines like:<br>
":MIX  :   PRATT  REGULARIZATION: 2.50E-04  GREED: 0.035"<br>
instead of what is listed in the mixer ReadMe, ie ":MIX :   MSR1a REGULARIZATION...."<br>
<br>
I'm not sure why I'm not using MSR1a in this case, maybe I'm missing something obvious?<br>
Thanks in advance for any help/hints/advice/insight,<br>
~Victoria<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
</div>
</div>
</div>

</blockquote></div>