As one would expect, neither the full electron nor ECP order ten variants of our model provide accurate estimations of the vertical electron affinities of the various small molecules as summarized in Table . Specifically, the full electron variant shows a much lesser accuracy, in terms of magnitude, of the two when compared to the equivalent ab initio calculated result with discrepancies with two or more orders of magnitude higher. These results come to no surprise given that our model is unable to model chemical bonds (refer to Figure of the supplementary material). Nevertheless, a detailed exploration and reporting of the strengths and weaknesses of our model is essential for those aiming to conduct accurate MD simulations.