Hi Anders, thanks very much for your help. I tried to use the ExternalPotential as advised to calculate the Hartree potential of a graphene ribbon with zero bias voltage using the DFTB method, the resulted plot is just a plane with uniform color (ie one value). This is very different from Hartree potential figures published in the literature that have may varying values across the ribbon. Any further advice? Thanks and regards, wanzhi.