Charge-dependent cavity radii for an accurate dielectric continuum model of solvation with emphasis on ions

Aqueous solutes with oxo, hydroxo, amino, methyl, chloro, bromo, and fluoro functionalities

Bojana Ginovska, Donald M. Camaioni, Michel Dupuis, Christine A. Schwerdtfeger, Quinn Gil

Research output: Contribution to journalArticle

25 Citations (Scopus)

Abstract

Dielectric continuum solvation models are widely used because they are a computationally efficacious way to simulate equilibrium properties of solutes. With advances that allow for molecular-shaped cavities, they have reached a high level of accuracy, in particular for neutral solutes. However, benchmark tests show that existing schemes for defining cavities are unable to consistently predict accurately the effects of solvation on ions, especially anions. This work involves the further development of a protocol put forth earlier for defining the cavities of aqueous solutes, with resulting advances that are most striking for anions. Molecular cavities are defined as interlocked spheres around atoms or groups of atoms in the solute, but the sphere radii are determined by simple empirically based expressions involving the effective atomic charges of the solute atoms (derived from molecular electrostatic potential) and base radii. Both of these terms are optimized for the different types of atoms or functional groups in a training set of neutral and charged solutes. Parameters in these expressions for radii were fitted by minimizing residuals between calculated and measured standard free energies of solvation (ΔG*), weighted by the uncertainty in the measured value. The calculations were performed using density functional theory with the B3LYP functional and the 6-311+G** basis set and the COnductor-like Screening MOdel (COSMO). The optimized radii definitions reproduce ΔG s* of neutral solutes and singly charged ions in the training set to within experimental uncertainty and, more importantly, accurately predict ΔGs* of compounds outside the training set, in particular anions (J. Phys. Chem. A 2003, 107, 5778). Inherent to this approach, the cavity definitions reflect the strength of specific solute-water interactions. We surmise that this feature underlies the success of the model, referred to as the CD-COSMO model for Charge-Dependent (also Camaioni-Dupuis) COSMO model. These findings offer encouragement that we can keep extending this scheme to other functional groups and obtain better accuracy in using continuum solvation models to predict equilibrium properties of aqueous ionic solutes. The approach is illustrated for a number of test cases, including the determination of acidities of an amine base, a study of the tautomerization equilibrium of a zwitterionic molecule (glycine), and calculating solvation energies of transition states toward a full characterization of reaction pathways in aqueous phase, here in SN2 exchange reactions. The calculated reaction barriers in aqueous solution are in excellent agreement with experimental values.

Original languageEnglish
Pages (from-to)10604-10613
Number of pages10
JournalJournal of Physical Chemistry A
Volume112
Issue number42
DOIs
Publication statusPublished - Oct 23 2008

Fingerprint

Solvation
solvation
solutes
Ions
continuums
cavities
radii
ions
Anions
Atoms
Screening
education
screening
conductors
Functional groups
anions
atoms
Acidity
Glycine
glycine

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

Charge-dependent cavity radii for an accurate dielectric continuum model of solvation with emphasis on ions : Aqueous solutes with oxo, hydroxo, amino, methyl, chloro, bromo, and fluoro functionalities. / Ginovska, Bojana; Camaioni, Donald M.; Dupuis, Michel; Schwerdtfeger, Christine A.; Gil, Quinn.

In: Journal of Physical Chemistry A, Vol. 112, No. 42, 23.10.2008, p. 10604-10613.

Research output: Contribution to journalArticle

@article{d089200fcb5a4f80b788d9a9211ab4fa,
title = "Charge-dependent cavity radii for an accurate dielectric continuum model of solvation with emphasis on ions: Aqueous solutes with oxo, hydroxo, amino, methyl, chloro, bromo, and fluoro functionalities",
abstract = "Dielectric continuum solvation models are widely used because they are a computationally efficacious way to simulate equilibrium properties of solutes. With advances that allow for molecular-shaped cavities, they have reached a high level of accuracy, in particular for neutral solutes. However, benchmark tests show that existing schemes for defining cavities are unable to consistently predict accurately the effects of solvation on ions, especially anions. This work involves the further development of a protocol put forth earlier for defining the cavities of aqueous solutes, with resulting advances that are most striking for anions. Molecular cavities are defined as interlocked spheres around atoms or groups of atoms in the solute, but the sphere radii are determined by simple empirically based expressions involving the effective atomic charges of the solute atoms (derived from molecular electrostatic potential) and base radii. Both of these terms are optimized for the different types of atoms or functional groups in a training set of neutral and charged solutes. Parameters in these expressions for radii were fitted by minimizing residuals between calculated and measured standard free energies of solvation (ΔG*), weighted by the uncertainty in the measured value. The calculations were performed using density functional theory with the B3LYP functional and the 6-311+G** basis set and the COnductor-like Screening MOdel (COSMO). The optimized radii definitions reproduce ΔG s* of neutral solutes and singly charged ions in the training set to within experimental uncertainty and, more importantly, accurately predict ΔGs* of compounds outside the training set, in particular anions (J. Phys. Chem. A 2003, 107, 5778). Inherent to this approach, the cavity definitions reflect the strength of specific solute-water interactions. We surmise that this feature underlies the success of the model, referred to as the CD-COSMO model for Charge-Dependent (also Camaioni-Dupuis) COSMO model. These findings offer encouragement that we can keep extending this scheme to other functional groups and obtain better accuracy in using continuum solvation models to predict equilibrium properties of aqueous ionic solutes. The approach is illustrated for a number of test cases, including the determination of acidities of an amine base, a study of the tautomerization equilibrium of a zwitterionic molecule (glycine), and calculating solvation energies of transition states toward a full characterization of reaction pathways in aqueous phase, here in SN2 exchange reactions. The calculated reaction barriers in aqueous solution are in excellent agreement with experimental values.",
author = "Bojana Ginovska and Camaioni, {Donald M.} and Michel Dupuis and Schwerdtfeger, {Christine A.} and Quinn Gil",
year = "2008",
month = "10",
day = "23",
doi = "10.1021/jp804092v",
language = "English",
volume = "112",
pages = "10604--10613",
journal = "Journal of Physical Chemistry A",
issn = "1089-5639",
publisher = "American Chemical Society",
number = "42",

}

