[Wien] QTL WARN differ between lapw2 and lapw2 -qtl
Peter Blaha
pblaha at theochem.tuwien.ac.at
Tue Nov 14 10:11:59 CET 2017
> 1) Why "x lapw2" and "x lapw2 -qtl" show a different warning?
because x lapw2 calculates only occupied states, while x lapw2
-qtl also includes all unoccupied states in you energy window.
The first :WAR is probably negligible. You don't show us your in1 nor
the relevant part of the scf1up/dn files, so I don't know where the 3d
energy parameter actually is. It comes because you used a rather large
sphere and maybe it comes from the -in1new switch.
The second problem is eventually a problem when you calculate a DOS for
higher energies.
Eventually, try to set a HCLO for Ni-d.
PS: why are you doing non-spinpolarized calculations for Ni ???
>
> Description:
> I caliculate fcc Ni( struct file is attached below) by "run_lapw -p -ec
> 0.0001 -cc 0.0001"and I got QTL-B warn from scf file:
> ```
> :WARN : QTL-B value eq. 2.82 in Band of energy 0.64306 ATOM= 1
> L= 2
> :WARN : You should change the E-parameter for this atom and L-value in
> case.in1 (or try the -in1new switch)
> ```
>
> OK, to research the problem, I caliculate QTL value from "x lapw2 -p
> -qtl." But there was only this warn( this is ok because I don' interest
> in energy ~ 2.0Ry ):
> ```
> QTL-B VALUE .EQ. 52.52583 in Band of energy 2.21312 ATOM= 1 L= 2
> Check for ghostbands or EIGENVALUES BELOW XX messages
> Adjust your Energy-parameters for this ATOM and L (or use -in1new
> switch), check RMTs !!!
> ```
>
> And QTL for Band energy 0.64306:
> BAND 9 E= 0.64306 WEIGHT=-0.0009021
> L= 0 L= 1 L= 2 D-EG: D-T2G: L= 3
> QINSID: 3.9478 19.1584 64.5465 28.0223 36.5242 0.2695
> Q(U) : 3.9410 19.1380 61.7270 26.8170 34.9100 0.2690
> Q(UE) : 0.0070 0.0190 2.8190 1.2050 1.6140 0.0000
> QOUT : 11.9931
>
> I think that there is no strange behavior(Q(U)>>Q(UE)). So, there should
> be no QTL-B warn.
> Please teach me where is my mistake.
>
> Note:
> 1) WIEN2k version is 17_1
> 2) I took Ni_RMT=2.32, RKMAX=9.0 and lvns=10.
> 3) I try to run_lapw with "-in1new 2" option, but this does not vanish
> this warning.
> 4) I try to change EPARAMETER=0.3(0.4, 0.5). I think this is reasonable
> for this situation(EF=0.568), but these do not vanish this warning also.
> 5) My struct file is following:
> ```
> F LATTICE,NONEQUIV.ATOMS: 1
> 225_Fm-3m
> MODE OF CALC=RELA
> unit=ang
> 6.614044 6.614044 6.614044 90.000000 90.000000
> 90.000000
> ATOM 1: X=0.00000000 Y=0.00000000 Z=0.00000000
> MULT= 1 ISPLIT= 2
> Ni NPT= 781 R0=0.00005000 RMT= 2.32 Z:
> 28.000
> LOCAL ROT MATRIX: 1.0000000 0.0000000 0.0000000
> 0.0000000 1.0000000 0.0000000
> 0.0000000 0.0000000 1.0000000
> ```
>
> With Best regards,
>
>
> _______________________________________________
> 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
>
--
P.Blaha
--------------------------------------------------------------------------
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
--------------------------------------------------------------------------
More information about the Wien
mailing list