In a broader analysis of the predicted partial charges from both our ECP and full electron fitted variants of our model, compared to their respective Mulliken charges computed at the CCSD(T)/CEP-31G and CCSD(T)/cc-pVTZ levels of theory, a consistent trend emerges. Both versions of our fitted model showcase acceptable accuracy in estimating Mulliken charges, particularly for configurations with lower energy (see Figure ). Notably, the ECP variant slightly outperforms in terms of accuracy, boasting a slightly larger coefficient of determination against its reference ab initio data. This marginal improvement in accuracy is likely also influenced by the inherent numerical error noise introduced due to the use of limited floating-point precision representation and the necessity to handle exponentiation of relatively large numbers, as evident in the decaying coefficients of the Gaussian basis functions used to approximate atomic electron contributions (refer to Table of the supplementary material), namely for fitted full electron fitted densities.