This version of our model is calibrated with the same reference data that was used to calibrate the previous nonpolarizable version. In our previous study, the valence electron densities of all the atoms, in the first three rows of the periodic table, are fitted, and tabulated, via least-squares minimization against the electronic density obtained from Density Functional Theory (DFT) calculations using the B3LYP 37 hybrid functional using Grimme's empirical dispersion formula GD3 using the Effective Core Potential (ECP) basisset CEP-31G for all atoms except hydrogen and helium where the full-electron basis-set cc-pVDZ is used instead. The water dimer interaction calculations used as reference data consist on a 25 × 200 uniform grid of calculations in an O • • • O distance coordinate, and a 25 × 200 uniform grid of calculations in an O • • • H distance coordinate, where the interatomic distances are sampled with 200 points whereas the angular component at which one of the molecules of water is rotated in an axis perpendicular to the plane of the molecule is sampled with 25 points, using CCSD(T) level of theory with the basis-set cc-pVTZ. To validate our model, comparisons between different energies predicted by our model and accurate ab initio calculations are made on a set of different complexes with small molecules, namely H 2 , N 2 , O 2 , and CO 2 . These comparisons include the binding energies between the two atoms, of all the diatomic molecules just mentioned, to determine its ability to predict covalent bonds. Comparisons between Vertical Electron Affinities (VEA), and interaction energies between dimer complexes of two identical molecules in their neutral and anion states are also provided. The comparisons of the dimer complexes consist in displacing one of the molecules by a distance ∆L ∈ [1, 5], in Å, and rotating it by an angle θ ∈ {0, π/4, π/2}.