The final limitation in our previous model that we intend to correct is its inability to polarize the molecules in a form that allows for charge redistribution. With this we mean that the atomic contributions to the total electron density can shift depending of the environment of the atoms. It is noted that our previous model accounts for dipole moments of the molecules to an extent, since the atomic density contributions are multiplied by scaling a scaling factor that ensures that these contributions agree with the Mulliken partial charges obtained from high accuracy ab initio calculations. This has the disadvantage of requiring ab initio calculations on every molecular species in the system that one aims to model.