Chelpg Gaussian, It can be done by using the keyword pop=chelpg (it means that charges from electrostatic .