<p dir="ltr">Me a culpa, I should have checked the mailing list first for the answers.</p>
<p dir="ltr">That said, this issue has come up enough times in the past that I think the UG should be tweaked so it is clearer. Let me try my interpretation, so I can be shot down if needed.</p>
<p dir="ltr">Within Wien2k magnetic effects can be approximately included in a number of ways. Some such as the spin-orbit coupling assume a direction for the spin vector (for all electrons actively considered), others such as Bext in orb specify a direction for an applied magnetic field (in Tesla) and use the same direction for the spin vector. (The two spin states are then either parallel or anti parallel to the specified direction.) When a direction is specified in case.inso or case.inorb this fixes the spin vector and (if used) the external magnetic field direction. Via the output files from lapwdm (case.scfdmXX) one can monitor how the angular momentum changes [1]. By using different directions for the spin vector (and field) one can probe how the energy changes and/orbital occupancies with assumed directions for the spin/external magnetic field. </p>
<p dir="ltr">To escape the assumption that the spin vectors all have one direction the Wienncm code has to be used.</p>
<p dir="ltr">[1] My addendum. Changes in the occupancies can be a soft electronic mode, i.e. very small changes in the energy for quite large changes in the density. The older mixing algorithms (MSEC1 or MSEC3) are not so good for soft modes and can stagnate. MSR1 is better and the next release (7.0) is much better. With onsite -eece &/or -orb it may help to "push" the mixer by either forcing a larger step (echo .2 > .msec or echo .1 > .pratt) or stopping, doing a force on the orbital potential (x orb -up; x orb -dn) then restarting with -NI. It is probably wise to check how the orbital momentum is converging (grep :ORB0 *.scf, perhaps other) and make sure that the mixing is not starving (grep GREED: *.scf and check the values are not small, e.g. 0.035).</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 May 4, 2015 6:22 PM, "Laurence Marks" <<a href="mailto:L-marks@northwestern.edu">L-marks@northwestern.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 dir="ltr">Typo:<div><br></div><div>"<span style="color:rgb(0,0,0);font-size:12.8000001907349px">although I remember don't symmetry operations being split into these two classes everywhere in the code"</span></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, May 4, 2015 at 6:04 PM, Laurence Marks <span dir="ltr"><<a href="mailto:L-marks@northwestern.edu" target="_blank">L-marks@northwestern.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">I am a newbie at -so, so a few simple questions.<div><br></div><div>a) What is the meaning of the orbital moment in case.scfdm* ? Is that the average direction projected to the global axis system?</div><div><br></div><div>b) What is the physical significance of the orbital moment being parallel (or not quite parallel) to the direction used in case.inso?</div><div><br></div><div>c) I understand that the results for different directions of B in case.inso reflect the magnetic anisotropy, but what are the units of field (if any)?</div><div><br></div><div>d) What else is worth looking at? The partial orbital moment (:POM) seems relevant, but what exactly is it?</div><div><br></div><div>e) I am "blindly" trusting that initso knows what it is doing, and have left the "B" symmetry operations in case.struct (although I remember symmetry operations being split into these two classes everywhere in the code). This seems to conflict with Pavel's notes, although those may be too old.</div><div><br></div><div>Thanks.<span><font color="#888888"><br clear="all"><div><br></div>-- <br><div><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>
</font></span></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div><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>
</blockquote></div>