On the representation of electrostatic fields around ab initio charge distributions.

We compare two methods (Mulliken charges and a distributed multipole analysis, DMA) of representing an ab initio charge distribution for calculating the electrostatic field and potential outside the molecule, using pyrimidine and the RNA base uracil as examples. This is done using a 3-D graphical display of the electrostatic fields, which, when used with real-time rotation, zooming and clipping, has many advantages for qualitatively assessing the electrostatic interactions of a molecule. The errors involved in using Mulliken point charges may be of similar magnitude to the total electrostatic field in regions which are important in recognition processes. The DMA representation automatically includes the anisotropic electrostatic effects of non-spherical features in the charge distribution of each atom, and yet the displayed electrostatic fields around the atoms which have lone-pair density do not show marked anisotropy.