<P><SPAN lang=EN-US style="FONT-SIZE: 10.5pt; FONT-FAMILY: 'Times New Roman'; mso-fareast-font-family: 宋体; mso-font-kerning: 1.0pt; mso-ansi-language: EN-US; mso-fareast-language: ZH-CN; mso-bidi-language: AR-SA">Dear prof. Blaha:</SPAN></P>

<P><SPAN lang=EN-US style="FONT-SIZE: 10.5pt; FONT-FAMILY: 'Times New Roman'; mso-fareast-font-family: 宋体; mso-font-kerning: 1.0pt; mso-ansi-language: EN-US; mso-fareast-language: ZH-CN; mso-bidi-language: AR-SA">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN lang=EN-US style="FONT-SIZE: 10.5pt; FONT-FAMILY: 'Times New Roman'; mso-fareast-font-family: 宋体; mso-font-kerning: 1.0pt; mso-ansi-language: EN-US; mso-fareast-language: ZH-CN; mso-bidi-language: AR-SA"><SPAN style="mso-spacerun: yes">&nbsp;</SPAN>I'm sorry for the wrong struct.in fact,&nbsp;I'm dealing with a &nbsp;2*2*2 Ni3Al supercell. I just want <SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">to find which site will be substituted by Re(Ni site or Al site). All goes well when i use the usual command run_lapw, except *.scf1 gives such information:</SPAN></SPAN></SPAN></P><SPAN lang=EN-US style="FONT-SIZE: 10.5pt; FONT-FAMILY: 'Times New Roman'; mso-fareast-font-family: 宋体; mso-font-kerning: 1.0pt;

 mso-ansi-language: EN-US; mso-fareast-language: ZH-CN; mso-bidi-language: AR-SA"><SPAN lang=EN-US style="FONT-SIZE: 10.5pt; FONT-FAMILY: 'Times New Roman'; mso-fareast-font-family: 宋体; mso-font-kerning: 1.0pt; mso-ansi-language: EN-US; mso-fareast-language: ZH-CN; mso-bidi-language: AR-SA"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt"><?xml:namespace prefix = o ns = "urn:schemas-microsoft-com:office:office" /><o:p>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">ATOMIC SPHERE DEPENDENT PARAMETERS FOR ATOM&nbsp; Re1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; OVERALL ENERGY PARAMETER IS&nbsp;&nbsp;&nbsp; 0.3000<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; OVERALL BASIS SET ON ATOM IS LAPW<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; E( 1)=&nbsp;&nbsp;&nbsp; 0.3000<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LAPW<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; E( 1)=&nbsp;&nbsp; -2.3950&nbsp;&nbsp; E(BOTTOM)=&nbsp;&nbsp; -2.750&nbsp;&nbsp; E(TOP)=&nbsp;&nbsp; -2.040<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LOCAL ORBITAL<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; E( 3)=&nbsp;&nbsp;&nbsp;

 0.3000<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LAPW<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; E( 3)=&nbsp;&nbsp; -1.9550&nbsp;&nbsp; E(BOTTOM)=&nbsp;&nbsp; -1.970&nbsp;&nbsp; E(TOP)=&nbsp;&nbsp; -1.940<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LOCAL ORBITAL<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; E( 2)=&nbsp;&nbsp;&nbsp; 0.3000&nbsp;&nbsp; E(BOTTOM)=&nbsp;&nbsp; -0.070&nbsp;&nbsp; E(TOP)= -200.000<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LAPW<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; E( 0)=&nbsp;&nbsp;&nbsp; 0.3000<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LAPW<o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">……<o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">so i think maybe there is something wrong with the linear energy of Re. The following&nbsp;is *.in1 and *.scf2:<o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">*.in1<o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">WFFIL&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (WFPRI, SUPWF) &nbsp; 7.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10&nbsp;&nbsp;&nbsp; 4 (R-MT*K-MAX; MAX L IN WF, V-NMT &nbsp; 0.30&nbsp;&nbsp;&nbsp; 6&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW) &nbsp;1&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 CONT 0 &nbsp;1&nbsp;&nbsp; -2.92&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.005 STOP 0 &nbsp;3&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 CONT 0 &nbsp;3&nbsp;&nbsp; -2.68&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.010 CONT 0 &nbsp;2&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.010 CONT 0 &nbsp;0&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 CONT 0 &nbsp; 0.30&nbsp;&nbsp;&nbsp; 4&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (GLOBAL E-PARAMETER WITH n OTHER CHOICES, glo!

bal

 APW/LAPW)<BR>&nbsp;1&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 CONT 0<BR>&nbsp;1&nbsp;&nbsp; -4.96&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.005 STOP 0<BR>&nbsp;2&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.010 CONT 0<BR>&nbsp;0&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 CONT 0<BR>&nbsp; 0.30&nbsp;&nbsp;&nbsp; 4&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<BR>&nbsp;1&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 CONT 0<BR>&nbsp;1&nbsp;&nbsp; -4.96&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.005 STOP 0<BR>&nbsp;2&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.010 CONT 0<BR>&nbsp;0&nbsp;&nbsp;&nbsp; 0.30&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 CONT 0<o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">……<o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">*.scf2<o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">:GMA&nbsp; : POTENTIAL AND CHARGE CUT-OFF&nbsp; 14.00 Ry**.5<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Energy to separate semicore and valencestates:&nbsp; -0.09763<o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt"><BR>:NOE&nbsp; : NUMBER OF ELECTRONS&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 467.000<o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">:FER&nbsp; : F E R M I - ENERGY(TETRAH.M.)=&nbsp;&nbsp; 0.75214<o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">&nbsp;<o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">Prof Blaha, as you suggusted, I should try the command run_lapw -in1new2, right?<o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt; LINE-HEIGHT: 130%"><SPAN lang=EN-US style="mso-bidi-font-size: 10.5pt">&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;&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;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <o:p></o:p></SPAN></P>

<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US>&nbsp;<o:p></o:p></SPAN></P>

<P></o:p></SPAN>&nbsp;</P></SPAN></SPAN><p><br><hr size=1><b>Do You Yahoo!?</b><br>

<a href="http://cn.rd.yahoo.com/mail_cn/tag/?http://cn.promo.yahoo.com/minisite/messenger1/" target=blank>暑期大片齐聚雅虎通 网络摄像头+雅虎通调频收音机等你来拿</a>