<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>&nbsp;&nbsp;&nbsp; 430&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(:)=label1(:) &nbsp;&nbsp;&nbsp; 431&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index_shift=0 &nbsp;&nbsp;&nbsp; 432&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; do jj1=1,3 &nbsp;&nbsp;&nbsp; 433&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; do ii1=1,39 &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 &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 &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)  &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' &nbsp;&nbsp;&nbsp; 438&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif &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 &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 &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)  &nbsp;&nbsp;&nbsp; 442&nbsp;&nbsp;&nb!

 sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&

nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmlabel1(ii1:ii1+5) ='p'//achar(92)//'sy'//achar(92)//'N' &nbsp;&nbsp;&nbsp; 443&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif &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 &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 &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)  &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' &n!

bsp;&nbsp;&nbsp; 448&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif &nbsp;&nbsp;&nbsp; 449&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; enddo !ii1 &nbsp;&nbsp;&nbsp; 450&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; do ii1=1,39 &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 &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 &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)  &nbsp;&nbsp;&nbsp; 454&nbsp;&nbsp;&nbsp;&nbsp;&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' &nbsp;&nbsp;&nbsp; 455&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif &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 &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 &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)  &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' &nbsp;&nbsp;&nbsp; 460&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif &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 &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 &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)  &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' &nbsp;&nbsp;&nbsp; 465&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;!

 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; end

if &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 &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 &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)  &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' &nbsp;&nbsp;&nbsp; 470&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif &nbsp;&nbsp;&nbsp; 471&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp!

;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; enddo !ii1 &nbsp;&nbsp;&nbsp; 472&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; do ii1=1,39 &nbsp;&nbsp;&nbsp; 473 if(xmlabel1(ii1:ii1+4).eq.'DX2Y2') then &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 &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)  &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' &nbsp;&nbsp;&nbsp; 477&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; endif ==============================!

 =======================================</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; hhguo@imr.ac.cn<BR></SPAN><br><br><br>