[Wien] MSR1a in case.inm but getting PRATT mixing

Kortan, Victoria R victoria-povilus at uiowa.edu
Thu Oct 16 01:01:21 CEST 2014


Dear Wien Users,

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)):
"runsp_lapw -orb -cc 0.1 -ec 0.25 -fc 20 -p"
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).

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:
":MIX  :   PRATT  REGULARIZATION: 2.50E-04  GREED: 0.035"
instead of what is listed in the mixer ReadMe, ie ":MIX :   MSR1a REGULARIZATION...."

I'm not sure why I'm not using MSR1a in this case, maybe I'm missing something obvious?
Thanks in advance for any help/hints/advice/insight,
~Victoria







-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20141015/75cd1956/attachment.html>


More information about the Wien mailing list