Abstract
Thole-style mutual induction models for molecular polarization have been adopted by several popular polarizable force fields (FFs) for their simplicity and transferability. The atomic polarizability parameters of these models are typically derived by fitting to ab initio or/and experimental molecular polarizabilities. In this work, we improve upon Thole polarizability parameters by employing both high-level quantum mechanics molecular polarizabilities and electrostatic potential (ESP) responses on three-dimensional grids. Our results indicate that the two approaches to derive atomic polarizability parameters are both effective, while the ESP approaches can also capture the polarization for the atoms with lone pair electrons. The resulting polarizability parameters have been validated on a set of over 7200 molecules covering the most common elements found in organic molecules (C, H, O, N, P, S, F, Cl, Br, and I). These parameters have also been tested on the experimentally measured molecular polarizabilities of 422 molecules. The final set of parameters derived in this work show notable improvement over the current AMOEBA set. The result is a highly transferable, expanded set of atomic polarizabilities defined by the local chemical environment in the form of SMARTS patterns. These parameters can be used directly in molecular mechanics polarizable potential energy functions such as AMOEBA, AMOEBA+, and other Thole-style models.