Review of force fields and intermolecular potentials used in atomistic computational materials research

Molecular simulation is a powerful computational tool for a broad range of applications including the examination of materials properties and accelerating drug discovery. At the heart of molecular simulation is the analytic potential energy function. These functions span the range of complexity from very simple functions used to model generic phenomena to complex functions designed to model chemical reactions. The complexity of the mathematical function impacts the computational speed and is typically linked to the accuracy of the results obtained from simulations that utilize the function. One approach to improving accuracy is to simply add more parameters and additional complexity to the analytic function. This approach is typically used in non-reactive force fields where the functional form is not derived from quantum mechanical principles. The form of other types of potentials, such as the bond-order potentials, is based on quantum mechanics and has led to varying levels of accuracy and transferability. When selecting a potential energy function for use in molecular simulations, the accuracy, transferability, and computational speed must all be considered. In this focused review, some of the more commonly used potential energy functions for molecular simulations are reviewed with an eye toward presenting their general forms, strengths, and weaknesses.