<p dir="ltr">I have only done a few -so calculations, but from the list they seem to lead to problems. A few suggestions:</p>
<p dir="ltr">a) Check for user error when you setup the in plane calculation. Have you used the right symmetry reduction? This seems to be the most common issue.</p>
<p dir="ltr">b) Check that you are using the patches for issues in initso that have been posted.</p>
<p dir="ltr">c) Look if your calculation is oscillating or just veeeeeery slowly converging. You can do a grep on some of the moments such as :ORB, :POM, :SPI in case.scf. I have seen these only slowly converge as the occupied states change.</p>
<p dir="ltr">d) Think carefully about the physics of your problem. Some time ago I did a test LDA+U calculation for a Ni (metal) surface. It did not converge, which was OK. In Wien2k you are solving a fixed-point problem, and always assuming that such a solution exists. I do not believe that anything guaruntees that a solution exists in all cases. ( This is different from doing a variational calculation.) Maybe with your distances your impurities want to have an AFM ordering in-plane, but you are preventing this.</p>
<p dir="ltr">These are just guesses, without more information it is hard to know.</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 Feb 12, 2016 06:19, "Islam, Md F" <<a href="mailto:islamm@uta.edu">islamm@uta.edu</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Wien2k users,<br>
<br>
I am trying to do a magnetic impurity calculation on a surface including spin-orbit coupling.<br>
I set the energy and charge convergence to 0.00001 and 0.001, respectively. The calculation is<br>
well converged if magnetization is in-plane but for out of plane magnetization it is not converging.<br>
While charge satisfies convergence criterion very well, the energy keeps oscillating near the 3rd/4th<br>
decimal places.. I tried with different mixing (MSR1: 0.05 to 0.01) but nothing seems to be working<br>
for about 2 months. For larger mixing, job crashes with QTL-B error. I would appreciate if you can<br>
give me some suggestion to resolve the problem.<br>
<br>
I am also wondering if anyone has done any anisotropy calculation with both self-consistent and<br>
non self-consistent methods and whether they give the same anisotropy energy.<br>
<br>
<br>
Thanks,<br>
Fhokrul<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>
</blockquote></div>