Electrostatic energy calculations to obtain optimized Dopant Site of a Slab


I am a new user of Pymatgen. I am trying to adsorb a molecule over a slab surface. Additionally I also need to calculate the optimized (preferred) dopant site with respect to the electrostatic energy. Can anyone guide me in this issue.

Thanks in Advance.