The results of molecular dynamics simulations are highly dependent on the interatomic potential (IP) used. Ab-initio molecular dynamics (AIMD) and neural network-based (NNIP) IPs have been very successful at reproducing structural properties of various material and molecular systems. However, these are too computationally expensive to access the timeand length-scales needed to ergodically investigate amorphous solid polymer electrolyte systems, which are in the hundreds of nanoseconds and tens of thousands of atoms. Thus, classical IPs which incorporate nonreactive covalent interactions, Lennard-Jones dispersion interactions, and nonpolarizable point charge Coulomb interactions are the most accessible approach for modeling these systems. However, the most prominent classical IPs such as OPLS, 8,9 PCFF+, and COMPASS were parameterized through decades of handtuning through a combination of bottom up fitting on computational chemistry and top-down fitting to reproduce experimental values. Additionally, these parameters are often obfuscated by commercial licenses. Furthermore, it is common practice to independently adjust these IP parameters (especially the partial charges assigned to each atom) for individual studies. For example, France-Lanord 12 and Molinari 4 both report report using the proprietary MedeA package from Materials Design to obtain the interatomic potential yet still use different partial charges on the TFSI atoms. A methodical and reproducible parameterization scheme to benchmark existing literature is needed.