mod_rhons Module

The charge density is developed in spherical harmonics \begin{eqaution} \rho(r) = \sum_{lm} \rho(lm,r) Y(r,lm) } \end{equation} \begin{eqaution} \rho(lm,r) = \int \rho(r) * Y(r,lm) \end{equation} in the case of spin-polarization : the spin density is developed in spherical harmonics : \begin{eqaution} sden(r) = \sum_{lm}{ sden(lm,r) Y(r,lm) } \end{equation} \begin{eqaution} sden(lm,r) = \int sden(r) T(r,lm) \end{equation} $n(r,e)$ is developed in \begin{eqaution} n(r,e) = Y(r,l'm') n(l'm',lm,r,e) Y(r,lm) \end{equation} Therefore a faltung of n(l'm',lm,r,e) with the gaunt coeffients has to be used to calculate the lm-contribution of the charge density. Calculate the valence density of states , in the spin-polarized case spin dependent. recognize that the density of states is always complex also in the case of real-energy-integation (ief>0) since in that case the energy integration is done parallel to the real energy axis but not on the real energy axis. In the last energy-spin loop the l-contribution of the valence charge is calculated.