Together these data are consistent with the characteristic overpolarization estimates and weighted residuals we obtained with the linear regression model fits described previously. In general, HF/6-31G* gas phase calculations do not overpolarize molecular charge distributions as strongly as a more sophisticated method using PCM benzene, a nonpolar solvent of low dielectric constant. The results for the HF/aug-cc-pV(T+d)Z gas phase calculations demonstrate that there is an inherent overpolarization due to the HF method, whereas the inconsistency in the strength of the overpolarization across different molecules is due to the 6-31G* basis set.