<DIV><FONT size=3>Dear Dr. Stefaan:</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; First, thank you very much for your helpful and enlightened reply, and your reply this time is very useful to me.</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; I think there was a misunderstanding about my first question:</FONT></DIV>
<DIV><FONT size=3>case.scfdmup<BR>&nbsp;*************************************************************<BR>&nbsp;Spin-polarized + s-o calculation, M||&nbsp; 0.000&nbsp; 0.000&nbsp; 1.000<BR>&nbsp;&nbsp; Calculation of &lt;X&gt;, X=c*Xr(r)*Xls(l,s)<BR>&nbsp;&nbsp; Xr(r)&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I<BR>&nbsp;&nbsp; Xls(l,s) = L(dzeta)<BR>&nbsp;&nbsp; c=&nbsp; 1.00000<BR>&nbsp;&nbsp; atom&nbsp;&nbsp; L&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; up&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; dn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; total<BR>&nbsp;:XOP&nbsp; 1&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp; 0.00332&nbsp;&nbsp;&nbsp;&nbsp; 1.25294&nbsp;&nbsp;&nbsp;&nbsp; 1.25626<BR>&nbsp;*************************************************************</FONT></DIV>
<DIV><FONT size=3></FONT>&nbsp;</DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; (1)My question is that<FONT color=#ff0066><STRONG> why</STRONG></FONT> the direction of orbital moment<FONT color=#ff0066><STRONG> is oppsite to</STRONG></FONT> the spin moment ? (spin moment is up, orbital moment is down, and I know the third rules of Hund's rules)</FONT></DIV>
<P><FONT size=3></FONT>&nbsp;</P>
<P><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;(2)Is setting or changing the occupational numbers of electrons in case.dmatdn reasonable ? ( I mean we are try to simulate the microscopic behavior by First principles, we should not change or set any parameters in order to get the result we want.)&nbsp; I don't know why my understanding is right or not ,please correct it if I am wrong, thank you.</FONT></P>
<P><FONT size=3></FONT>&nbsp;</P>
<P><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;(3)I checked the DOS by GGA or GGA+SO, the f_up states are always localized with or without plus U. So my understanding is that plus U can only redistribute the electrons of f_dn states, or&nbsp;localize the f_dn states electrons&nbsp;? Am I right ?</FONT></P>
<P><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;</FONT></P>
<P><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;(4)I have read a book "Condensed Matter Physics" before, it said spin_orbit spliting energy is about 300 meV, so&nbsp;what the effect of SO on DOS ? </FONT></P>
<P><FONT size=3>&nbsp; </FONT></P>
<P><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp; Thank you very much, Stefaan. (Anyway, I will try to change the dmatup files and recalculate again first)</FONT></P>
<P><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp; Cheers.</FONT></P>
<P><FONT size=3></FONT>&nbsp;</P>
<P><FONT size=3>Your sincerely:</FONT></P>
<P><FONT size=3>Hui Wang</P>
<DIV><BR></DIV></FONT> <br><!-- urlfiles --><br><br><!-- footer --><br>