[Wien] problem with in1new in spin polarized calculations

[Ilya Elfimov]

I think, there is a problem with -in1new option in some materials  
with half filled d of f orbitals. Take EuO, for example.
Eu^2 has seven 4f electrons which are all spin parallel (S=7/2 high  
spin state).
Now, if all f-electrons are spin up then write_in1_lapw -up will  
produce case.in1new which has f orbitals in it. However,  
write_in1_lapw -dn will drop 4f's simply because their partial  
charges are almost zero. The final case.in1 is, of course, spin down  
one.
I would think that this is easy to fix by looking at the total charge  
in some orbital rather than spin up and spin down separately. I am  
wondering if this is on the list for the next version.
Related question: is there any strong reason to keep case.in1 spin  
independent?
It seems a bit strange to use energies obtained from spin down for  
spin up too.
In case if someone is interested, here is case.struct for EuO.

Thanks,
Ilya

-------------- next part --------------
A non-text attachment was scrubbed...
Name: case.struct
Type: application/octet-stream
Size: 4022 bytes
Desc: not available
Url : http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20060928/5be16361/case.obj
-------------- next part --------------