The ECP order ten variant of our model replicated the ab initio calculated energies with a fair degree of accuracy, except for the N 2 + N 2 complexes. For the H 2 + H 2 complexes the EFP results displayed higher accuracy compared to our model, even at the higher energy regimes. Nevertheless, our model also displayed a fair degree of accuracy for this complex. For the N 2 + N 2 complexes both our model and MMFF9S displayed similar levels of accuracy, both differing to the ideal ab initio result by a factor of approximately 3.0 (Figure ). For the O 2 + O 2 and CO 2 + CO 2 complexes our model displayed high levels accuracy outperforming all the other molecular force fields considered in the high energy regions. Notably, for the CO 2 + CO 2 complexes MB-Pol displayed a very close to linear relation against the ab initio calculated data, however, a constant multiplicative factor of roughly 6.3 between the ab initio and MB-Pol calculated results is observed as shown in the bottom right panel of Figure . All the comparisons considered against other established force field models are constrained to instances wherein the molecules of the dimer complexes are encompassed within their default fragment libraries, and/or those cases where the models do not fail to run. in the supplementary material, despite the its many improvements. Therefore, this version of our model is inadequate for intramolecular geometry optimization, hence, the geometries of the molecules in the system must be provided either from a reliable database, or a theoretical calculation capable of making such kind molecular geometry optimization.