Resolution of a Challenge for Solvation Modeling: Calculation of Dicarboxylic Acid Dissociation Constants Using Mixed Discrete-Continuum Solvation Models
First and second dissociation constants (pKa values) of oxalic acid, malonic acid, and adipic acid were computed by using a number of theoretical protocols based on density functional theory and using both continuum solvation models and mixed discrete-continuum solvation models. We show that fully implicit solvation models (in which the entire solvent is represented by a dielectric continuum) fail badly for dicarboxylic acids with mean unsigned errors averaged over six pKa values) of 2.4-9.0 log units, depending on the particular implicit model used. The use of water-solute clusters and accounting for multiple conformations in solution significantly improve the performance of both generalized Born solvation models and models that solve the nonhomogeneous dielectric Poisson equation for bulk electrostatics. The four most successful models have mean unsigned errors of only 0.6-0.8 log units.
- Research Organization:
- Pacific Northwest National Laboratory (PNNL), Richland, WA (US), Environmental Molecular Sciences Laboratory (EMSL)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC05-76RL01830
- OSTI ID:
- 1120103
- Journal Information:
- Journal of Physical Chemistry Letters, Journal Name: Journal of Physical Chemistry Letters Journal Issue: 11 Vol. 3; ISSN 1948-7185
- Publisher:
- American Chemical Society
- Country of Publication:
- United States
- Language:
- English
Similar Records
Calculating Free Energy Changes in Continuum Solvation Models
Thermodynamics of complexation of lanthanides by dicarboxylate ligands