<div>Dear all,<br><br>I intend to calculate the orbital moment and MAE by using Force theorem. First I did a self-consistent calculation then did the non-scf calculation with different magnetization directions by using the commands as following.<br>

<br>x lapwso -up -p -orb -c<br>x lapw2 -up -p -so -c<br>x lapw2 -dn -p -so -c<br>x lapwdm -up -p -so -c<br><br>(-orb for using GGA+U in the system)<br><br>In this way, I obtain the orbital moment in the case.scfdmup file, for example<br>

<br>----------------------------------------------------------------------------------------------------------------------------<br>:ORB001:  ORBITAL MOMENT: -0.00187  0.00855  0.49744 PROJECTION ON M  0.49744<br>:SPI001:  SPIN MOMENT:  -0.00002   0.00008   1.90505 PROJECTION ON M  1.90505<br>

----------------------------------------------------------------------------------------------------------------------------<br><br>however, I found there is a difference (about 0.5 muB) between the spin moment in case.scfdmup and case.scf produced by the scf calculation.<br>

<br>----------------------------------------------------------------------------------------------------------------------------<br>:MMI001: MAGNETIC MOMENT IN SPHERE   1    =    2.45641<br>----------------------------------------------------------------------------------------------------------------------------<br>

<br>I have searched the maillist and found the answer to the difference provided by Prof. Novak<br><div></div><div></div><blockquote><blockquote style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote">

 in <span style class="J-JK9eJ-PJVNOc">LAPWDM</span> you calculate the spin moment from selected electrons only (usually d or f),<br>while moment in the sphere is the sum from all electrons: you can check it by running<span style class="J-JK9eJ-PJVNOc"><br>

LAPWDM</span> for all s,p, d and f states.<br></blockquote></blockquote></div><div>But I think the spin moment from s, p states would not larger than 0.5 muB. So how the difference comes from? spin orbit coupling?<br><br>

Any suggestion will be appreciated, thank you in advanced.<br><br>Best,<br></div><div><br></div><div><br></div>-- <br>Bin Shao, Ph.D. Candidate<br>College of Information Technical Science, Nankai University<br>94 Weijin Rd. Nankai Dist. Tianjin 300071, China<br>

Email: <a href="mailto:bshao@mail.nankai.edu.cn" target="_blank">bshao@mail.nankai.edu.cn</a><br>
<br>