So the final unit should then be bohr*e, great. I am actually using electron difference density with a script adapted from your tutorial, but thanks for the heads up. I have a final curious question, would this dipole approach actually be applicable for periodic boundary conditions or is this purely molecular?
I am trying to calculate both polarization using the dipole for a monolayer with a bias. The bias I create with the metal plates but the dipole and seem to respond rather drastically yielding too high values, so I just wondered if maybe I shouldn't do dipole with periodic boundaries?
Hope you have some insight