TY - JOUR

T1 - Charge-dependent cavity radii for an accurate dielectric continuum model of solvation with emphasis on ions

T2 - Aqueous solutes with oxo, hydroxo, amino, methyl, chloro, bromo, and fluoro functionalities

AU - Ginovska, Bojana

AU - Camaioni, Donald M.

AU - Dupuis, Michel

AU - Schwerdtfeger, Christine A.

AU - Gil, Quinn

PY - 2008/10/23

Y1 - 2008/10/23

N2 - Dielectric continuum solvation models are widely used because they are a computationally efficacious way to simulate equilibrium properties of solutes. With advances that allow for molecular-shaped cavities, they have reached a high level of accuracy, in particular for neutral solutes. However, benchmark tests show that existing schemes for defining cavities are unable to consistently predict accurately the effects of solvation on ions, especially anions. This work involves the further development of a protocol put forth earlier for defining the cavities of aqueous solutes, with resulting advances that are most striking for anions. Molecular cavities are defined as interlocked spheres around atoms or groups of atoms in the solute, but the sphere radii are determined by simple empirically based expressions involving the effective atomic charges of the solute atoms (derived from molecular electrostatic potential) and base radii. Both of these terms are optimized for the different types of atoms or functional groups in a training set of neutral and charged solutes. Parameters in these expressions for radii were fitted by minimizing residuals between calculated and measured standard free energies of solvation (ΔG*), weighted by the uncertainty in the measured value. The calculations were performed using density functional theory with the B3LYP functional and the 6-311+G** basis set and the COnductor-like Screening MOdel (COSMO). The optimized radii definitions reproduce ΔG s* of neutral solutes and singly charged ions in the training set to within experimental uncertainty and, more importantly, accurately predict ΔGs* of compounds outside the training set, in particular anions (J. Phys. Chem. A 2003, 107, 5778). Inherent to this approach, the cavity definitions reflect the strength of specific solute-water interactions. We surmise that this feature underlies the success of the model, referred to as the CD-COSMO model for Charge-Dependent (also Camaioni-Dupuis) COSMO model. These findings offer encouragement that we can keep extending this scheme to other functional groups and obtain better accuracy in using continuum solvation models to predict equilibrium properties of aqueous ionic solutes. The approach is illustrated for a number of test cases, including the determination of acidities of an amine base, a study of the tautomerization equilibrium of a zwitterionic molecule (glycine), and calculating solvation energies of transition states toward a full characterization of reaction pathways in aqueous phase, here in SN2 exchange reactions. The calculated reaction barriers in aqueous solution are in excellent agreement with experimental values.

AB - Dielectric continuum solvation models are widely used because they are a computationally efficacious way to simulate equilibrium properties of solutes. With advances that allow for molecular-shaped cavities, they have reached a high level of accuracy, in particular for neutral solutes. However, benchmark tests show that existing schemes for defining cavities are unable to consistently predict accurately the effects of solvation on ions, especially anions. This work involves the further development of a protocol put forth earlier for defining the cavities of aqueous solutes, with resulting advances that are most striking for anions. Molecular cavities are defined as interlocked spheres around atoms or groups of atoms in the solute, but the sphere radii are determined by simple empirically based expressions involving the effective atomic charges of the solute atoms (derived from molecular electrostatic potential) and base radii. Both of these terms are optimized for the different types of atoms or functional groups in a training set of neutral and charged solutes. Parameters in these expressions for radii were fitted by minimizing residuals between calculated and measured standard free energies of solvation (ΔG*), weighted by the uncertainty in the measured value. The calculations were performed using density functional theory with the B3LYP functional and the 6-311+G** basis set and the COnductor-like Screening MOdel (COSMO). The optimized radii definitions reproduce ΔG s* of neutral solutes and singly charged ions in the training set to within experimental uncertainty and, more importantly, accurately predict ΔGs* of compounds outside the training set, in particular anions (J. Phys. Chem. A 2003, 107, 5778). Inherent to this approach, the cavity definitions reflect the strength of specific solute-water interactions. We surmise that this feature underlies the success of the model, referred to as the CD-COSMO model for Charge-Dependent (also Camaioni-Dupuis) COSMO model. These findings offer encouragement that we can keep extending this scheme to other functional groups and obtain better accuracy in using continuum solvation models to predict equilibrium properties of aqueous ionic solutes. The approach is illustrated for a number of test cases, including the determination of acidities of an amine base, a study of the tautomerization equilibrium of a zwitterionic molecule (glycine), and calculating solvation energies of transition states toward a full characterization of reaction pathways in aqueous phase, here in SN2 exchange reactions. The calculated reaction barriers in aqueous solution are in excellent agreement with experimental values.

UR - http://www.scopus.com/inward/record.url?scp=55749093054&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=55749093054&partnerID=8YFLogxK

U2 - 10.1021/jp804092v

DO - 10.1021/jp804092v

M3 - Article

VL - 112

SP - 10604

EP - 10613

JO - Journal of Physical Chemistry A

JF - Journal of Physical Chemistry A

SN - 1089-5639

IS - 42

ER -