<div dir="ltr"><div><div><div>Dear Peter <br><span style="color:rgb(255,0,0)">> You should probably repeat this.<br></span></div><div><span style="color:rgb(255,0,0)"><span style="color:rgb(0,0,0)">what's your</span> </span>purpose?</div>
I must repeat which calculations? scf or lapwdm?<br><br>> <span style="color:rgb(255,0,0)">I can hardly believe that changing alpha to 0.1 makes the orbital moment
exactly zero. Even plain GGA should give you some moment </span><br></div><div>why do you have this believe?<br></div> base on my results? or other things?<br></div>thanks<br><div><br></div></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Mon, Dec 9, 2013 at 12:30 PM, Peter Blaha <span dir="ltr"><<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
You should probably repeat this.<br>
<br>
I can hardly believe that changing alpha to 0.1 makes the orbital moment exactly zero. Even plain GGA should give you some moment (at least as some spin-moment is still present)....<div class="im HOEnZb"><br>
<br>
<br>
Dear Peter<br>
thanks for your reply<br>
> Maybe you forgot to add -so for the second calculation ?<br>
no I didn't forget -so.<br>
I check the log file. this file shows that there is -so switch in SCF calculations (I putted a section of this file in my first email).<br>
I run the lapwdm program similar to log that is:<br>
x lapwdm -up -p -so -c<br>
program is done with any error.(i check it again)<br>
but when I check the case.scfdmup see that orb contribution is zero<br>
[yazdani@cm6 case]$ cat case.scfdmup<br>
Spin-polarized + s-o calculation, M|| 1.000 1.000 1.000<br>
Calculation of <X>, X=c*Xr(r)*Xls(l,s)<br>
Xr(r) = I<br>
Xls(l,s) = L(dzeta)<br>
c= 1.00000<br>
atom L up dn total<br>
:XOP001 3 0.00000 0.00000 0.00000 0.00000<br>
:XOP002 3 0.00000 0.00000 0.00000 0.00000<br>
[yazdani@cm6case]$<br>
while when I grep the ORB001 in the scf file see that orbital magnetic moment isn't zero.<br>
I also follow similar process for the alpha=0.2 and see non zero orb contribution in the case.scfdmup.<br>
As I said in the second email I run this program for the Ni as sample and drive non zero orb contribution.<br>
I see this problem for my case with B3Pw91 with alpha=0.1 functional.<br>
<br>
<br>
<br>
<br></div><span class="HOEnZb"><font color="#888888">
-- <br>
<br>
P.Blaha<br>
------------------------------<u></u>------------------------------<u></u>--------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300 FAX: +43-1-58801-165982<br>
Email: <a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a> WWW: <a href="http://info.tuwien.ac.at/theochem/" target="_blank">http://info.tuwien.ac.at/<u></u>theochem/</a><br>
------------------------------<u></u>------------------------------<u></u>--------------</font></span><div class="HOEnZb"><div class="h5"><br>
______________________________<u></u>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.<u></u>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.<u></u>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" target="_blank">http://www.mail-archive.com/<u></u>wien@zeus.theochem.tuwien.ac.<u></u>at/index.html</a><br>
</div></div></blockquote></div><br></div>