<P>郭怀红你好，</P>
<P>我也是WIEN2k的用户，目前从事“百人计划”的研究，我正在做硼酸盐矿物中B原子位置p电子轨道的px, py和pz的分解，就此我想请教你：</P>
<P>1、我的矿物是三斜晶系；</P>
<P>2、px, py和pz的分解和projected时候，是否应该在PAS坐标系？是否需要在local rotation matrix和global coordinate system间转换？</P>
<P>3、如需要这样的转换，应如何编辑相应的输入文件如case.inq?</P>
<P>谢谢！</P>
<P>周兵<BR><BR><BR></P>
<BLOCKQUOTE style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #b6b6b6 2px solid; MARGIN-RIGHT: 0px" name="replyContent">-----原始邮件-----<BR><B>发件人:</B> guohuaihong &lt;hhguo@imr.ac.cn&gt;<BR><B>发送时间:</B> 2012年7月8日 星期日<BR><B>收件人:</B> "A Mailing list for WIEN2k users" &lt;wien@zeus.theochem.tuwien.ac.at&gt;<BR><B>抄送:</B> <BR><B>主题:</B> [Wien] "case.insp" to define band character for plotting<BR><BR><SPAN id="yiv1110133197yui_3_2_0_18_134154424361457">Dear Wien2K users and Prof. Blaha:<BR><BR></SPAN><SPAN lang="">
<DIV><SPAN id="yiv1110133197yui_3_2_0_18_134154424361464"></SPAN></DIV>
<DIV><SPAN id="yiv1110133197yui_3_2_0_18_134154424361467">When calculating band structure with band character for plotting, we need to change <BR>"case.insp" file accordingly. I have two questions concerning with 'jtype':<BR><BR></SPAN></DIV>
<DIV><SPAN>1.&nbsp;&nbsp; Is jtype (0,1,2,3) correspoding to orbital (s,p,d,f) for each atom speficified ?<BR><BR></SPAN></DIV>
<DIV><SPAN>2. &nbsp; How to define 'jtype' if I want to calculate band character of projected orbitals, such as Px, Py<BR>and dx2-y2 and so on? <BR><BR>In the </SPAN><SPAN><SPAN>SRC_spaghetti/spag.f, PX, PY, PZ, DXY, DYZ, DXZ, DX2Y2 and DZ2 as listed below may give<BR>some hint, but </SPAN>I am still confused with it.<BR><BR></SPAN></DIV>
<DIV><SPAN>===============================================</SPAN></DIV>
<DIV><SPAN>&nbsp;&nbsp;&nbsp; 430&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(:)=label1(:)<BR>&nbsp;&nbsp;&nbsp; 431&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index_shift=0<BR>&nbsp;&nbsp;&nbsp; 432&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; do jj1=1,3<BR>&nbsp;&nbsp;&nbsp; 433&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; do ii1=1,39<BR>&nbsp;&nbsp;&nbsp; 434&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if(xmlabel1(ii1:ii1+1).eq.'PX') then<BR>&nbsp;&nbsp;&nbsp; 435&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index_shift = index_shift+4<BR>&nbsp;&nbsp;&nbsp; 436&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(! ii1+6:12+index_shift) = xmlabel1(ii1+2:12+index_shift-4) <BR>&nbsp;&nbsp;&nbsp; 437&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1:ii1+5) ='p'//achar(92)//'sx'//achar(92)//'N'<BR>&nbsp;&nbsp;&nbsp; 438&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif<BR>&nbsp;&nbsp;&nbsp; 439&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if(xmlabel1(ii1:ii1+1).eq.'PY') then<BR>&nbsp;&nbsp;&nbsp; 440&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index_shift = index_shift+4<BR>&nbsp;&nbsp;&nbsp; 441&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1+6:12+index_shift) = xmlabel1(ii1+2:12+index_shift-4) <BR>&nbsp;&nbsp;&nbsp; 442&nbsp;&nbsp;&amp;nb! sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&amp; nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1:ii1+5) ='p'//achar(92)//'sy'//achar(92)//'N'<BR>&nbsp;&nbsp;&nbsp; 443&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif<BR>&nbsp;&nbsp;&nbsp; 444&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if(xmlabel1(ii1:ii1+1).eq.'PZ') then<BR>&nbsp;&nbsp;&nbsp; 445&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index_shift = index_shift+4<BR>&nbsp;&nbsp;&nbsp; 446&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1+6:12+index_shift) = xmlabel1(ii1+2:12+index_shift-4) <BR>&nbsp;&nbsp;&nbsp; 447&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1:ii1+5) ='p'//achar(92)//'sz'//achar(92)//'N'<BR>&amp;n! bsp;&nbsp;&nbsp; 448&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif<BR>&nbsp;&nbsp;&nbsp; 449&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; enddo !ii1<BR>&nbsp;&nbsp;&nbsp; 450&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; do ii1=1,39<BR>&nbsp;&nbsp;&nbsp; 451&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if(xmlabel1(ii1:ii1+2).eq.'DZ2') then<BR>&nbsp;&nbsp;&nbsp; 452&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index_shift = index_shift+6<BR>&nbsp;&nbsp;&nbsp; 453&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1+9:12+index_shift) = xmlabel1(ii1+3:12+index_shift-6) <BR>&nbsp;&nbsp;&nbsp; 454&nbsp;&nbsp;&nbsp;&nbsp;&amp;n! bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1:ii1+8) ='d'//achar(92)//'sz'//achar(92)//'S2'//achar(92)//'N'<BR>&nbsp;&nbsp;&nbsp; 455&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif<BR>&nbsp;&nbsp;&nbsp; 456&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if(xmlabel1(ii1:ii1+2).eq.'DXY') then<BR>&nbsp;&nbsp;&nbsp; 457&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index_shift = index_shift+4<BR>&nbsp;&nbsp;&nbsp; 458&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1+7:12+index_shift) = xmlabel1(ii1+3:12+index_shift-4) <BR>&nbsp;&nbsp;&nbsp; 459&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1:ii1+6) ='d'//achar(92)//'sxy'//achar(92)//'! N'//achar(92)//'N'<BR>&nbsp;&nbsp;&nbsp; 460&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif<BR>&nbsp;&nbsp;&nbsp; 461&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if(xmlabel1(ii1:ii1+2).eq.'DXZ') then<BR>&nbsp;&nbsp;&nbsp; 462&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index_shift = index_shift+4<BR>&nbsp;&nbsp;&nbsp; 463&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1+7:12+index_shift) = xmlabel1(ii1+3:12+index_shift-4) <BR>&nbsp;&nbsp;&nbsp; 464&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1:ii1+6) ='d'//achar(92)//'sxz'//achar(92)//'N'//achar(92)//'N'<BR>&nbsp;&nbsp;&nbsp; 465&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;! &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; end if<BR>&nbsp;&nbsp;&nbsp; 466&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if(xmlabel1(ii1:ii1+2).eq.'DYZ') then<BR>&nbsp;&nbsp;&nbsp; 467&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index_shift = index_shift+4<BR>&nbsp;&nbsp;&nbsp; 468&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1+7:12+index_shift) = xmlabel1(ii1+3:12+index_shift-4) <BR>&nbsp;&nbsp;&nbsp; 469&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1:ii1+6) ='d'//achar(92)//'syz'//achar(92)//'N'//achar(92)//'N'<BR>&nbsp;&nbsp;&nbsp; 470&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif<BR>&nbsp;&nbsp;&nbsp; 471&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&amp;nbsp! ;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; enddo !ii1<BR>&nbsp;&nbsp;&nbsp; 472&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; do ii1=1,39<BR>&nbsp;&nbsp;&nbsp; 473 if(xmlabel1(ii1:ii1+4).eq.'DX2Y2') then<BR>&nbsp;&nbsp;&nbsp; 474&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index_shift = index_shift+13<BR>&nbsp;&nbsp;&nbsp; 475&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1+18:12+index_shift) = xmlabel1(ii1+5:12+index_shift-13) <BR>&nbsp;&nbsp;&nbsp; 476&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1:ii1+17)='d'//achar(92)//'sx'//achar(92)//'S2'//achar(92)//'N'//achar(92)//'s-y'//achar(92)//'S2'//achar(92)//'N'<BR>&nbsp;&nbsp;&nbsp; 477&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif<BR>==============================! =======================================</SPAN></DIV>
<DIV><SPAN><BR><BR>Can you please give me some tips about it, thanks in advance.</SPAN></DIV>
<DIV><SPAN></SPAN>&nbsp;</DIV>
<DIV><SPAN>Yours,</SPAN></DIV>
<DIV><SPAN></SPAN><SPAN><BR>H.H. GUO</SPAN></DIV>
<DIV><SPAN></SPAN></DIV></SPAN><BR><SPAN>Magnetism and Magnetic Materials Division<BR>Shenyang Materials Science National Laboratory<BR>Institute of Metal Research <BR>Chinese Academy of Sciences<BR>72 Wenhua Road,Shenyang 110016, China<BR><BR><BR>+86-15140243901 (mobile)<BR>work: &nbsp; &nbsp; <A href="mailto:hhguo@imr.ac.cn" target="_blank">hhguo@imr.ac.cn</A><BR></SPAN><BR><BR><BR></BLOCKQUOTE><BR><SPAN></SPAN><br><br><br>