[Wien] What is p-1/2 LOs ?

Stefaan Cottenier stefaan.cottenier at fys.kuleuven.ac.be
Sun Mar 13 16:17:41 CET 2005


>     I am doing a calculation with spin orbit coupling. I need to define 
> the case.inso files, but I
> dont know how to add the p-1/2 LOs,ie. the El to increase accuracy as 
> the instruction said.I
> use the El in case.in1 files, but the result is poor convergence.I want 
> to know what is the p-1/2
> LOs and how to define them.I appreciate your help.

In a relativistic treatment, p-electrons are indexed not by l=1 and 
their spin s, but by the quantum number j, which is either l-s=1/2 or 
l+s=3/2. In a density of states, you see this as a split of the p-peaks 
in two parts. This splitting is reproduced by LDA+spin-orbit poorly, 
unless you go to very high Emax (>10). In order to get the correct 
splitting using smaller Emax (3), you can 'help' the basis set a bit by 
adding LO's that have the proper shape for these p-1/2 orbitals.

This is useful only for heavy atoms (5p, 6p). If this is not your case, 
then do not add p-1/2 LO's.

Stefaan




More information about the Wien mailing list