<DIV>
<DIV><FONT size=3>Dear Wien users:<BR>&nbsp;&nbsp;&nbsp; I am running wien version wien2k 08_03 on a machine of type Xeon with operating system suse linux, fortran compiler ifort 10.1.021 and math libraries mkl 10.0.2.018.</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;I have been trying to&nbsp;simulate hcp Tb for weeks, no matter&nbsp;whatever I did (including change the parameters&nbsp;and use different method such as GGA+U or SO), <FONT color=#ff0066><STRONG>I still couldn't get the right magnetic moment of Tb atom, and the convergence criterion was met.</STRONG></FONT> The experimental magnetic moment of Tb is around 10 uB, but the largest magnetic moment I got by wien2k_08.3 was around 7.6 uB (GGA+U+SO).</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; My computational details are as follows:</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; here are all the import input files are prepared before runsp_lapw:</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; <FONT color=#0000ff><STRONG>case.struct</STRONG></FONT></FONT></DIV>
<DIV><FONT color=#000000 size=3>*************************************************************</FONT></DIV>
<DIV><FONT color=#0000ff size=3><FONT color=#000000 size=2>HCP_Tb<BR>H&nbsp;&nbsp; LATTICE,NONEQUIV.ATOMS:&nbsp; 1 194 P63/mmc<BR>MODE OF CALC=RELA unit=bohr<BR>&nbsp; 6.803402&nbsp; 6.803402 10.767485 90.000000 90.000000120.000000<BR>ATOM&nbsp; -1: X=0.33333333 Y=0.66666667 Z=0.75000000<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; MULT= 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ISPLIT= 4<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -1: X=0.66666667 Y=0.33333333 Z=0.25000000<BR>Tb1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NPT=&nbsp; 781&nbsp; R0=0.00010000 RMT=&nbsp;&nbsp; 2.50000&nbsp;&nbsp; Z: 65.0<BR>LOCAL ROT MATRIX:&nbsp;&nbsp;&nbsp; 1.0000000 0.0000000 0.0000000<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000000 1.0000000 0.0000000<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000000 0.0000000 1.0000000<BR>&nbsp; 24&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NUMBER OF SYMMETRY OPERATIONS</FONT><BR><FONT color=#000000>*************************************************************</FONT></FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; <FONT color=#0000ff><STRONG>case.in1</STRONG></FONT></FONT></DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV>WFFIL&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (WFPRI, SUPWF)<BR>&nbsp; 7.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 12&nbsp;&nbsp;&nbsp; 4 (R-MT*K-MAX; MAX L IN WF, V-NMT<BR>&nbsp; 0.30&nbsp;&nbsp;&nbsp; 6&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR>&nbsp;0&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 CONT 1<BR>&nbsp;0&nbsp;&nbsp; -3.53&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.005 STOP 1<BR>&nbsp;1&nbsp;&nbsp; -1.70&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.010 CONT 1<BR>&nbsp;1&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 CONT 1<BR>&nbsp;3&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.010 CONT 1<BR>&nbsp;2&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.010 CONT 1<BR>K-VECTORS FROM UNIT:4&nbsp;&nbsp; -15.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5.0&nbsp;&nbsp;&nbsp; 48&nbsp;&nbsp; emin/emax/nband</DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; <FONT color=#0000ff><STRONG>case.in2</STRONG></FONT></FONT></DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV><FONT size=2>TOT&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (TOT,FOR,QTL,EFG,FERMI)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -9.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 38.0 0.50 0.05&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; EMIN, NE, ESEPERMIN, ESEPER0<BR>TETRA&nbsp;&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (GAUSS,ROOT,TEMP,TETRA,ALL&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; eval)<BR>&nbsp; 0 0&nbsp; 2 0&nbsp; 3 3&nbsp; 4 0&nbsp; 5 3&nbsp; 6 0&nbsp; 6 6<BR>&nbsp;12.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; GMAX<BR>NOFILE&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; FILE/NOFILE&nbsp; write recprlist</FONT></DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; <STRONG><FONT color=#0000ff>case.inm</FONT></STRONG></FONT></DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV>MSEC1&nbsp; 0.0&nbsp;&nbsp; YES&nbsp; (BROYD/PRATT, extra charge (+1 for additional e), norm)<BR>0.02&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing FACTOR for BROYD/PRATT scheme<BR>1.00&nbsp; 1.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; PW and CLM-scaling factors<BR>9999&nbsp; 8&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; idum, HISTORY</DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; <STRONG><FONT color=#0000ff>case.inst</FONT></STRONG></FONT></DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV>Tb<BR>Xe 4<BR>4, 3,3.0&nbsp; N<BR>4, 3,1.0&nbsp; N<BR>4,-4,4.0&nbsp; N<BR>4,-4,0.0&nbsp; N<BR>5, 2,1.0&nbsp; N<BR>5, 2,0.0&nbsp; N<BR>6,-1,1.0&nbsp; N<BR>6,-1,1.0&nbsp; N<BR>****<BR>****&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; END of input (instgen_lapw)</DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; <STRONG><FONT color=#0000ff>case.indm</FONT></STRONG></FONT></DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV>-15.&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</DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; <FONT color=#0000ff><STRONG>case.indmc</STRONG></FONT></FONT></DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV>-15.&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</DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;<FONT color=#0000ff><STRONG>case.inorb&nbsp; <FONT color=#cc0099>(here, the value of u varies from 0.0 to 0.7 Ry)</FONT></STRONG></FONT></FONT></DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV>&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; 0.10 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J (Ry)&nbsp;&nbsp; Note: we recommend to use U_eff = U-J and J=0</DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; <STRONG><FONT color=#0000ff>case.inso</FONT></STRONG></FONT></DIV>
<DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV>WFFIL<BR>&nbsp;4&nbsp; 1&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; llmax,ipr,kpot<BR>&nbsp;-15.0000&nbsp;&nbsp; 5.0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; emin,emax (output energy window)<BR>&nbsp;&nbsp; 0.&nbsp; 0.&nbsp; 1.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; direction of magnetization (lattice vectors)<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; number of atoms for which RLO is added<BR>&nbsp;1&nbsp;&nbsp; -1.7&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.01&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; atom number,e-lo,de (case.in1), repeat NX times<BR>&nbsp;0 0 0 0 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; number of atoms for which SO is switch off; atoms</DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV></DIV>
<DIV><FONT size=3>other parameters: RMT = 2.5, Kpoints: 1000, IBZ: 76</FONT></DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV><FONT color=#ff0066 size=3><STRONG></STRONG></FONT>&nbsp;</DIV>
<DIV><FONT color=#ff0066 size=3><STRONG>Result: </STRONG></FONT></DIV>
<DIV><STRONG><FONT color=#009933 size=3>by GGA+U:</FONT></STRONG></DIV>
<DIV><FONT color=#0000cc size=3>command line: runsp -ec 0.0001 -cc 0.0001 -orb -i 1000</FONT></DIV>
<DIV>&nbsp;U(Ry)&nbsp;&nbsp;&nbsp;&nbsp;E(Ry/cell)&nbsp;&nbsp;&nbsp;&nbsp; M(uB/atom)<BR>&nbsp; 0.0&nbsp;&nbsp;&nbsp; -46876.362489&nbsp;&nbsp;&nbsp; 5.86499<BR>&nbsp; 0.1&nbsp;&nbsp;&nbsp; -46876.362491&nbsp;&nbsp;&nbsp; 5.86434 ( whe Ueff = 0.1 Ry, using -orb doesn't converge, so i&nbsp;use -orbc)<BR>&nbsp; 0.2&nbsp;&nbsp;&nbsp; -46876.263915&nbsp;&nbsp;&nbsp; 6.20827<BR>&nbsp; 0.3&nbsp;&nbsp;&nbsp; -46876.270606&nbsp;&nbsp;&nbsp; 6.25566<BR>&nbsp; 0.4&nbsp;&nbsp;&nbsp; -46876.260620&nbsp;&nbsp;&nbsp; 6.26501<BR>&nbsp; 0.5&nbsp;&nbsp;&nbsp; -46876.232366&nbsp;&nbsp;&nbsp; 6.26697<BR>&nbsp; 0.6&nbsp;&nbsp;&nbsp; -46876.243298&nbsp;&nbsp;&nbsp; 6.27250<BR>&nbsp; 0.7&nbsp;&nbsp;&nbsp; -46876.216162&nbsp;&nbsp;&nbsp; 6.27217</DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT color=#009933 size=3><STRONG>by GGA+U+SO</STRONG> <FONT color=#cc0099>(based on the result of GGA+U, and rm *.broy*)</FONT></FONT><BR><FONT color=#0000cc size=3>command line: runsp -ec 0.0001 -cc 0.0001 -orb -so -i 1000</FONT></DIV>
<DIV>U(Ry)&nbsp;&nbsp;&nbsp; E(Ry/cell)&nbsp;&nbsp;&nbsp;M_spin(uB/atom)&nbsp; M_orbit(uB/atom)<BR>0.0&nbsp;&nbsp;&nbsp;&nbsp; -46876.609319&nbsp;&nbsp;&nbsp;&nbsp; 5.81268&nbsp;&nbsp;&nbsp;&nbsp; 1.25824<BR>0.1&nbsp;&nbsp;&nbsp;&nbsp; -46876.335016&nbsp;&nbsp;&nbsp;&nbsp; 6.04722&nbsp;&nbsp;&nbsp;&nbsp; 0.00000 (GGA+U+SO doesn't converge, i don't know why ?)<BR>0.2&nbsp;&nbsp;&nbsp;&nbsp; -46876.507404&nbsp;&nbsp;&nbsp;&nbsp; 6.19715&nbsp;&nbsp;&nbsp;&nbsp; 1.40041<BR>0.3&nbsp;&nbsp;&nbsp;&nbsp; -46876.498692&nbsp;&nbsp;&nbsp;&nbsp; 6.23715&nbsp;&nbsp;&nbsp;&nbsp; 0.01548<BR>0.4&nbsp;&nbsp;&nbsp;&nbsp; -46876.487014&nbsp;&nbsp;&nbsp;&nbsp; 6.25003&nbsp;&nbsp;&nbsp;&nbsp; 0.00605<BR>0.5&nbsp;&nbsp;&nbsp;&nbsp; -46876.470927&nbsp;&nbsp;&nbsp;&nbsp; 6.25129&nbsp;&nbsp;&nbsp;&nbsp; 1.33299<BR>0.6&nbsp;&nbsp;&nbsp;&nbsp; -46876.467519&nbsp;&nbsp;&nbsp;&nbsp; 6.26165&nbsp;&nbsp;&nbsp;&nbsp; 0.00113<BR>0.7&nbsp;&nbsp;&nbsp;&nbsp; -46876.452793&nbsp;&nbsp;&nbsp;&nbsp; 6.25998&nbsp;&nbsp;&nbsp;&nbsp; 1.36656</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; </DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; <FONT color=#ff0000 size=4><STRONG>Two questions:</STRONG></FONT></DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; <FONT size=3>(1)From the data above, the largest magnetic moment of atom Tb is about 7.6 uB, which is still small than the experimental magnetic moment (10 uB). How can I get the right magnetic moment of Tb?</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; (2)About GGA+U, it need case.indm or case.indmc ? runsp -orb or runsp -orbc ?</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; About GGA+U+SO, there are <FONT color=#ff00ff><STRONG>two Tb atoms which are at the unequivalent position</STRONG></FONT>, how to change the case.inso files ?</FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<STRONG><FONT color=#0000ff>case.inso</FONT></STRONG>
<DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV><FONT size=2>WFFIL<BR>&nbsp;4&nbsp; 1&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; llmax,ipr,kpot<BR>&nbsp;-15.0000&nbsp;&nbsp; 5.0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; emin,emax (output energy window)<BR>&nbsp;&nbsp; 0.&nbsp; 0.&nbsp; 1.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; direction of magnetization (lattice vectors)<BR>&nbsp;</FONT><FONT size=2><FONT color=#ff00ff>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;number of atoms for which RLO is added<BR>&nbsp;1&nbsp;&nbsp; -1.7&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.01&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; atom number,e-lo,de (case.in1), repeat NX times </FONT><FONT color=#0000ff>( if i wanna add&nbsp;RLO to the second Tb atom, what should I do ?)<BR></FONT>&nbsp;0 0 0 0 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; number of atoms for which SO is switch off; atoms</FONT></DIV>
<DIV><FONT size=3>*************************************************************</FONT></DIV>
<DIV>&nbsp;</DIV></DIV></FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; I know havey rare earth elements are hard to deal with , and I really need someone who can give me some enlightenment. Any suggestion will be greatly appreciated.</FONT></DIV><FONT size=3></FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp;&nbsp;I am looking forward to your reply.<IMG src="http://mail.cstnet.cn/funcs/images/smile01.gif"></FONT></DIV>
<DIV><FONT size=3>&nbsp;&nbsp;&nbsp; Cheers.</FONT></DIV>
<DIV><FONT size=3></FONT>&nbsp;</DIV>
<DIV><FONT size=3>Yours sincerely</FONT></DIV>
<DIV><FONT size=3>Hui Wang</DIV>
<DIV><BR></DIV></FONT>
<DIV><FONT color=#0000cc size=3></FONT>&nbsp;</DIV>
<DIV><STRONG><FONT color=#009933 size=3></FONT></STRONG>&nbsp;</DIV>
<DIV><FONT size=3></FONT>&nbsp;</DIV><br><!-- sign --><br>=========================================================<br>Magnetism&nbsp;and&nbsp;Magnetic&nbsp;Materials&nbsp;Division<br>Shenyang&nbsp;Materials&nbsp;Science&nbsp;National&nbsp;Laboratory<br>Institute&nbsp;of&nbsp;Metal&nbsp;Research<br>Chinese&nbsp;Academy&nbsp;of&nbsp;Sciences<br>72&nbsp;Wenhua&nbsp;Road,&nbsp;Shenyang&nbsp;110016,&nbsp;P.&nbsp;R.&nbsp;China<br>Tel:&nbsp;+86-24-83978845<br>PHD:&nbsp;Wanghui<br>Email:&nbsp;hwang@imr.ac.cn&nbsp;<br>========================================================= <br><!-- urlfiles --><br><br><!-- footer --><br>