Molecular dynamics and Monte Carlo simulations are increasingly being used in a variety of fields from the study of materials and biopolymers to computer aided drug design . The accuracy of these simulations depends on the accuracy of the their potential function for calculation the energy and forces in the system. For these large systems, that potential function is typically a fixed charge, all atom force field which uses the positions of all atoms in the system to calculate the potential energy . While this is a simplistic representation, molecular simulations using force fields have shown promise in a variety of applications such as calculating hydration and solvation free energies as well as protein-ligand binding . That being said, the literature is littered with examples of force field limitations . Automating parameterization would allow tests of which functional forms result in the best and most general force fields .