Further refinements of next-generation force fields — Nonempirical localization of off-centered points in molecules

Abstract

In the present work, we investigate different possibilities for the nonempirical localization of nonatomic centers within the context of the design of new generation polarizable force fields. To do so, the positions of electron localization function (ELF) and electron pair localization function (EPLF) attractors and of Boys orbital centroids are determined for a set of thirteen saturated and conjugated biologically relevant molecules. We consider the similarities and differences in the representations of localized lone pairs and π delocalized systems by these approaches, as well as the effects of the basis sets and of the level of quantum chemistry (QC). All three methods give consistent results upon dealing with the localized lone pairs. Concerning aromatic systems, whereas ELF and EPLF methods give mutually consistent results, the Boys scheme breaks the symmetry by alternating the electron distributions along the C–C bonds providing a different distribution of off-centered points. We then investigate the influence of lone pair localization in the water model of the Gaussian electrostatic model/sum of interactions between fragments ab initio computed (GEM/SIBFA) polarizable force field, which embodies an explicit representation of the lone pairs. This is done, in a series of mono- and poly-hydrated Zn(II) complexes, by comparing the overlap-dependent repulsion (Erep) and charge-transfer (Ect) contributions to their QC counterparts. We resort to either the current GEM/SIBFA water lone pair internal coordinates or to the ones derived from the previously mentioned localization procedures. The latter enables closer reproductions by Erep and Ect (GEM/SIBFA) of their exchange repulsion (Eexch) and Ect QC counterparts. The present preliminary results show that QC localization procedures can be used to derive accurate, nonempirical positions for off-centered points intervening in the formulation of the overlap-dependent contributions of next-generation polarizable force fields.