Adding Explicit Solvent Molecules to Continuum Solvent Calculations for the Calculation of Aqueous Acid Dissociation Constants

• Published in: The journal of physical chemistry. A, Molecules, spectroscopy, kinetics, environment, & general theory Molecules, spectroscopy, kinetics, environment, & general theory, 2006-02, Vol.110 (7), p.2493-2499
• Main Authors: Kelly, Casey P, Cramer, Christopher J, Truhlar, Donald G
• Format: Article
• Language: English
• Subjects: Acids - chemistry; Algorithms; Anions; Bicarbonates - chemistry; Hydrogen Bonding; Solvents - chemistry; Thermodynamics; Water - chemistry
• ISSN: 1089-5639; 1520-5215
• DOI: 10.1021/jp055336f

Aqueous acid dissociation free energies for a diverse set of 57 monoprotic acids have been calculated using a combination of experimental and calculated gas and liquid-phase free energies. For ionic species, aqueous solvation free energies were calculated using the recently developed SM6 continuum solvation model (Kelly, C. P.; Cramer, C. J.; Truhlar, D. G. J. Chem. Theory Comput. 2005, 1, 1133). This model combines a dielectric continuum with atomic surface tensions to account for bulk solvent effects. For some of the acids studied, a combined approach that involves attaching a single explicit water molecule to the conjugate base (anion), and then surrounding the resulting anion−water cluster by a dielectric continuum, significantly improves the agreement between the calculated pK a value and experiment. This suggests that for some anions, particularly those concentrating charge on a single exposed heteroatom, augmenting implicit solvent calculations with a single explicit water molecule is required, and adequate, to account for strong short-range hydrogen bonding interactions between the anion and the solvent. We also demonstrate the effect of adding several explicit waters by calculating the pK a of bicarbonate (HCO3 -) using as the conjugate base carbonate (CO3 2-) bound by up to three explicit water molecules.