Abstract:Molecular mechanical force field (FF) methods are useful in studying condensed phase properties. They are complementary to experiments and can often go beyond experiments in atomic details. Even if a FF is specific for studying structures, dynamics, and functions of biomolecules, it is still important for the FF to accurately reproduce the experimental liquid properties of small molecules that represent the chemical moieties of biomolecules. Otherwise, the force field may not describe the structures and energies of macromolecules in aqueous solutions properly. In this work, we have carried out a systematic study to evaluate the General AMBER Force Field (GAFF) in studying densities and heats of vaporization for a large set of organic molecules that covers the most common chemical functional groups. The latest techniques, such as the particle mesh Ewald (PME) for calculating electrostatic energies and Langevin dynamics for scaling temperatures, have been applied in the molecular dynamics (MD) simulations. For density, the average percent error (APE) of 71 organic compounds is 4.43% when compared to the experimental values. More encouragingly, the APE drops to 3.43% after the exclusion of two outliers and four other compounds for which the experimental densities have been measured with pressures higher than 1.0 atm. For the heat of vaporization, several protocols have been investigated, and the best one, P4/ntt0, achieves an average unsigned error (AUE) and a root-mean-square error (RMSE) of 0.93 and 1.20 kcal/mol, respectively. How to reduce the prediction errors through proper van der Waals (vdW) parametrization has been discussed. An encouraging finding in vdW parametrization is that both densities and heats of vaporization approach their "ideal" values in a synchronous fashion when vdW parameters are tuned. The following hydration free energy calculation using thermodynamic integration further justifies the vdW refinement. We conclude that simple vdW parametrization can significantly reduce the prediction errors. We believe that GAFF can greatly improve its performance in predicting liquid properties of organic molecules after a systematic vdW parametrization, which will be reported in a separate paper.

Rapid Prediction of Solvation Free Energy and Vapor Pressure of Liquid and Solid from Molecular Dynamics Simulation

A Novel Multiscale Approach for Rapid Prediction of Phase Behaviors with Consideration of Molecular Conformations

A molecular dynamics method for calculating molecular volume changes appropriate for biomolecular simulation

An Efficient Method for Computing Excess Free Energy of Liquid

A generalized approach for rapid entropy calculation of liquids and solids

Fast Computation of Solvation Free Energies with Molecular Density Functional Theory: Thermodynamic-Ensemble Partial Molar Volume Corrections

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman-Rubin Convergence Diagnostics

A High-Throughput Computing Procedure for Predicting Vapor-Liquid Equilibria of Binary Mixtures – Using Carbon Dioxide and N-Alkanes As Examples

A Novel Method for Fluid Pour-Point Prediction by Molecular Dynamics Simulations

High-Throughput Prediction of the Hydration Free Energies of Small Molecules from A Classical Density Functional Theory

Comparison of Two Simulation Methods to Compute Solvation Free Energies and Partition Coefficients

Liquid-Vapor Coexistence and Spontaneous Evaporation at Atmospheric Pressure of Common Rigid Three-Point Water Models in Molecular Simulations

Application of Molecular Dynamics Simulations in Molecular Property Prediction II: Diffusion Coefficient.

First Principles Free-Energy Theory of Solvation with Atomic Scale Liquid Structure

Prediction of thermodynamic properties of fluid by molecular simulation method

Correction to Application of Molecular Dynamics Simulations in Molecular Property Prediction. 1. Density and Heat of Vaporization.

Implementation of Solvation Free Energy Framework to Predict the Vapor-Liquid Equilibrium Behaviors for the Water–hydrazine and Ethanol–water Azeotropic Systems

Combining Molecular Dynamics and Machine Learning to Predict Self-Solvation Free Energies and Limiting Activity Coefficients

Solvation Free Energy Calculations with Quantum Mechanics/Molecular Mechanics and Machine Learning Models.

Molecular Dynamics Investigation on the Vapor–liquid Interface Behavior of Long-Chain Alkanes, Alcohols, and Their Mixtures