In Molecular Dynamics (MD) simulations the selection of force field(s) is as a central factor one needs to consider to obtain accurate results. The rationale behind employing force fields stems in their ability to transition from having to solve a difficult high-dimensional partial differential equation in the form of Schrödinger's equation, to a tractable many-body problem in the form of a system of coupled second-order ordinary differential equations (ODEs) embodied in Newton's laws of motion.