<div dir="ltr">Yup, this can happen. A few ideas:<div><br></div><div>a) Turn off MSR1a for a bit and converge somewhat the density. It can be dangerous to run MSR1a when :DIS is large (e.g. 2.0). Don't overconverge, e.g. just bring it down to something like -cc 0.1.</div><div><br></div><div>b) Add an f-state in case.in1 high above E_F, e.g. at 1.5 to try and avoid getting the f states occupied.</div><div><br></div><div>A "deep" problem is whether you have a reasonable model for the oxide surface, or if it is physically unrealistic. A structure which is unreasonable, for instance would require La2+ for valence neutrality will never converge to anything reasonable in my experience.</div><div><br></div><div>And, as I mentioned before, structures with inversion symmetry in my experience are much better behaved than those without. While you may have more atoms because you only have real components of the PW's they can be smaller calculations than ones without inversion symmetry.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Nov 19, 2015 at 11:27 AM, Elias Assmann <span dir="ltr"><<a href="mailto:elias.assmann@gmail.com" target="_blank">elias.assmann@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi List,<br>
<br>
I have a big calculation running with MSR1a (one of those oxide surfaces<br>
I asked about earlier).  After a while, it stopped due to touching<br>
spheres, so I used clminter with `setrmt -r 3' and restarted.  After 130<br>
more iterations, the calculation stops with a QTL-B error<br>
<br>
   QTL-B VALUE .EQ.    9.42627   in Band of energy   -0.21251   ATOM=<br>
129   L=  3<br>
:WARN : QTL-B value eq.  29.45 in Band of energy  -0.35946  ATOM=  134<br>
L=  3<br>
<br>
Now, atom 134 is a La and its f-states should be *empty*.  And so they<br>
were for most of the calculation, but ~30 iterations after clminter,<br>
they start to be full (they pretty much jump from 0 to 7/spin, here is a<br>
plot: <<a href="http://i.imgur.com/DZavhiq.png" rel="noreferrer" target="_blank">http://i.imgur.com/DZavhiq.png</a>>; the scale goes to 7).<br>
<br>
This actually runs (with some oscillations) for ~100 iterations, then<br>
the QTL-B error happens.<br>
<br>
What could be the cause of this?<br>
<br>
The R-MT for La is 2.28, before the clminter it was 2.35.  ‘scf1’ section:<br>
<br>
          ATOMIC SPHERE DEPENDENT PARAMETERS FOR ATOM  La lo<br>
:e__0134: OVERALL ENERGY PARAMETER IS    0.2434<br>
          OVERALL BASIS SET ON ATOM IS LAPW<br>
:E0_0134: E( 0)=   -2.9884   E(BOTTOM)=   -3.787   E(TOP)=   -2.190  4<br>
5   177<br>
             APW+lo<br>
:E0_0134: E( 0)=    0.6434<br>
             LOCAL ORBITAL<br>
:E1_0134: E( 1)=   -1.5836   E(BOTTOM)=   -2.710   E(TOP)=   -0.457  3<br>
4   218<br>
             APW+lo<br>
:E1_0134: E( 1)=    0.6434<br>
             LOCAL ORBITAL<br>
:E2_0134: E( 2)=    0.2434   E(BOTTOM)=   -0.980   E(TOP)= -200.000  2<br>
-1   175<br>
             APW+lo<br>
<br>
<br>
Thanks,<br>
<br>
        Elias<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
Elias Assmann<br>
Institute of Theoretical and Computational Physics<br>
TU Graz                   ⟨<a href="https://itp.tugraz.at/" rel="noreferrer" target="_blank">https://itp.tugraz.at/</a>⟩<br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr">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><div>Corrosion in 4D: <a href="http://MURI4D.numis.northwestern.edu" target="_blank">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</div></div></div>
</div>