Molecular dynamics force fields (FFs) have been used for decades to theoretically study the dynamics and variable structure of ribonucleic acids (RNA). Currently, there exist a handful of force fields with specific parameter sets used to study nucleic acids, including CHARMM, 1,2 OPLS, and AMBER. Since their initial introduction, many of these models have continuously been improved upon to better describe the unique and defining characteristics of RNA, such as backbone flexibility, 2,7 base stacking, phosphate electrostatics, and hydrogen bonding. While each of these models may differ in their parameterization strategy, what most of them have in common is the formulation of the energy functions that describe nonbonded interactions. Traditionally, these interactions have been represented by three separate empirical functions that approximate the non-covalent interactions between base pairs and other parts of RNA. Among these interactions are atom centered chargecharge Coulomb energies and a combination of London dispersion and Pauli repulsion forces into a single Lennard-Jones term. Although there is a long history of using these force fields to assist experimental structure prediction, the formulation of these models over two decades ago was mostly intended to circumvent the lack of computational power rather than to replicate physically realistic interactions.