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.

Evaluation of Three Heat of Formation Calculation Methods for High-throughput Energetic Molecule Design and Screening

A Method for Calculating the Heats of Formation of Medium-Sized and Large-Sized Molecules

A Comparison of the Accuracy of Semi-empirical PM3, PDDG and PM6 Methods in Predicting Heats of Formation for Organic Compounds

Heats of Formation for Energetic Compounds Calculated Using Atomization Reactions

Semi-empirical MO Investigation on Heats of Formation for Energetic Compounds

High-Throughput Predictions of Accurate Enthalpies of Formation for Larger Molecules Utilizing the Bond Difference Correction Method.

High-Throughput Design of Energetic Molecules

Accurate Prediction of Heat of Formation by Combining Hartree-Fock/Density Functional Theory Calculation with Linear Regression Correction Approach

The comparative study for the heat of formation of organic molecules calculated by different methods

First principles reactive simulation for equation of state prediction

A Simple Method for the Prediction of the Detonation Performances of Metal-Containing Explosives

Heats of formation of medium-sized organic compounds from contemporary electronic structure methods

The X1 Family of Methods That Combines B3LYP with Neural Network Corrections for an Accurate Yet Efficient Prediction of Thermochemistry

Method and automatic program for accurate thermodynamic data of reaction mechanisms for combustion modeling

High-precision Standard Enthalpy of Formation for Polycyclic Aromatic Hydrocarbons Predicting from General Connectivity Based Hierarchy with Discrete Correction of Atomization Energy

The X1 Method For Accurate And Efficient Prediction Of Heats Of Formation

First-Principles Performance Prediction of High Explosives Enabled by Machine Learning

Corrections of Molecular Morphology and Hydrogen Bond for Improved Crystal Density Prediction

A Comparison of High-Level Theoretical Methods to Predict the Heats of Formation of Azo Compounds

Assessment of hybrid, meta-hybrid-GGA, and long-range corrected density functionals for the estimation of enthalpies of formation, barrier heights, and ionisation potentials of selected C1–C5 oxygenates

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