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

The accuracy of predicting thermophysical properties through molecular dynamics simulations is constrained by the precision of models used to describe molecular interactions. The Open Force Field initiative has established a computational structure to develop new models and introduced two non-polarizable force fields, Parsley and Sage. Sage version 2.0.0 focused on refining Lennard-Jones parameters to reflect thermophysical properties accurately. In this context, we evaluate the ability of our introduced D-MBIS non-bonded force field parameters to replicate liquid densities and enthalpies of evaporation of 49 neutral compounds from the ThermoML database using the openff-evaluator package. Our findings confirm that our ab-initio derived non-bonded force field parameters accurately mirror both thermophysical properties with a high degree of precision.