TY - JOUR
T1 - Predicting hydration Gibbs energies of alkyl-aromatics using molecular simulation
T2 - a comparison of current force fields and the development of a new parameter set for accurate solvation data
AU - Garrido, Nuno M.
AU - Jorge, Miguel
AU - Queimada, Antonio J.
AU - Gomes, Jose R. B.
AU - Economou, Ioannis G.
AU - Macedo, Eugenia A.
PY - 2011
Y1 - 2011
N2 - The Gibbs energy of hydration is an important quantity to understand the molecular behavior in aqueous systems at constant temperature and pressure. In this work we review the performance of some popular force fields, namely TraPPE, OPLS-AA and Gromos, in reproducing the experimental Gibbs energies of hydration of several alkyl-aromatic compounds-benzene, mono-, di- and tri-substituted alkylbenzenes-using molecular simulation techniques. In the second part of the paper, we report a new model that is able to improve such hydration energy predictions, based on Lennard Jones parameters from the recent TraPPE-EH force field and atomic partial charges obtained from natural population analysis of density functional theory calculations. We apply a scaling factor determined by fitting the experimental hydration energy of only two solutes, and then present a simple rule to generate atomic partial charges for different substituted alkyl-aromatics. This rule has the added advantages of eliminating the unnecessary assumption of fixed charge on every substituted carbon atom and providing a simple guideline for extrapolating the charge assignment to any multi-substituted alkyl-aromatic molecule. The point charges derived here yield excellent predictions of experimental Gibbs energies of hydration, with an overall absolute average deviation of less than 0.6 kJ mol(-1). This new parameter set can also give good predictive performance for other thermodynamic properties and liquid structural information.
AB - The Gibbs energy of hydration is an important quantity to understand the molecular behavior in aqueous systems at constant temperature and pressure. In this work we review the performance of some popular force fields, namely TraPPE, OPLS-AA and Gromos, in reproducing the experimental Gibbs energies of hydration of several alkyl-aromatic compounds-benzene, mono-, di- and tri-substituted alkylbenzenes-using molecular simulation techniques. In the second part of the paper, we report a new model that is able to improve such hydration energy predictions, based on Lennard Jones parameters from the recent TraPPE-EH force field and atomic partial charges obtained from natural population analysis of density functional theory calculations. We apply a scaling factor determined by fitting the experimental hydration energy of only two solutes, and then present a simple rule to generate atomic partial charges for different substituted alkyl-aromatics. This rule has the added advantages of eliminating the unnecessary assumption of fixed charge on every substituted carbon atom and providing a simple guideline for extrapolating the charge assignment to any multi-substituted alkyl-aromatic molecule. The point charges derived here yield excellent predictions of experimental Gibbs energies of hydration, with an overall absolute average deviation of less than 0.6 kJ mol(-1). This new parameter set can also give good predictive performance for other thermodynamic properties and liquid structural information.
KW - hydration Gibbs energies
KW - alkyl-aromatics
KW - molecular simulation
KW - force fields
KW - parameter set
KW - accurate solvation data
UR - http://pubs.rsc.org/en/journals/journalissues/cp
U2 - 10.1039/c1cp21245a
DO - 10.1039/c1cp21245a
M3 - Article
SN - 1463-9076
VL - 13
SP - 17384
EP - 17394
JO - Physical Chemistry Chemical Physics
JF - Physical Chemistry Chemical Physics
IS - 38
ER -