title: Force Field Benchmark of Organic Liquids. 2. Gibbs Energy of Solvation.
abstract: Quantitative prediction of physical properties of liquids is a longstanding goal of molecular simulation. Here, we evaluate the predictive power of the Generalized Amber Force Field (Wang et al. J. Comput. Chem. 2004, 25, 1157-1174) for the Gibbs energy of solvation of organic molecules in organic solvents using the thermodynamics integration (TI) method. The results are compared to experimental data, to a model based on quantitative structure property relations (QSPR), and to the conductor-like screening models for realistic solvation (COSMO-RS) model. Although the TI calculations yield slightly better correlation to experimental results than the other models, in all fairness we should conclude that the difference between the models is minor since both QSPR and COSMO-RS yield a slightly lower RMSD from that of the experiment (<3.5 kJ/mol). By analyzing which molecules (either as solvents or solutes) are outliers in the TI calculations, we can pinpoint where additional parametrization efforts are needed. For the force field based TI calculations, deviations from the experiment occur in particular when compounds containing nitro or ester groups are solvated into other liquids, suggesting that the interaction between these groups and solvents may be too strong. In the COSMO-RS calculations, outliers mainly occur when compounds containing (in particular aromatic) rings are solvated despite using a ring correction term in the calculations.
authors.affiliation: †Department of Chemistry, Zhejiang University, Hangzhou 310027, China.|‡Uppsala Center for Computational Chemistry, Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden.
authors.forename: Jin
authors.initials: J
authors.lastname: Zhang
authors.affiliation: ‡Uppsala Center for Computational Chemistry, Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden.|§Department of Biology, School of Arts and Sciences, National University of Mongolia, Ulaanbaatar 14200, Mongolia.
authors.forename: Badamkhatan
authors.initials: B
authors.lastname: Tuguldur
authors.affiliation: ‡Uppsala Center for Computational Chemistry, Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden.
authors.forename: David
authors.initials: D
authors.lastname: van der Spoel
mesh_terms: D019985:Benchmarking; D048088:Informatics / Q000379:methods*; D008958:Models, Molecular; D008968:Molecular Conformation; D009930:Organic Chemicals / Q000737:chemistry*; D021281:Quantitative Structure-Activity Relationship; D012997:Solvents / Q000737:chemistry*; D013816:Thermodynamics
chemical_list: D009930:Organic Chemicals; D012997:Solvents