Electrostatic potentials and fields from density expansions of deformed atoms in molecules

The exact representation of the molecular density by means of atomic expansions, consisting in spherical harmonics times analytical radial factors, is employed for the calculation of electrostatic potentials, fields, and forces. The resulting procedure is equivalent to an atomic multipolar expansion in the long-range regions, but works with similar efficiency and accuracy in the short-range region, where multipolar expansions are not valid. The performances of this procedure are tested on the calculation of the electrostatic potential contour maps and electrostatic field flux lines of water and nitrobenzene, computed from high-quality molecular electron densities obtained with Slater basis sets.