How can I plot the Vxc from taking the difference between Vtot and Vhar from LOCPOT by using pymatgen?