The results of the sixth blind test of organic crystal structure prediction methods are presented and discussed, highlighting progress for salts, hydrates and bulky flexible molecules, as well as on-going challenges.

A new induced dipole polarization model based on interacting Gaussian charge densities is presented. In contrast to the original induced point dipole model, the Gaussian polarization model is capable of finite interactions at short distances. Aspects of convergence related to the Gaussian model will be explored. The Gaussian polarization model is compared with the damped Thole-induced dipole model and the point dipole model. It will be shown that the Gaussian polarization model performs slightly better than the Thole model in terms of fitting to molecular polarizability tensors. An advantage of the model based on Gaussian charge distribution is that it can be easily generalized to other multipole moments and provide effective damping for both permanent electrostatic and polarization models. Finally, a method of parameterizing polarizabilities is presented. This method is based on probing a molecule with point charges and fitting polarizabilities to electrostatic potential. In contrast to the generic atom type polarizabilities fit to molecular polarizability tensors, probed polarizabilities are significantly more accurate in terms of reproducing molecular polarizability tensors and electrostatic potential, while retaining conformational transferability.

An electrostatic model based on charge density is proposed as a model for future force fields. The model is composed of a nucleus and a single Slater-type contracted Gaussian multipole charge density on each atom. The Gaussian multipoles are fit to the electrostatic potential (ESP) calculated at the B3LYP/6-31G* and HF/aug-cc-pVTZ levels of theory and tested by comparing electrostatic dimer energies, inter-molecular density overlap integrals, and permanent molecular multipole moments with their respective ab initio values. For the case of water, the atomic Gaussian multipole moments Q lm are shown to be a smooth function of internal geometry (bond length and bond angle), which can be approximated by a truncated linear Taylor series. In addition, results are given when the Gaussian multipole charge density is applied to a model for exchange-repulsion energy based on the inter-molecular density overlap.