A fixed-charge model for alcohol polarization in the condensed phase, and its role in small molecule hydration

We present a simple optimization strategy for incorporating experimental dielectric response information on neat liquids in classical molecular models of alcohol. Using this strategy, we determine simple and transferable hydroxyl modulation rules that, when applied to an existing molecular parameter set, result in a newly dielectric corrected (DC) parameter set. We applied these rules to the general Amber force field (GAFF) to form an initial set of GAFF-DC parameters, and we found this to lead to significant improvement in the calculated dielectric constant and hydration free energy values for a wide variety of small molecule alcohol models. Tests of the GAFF-DC parameters in the SAMPL4 blind prediction event for hydration show these changes improve agreement with experiment. Surprisingly, these simple modifications also outperform detailed quantum mechanical electric field calculations using a self-consistent reaction field environment coupling term. This work provides a potential benchmark for future developments in methods for representing condensed-phase environments in electronic structure calculations.