[Wien] Electron density differences

[MoonfalL]

How to get the Electron density differences from Wien2K.
For example, GaN with hexagonal structure,
               *.in5 file:
                     0 0 0 4
                     4 8 0 4
                     0 0 8 4
                     6 6 6
                    100 100
                    DIFF SUB
                    ANG VAL NODEBUG
                    ORTHO
But I get the electron density difference all be plus; I think I should
subtract the total electron densities of the N and Ga atoms from the bulk,
but how can I do; It is something wrong with my input?

Thanks!!




                                                 Guan Pengfei
                                                 Tsinghua Universetiry
                                                 Department od Physics
                                                 Beijing   China
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20060628/8df60cd6/attachment.html