An improved method for computing potential-derived charges is described which is based upon the CHELP program available from QCPE.' This approach (CHELPG) is shown to be considerably less dependent upon molecular orientation than the original CHELP program. In the second part of this work, the CHELPG point selection algorithm was used to analyze the changes in the potential-derived charges in formamide during rotation about the C-N bond. In order to achieve a level of rotational invariance less than 10% of the magnitude of the electronic effects studied, a n equally-spaced array of points 0.3 A apart was required. Points found to be greater than 2.8 A from any nucleus were eliminated, along with all points contained within the defined VDW distances from each of the atoms. The results are compared to those obtained by using CHELP. Even when large numbers of points (ca. 3000) were sampled using the CHELP selection routine, the results did not indicate a satisfactory level of rotational invariance. On the basis of these results, the original CHELP program was found to be inadequate for analyzing internal rotations.

No abstract

We describe several improvements to the reaction field model for the ab initio determination of solvation effects. First, the simple spherical cavity model is expanded to include higher-order electrostatic interactions. Second, two new and efficient implementations of the polarizable continuum model (PCM) are described, which allow a more realistic specification of the solute cavity as well as infinite-order electrostatics. Electron correlation effects are evaluated using the B3LYP density functional and Möller−Plesset perturbation theory to second order. An assessment of the importance of these various factors is made by comparing theoretical results to the experimentally known conformational equilibrium between syn and anti forms of furfuraldehyde and the C−C rotational barrier of (2-nitrovinyl)amine. Comparisons are also made with calculations that employ an ellipsoidal cavity with sixth-order electrostatics. Optimization using a simple Onsager model appears to be sufficient to evaluate the important geometry changes in solution. Energies obtained from the spherical and ellipsoidal cavity models often exhibit poor convergence in the truncated electrostatic series. Correlation to experiment is much improved when an infinite-order PCM method is used.

No abstract

Atomic charges were obtained from ab initio molecular orbital calculations using a variety of procedures to compare them and assess their utility. Two procedures based on the molecular orbitals were examined, the Mulliken population analysis and the Weinhold-Reed Natural Population Analysis. Two procedures using the charge density distribution were included the Hirshfeld procedure and Bader's Atoms in Molecules method. Charges also were derived by fitting the electrostatic potential (CHELPG) and making use of the atomic polar tensors (GAIT). The procedures were f i s t examined for basis set independence, and then applied to a group of hydrocarbons. The dipole moments for these molecules were computed from the various atomic charges and compared to the total SCF dipole moments. This was followed by an examination of a series of substituted methanes, simple hydrides, and a group of typical organic compounds such as carbonyl derivatives, nitriles, and nitro compounds. In some cases, the ability of the charges to reproduce electrostatic potentials was examined.