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

Peter Blaha pblaha at theochem.tuwien.ac.at
Sun Apr 5 13:00:24 CEST 2020


As expected, your struct file is the problem.  Why are you using P1 
symmetry ??

In this paper they did a NON-SPINPOLARIZED calculation in LDA!!!

So when you want to reproduce their results, remove from your struct 
file all the labels, i.e.

"V 1"  --> "V  "   and so on.

Please remember, the struct file is sensitive to the colums. So you have 
to overwrite the 1,2,.. by blanks, you must NOT delete the 1, because 
then the rest of the line moves to the left!

Then follow the recommendations of nn and sgroup.

You should obtain a struct file with 2 non-equivalent positions (one V, 
one O)  and ISPLI=4 instead if ISPLIT=8



Am 05.04.2020 um 08:05 schrieb Gavin Abo:
> If I interpret Prof. Blaha's post at [1] correctly, "x lapw2 -qtl" 
> should automatically give what you need such that you would not have to 
> use "x qtl" (thus, you would not need case.inq if that is the case).
> 
> Referring to the description for line 6 in section '4.3 The “master 
> input“ file case.struct' of the WIEN2k 19.1 usersguide [2] on page 41, 
> you will probably just want to check to insure that isplit was set to 4 
> in case.struct after the "x symmetry" step during initialization with 
> init_lapw.
> 
> You might recall the TiC getting started example of section "3.11.2 
> Density of States (DOS)" on page 24 of the usersguide [2].  Your 
> description below kind of sounds like you might have only plotted 'tot' 
> for the total DOS rather than the PDOS for an individual atom.
> 
> For your V2O3 calculation, I think you will want to make sure you 
> selected the d-z2, (d-xy,dx2-y2), and (d-xz,dyz) by setting case.int 
> appropriately.
> 
> [1] 
> https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg19805.html
> [2] http://susi.theochem.tuwien.ac.at/reg_user/textbooks/usersguide.pdf
> 
> On 4/4/2020 11:26 PM, Wasim Raja Mondal wrote:
>> Gavin,
>>          Thank you for your reply. Now, I understand the group table. 
>> I am getting a1g, eg\pi and \eg \sigma as reported in literature for 
>> example in Fig5 of the paper I provided in the link. But, in my DOS 
>> a1g, eg\pi and eg\sigma are still mixed up.
>>
>> For that, I tried in case.inq file different new z and x axis. But I 
>> am not getting the DOS as reported in the Fig 5 of the paper (link I 
>> provided.)
>>
>> I know I should not ask this type of question here. But I am 
>> struggling bit for many days with this structure. My question is : 
>> starting with cif file , running LAPW, calculating DOS using qtl with 
>> option 4, what new z and x axis should I define in case.inq file to 
>> reproduce the LDA dos for this structure. It is just simple LDA 
>> calculation.
>>
>>
>> Please let me know if you find with just one look.
>>
>>
>> Thanks
>> Wasim
> _______________________________________________
> 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.Materials Chemistry
TU Vienna
Getreidemarkt 9
A-1060 Vienna
Austria
+43-1-5880115671


More information about the Wien mailing list