<div dir="ltr"><div dir="ltr"><div dir="ltr">Dear Prof. Blaha,<div>                            Thank you for your reply. I understand now the space group and also no the unitary transformation for the V2O3 which I am giving here:</div><div><br></div><div>                               a1g =  d3z2-1</div><div><br></div><div>                            eg\pi 1= \sqrt(2/3) dxy + 1/sqrt(3) dxz</div><div><br></div><div>                           eq\pi2 = -\sqrt(2/3) dx2-y2 -1/sqrt(3) dyz</div><div><br></div><div>I want to use this unitary transformation. For that My case.inq file is given below:</div><div><br></div><div><br></div><div><pre style="color:rgb(0,0,0)">-9.  3.          Emin Emax
 2               number of selected atoms
 1  -2  0  0      iatom1  qsplit1 symmetrize loro
 2  0  1         nL1 p  d
 2  6  0  0      iatom2  qsplit2 symmetrize loro
 4  0  1  2  3   nL2 s  p  d  f</pre><pre style="color:rgb(0,0,0)"><br></pre><pre style="color:rgb(0,0,0)">In the above file, iatom1 is oxygen and iatom2 is vanadium.
<br></pre><pre style="color:rgb(0,0,0)">From wien2kroot, I have copied case.cf_eg_t2g and renamed it as case.cf2. Next, when I run x qtl, I am getting this error:</pre><pre style="color:rgb(0,0,0)"><br></pre><pre style="color:rgb(0,0,0)"><p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">forrtl: severe (59): list-directed I/O syntax error, unit -5, file Internal List-Directed Read</span></p>
<p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">Image              PC                Routine            Line        Source             </span></p>
<p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">qtl                000000000043A2BB  Unknown               Unknown  Unknown</span></p>
<p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">qtl                00000000004577E1  Unknown               Unknown  Unknown</span></p>
<p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">qtl                000000000041D392  readc_                     44  readc.f</span></p>
<p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">qtl                000000000041969A  MAIN__                    330  qtlmain.f</span></p>
<p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">qtl                00000000004046E2  Unknown               Unknown  Unknown</span></p>
<p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"><a href="http://libc-2.17.so">libc-2.17.so</a>       00002B7E158E2505  __libc_start_main     Unknown  Unknown</span></p>
<p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">qtl                00000000004045E9  Unknown               Unknown  Unknown</span></p>
<p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">0.011u 0.013s 0:00.02 100.0%<span class="gmail-Apple-tab-span">  </span>0+0k 0+32io 0pf+0w</span></p>
<p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">error: command   /usr/local/wien2k-19.1/qtl qtl.def   failed</span></p><p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">Should I have to edit new2.cif file according to my unitary matrix as I described above my a1g,eg\pi,eg\sigma? I am using wien2k 19.1</span></p><p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">Thanking you</span></p><p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-stretch:normal;font-size:15px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">Wasim</span></p></pre><pre style="color:rgb(0,0,0)"><br></pre></div><div><br></div><div>  </div><div> </div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Mar 26, 2020 at 3:59 AM Peter Blaha <<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">Wien2k will automatically calculate the symmetry related splitting of <br>
the d-states.<br>
<br>
Example 1: Octahedral (cubic) symmetry:  WIEN2k sets ISPLIT=2 in <br>
case.struct during initialization. This symmetry was detected by x <br>
symmetry and you can find the proper point group of your atoms in <br>
case.outputs.<br>
Because of the ISPLIT=2, when you run   x lapw2 -qtl, the case.qtl file <br>
will automatically contain for an atom with eg. Oh symmetry:<br>
s,p,d,d-eg,d-t2g,f<br>
<br>
Please look at the header of case.qtl, where it tells you the <br>
decomposition of the states<br>
<br>
Example 2.: hexagonal symmetry (as in hcp Zn,Be or Mg) In case.outputs <br>
you find:<br>
pointgroup is -6m2, leading to ISPLIT=4 (see UG for definitions of ISPLIT)<br>
and the header of case.qtl givs:<br>
s,p,pz,pxy,d,d-z2,(d-xy,dx2-y2),(d-xz,dyz),f<br>
So the p-states have been automatically splitted into pz and px+py, <br>
while the 5 dstates are splitted into 3 groups, one 1-dimensional and 2 <br>
2-dimensional irreducible representations. Wien2k does not give you the <br>
"label" of the irrep (like A1g or Eg) as you are looking for, but use a <br>
good "group-theory character table" for pointgroup -6m2 and you will <br>
find the "name of the representation" (A1g,..) and the corresponding <br>
"basis functions" (z, x+y; z**2, xz,xz, ..) so that you can make the <br>
corresponding asignment between or splitted partial charges and the <br>
"labels".<br>
<br>
More complicated is the calculation of "approximate symmetries" (like <br>
eg-t2g) in a distorted octahedron), because in that case WIEN2k cannot <br>
do it automatically for you, but you can combine the single orbitals in <br>
the PROPER coordinate system yourself. In this case the x qtl program <br>
can be useful since it allows to calculate the PDOS in a rotated <br>
coordinate frame pointing into your "approximate octahedron".<br>
<br>
PS: In corundum structure, usually it also has ISPLIT=4 and the d-z2 <br>
orbital corresponds to A1g....<br>
<br>
<br>
Am 26.03.2020 um 01:31 schrieb Wasim Raja Mondal:<br>
> Dear Wien2k experts,<br>
>                                    I am doing DOS calculation to <br>
> reproduce some data for V2O3 corundum trigonal structure. In this <br>
> distorted structure octahedra of oxygen atoms, d orbital should be split <br>
> into single degenerate a1g and double degenerate eg. I have calculated <br>
> projected density of states and I am not finding that. So, I have <br>
> realized that I have to use rotation and for that I want to use QTL <br>
> program. So far, I have done following steps:<br>
> <br>
>        1. x cif2struct<br>
> <br>
>        2. init_lapw<br>
> <br>
>       3. run_lapw<br>
> <br>
>      4.  edited   case.inq  and used the option qsplit=6<br>
> <br>
>     5. copied templet for case_cf_eg_t2g in my directory.<br>
>     6. cp case_cf_eg_t2g to case_cf2 since my atom number V is assigned as 2<br>
> <br>
>    7. Getting error.<br>
> <br>
> <br>
> Can you suggest any solution? Is this the way to project d to a1g and eg?<br>
> <br>
> Thanks in advance<br>
> Wasim<br>
> <br>
> <br>
> On Wed, Mar 25, 2020 at 1:49 PM Wasim Raja Mondal <br>
> <<a href="mailto:wasimr.mondal@gmail.com" target="_blank">wasimr.mondal@gmail.com</a> <mailto:<a href="mailto:wasimr.mondal@gmail.com" target="_blank">wasimr.mondal@gmail.com</a>>> wrote:<br>
> <br>
>     Dear expert,<br>
>                         I am trying to project out t2g and a1g states of<br>
>     V d orbitals. For that I want to use qtl program for DOS<br>
>     calculation. I am chosing QSPLIT=6 for unitary transformation. For<br>
>     that I copied templet as case.cf_d_eg_t2g .<br>
> <br>
>     Can you please suggest me how can I proceed next?<br>
> <br>
> <br>
> <br>
>     Thanks<br>
> <br>
> <br>
> _______________________________________________<br>
> Wien mailing list<br>
> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <br>
<br>
-- <br>
--------------------------------------------------------------------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300             FAX: +43-1-58801-165982<br>
Email: <a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a>    WIEN2k: <a href="http://www.wien2k.at" rel="noreferrer" target="_blank">http://www.wien2k.at</a><br>
WWW: <br>
<a href="http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------" rel="noreferrer" target="_blank">http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------</a> <br>
<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>