Accurate molecular interaction potentials, usually referred to as force fields, are an essential part of modern molecular dynamics workflows. For common applications such as simulations of proteins and computer aided drug design (CADD), the simple fixed-charge force field functional form is generally used. This formulation splits the potential energy of molecules into discrete components, with separate energy terms for each component . Broadly, these can be divided into the bonded (or valence) components, which comprise the energies of bond lengths, bond angles, and torsional angles, and the non-bonded components, representing short-range dispersion repulsion interactions and longer-range Coulombic interactions. This type of force field has been successful in many applications because of its simplicity, interpretability, and computational efficiency. Many studies have used these force fields to probe the mechanisms of protein dynamics , and they have become widely adopted in the pharmaceutical industry as a means of screening drug candidate molecules in silico . While these force fields are quite simple in their functional form, their accuracy is dependent on hundreds to thousands of parameters, which dictate the strength of interactions in different chemical environments. Decades of effort from the computational chemistry community has produced many different parameter sets to cover a wide range of chemistries , largely by fitting parameters to quantum mechanics (QM) calculations and experimental physical properties .