<DIV>
<DIV><FONT color=#000000 size=4>Dear&nbsp;Dr.Stefaan:</FONT></DIV>
<DIV><FONT color=#000000 size=4>&nbsp;&nbsp;&nbsp; First,thank you so much for your kindly reply, after reading some mailing list I think I made some mistake during the calculation. After updating some parameters,&nbsp;my computaional detail are as follows(the same as last time are omited:</FONT></DIV>
<DIV><FONT size=4></FONT>&nbsp;</DIV>
<DIV><FONT color=#000000 size=5>order: GGA==&gt;&gt;GGA+SO==&gt;&gt;GGA+SO+U</FONT></DIV>
<DIV>
<DIV>
<DIV><FONT size=3></FONT>&nbsp;</DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;After a&nbsp;regular GGA+SO scf, commamd line: <FONT color=#0000ff><STRONG>runsp -ec 0.0001 -cc 0.0001 -so -i 1000</STRONG></FONT> 
<DIV><FONT color=#ff00ff size=3><FONT color=#000000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; then i got some results as follows:</FONT></FONT></DIV>
<DIV><FONT color=#ff00ff size=3><FONT color=#000000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :ENE:&nbsp;&nbsp;&nbsp; -46876.608343 Ry</FONT></FONT></DIV>
<DIV><FONT color=#ff00ff size=3><FONT color=#000000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;:MMTOT:&nbsp;&nbsp;11.92563&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; uB</FONT></FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :MMI001: 5.81891&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; uB</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;:MMINT:&nbsp; 0.28781 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; uB</FONT></DIV></FONT></DIV></DIV></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<FONT color=#ff0000>&nbsp;<STRONG></STRONG></FONT></FONT></DIV>
<DIV><STRONG><FONT color=#ff0000 size=3>&nbsp;&nbsp;&nbsp;&nbsp; <FONT color=#0000ff>I changed the case.indmc with '1 3' and x lapwdm -c -so -up, and I got the following result:</FONT></FONT></STRONG></DIV>
<DIV><STRONG><FONT color=#0000ff size=3></FONT></STRONG>&nbsp;</DIV>
<DIV><FONT color=#000000 size=3>case.scfdmup</FONT></DIV>
<DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV></DIV>
<DIV><FONT color=#ff0000><FONT color=#000000 size=2>Spin-polarized + s-o calculation, M||&nbsp; 0.000&nbsp; 0.000&nbsp; 1.000<BR>&nbsp; Calculation of &lt;X&gt;, X=c*Xr(r)*Xls(l,s)<BR>&nbsp; Xr(r)&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I<BR>&nbsp; Xls(l,s) = L(dzeta)<BR>&nbsp; c=&nbsp; 1.00000<BR>&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>:XOP&nbsp; 1&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp; 0.00332&nbsp;&nbsp;&nbsp;&nbsp; 1.25294&nbsp;&nbsp;&nbsp;&nbsp; 1.25626</FONT></FONT></DIV><FONT color=#ff0000><FONT color=#000000 size=2>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV><BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; This is the orbital moment with SO only, and it is 1.25294 uB, direction is down.</DIV>
<DIV></FONT></FONT>&nbsp;</DIV>
<DIV><FONT color=#ff0000><FONT color=#000000 size=2>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT color=#ff0066 size=3><STRONG>Here, I am confused: </STRONG></FONT></FONT></FONT></DIV>
<DIV><FONT color=#ff0000><FONT color=#000000><STRONG><FONT color=#ff0066 size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;1&gt; The orbital moment is still too small, it should be around 5 uB. I don't know why ?</FONT></STRONG></FONT></FONT></DIV>
<DIV><FONT color=#ff0000><FONT color=#000000><STRONG><FONT color=#ff0066 size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;2&gt; In case.inst,&nbsp;the initial&nbsp;spin is up, I think the orbital moment should be in the same</FONT></STRONG></FONT></FONT></DIV>
<DIV><FONT color=#ff0000><FONT color=#000000><STRONG><FONT color=#ff0066 size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; direction as spin. Or </FONT></STRONG></FONT></FONT><FONT color=#ff0000><FONT color=#000000><STRONG><FONT color=#ff0066 size=3>how can the total moment be around 10 uB ??</FONT></STRONG></DIV></FONT></FONT>
<DIV><STRONG><FONT color=#ff0000 size=3></FONT></STRONG>&nbsp;</DIV>
<DIV><STRONG><FONT color=#ff0000 size=3></FONT></STRONG>&nbsp;</DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<STRONG><FONT color=#ff00ff>(3)</FONT></STRONG>based on the GGA+SO scf, then i added the U parameter.</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; case.inorb and case.indm were prepared as follows:</FONT></DIV>
<DIV><FONT size=3></FONT>&nbsp;</DIV>
<DIV><FONT size=3>
<DIV><FONT color=#0000ff size=2><STRONG>case.indm</STRONG></FONT></DIV>
<DIV><FONT size=2>*************************************************************</FONT></DIV>
<DIV><FONT size=2>-9.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Emin cutoff energy<BR>&nbsp;1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; number of atoms for which density matrix is calculated<BR>&nbsp;1&nbsp; 1&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index of 1st atom, number of L's, L1<BR>&nbsp;0 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; r-index, (l,s)index</FONT></DIV>
<DIV><FONT size=2>*************************************************************</FONT></DIV>
<DIV><FONT color=#0000ff size=2><STRONG>case.inorb&nbsp; </STRONG></FONT></DIV>
<DIV><FONT size=2>*************************************************************</FONT></DIV>
<DIV><FONT size=2>&nbsp; 1&nbsp; 1&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nmod, natorb, ipr<BR>PRATT&nbsp; 1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BROYD/PRATT, mixing<BR>&nbsp; 1 1 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb<BR>&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nsic 0..AFM, 1..SIC, 2..HFM<BR>&nbsp;&nbsp;<FONT color=#ff00ff>0.50 0.00</FONT>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J (Ry)&nbsp;&nbsp; Note: we recommend to use U_eff = U-J and J=0</FONT></DIV>
<DIV><FONT size=2>*************************************************************</FONT></DIV>
<DIV><FONT size=2></FONT>&nbsp;</DIV></FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; run a regular GGA+SO scf, commamd line: <FONT color=#0000ff><STRONG>runsp -ec 0.0001 -cc 0.0001 -so -orb&nbsp;-i 1000</STRONG></FONT> 
<DIV><FONT color=#ff00ff size=3><FONT color=#000000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; then i got some results as follows:</FONT></FONT></DIV>
<DIV><FONT color=#ff00ff size=3><FONT color=#000000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :ENE:&nbsp;&nbsp;&nbsp; -46876.47444 Ry</FONT></FONT></DIV>
<DIV><FONT color=#ff00ff size=3><FONT color=#000000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;:MMTOT:&nbsp;&nbsp;13.45593&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; uB</FONT></FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :MMI001:&nbsp;6.25394&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; uB</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;:MMINT:&nbsp;&nbsp;0.94806&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; uB</FONT></DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <FONT color=#ff0000><STRONG>:ORB001:&nbsp;0.00439&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; uB </STRONG></FONT></DIV>
<DIV><STRONG><FONT color=#ff0000></FONT></STRONG>&nbsp;</DIV>
<DIV><STRONG><FONT color=#ff0000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </FONT><FONT color=#ff0066>After the three steps above, I found orbital moment in case.scf which is 0.00439 uB. comparing with the "&nbsp;1.25294 uB" generated last step,&nbsp;why it's&nbsp;so small after plus U( Ueff =&nbsp;0.5 Ry) ??</FONT></STRONG></DIV>
<DIV>&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Then, I changed the last line of case.indmc from '0 0'&nbsp;to '1 3' after convergence, </DIV>
<DIV><FONT color=#000000><FONT color=#0000ff size=2><STRONG>case.indmc</STRONG></FONT></DIV>
<DIV>
<DIV><FONT size=2>*************************************************************</FONT></DIV>
<DIV><FONT size=2>-9.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Emin cutoff energy<BR>&nbsp;1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; number of atoms for which density matrix is calculated<BR>&nbsp;1&nbsp; 1&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index of 1st atom, number of L's, L1<BR>&nbsp;<FONT color=#ff0066>1 3</FONT>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; r-index, (l,s)index</FONT></DIV>
<DIV><FONT size=2>*************************************************************</FONT></DIV>
<DIV><FONT size=2></FONT>&nbsp;</DIV></FONT></DIV>
<DIV><FONT color=#000000>&nbsp;&nbsp;&nbsp; then I run <FONT color=#0000ff><STRONG>'x lapwdm -c -so -up'</STRONG></FONT>,&nbsp;and&nbsp;I&nbsp;find the orbital moment in case.scfdmup is as follows:</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><STRONG><FONT color=#0000ff size=2>case.scfdmup</FONT></STRONG></DIV>
<DIV><FONT size=2>*************************************************************</FONT></DIV>
<DIV><FONT size=2>&nbsp;Spin-polarized + s-o calculation, M||&nbsp; 0.000&nbsp; 0.000&nbsp; 1.000<BR>&nbsp; Calculation of &lt;X&gt;, X=c*Xr(r)*Xls(l,s)<BR>&nbsp; Xr(r)&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I<BR>&nbsp; Xls(l,s) = L(dzeta)<BR>&nbsp; c=&nbsp; 1.00000<BR>&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>:XOP&nbsp; 1&nbsp; 3&nbsp;&nbsp;&nbsp; -0.00187&nbsp;&nbsp;&nbsp;&nbsp; 0.00626&nbsp;&nbsp;&nbsp;&nbsp; 0.00439<BR>*************************************************************</FONT></DIV>
<DIV><FONT size=2></FONT>&nbsp;</DIV>
<DIV><FONT size=2>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</FONT><FONT color=#ff0033 size=4><STRONG> Orbital moment is 0.00439 uB, I was completely confused ??????</STRONG></FONT></DIV>
<DIV><STRONG><FONT color=#ff0033 size=4></FONT></STRONG>&nbsp;</DIV>
<DIV><FONT color=#ff0033 size=4><STRONG></STRONG></FONT>&nbsp;</DIV>
<DIV><FONT color=#ff0033 size=4><STRONG>&nbsp;&nbsp;&nbsp;</STRONG><FONT color=#000000 size=3>The following files are&nbsp;the dmatup/dmatdn files: ( I didn't get it)</FONT></FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV>case.dmatup</DIV>
<DIV><FONT size=2>*************************************************************</FONT></DIV>
<DIV><FONT size=2>&nbsp;&nbsp;&nbsp; 1 atom density matrix<BR>&nbsp;&nbsp;&nbsp; 3&nbsp; 0.000000&nbsp; 0.000000 -0.001873 L, Lx,Ly,Lz in global orthogonal system<BR>&nbsp; 0.99122776E+00 -0.11754944E-36&nbsp;&nbsp;&nbsp; 0.15671123E-12 -0.12705494E-20<BR>&nbsp;-0.33034285E-19&nbsp; 0.84703295E-21&nbsp;&nbsp;&nbsp; 0.12447650E-12&nbsp; 0.13764285E-20<BR>&nbsp;-0.29646153E-19&nbsp; 0.00000000E+00&nbsp;&nbsp; -0.58634492E-12&nbsp; 0.84703295E-21<BR>&nbsp;-0.51769973E-03&nbsp; 0.19058241E-20<BR>&nbsp; 0.15671123E-12&nbsp; 0.12705494E-20&nbsp;&nbsp;&nbsp; 0.98989657E+00&nbsp; 0.42317797E-36<BR>&nbsp;-0.16740146E-12 -0.16940659E-20&nbsp;&nbsp; -0.72844833E-19&nbsp; 0.00000000E+00<BR>&nbsp; 0.78433837E-12&nbsp; 0.42351647E-21&nbsp;&nbsp; -0.12705494E-19&nbsp; 0.00000000E+00<BR>&nbsp; 0.51044993E-12 -0.25410988E-20<BR>&nbsp;-0.33034285E-19 -0.84703295E-21&nbsp;&nbsp; -0.16740146E-12&nbsp; 0.16940659E-20<BR>&nbsp; 0.98949500E+00&nbsp; 0.56423729E-36&nbsp;&nbsp; -0.11551635E-13&nbsp; 0.21175824E-21<BR>&nbsp; 0.51669010E-19&nbsp; 0.33881318E-20&nbsp;&nbsp; -0.97451913E-12 -0.21175824E-20<BR>&nbsp;-0.38963516E-19 -0.33881318E-20<BR>&nbsp; 0.12447650E-12 -0.13764285E-20&nbsp;&nbsp; -0.72844833E-19&nbsp; 0.00000000E+00<BR>&nbsp;-0.11551635E-13 -0.21175824E-21&nbsp;&nbsp;&nbsp; 0.98721653E+00 -0.29779190E-36<BR>&nbsp; 0.68056484E-14&nbsp; 0.21175824E-21&nbsp;&nbsp; -0.78668185E-19&nbsp; 0.33881318E-20<BR>&nbsp;-0.10434050E-12&nbsp; 0.21175824E-21<BR>&nbsp;-0.29646153E-19&nbsp; 0.00000000E+00&nbsp;&nbsp;&nbsp; 0.78433837E-12 -0.42351647E-21<BR>&nbsp; 0.51669010E-19 -0.33881318E-20&nbsp;&nbsp;&nbsp; 0.68056484E-14 -0.21175824E-21<BR>&nbsp; 0.98801066E+00&nbsp; 0.31346516E-36&nbsp;&nbsp;&nbsp; 0.18533121E-12 -0.42351647E-21<BR>&nbsp;-0.39069395E-19&nbsp; 0.16940659E-20<BR>&nbsp;-0.58634492E-12 -0.84703295E-21&nbsp;&nbsp; -0.12705494E-19&nbsp; 0.00000000E+00<BR>&nbsp;-0.97451913E-12&nbsp; 0.21175824E-20&nbsp;&nbsp; -0.78668185E-19 -0.33881318E-20<BR>&nbsp; 0.18533121E-12&nbsp; 0.42351647E-21&nbsp;&nbsp;&nbsp; 0.98954264E+00 -0.20375235E-36<BR>&nbsp;-0.15970123E-12 -0.21175824E-21<BR>&nbsp;-0.51769973E-03 -0.19058241E-20&nbsp;&nbsp;&nbsp; 0.51044993E-12&nbsp; 0.25410988E-20<BR>&nbsp;-0.38963516E-19&nbsp; 0.33881318E-20&nbsp;&nbsp; -0.10434050E-12 -0.21175824E-21<BR>&nbsp;-0.39069395E-19 -0.16940659E-20&nbsp;&nbsp; -0.15970123E-12&nbsp; 0.21175824E-21<BR>&nbsp; 0.99133429E+00&nbsp; 0.10971281E-36</FONT><BR><FONT size=2>*************************************************************</FONT></DIV>
<DIV><FONT size=2></FONT>&nbsp;</DIV>
<DIV><FONT size=2>case.dmatdn</FONT></DIV>
<DIV><FONT size=2>*************************************************************</FONT></DIV>
<DIV>&nbsp;<FONT size=2> 1 atom density matrix<BR>&nbsp;&nbsp;&nbsp; 3&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.006263 L, Lx,Ly,Lz in global orthogonal system<BR>&nbsp; 0.76999012E-03&nbsp; 0.14497764E-36&nbsp;&nbsp; -0.82982858E-13&nbsp; 0.44998625E-21<BR>&nbsp;-0.72844833E-19&nbsp; 0.35998900E-20&nbsp;&nbsp; -0.86064816E-13 -0.16940659E-20<BR>&nbsp;-0.30969642E-19 -0.50821977E-20&nbsp;&nbsp;&nbsp; 0.19253073E-11&nbsp; 0.19905274E-19<BR>&nbsp;-0.23424275E-03 -0.13552527E-19<BR>&nbsp;-0.82982858E-13 -0.44998625E-21&nbsp;&nbsp;&nbsp; 0.16025310E-02&nbsp; 0.38399482E-36<BR>&nbsp; 0.21974289E-12&nbsp; 0.84703295E-21&nbsp;&nbsp;&nbsp; 0.18973538E-18&nbsp; 0.00000000E+00<BR>&nbsp; 0.48535084E-12 -0.16940659E-20&nbsp;&nbsp; -0.27952087E-19 -0.25410988E-20<BR>&nbsp;-0.18984536E-11 -0.15246593E-19<BR>&nbsp;-0.72844833E-19 -0.35998900E-20&nbsp;&nbsp;&nbsp; 0.21974289E-12 -0.84703295E-21<BR>&nbsp; 0.13365436E-02 -0.13322269E-36&nbsp;&nbsp;&nbsp; 0.21297389E-12&nbsp; 0.14823077E-20<BR>&nbsp;-0.75385932E-19&nbsp; 0.33881318E-20&nbsp;&nbsp; -0.17868594E-12&nbsp; 0.78085850E-21<BR>&nbsp;-0.32187252E-19 -0.16940659E-20<BR>&nbsp;-0.86064816E-13&nbsp; 0.16940659E-20&nbsp;&nbsp;&nbsp; 0.18973538E-18&nbsp; 0.00000000E+00<BR>&nbsp; 0.21297389E-12 -0.14823077E-20&nbsp;&nbsp;&nbsp; 0.98011221E+00&nbsp; 0.51721751E-36<BR>&nbsp;-0.20662598E-12&nbsp; 0.00000000E+00&nbsp;&nbsp;&nbsp; 0.19989978E-18&nbsp; 0.13552527E-19<BR>&nbsp; 0.23065083E-13 -0.16940659E-20<BR>&nbsp;-0.30969642E-19&nbsp; 0.50821977E-20&nbsp;&nbsp;&nbsp; 0.48535084E-12&nbsp; 0.16940659E-20<BR>&nbsp;-0.75385932E-19 -0.33881318E-20&nbsp;&nbsp; -0.20662598E-12&nbsp; 0.00000000E+00<BR>&nbsp; 0.20209675E-02 -0.76798964E-36&nbsp;&nbsp; -0.22783686E-12 -0.15881868E-20<BR>&nbsp;-0.92008954E-19 -0.16940659E-20<BR>&nbsp; 0.19253073E-11 -0.19905274E-19&nbsp;&nbsp; -0.27952087E-19&nbsp; 0.25410988E-20<BR>&nbsp;-0.17868594E-12 -0.78085850E-21&nbsp;&nbsp;&nbsp; 0.19989978E-18 -0.13552527E-19<BR>&nbsp;-0.22783686E-12&nbsp; 0.15881868E-20&nbsp;&nbsp;&nbsp; 0.26211585E-02 -0.65827684E-36<BR>&nbsp;-0.67107735E-13&nbsp; 0.21175824E-20<BR>&nbsp;-0.23424275E-03&nbsp; 0.13552527E-19&nbsp;&nbsp; -0.18984536E-11&nbsp; 0.15246593E-19<BR>&nbsp;-0.32187252E-19&nbsp; 0.16940659E-20&nbsp;&nbsp;&nbsp; 0.23065083E-13&nbsp; 0.16940659E-20<BR>&nbsp;-0.92008954E-19&nbsp; 0.16940659E-20&nbsp;&nbsp; -0.67107735E-13 -0.21175824E-20<BR>&nbsp; 0.19503738E-02 -0.34481168E-36</FONT></DIV>
<DIV><FONT size=2>*************************************************************</FONT></DIV>
<DIV><FONT size=2></FONT>&nbsp;</DIV>
<DIV><STRONG><FONT color=#ff0066>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Dr. Stefaan, above are the dmatup/dn files of this&nbsp;case. Accutually, I don't get this dmatup/dn files. </FONT></STRONG></DIV>
<DIV><STRONG><FONT color=#ff0066>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;1&gt;Please tell me what these numbers mean? </FONT></STRONG></DIV>
<DIV><STRONG><FONT color=#ff0066>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;2&gt;And how can I write an occupation in the dmat files that gives the occupation&nbsp;I want ? </FONT></STRONG></DIV>
<DIV><STRONG><FONT color=#ff0066>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &lt;3&gt;After that, what should I do next to get an correct total magnetic moment?</FONT></STRONG></DIV>
<DIV><STRONG><FONT color=#ff0066></FONT></STRONG>&nbsp;</DIV>
<DIV><FONT color=#ff0066><STRONG>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</STRONG></FONT><FONT color=#000000>Aagin, thank you so much for your kindly and detailed reply, Dr.Stefaan, I am looking forward to your coming enlightenment.<IMG src="http://mail.cstnet.cn/funcs/images/smile01.gif"></FONT></DIV>
<DIV><FONT color=#000000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Cheers.</FONT></DIV>
<DIV><FONT color=#000000></FONT>&nbsp;</DIV>
<DIV><FONT color=#000000>Yours sincerely:</FONT></DIV>
<DIV><FONT color=#000000>Hui Wang</FONT></DIV></FONT></DIV></DIV>
<DIV>&nbsp;</DIV> <br><!-- urlfiles --><br><br><!-- footer --><br>