Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published August 14, 2008 | Supplemental Material
Journal Article Open

Calculation of Solvation Free Energies of Charged Solutes Using Mixed Cluster/Continuum Models

Abstract

We derive a consistent approach for predicting the solvation free energies of charged solutes in the presence of implicit and explicit solvents. We find that some published methodologies make systematic errors in the computed free energies because of the incorrect accounting of the standard state corrections for water molecules or water clusters present in the thermodynamic cycle. This problem can be avoided by using the same standard state for each species involved in the reaction under consideration. We analyze two different thermodynamic cycles for calculating the solvation free energies of ionic solutes: (1) the cluster cycle with an n water cluster as a reagent and (2) the monomer cycle with n distinct water molecules as reagents. The use of the cluster cycle gives solvation free energies that are in excellent agreement with the experimental values obtained from studies of ion−water clusters. The mean absolute errors are 0.8 kcal/mol for H^+ and 2.0 kcal/mol for Cu2^+. Conversely, calculations using the monomer cycle lead to mean absolute errors that are >10 kcal/mol for H^+ and >30 kcal/mol for Cu2^+. The presence of hydrogen-bonded clusters of similar size on the left- and right-hand sides of the reaction cycle results in the cancelation of the systematic errors in the calculated free energies. Using the cluster cycle with 1 solvation shell leads to errors of 5 kcal/mol for H^+ (6 waters) and 27 kcal/mol for Cu2^+ (6 waters), whereas using 2 solvation shells leads to accuracies of 2 kcal/mol for Cu2^+ (18 waters) and 1 kcal/mol for H^+ (10 waters).

Additional Information

© 2008 American Chemical Society. Received 9 October 2007. Published online 23 July 2008. Published in print 1 August 2008. Funding for this work was provided by the National Science Foundation (NIRT CTS award no. 0506951) and by the U.S. Environmental Protection Agency (STAR grant RD-83252501). The computational facilities used in these studies were funded by grants from ARO-DURIP, ONR-DURIP, and NSF-MRI.

Attached Files

Supplemental Material - jp802665d-file003.txt

Files

jp802665d-file003.txt
Files (94.5 kB)
Name Size Download all
md5:435c50f6f4eb15aa409bada755929111
94.5 kB Preview Download

Additional details

Created:
August 19, 2023
Modified:
October 26, 2023