Premium
A reliable and efficient first principles‐based method for predicting p K a values. III. Adding explicit water molecules: Can the theoretical slope be reproduced and p K a values predicted more accurately?
Author(s) -
Zhang Shuming
Publication year - 2011
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/jcc.22886
Subject(s) - solvation , chemistry , molecule , aqueous solution , thermodynamics , deprotonation , ion , implicit solvation , computational chemistry , physics , organic chemistry
Abstract A popular method for predicting p K a values for organic molecules in aqueous solution is to establish empirical linear least‐squares fits between calculated deprotonation energies and known experimental p K a values. In virtually all such calculations, the empirically observed slope of the p K a vs. Δ E fit is significantly less than the theoretical value, 1/(2.303 RT ) (which is 0.73 mol/kcal at room temperature). In our own continuum solvation calculations (Zhang et al., J Phys Chem A 2010, 114, 432), the empirical slope for carboxylic acids was only 0.23 mol/kcal, despite the excellent fit to the experimental p K a values. There has been much speculation about the reason for this phenomenon. Although the Δ E – p K a relation neglects entropic effects, these are expected to largely cancel. The most likely cause for the strange behavior of the fitted slope is explicit solute–solvent (water) interactions, especially involving the ions, which cannot be described accurately by continuum solvation models. We used our previously developed p K a protocol (OLYP/6‐311+G(d,p)//3‐21G(d) with the COSMO solvation model) to investigate the effect of adding one or two explicit water molecules to the system. The slopes for organic acids (especially carboxylic acids) are much closer to the theoretical value when explicit water molecules are added to both the neutral molecule and the anion. However, explicit water molecules have almost no effect on the slopes for organic bases. Adding explicit water molecules to the ions only produces intermediate results. Unfortunately, linear fits involving explicit water molecules have much larger errors than with continuum solvation models alone and are also much more expensive. Consequently, they are not suitable for large‐scale p K a calculations. The results compared with literature values showed that our predicted p K a s are more accurate. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011