To enhance the accuracy of the full electron variant of our model, for high-energy interactions, the incorporation of similar polarization coefficients into the core shell contributions is needed. That is, no longer treating the inner shell contributions as frozen electron clouds. Incorporating these inner shell coefficients would require significant changes to the fitting of the parameters used to model partial charges since relying solely on Mulliken charges as reference would be insufficient. Instead, this would require the consideration of shell-specific contributions to the electron density for each configuration within our calibration dataset. This more nuanced approach is beyond the scope of this study. Moreover, the inclusion of post Born-Oppenheimer corrections for high-energy configurations is another factor to consider.