Whether these models are designed to act as a more generalized force field, like MASTIFF, MEDFF, or GMM/GEM, or tailored towards specific systems such as as the HIPPO water model, their use of atom-centered densities succeed in accounting for electron overlap effects that traditional MM models exclude from their Pauli repulsion or frozen electrostatics terms. What makes these models unique is that the interactions between monomers are explicitly formulated and parameterized to replicate energies from ab initio decomposition schemes, such as SAPT or CSOV, that map ab initio energies onto molecular mechanics counterparts. Following in their lead, we chose to both parameterize AmberFD using a similar decomposition scheme, as well as formulate the model energies themselves off of the appropriately termed Energy Decomposition Analysis (EDA) based on the use of absolutely localized molecular orbitals (ALMO). We note here that EDA is not the only method available, and there is not necessarily a unique prescription to correlate QM energies with a classical representation, as seen by some of the other and sometimes conflicting decomposition schemes of DEDA, NEDA, CSOV, and SAPT; we simply choose to use EDA over its ease of use and not being dependent on a perturbation theory. However, choosing one of these decomposition methods (or another by preference) as a baseline provides a set of energies to fit against that are ab initio in origin and allows for a systematic way to further improve or assess force field predictions in the future. In the remainder of the paper, we will discuss our efforts to parameterize such a model specifically tailored to RNA nucleobasenucleobase interactions using EDA. We will also compare the results of both a training and testing data set of base pair energies and other small molecule properties, and how AmberFD compares to two popular polarizable force fields for nucleic acids, the CHARMMM Drude oscillator and the AMOEBA multipole models.