[Wien] Using option Qsplit=6 in qtl to rotate to get t2g and eg

Wasim Raja Mondal wasimr.mondal at gmail.com
Sat Mar 28 01:11:08 CET 2020


Dear Prof. Blaha,
                            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:

                               a1g =  d3z2-1

                            eg\pi 1= \sqrt(2/3) dxy + 1/sqrt(3) dxz

                           eq\pi2 = -\sqrt(2/3) dx2-y2 -1/sqrt(3) dyz

I want to use this unitary transformation. For that My case.inq file is
given below:


-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


In the above file, iatom1 is oxygen and iatom2 is vanadium.

>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:


forrtl: severe (59): list-directed I/O syntax error, unit -5, file
Internal List-Directed Read

Image              PC                Routine            Line
Source

qtl                000000000043A2BB  Unknown               Unknown  Unknown

qtl                00000000004577E1  Unknown               Unknown  Unknown

qtl                000000000041D392  readc_                     44  readc.f

qtl                000000000041969A  MAIN__                    330  qtlmain.f

qtl                00000000004046E2  Unknown               Unknown  Unknown

libc-2.17.so       00002B7E158E2505  __libc_start_main     Unknown  Unknown

qtl                00000000004045E9  Unknown               Unknown  Unknown

0.011u 0.013s 0:00.02 100.0%	0+0k 0+32io 0pf+0w

error: command   /usr/local/wien2k-19.1/qtl qtl.def   failed



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


Thanking you


Wasim






On Thu, Mar 26, 2020 at 3:59 AM Peter Blaha <pblaha at theochem.tuwien.ac.at>
wrote:

> Wien2k will automatically calculate the symmetry related splitting of
> the d-states.
>
> Example 1: Octahedral (cubic) symmetry:  WIEN2k sets ISPLIT=2 in
> case.struct during initialization. This symmetry was detected by x
> symmetry and you can find the proper point group of your atoms in
> case.outputs.
> Because of the ISPLIT=2, when you run   x lapw2 -qtl, the case.qtl file
> will automatically contain for an atom with eg. Oh symmetry:
> s,p,d,d-eg,d-t2g,f
>
> Please look at the header of case.qtl, where it tells you the
> decomposition of the states
>
> Example 2.: hexagonal symmetry (as in hcp Zn,Be or Mg) In case.outputs
> you find:
> pointgroup is -6m2, leading to ISPLIT=4 (see UG for definitions of ISPLIT)
> and the header of case.qtl givs:
> s,p,pz,pxy,d,d-z2,(d-xy,dx2-y2),(d-xz,dyz),f
> So the p-states have been automatically splitted into pz and px+py,
> while the 5 dstates are splitted into 3 groups, one 1-dimensional and 2
> 2-dimensional irreducible representations. Wien2k does not give you the
> "label" of the irrep (like A1g or Eg) as you are looking for, but use a
> good "group-theory character table" for pointgroup -6m2 and you will
> find the "name of the representation" (A1g,..) and the corresponding
> "basis functions" (z, x+y; z**2, xz,xz, ..) so that you can make the
> corresponding asignment between or splitted partial charges and the
> "labels".
>
> More complicated is the calculation of "approximate symmetries" (like
> eg-t2g) in a distorted octahedron), because in that case WIEN2k cannot
> do it automatically for you, but you can combine the single orbitals in
> the PROPER coordinate system yourself. In this case the x qtl program
> can be useful since it allows to calculate the PDOS in a rotated
> coordinate frame pointing into your "approximate octahedron".
>
> PS: In corundum structure, usually it also has ISPLIT=4 and the d-z2
> orbital corresponds to A1g....
>
>
> Am 26.03.2020 um 01:31 schrieb Wasim Raja Mondal:
> > Dear Wien2k experts,
> >                                    I am doing DOS calculation to
> > reproduce some data for V2O3 corundum trigonal structure. In this
> > distorted structure octahedra of oxygen atoms, d orbital should be split
> > into single degenerate a1g and double degenerate eg. I have calculated
> > projected density of states and I am not finding that. So, I have
> > realized that I have to use rotation and for that I want to use QTL
> > program. So far, I have done following steps:
> >
> >        1. x cif2struct
> >
> >        2. init_lapw
> >
> >       3. run_lapw
> >
> >      4.  edited   case.inq  and used the option qsplit=6
> >
> >     5. copied templet for case_cf_eg_t2g in my directory.
> >     6. cp case_cf_eg_t2g to case_cf2 since my atom number V is assigned
> as 2
> >
> >    7. Getting error.
> >
> >
> > Can you suggest any solution? Is this the way to project d to a1g and eg?
> >
> > Thanks in advance
> > Wasim
> >
> >
> > On Wed, Mar 25, 2020 at 1:49 PM Wasim Raja Mondal
> > <wasimr.mondal at gmail.com <mailto:wasimr.mondal at gmail.com>> wrote:
> >
> >     Dear expert,
> >                         I am trying to project out t2g and a1g states of
> >     V d orbitals. For that I want to use qtl program for DOS
> >     calculation. I am chosing QSPLIT=6 for unitary transformation. For
> >     that I copied templet as case.cf_d_eg_t2g .
> >
> >     Can you please suggest me how can I proceed next?
> >
> >
> >
> >     Thanks
> >
> >
> > _______________________________________________
> > Wien mailing list
> > Wien at zeus.theochem.tuwien.ac.at
> > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
> > SEARCH the MAILING-LIST at:
> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
> >
>
> --
> --------------------------------------------------------------------------
> Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
> Phone: +43-1-58801-165300             FAX: +43-1-58801-165982
> Email: blaha at theochem.tuwien.ac.at    WIEN2k: http://www.wien2k.at
> WWW:
>
> http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------
>
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
> SEARCH the MAILING-LIST at:
> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20200327/59ac64fd/attachment.html>


More information about the Wien mailing list