[Wien] LSDA+U on Yb compounds

[Stefaan Cottenier]

Can you verify explicitely whether you have 13 or 14 4f-electrons? Do  
this by looking at the partial charges for f (:QTL), or by summing the  
diagonal elements of the up and dn density matrices, or by looking at  
the value of the integrated partial f-dos at the fermi level in  
case.outputt (and preferably do it by all three methods, they should  
give the same answer).

If you really have 14 electrons, you might have ended up with divalent  
Yb. If you have 13 electrons, perhaps your unoccupied state is diffuse  
and not at first sight recognizable as a peak.

Stefaan

> I am performing LSDA+U calculations on a series of Yb compounds in which the
> valence of Yb is close to 3+. I was expecting to get a splitting between the
> occupied and unoccupied 4f states, but my DOS only shows the occupied
> 4f-states located well below the Fermi level, without any empty 4f-states.
> I'm doing the following:
>
> runsp_lapw  -cc 0.001;
> runsp_lapw -so -cc 0.001;
> runsp_lapw -so -orb -cc 0.001.
>
> for DOS I do:
> x lapw1 -orb -up/dn
> x lapwso -up -orb
> x lapw2 -qtl -so -up/dn.
>
>
> It also seems that with LSDA+U there is no difference between -up and -dn
> states. (and I'm not using runsp_c_lapw).
> Is there any wien2k user who was able to observe the empty 4f states of Yb,
> located at or above Fermi level?
>
> Thanks,
> Zsolt Rak
>
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
>
>



-- 
Stefaan Cottenier
Instituut voor Kern- en Stralingsfysica
K.U.Leuven
Celestijnenlaan 200 D
B-3001 Leuven (Belgium)

tel: + 32 16 32 71 45
fax: + 32 16 32 79 85
e-mail: stefaan.cottenier at fys.kuleuven.be


Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm