Premium
Incorporating solvation effects into density functional theory: Calculation of absolute acidities
Author(s) -
Richardson W. H.,
Peng C.,
Bashford D.,
Noodleman L.,
Case D. A.
Publication year - 1997
Publication title -
international journal of quantum chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.484
H-Index - 105
eISSN - 1097-461X
pISSN - 0020-7608
DOI - 10.1002/(sici)1097-461x(1997)61:2<207::aid-qua3>3.0.co;2-#
Subject(s) - solvation , density functional theory , computational chemistry , chemistry , molecule , organic chemistry
An approach to the calculation of molecular electronic structures, solvation energies, and pK a values in condensed phases is described. The electronic structure of the solute is described by density functional quantum mechanics, and electrostatic features of environmental effects are modeled through external charge distributions and continuum dielectrics. The reaction potential produced by a mode of the molecular charge distribution is computed via finite‐difference solutions to the Poisson‐Boltzmann equation and incorporated into the self‐consistent field procedure. Here we report results on three sets of organic acids, whose p K a values range over 16 pH units. The first set provides models for ionizable side chains in proteins; the second set considers the effects of substituting one to three chlorine atoms for hydrogens in acetic acid; and the final set consists of 4‐substituted‐bicyclo‐[2.2.2]‐octanecarboxylic acids. Successful prediction of “absolute” p K a values places stringent requirements on the computation of gas‐phase proton affinities and on the response to solvation. In some cases the current model shows substantial errors, but overall the results and trends are in good agreement with experiment. Prospects for extending this approach to more complex systems such as proteins are briefly discussed. © 1997 John Wiley & Sons, Inc.