<p dir="ltr">Sorry Peter, but if every cycle is a PRATT step with 0.035 something is wrong. One way this can occur is if the mixer thinks there was no vorb in the previous cycle but there is one now (or the opposite). Then a failsafe Pratt step is used.</p>
<p dir="ltr">I am not certain that a clean user message is printed out, particularly with a older version. There might be something in case.scfm or case.outputm. If nothing is apparent please send me those files and I can probably work it out.</p>
<div class="gmail_quot<blockquote class=" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">case.inorb has nothing to do with the mixer (and is ok).<br>
<br>
Ceckout in more detail what case.scfm says about mixing, but when your<br>
positions change, FRMS gets smaller and finally the forces are close to<br>
zero and the :ENE is lower than in the unrelaxed case, it should be fine.<br>
<br>
On 10/17/2014 02:43 AM, Kortan, Victoria R wrote:<br>
> Thank you very much for the replies.<br>
><br>
> I have deleted the *bro* files and moved the case.scf to case.scf_old,<br>
> and I still get the same result ie all 18 cycles still print:<br>
> ":MIX : PRATT REGULARIZATION: 2.50E-04 GREED: 0.035"<br>
> Maybe there is something wrong with my case.inorb file (I'll add the<br>
> text to the bottom of the e-mail, just in case the error is there).<br>
><br>
> But, as has been pointed out, since :FRMS is being printed in case.scf<br>
> MSR1a must be in use - so is it safe to assume this is sort of a moot point?<br>
><br>
> Thanks,<br>
> ~Victoria<br>
><br>
><br>
> (single Er, atom #65, in Si64 supercell)<br>
> -------------------- case.inorb --------------------<br>
> 1 1 0 nmod, natorb, ipr<br>
> PRATT 1.0 BROYD/PRATT, mixing<br>
> 65 1 3 iatom nlorb, lorb<br>
> 1 nsic 0..AFM, 1..SIC, 2..HFM<br>
> 0.63 0.00 U J (Ry) Note: we recommend to use U_eff = U-J<br>
> and J=0<br>
> ------------------------------------------------------<br>
><br>
> ------------------------------------------------------------------------<br>
> *From:* <a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a><br>
> [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>] on behalf of Laurence Marks<br>
> [<a href="mailto:L-marks@northwestern.edu">L-marks@northwestern.edu</a>]<br>
> *Sent:* Wednesday, October 15, 2014 6:28 PM<br>
> *To:* A Mailing list for WIEN2k users<br>
> *Subject:* Re: [Wien] MSR1a in case.inm but getting PRATT mixing<br>
><br>
> First, MSR1a is a bit better in version 13, so updating would help. (13<br>
> & 14 are about the same, the next one should be better, fingers crossed.)<br>
><br>
> 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<br>
> problem or OS crash. Do mv case.scf case.scf_old<br>
> d) Something else.<br>
><br>
> The mixer thinks there is something inconsistent between the previous<br>
> history files and the current cycle, so is reverting to a failsafe<br>
> Pratt. If none of the above work please send case.scf to my private email.<br>
><br>
> ___________________________<br>
> Professor Laurence Marks<br>
> Department of Materials Science and Engineering<br>
> Northwestern University<br>
> <a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a> <<a href="http://www.numis.northwestern.edu" target="_blank">http://www.numis.northwestern.edu</a>><br>
> <a href="http://MURI4D.numis.northwestern.edu" target="_blank">MURI4D.numis.northwestern.edu</a> <<a href="http://MURI4D.numis.northwestern.edu" target="_blank">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<br>
> nobody else has thought"<br>
> Albert Szent-Gyorgi<br>
><br>
> On Oct 15, 2014 6:01 PM, "Kortan, Victoria R"<br>
> <<a href="mailto:victoria-povilus@uiowa.edu">victoria-povilus@uiowa.edu</a> <mailto:<a href="mailto:victoria-povilus@uiowa.edu">victoria-povilus@uiowa.edu</a>>> wrote:<br>
><br>
> Dear Wien Users,<br>
><br>
> I am running WIEN2k_12.1 and attempting to calculate the total<br>
> energy of an Er atom at an interstitial position in a 2x2x2 Si<br>
> supercell with relaxed atomic positions. This is an LDA+U(for the f<br>
> orbitals of the Er) spin polarized calculation. I have run (as<br>
> 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<br>
> calculation) and it looks like the calculation is headed in the<br>
> correct direction (FRMS decreasing, DIS decreasing).<br>
><br>
> I would like to use the MSR1a mixing method, but am having trouble<br>
> making that happen. I have changed MSR1 to MSR1a in case.inm, but<br>
> 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<br>
> REGULARIZATION...."<br>
><br>
> I'm not sure why I'm not using MSR1a in this case, maybe I'm missing<br>
> something obvious?<br>
> Thanks in advance for any help/hints/advice/insight,<br>
> ~Victoria<br>
><br>
><br>
><br>
><br>
><br>
><br>
><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" 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" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
><br>
<br>
--<br>
<br>
P.Blaha<br>
--------------------------------------------------------------------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300 FAX: +43-1-58801-165982<br>
Email: <a href="mailto:blaha@theochem.tuwien.ac.at">blaha@theochem.tuwien.ac.at</a> WIEN2k: <a href="http://www.wien2k.at" target="_blank">http://www.wien2k.at</a><br>
WWW: <a href="http://www.imc.tuwien.ac.at/staff/tc_group_e.php" target="_blank">http://www.imc.tuwien.ac.at/staff/tc_group_e.php</a><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" 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" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</div>