A Simple Method for Including Polarization Effects in Solvation Free Energy Calculations When Using Fixed-Charge Force Fields: Alchemically Polarized Charges.

11 June 2020, Version 2
This content is a preprint and has not undergone peer review at the time of posting.


The incorporation of polarizability in classical force-field molecular simulations is an ongoing area of research. We focus here on its application to hydration free energy simulations of organic molecules. In contrast to computationally complex approaches involving the development of explicitly polarizable force fields, we present herein a simple methodology for incorporating polarization into such simulations using standard fixed-charge force-fields, which we call the Alchemically Polarized Charges (APolQ) method. APolQ employs a standard classical alchemical free energy change simulation to calculate the free energy difference between a fully polarized solute particle in a condensed phase and its unpolarized state in a vacuum. One electronic structure (ES) calculation to of the electron densities is required for each state: for the former, we use a Polarizable Continuum Model (PCM), and for the latter we use vacuum-phase electronic structure calculations.

We applied APolQ to hydration free energy data for a test set of 45 neutral solute molecules in the FreeSolv database, and compared results obtained using three different water models (SPC/E, TIP3P, OPC3) and using MBIS and RESP partial charge methodologies. ES calculations were carried out at the MP2 level of theory and with cc-pVTZ and aug-cc-pVTZ basis sets. In comparison with AM1-BCC, we found that APolQ outperforms it for the test set. Despite our method using default GAFF parameters, the MBIS partial charges yield Absolute Average Deviations (AAD) 1.5 to 1.9 kJ·mol−1 lower than AM1-BCC.

We conjecture that this method can be further improved by fitting the Lennard-Jones and torsional parameters to partial charges derived using MBIS or RESP methodologies.


chemical potential
Free energy calculations
Partial Charge Models
Hydration Free Energy Calculations
Molecular Simulation Studies
molecular dynamics

Supplementary materials

Supporting Information APolQ

Supplementary weblinks


Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.