Dear Prof. Blaha,<div><br></div><div>Thanks for the answers. <br><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Either you change the x_lapw script, or, more simple,<br>


create the uplapw1.def and dnlapw1.def files once   (x lapw1 -up/-dn -d)<br>
edit up/dnlapw1.def ach change the path to the output1up/dn file<br>
run:  lapw1(c) uplapw1.def</blockquote><div><br></div><div>I found the line where the def file was created, it was a simple matter to change the line to make the output stay in the scratch directory. However, I never noticed that one could only create the def file - my fault here...</div>

<div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div class="im"><br></div>
You don&#39;t need to run initso for k-convergence tests again. Just do it once, and afterwards do:  x kgen -so</blockquote><div><br></div><div>I actually tried this after running a complete initso cycle, but it gave me an error, saying it could not open the file case.ksym .</div>

<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div class="im"></div>
If you know the magnetization direction from exp, then use it in your calc. Otherwise you should in principle test various possibilities.<br>
Fortunately, when the elements are not too heavy and the structure is not very anisotropic, you may not see much difference in a &quot;bandstructure&quot;, since the magnetocristalline anisotropy usually ranges from a few  mRy to even only a few microRy !!<br>

</blockquote><div><br></div><div>Just a bit more of bothering on this topic: are there any guidelines or principles for choosing it, so as to restrict the number of possibilities? In my case, I don&#39;t have any experimental input, since it&#39;s a material we are proposing theoretically. Is there any way of finding out the direction of the magnetization in the Wien2K output files when spin-orbit is not used, so that I could use it as a first approximation to the magnetization axis?</div>

<div><br></div><div>Thanks,</div><div><br></div><div>Marcos</div><div><div><br></div><div>---------------------------</div><div>Marcos Verissimo Alves</div><div>Post-doctoral Researcher</div><div>Universidad de Cantabria</div>

<div>Santander, Spain</div></div></div></div>