  We investigate regularity properties of molecular one-electron densities rho
near the nuclei. In particular we derive a representation rho(x)=mu(x)*(e^F(x))
with an explicit function F, only depending on the nuclear charges and the
positions of the nuclei, such that mu belongs to C^{1,1}(R^3), i.e., mu has
locally essentially bounded second derivatives. An example constructed using
Hydrogenic eigenfunctions shows that this regularity result is sharp. For
atomic eigenfunctions which are either even or odd with respect to inversion in
the origin, we prove that mu is even C^{2,\alpha}(R^3) for all alpha in (0,1).
Placing one nucleus at the origin we study rho in polar coordinates x=r*omega
and investigate rho'(r,omega) and rho''(r,omega) for fixed omega as r tends to
zero. We prove non-isotropic cusp conditions of first and second order, which
generalize Kato's classical result.
