Despite the progress made, there are aspects that warrant further exploration, such as the incorporation of polarization effects in the inner core shells, thereby avoiding the treatment of them as frozen. This probably is important in high-energy molecular dynamics simulations. Additionally, efforts could be directed towards deriving exact or approximate expressions when infinite expansions are employed to mitigate the numerical noise associated with low-energy interactions. Nevertheless, for up to a tenth order truncated series expansion the numerical noise for the ECP calibrated variant is negligible for the various dimer configurations explored in this study, with reasonably accurate predictions for both the interaction energies and partial charges for neutral low energy configurations.