Dear Prof.<span title="Gavin Abo <gsabo@crimson.ua.edu>"> Gavin Abo</span>:<br>
<br>
Thanks for your reply. The column counting is more clear to me. For instance, <br>
if PZ and DZ2 is concerned, according to the 'case.qtl file',<br>
<br>
--------------------<br>
LATTICE CONST.= 5.5755 5.5755 20.2062 FERMI ENERGY= 0.99511<br>
425 < NMAT < 463 SPIN=1 NAT= 6 SO 0<br>
JATOM 1 MULT= 1 ISPLIT= 8 tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3<br>
JATOM 2 MULT= 1 ISPLIT= 8 tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3<br>
-----------------------<br>
<br>
I use "6" and "8" for the "jtype", respectively. <br><br>Now I have one new question, do you think if it is possible to mess up with cathesian direction? <br>I mean, I expect to get a big weight from Pz orbital, but result gives out Px, similar problem <br>occurs between the weight of DZ2 and that of DX2Y2. I have double checked the Catesian <br>
coordinate for the structure, but found no problem. Do you think if there exists some bugs in <br>the process of orbital decomposition? My version is WIEN2K_09.1.<br>
<br>
Thanks again!
<p>Yours,</p>
<p>H.H.GUO</p>
<span>--<br>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>+86-15140243901 (mobile)<br>work: hhguo@imr.ac.cn</span><blockquote name="replyContent" style="padding-left:5px;margin-left:5px;border-left:#b6b6b6 2px solid;margin-right:0">-----原始邮件-----<br>
<b>发件人:</b> "Gavin Abo" <gsabo@crimson.ua.edu><br>
<b>发送时间:</b> 2012年7月8日 星期日<br>
<b>收件人:</b> "A Mailing list for WIEN2k users" <wien@zeus.theochem.tuwien.ac.at><br>
<b>抄送:</b> <br>
<b>主题:</b> Re: [Wien] "case.insp" to define band character for plotting<br><br>
<div class="moz-cite-prefix">Click "Edit case.insp" in w2web under
Utils. -> Bandstructure. You should see on the screen
something like:<br>
<br>
ATOM atom_name: jatom tot, s, p, d, ...<br>
<br>
If it is not seen in a version older than Wien2k 11.1 or you
cannot use w2web, you can see around line 5 in case.qtl, where one
line should look like:<br>
<br>
JATOM 1 MULT= 1 ISPLIT = 2 tot,0,1,2, ...<br>
<br>
The (0,1,2,3) in case.qtl correspond to (s,p,d,f).<br>
<br>
The jtype (1,2,3,4,...) usually corresponds to (tot,s,p,d,...).
However, look at your case.qtl as mentioned above and count the
columns to determine the number for jtype.<br>
<br>
If you use "plot bandstructure" in w2web, it should also label on
the plot such as "tot" if jtype 1 was used in case.insp.<br>
<br>
On 7/8/2012 7:17 AM, guohuaihong wrote:<br>
</div>
<blockquote cite="mid:1f0b2c9.4781.13866bd0d71.Coremail.hhguo@imr.ac.cn" type="cite"><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. Is jtype (0,1,2,3) correspoding to orbital
(s,p,d,f) for each atom speficified ?<br>
<br>
</span></div>
<div><span>2. 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> 430 xmlabel1(:)=label1(:)<br>
431 index_shift=0<br>
432 do jj1=1,3<br>
433 do ii1=1,39<br>
434 if(xmlabel1(ii1:ii1+1).eq.'PX')
then<br>
435 index_shift = index_shift+4<br>
436 xmlabel1(! ii1+6:12+index_shift)
= xmlabel1(ii1+2:12+index_shift-4) <br>
437 xmlabel1(ii1:ii1+5)
='p'//achar(92)//'sx'//achar(92)//'N'<br>
438 endif<br>
439 if(xmlabel1(ii1:ii1+1).eq.'PY')
then<br>
440 index_shift = index_shift+4<br>
441 xmlabel1(ii1+6:12+index_shift) =
xmlabel1(ii1+2:12+index_shift-4) <br>
442 &nb! sp; &
nbsp; xmlabel1(ii1:ii1+5)
='p'//achar(92)//'sy'//achar(92)//'N'<br>
443 endif<br>
444 if(xmlabel1(ii1:ii1+1).eq.'PZ')
then<br>
445 index_shift = index_shift+4<br>
446 xmlabel1(ii1+6:12+index_shift) =
xmlabel1(ii1+2:12+index_shift-4) <br>
447 xmlabel1(ii1:ii1+5)
='p'//achar(92)//'sz'//achar(92)//'N'<br>
&n! bsp; 448 endif<br>
449 enddo !ii1<br>
450 do ii1=1,39<br>
451 if(xmlabel1(ii1:ii1+2).eq.'DZ2')
then<br>
452 index_shift = index_shift+6<br>
453 xmlabel1(ii1+9:12+index_shift) =
xmlabel1(ii1+3:12+index_shift-6) <br>
454 &n! bsp;
xmlabel1(ii1:ii1+8)
='d'//achar(92)//'sz'//achar(92)//'S2'//achar(92)//'N'<br>
455 endif<br>
456 if(xmlabel1(ii1:ii1+2).eq.'DXY')
then<br>
457 index_shift = index_shift+4<br>
458 xmlabel1(ii1+7:12+index_shift) =
xmlabel1(ii1+3:12+index_shift-4) <br>
459 xmlabel1(ii1:ii1+6)
='d'//achar(92)//'sxy'//achar(92)//'! N'//achar(92)//'N'<br>
460 endif<br>
461 if(xmlabel1(ii1:ii1+2).eq.'DXZ')
then<br>
462 index_shift = index_shift+4<br>
463 xmlabel1(ii1+7:12+index_shift) =
xmlabel1(ii1+3:12+index_shift-4) <br>
464 xmlabel1(ii1:ii1+6)
='d'//achar(92)//'sxz'//achar(92)//'N'//achar(92)//'N'<br>
465 ! end
if<br>
466 if(xmlabel1(ii1:ii1+2).eq.'DYZ')
then<br>
467 index_shift = index_shift+4<br>
468 xmlabel1(ii1+7:12+index_shift) =
xmlabel1(ii1+3:12+index_shift-4) <br>
469 xmlabel1(ii1:ii1+6)
='d'//achar(92)//'syz'//achar(92)//'N'//achar(92)//'N'<br>
470 endif<br>
471 ! ; enddo !ii1<br>
472 do ii1=1,39<br>
473 if(xmlabel1(ii1:ii1+4).eq.'DX2Y2') then<br>
474 index_shift = index_shift+13<br>
475 xmlabel1(ii1+18:12+index_shift)
= xmlabel1(ii1+5:12+index_shift-13) <br>
476
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>
477 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> </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: <a target="_blank" class="moz-txt-link-abbreviated" href="mailto:hhguo@imr.ac.cn">hhguo@imr.ac.cn</a><br>
</span><br>
<br>
<br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Wien mailing list
<a target="_blank" class="moz-txt-link-abbreviated" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>
<a target="_blank" class="moz-txt-link-freetext" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>
</pre>
</blockquote>
<br>
<br>
</blockquote><span><br></span><br><br><br>