skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Parameterizing Potential-Derived Charge-Dependent Cavity Radii for Continuum Solvation Models: Aqueous Solutes with Oxo, Hydroxo, Amino, and Methyl Functionalities

Abstract

In recent years, computer processing capabilities have improved, allowing computationally demanding solvation models to become practically implemented. Since the determination of experimental solvation free energies for short-lived intermediates such as radicals may be difficult, expensive, and time consuming, developing computational solvation models has become an effective alternative approach. In this work, the COSMO (COnductor like Screening MOdel) continuum solvation model was used to predict solvation free energies of a set of solutes with oxo, hydroxyl, amino, and methyl functional groups in the solvent water (ε=78.39). The continuum solvation model requires a cavity to be defined around a solute to predict solvation free energy (ΔGs*). To reproduce experimental values, the size and shape of the cavity must reflect the strength of solute/solvent interactions. The cavity is defined as interlocking spheres around atoms or groups of atoms in the solute. The sphere radii are defined by linear functions parameterized in terms of potential derived charges known as CHELPG charges. A training set of neutral and ionic solutes was created, and coefficients in the radii definitions were fitted to reproduce experimental ΔGs* values by minimizing residuals between experimental and calculated ΔGs* values for compounds in the training set. Data used in fitting themore » radii definitions’ coefficients was generated by performing electronic structure and solvation computations, using varying cavity sizes. The calculations were done using Density Functional Theory with the B3LYP functional and 6-311+G** basis set in the Gaussian98 package of programs. Radii definitions reproduce ΔGs* of neutrals and singly-charged ions in training set to within experimental uncertainty and accurately predict ΔGs* of some compounds outside the training set. The mean unsigned error of this training set is 0.18 kcal/mol. These findings suggest that the protocol described here for developing potential derived charge dependent cavity definitions may successfully be extended to more functional groups to increase the applicability of the scheme and obtain greater accuracy in using continuum solvation models in predicting equilibrium properties of aqueous solutes.« less

Authors:
;
Publication Date:
Research Org.:
Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
993378
Report Number(s):
PNNL-SA-46271
KL0101000; TRN: US201023%%221
DOE Contract Number:  
AC05-76RL01830
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Undergraduate Research, VI:147; Journal Volume: 6
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; SOLVATION; MATHEMATICAL MODELS; FREE ENERGY; C CODES; HYDROXY COMPOUNDS; AMINES; OXYGEN COMPOUNDS; METHYL RADICALS; DENSITY FUNCTIONAL METHOD; WATER; Continuum; solvation; models

Citation Formats

Schwerdtfeger, Christine A, and Camaioni, Donald M. Parameterizing Potential-Derived Charge-Dependent Cavity Radii for Continuum Solvation Models: Aqueous Solutes with Oxo, Hydroxo, Amino, and Methyl Functionalities. United States: N. p., 2006. Web.
Schwerdtfeger, Christine A, & Camaioni, Donald M. Parameterizing Potential-Derived Charge-Dependent Cavity Radii for Continuum Solvation Models: Aqueous Solutes with Oxo, Hydroxo, Amino, and Methyl Functionalities. United States.
Schwerdtfeger, Christine A, and Camaioni, Donald M. Tue . "Parameterizing Potential-Derived Charge-Dependent Cavity Radii for Continuum Solvation Models: Aqueous Solutes with Oxo, Hydroxo, Amino, and Methyl Functionalities". United States. doi:.
@article{osti_993378,
title = {Parameterizing Potential-Derived Charge-Dependent Cavity Radii for Continuum Solvation Models: Aqueous Solutes with Oxo, Hydroxo, Amino, and Methyl Functionalities},
author = {Schwerdtfeger, Christine A and Camaioni, Donald M},
abstractNote = {In recent years, computer processing capabilities have improved, allowing computationally demanding solvation models to become practically implemented. Since the determination of experimental solvation free energies for short-lived intermediates such as radicals may be difficult, expensive, and time consuming, developing computational solvation models has become an effective alternative approach. In this work, the COSMO (COnductor like Screening MOdel) continuum solvation model was used to predict solvation free energies of a set of solutes with oxo, hydroxyl, amino, and methyl functional groups in the solvent water (ε=78.39). The continuum solvation model requires a cavity to be defined around a solute to predict solvation free energy (ΔGs*). To reproduce experimental values, the size and shape of the cavity must reflect the strength of solute/solvent interactions. The cavity is defined as interlocking spheres around atoms or groups of atoms in the solute. The sphere radii are defined by linear functions parameterized in terms of potential derived charges known as CHELPG charges. A training set of neutral and ionic solutes was created, and coefficients in the radii definitions were fitted to reproduce experimental ΔGs* values by minimizing residuals between experimental and calculated ΔGs* values for compounds in the training set. Data used in fitting the radii definitions’ coefficients was generated by performing electronic structure and solvation computations, using varying cavity sizes. The calculations were done using Density Functional Theory with the B3LYP functional and 6-311+G** basis set in the Gaussian98 package of programs. Radii definitions reproduce ΔGs* of neutrals and singly-charged ions in training set to within experimental uncertainty and accurately predict ΔGs* of some compounds outside the training set. The mean unsigned error of this training set is 0.18 kcal/mol. These findings suggest that the protocol described here for developing potential derived charge dependent cavity definitions may successfully be extended to more functional groups to increase the applicability of the scheme and obtain greater accuracy in using continuum solvation models in predicting equilibrium properties of aqueous solutes.},
doi = {},
journal = {Journal of Undergraduate Research, VI:147},
number = ,
volume = 6,
place = {United States},
year = {Tue Jan 31 00:00:00 EST 2006},
month = {Tue Jan 31 00:00:00 EST 2006}
}