What is your number of :ITE? You would be more patient for AFM calculations. I suggest to run:<br>runsp_lapw -cc 0.00001 -i 500 ...<br>Sometimes after 200 iterations the MMTO starts to approach to zero. In this case, MMTOT~0, the MMI(Cr1) = -MMI(Cr2).<br><br><b><i>明文美 <iphyboy@hotmail.com></i></b> wrote:<blockquote class="replbq" style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"> <style> .hmmessage P { margin:0px; padding:0px } body.hmmessage { FONT-SIZE: 9pt; FONT-FAMILY:Tahoma } </style> Dear Wien2k users,<br> I have a very confusing a problem. I want to do a AFM calculation of LaOCrAs, and also I inverse the spin up and spin down. <br> in LaOCrAs.inst. But the converged result shows that the magnetic moments in the two Cr atoms are spin polarized in the parallel direction,<br> rather than in anti-parallel direction(see below), which is totally equivlent to a FM
calculation. The similar problems also occured when I try to calculate LaOMnAs, LaOFeAs, etc.The AFM calculations are just carried out by commanding runsp_lapw after I filps the two Cr atoms spins <br> and set magnetic moments of non-magnetic atoms to be zero. Did I make some mistakes in these AFM calculation<br> <br> Any suggestion will be greatly appreciated !<br> <br> Ming Wenmei <br> #########################################################################################<br> magnetic moments of all atom in MT spheres<br> :MMI001: MAGNETIC MOMENT IN SPHERE 1 = -0.01026<br>:MMI002: MAGNETIC MOMENT IN SPHERE 2 = -0.00031<br>:MMI003: MAGNETIC MOMENT IN SPHERE 3 = -0.00032<br>:MMI004: MAGNETIC MOMENT IN SPHERE
4 = 1.94370 ++++++ Cr1<br>:MMI005: MAGNETIC MOMENT IN SPHERE 5 = 1.95533 ++++++ Cr2<br>:MMI006: MAGNETIC MOMENT IN SPHERE 6 = -0.07087<br> ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++<br> ###################################################################################################<br> <br> Here is my Struct file and the LaOCrAs.inst file are as the followings, and Version of My Wien2k is 01_08<br> <br> ###########################################################################################<br> LaOCrAs P
6 <br> RELA <br> 7.437019 7.437019 17.662621 90.0! 00000 90.000000 90.000000 <br>ATOM -1: X=0.25000000 Y=0.25000000
Z=0.13978483<br> MULT= 2 ISPLIT= 8<br> -1: X=0.75000000 Y=0.75000000 Z=0.86021517<br>La1 NPT= 781 R0=0.00010000 RMT= 2.35 Z: 57.0 <br>LOCAL ROT MATRIX: 1.0000000 0.0000000 0.0000000<br> 0.0000000 1.0000000 0.0000000<br> &n! bsp; 0.0000000 0.0000000 1.0000000<br>ATOM -2: X=0.7500000 0 Y=0.25000000 Z=0.00000000<br> MULT=
1 ISPLIT=-2<br>O 2 NPT= 781 R0=0.00010000 RMT= 2.09 Z: 8.0 <br>LOCAL ROT MATRIX: 0.7071068-0.7071068 0.0000000<br> 0.7071068 0.7071068 0.0000000<br> 0.0000000 0.0000000 1.0000000<br>ATOM -3: X=0.25000000 Y=0.75000000 Z=0.00000000<br> MULT= 1 ISPLIT=-2<br>O 3&nb! sp; NPT= 781 R0=0.00010000
RMT= 2.09 Z: 8.0 <br>LOCAL ROT MATRIX: 0.7071068-0.7071068 0.0000000<br> 0.7071068 0.7071068 0.0000000<br> 0.0000000 0.0000000 1.0000000<br>ATOM -4: X=0.75000000 Y=0.25000000 Z=0.50000000<br> MULT= 1 ISPLIT=-2<br>Cr1 NPT= 781 R0=0.00010000 RMT= 2.46 Z: 24.0 &n!
bsp; <br>LOCAL ROT MATRIX: 0.7071068-0.7071068 0.0000000<br> 0.7071068 0.7071068 0.0000000<br> 0.0000000 0.0000000 1.0000000<br>ATOM -5: X=0.25000000 Y=0.75000000 Z=0.50000000<br> MULT= 1 ISPLIT=-2<br>Cr2 NPT= 781 R0=0.00010000 RMT= 2.46 Z: 24.0 <br>LOCAL ROT MATRIX: 0.7071068-0.7071068
0.0000000<br> &n! bsp; 0.7071068 0.7071068 0.0000000<br> 0.0000000 0.0000000 1.0000000<br>ATOM -6: X=0.25000000 Y=0.25000000 Z=0.65917688<br> MULT= 2 ISPLIT= 8<br> -6: X=0.75000000 Y=0.75000000 Z=0.34082312<br>As6 NPT= 781 R0=0.00010000 RMT= 2.18 Z: 33.0 <br>LOCAL ROT MATRIX: 1.0000000 0.0000000
0.0000000<br> 0.0000000 1.0000000 0.0000000<br> &n! bsp; 0.0000000 0.0000000 1.0000000<br> 8 NUMBER OF SYMMETRY OPERATIONS<br> 0-1 0 0.00000000<br>-1 0 0 0.00000000<br> 0 0-1 0.00000000<br> 1<br> 1 0 0 0.00000000<br> 0 1 0 0.00000000<br> 0 0 1 0.00000000<br> 2<br>-1 0 0 0.50000000<br> 0 1 0 0.00000000<br> 0 0 1 0.00000000<br> 3<br>-1 0 0 0.50000000<br> 0-1 0 0.50000000<br> 0 0 1 0.00000000<br> 4<br> 0 1 0 0.50000000<br>-1 0 0 0.00000000<br> 0 0-1 0.00000000<br>
5<br> 0-1 0 0.00000000<br> 1 0 0 0.50000000<br> 0 0-1 0.00000000<br> 6<br> 0 1 0 0.50000000<br> 1 0 0 0.50000000<br> 0 0-1 0.00000000<br> 7<br> 1 0 0 0.00000000<br> 0-1 0 0! .50000000<br> 0 0 1 0.00000000<br> 8<br> +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++<br> +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++<br> La<br>Xe 2 <br>5, 2,0.5 N<br>5, 2,0.5 N<br>6,-1,1.0 N<br>6,-1,1.0 N<br>O<br>He 3 <br>2,-1,1.0 N<br>2,-1,1.0 N<br>2, 1,1.0 N<br>2, 1,1.0 N<br>2,-2,1.0 N<br>2,-2,1.0 N<br>O<br>He 3 <br>2,-1,1.0 N<br>2,-1,1.0 N<br>2, 1,1.0 N<br>2, 1,1.0 N<br>2,-2,1.0 N<br>2,-2,1.0 N<br>Cr<br>Ar 3 <br>3, 2,2.0 N<br>3, 2,2.0
N<br>3,-3,1.0 N<br>3,-3,0.0 N<br>4,-1,1.0 N<br>4,-1,0.0 N<br>Cr<br>Ar 3 <br>3, 2,2.0 N<br>3, 2,2.0 N<br>3,-3,0.0 N<br>3,-3,1.0 N<br>4,-1,0.0 N<br>4,-1,1.0 N<br>As<br>Ar 5 <br>3, 2,2.0 N<br>3, 2,2.0 N<br>3,-3,3.0 N<br>3,-3,3.0 N<br>4,-1,1.0 N<br>4,-1,1.0 N<br>4, 1,0.5 N<br>4, 1,0.5 N<br>4,-2,1.0 N<br>4,-2,1.0 N<br>****<br>**** END of input (instgen_lapw)<br> <br> <br> <br><br><hr>轻松把Hotmail下载到本地,试试 Windows Live Mail。 <a href="http://get.live.com/wl/all" target="_new">立即尝试!</a>_______________________________________________<br>Wien mailing list<br>Wien@zeus.theochem.tuwien.ac.at<br>http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien<br></blockquote><br><p> 
<hr size=1>Be a better friend, newshound, and
know-it-all with Yahoo! Mobile. <a href="http://us.rd.yahoo.com/evt=51733/*http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ "> Try it now.</a>