[Wien] Basic question about DOS of onsite-hybrid functional calculation
tran at theochem.tuwien.ac.at
tran at theochem.tuwien.ac.at
Thu Jan 28 09:55:44 CET 2016
Hi,
Yes, as you said, NE in case.in2 was not properly copied in
case.in2eece. This is due to a bug in the script runeece_lapw.
In runeece_lapw, replace these two lines
set nocut=`head -2 $file.in2${complex} |tail -1 | cut -c 12-80`
echo $emin $nocut >> $file.in2${complex}eece
by these two
set nocut=`head -2 $file.in2${complex} |tail -1 `
echo $emin $nocut[2-] >> $file.in2${complex}eece
Then it should be fine.
F. Tran
>Hi, all:
>
>I have a basic question on the DOS calculation of onsite-hybrid
functional
>(ONF) calculation. When I check the UG, I can’t find any notes about the
>DOS plotting of ONF, so I assume the program flow is the same in ONF as
>normal (x lapw2 -p -c -up/dn -qtl, then x tetra -p -up/dn). But when I
>compare the DOS from ONF and pure PBE-GGA calculations, I find that
>there’s no any change in the selected orbitals with ONF. This stays even
>when I change the mixing factor alpha from 0.25 to 1.0, which confuses me.
>Can anyone tell me whether the way I use lapw2 is wrong?
>
>Some details about my calculations:
>I’m calculating a Cr dopant in 4H-SiC (Si31 C32 Cr). Only the d orbitals
>of Cr is treated with ONF. HYBR mode is used.
>
>If any further information is needed, please let me know.
>
>Thanks,
>Wenhao
More information about the Wien
mailing list