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.

Molecular Simulation of Caloric Properties of Fluids Modelled by Force Fields with Intramolecular Contributions: Application to Heat Capacities

Contribution of the Hydroxyl Group on Interfacial Heat Conduction of Monohydric Alcohols: A Molecular Dynamics Study

Force and Heat Current Formulas for Many-Body Potentials in Molecular Dynamics Simulations with Applications to Thermal Conductivity Calculations

Heat Capacities of Fluids: the Performance of Various Equations of State

Modeling Heat Capacity of Supercritical Fluids: A Phonon Theory Approach

MOLECULAR SIMULATION OF SOME THERMOPHYSICAL PROPERTIES AND PHENOMENA

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

Atomistic simulation of barocaloric effects

A general heat capacity model of simple liquids using phonon theory

Phonon gas model for thermal conductivity of dense, strongly interacting liquids

Modeling the Viscosity of Hydrofluorocarbons, Hydrofluoroolefins and Their Binary Mixtures Using Residual Entropy Scaling and Cubic-Plus-association Equation of State

Thermodynamic response functions of fluids: a microscopic approach based on NpT Monte Carlo

Non-bonded force field parameters from MBIS partitioning of the molecular electron density improve thermophysical properties prediction of organic liquids

Accurate Thermodynamic Properties of the Neon Fluid from Ab Initio Potentials and Uncertainty Assessments of Simulation Data

A New Model for Specific Heat Capacity of Real Gas

Modeling heat capacity of liquid hydrofluorocarbons and hydrofluoroolefins

Molecular Thermodynamic Model for Fluids Containing Associated Molecules

Thermal conductivity calculations of binary liquid organic mixtures by molecular dynamics simulation and its interpretation of microscopic heat transfer mechanism

Round Robin Study: Molecular Simulation of Thermodynamic Properties from Models with Internal Degrees of Freedom

Non-destructive molecular-dynamics simulation of the chemical potential of a fluid

Transport properties of carbon dioxide and methane from molecular dynamics simulations