External potential at atomic sites

Hello All,

I am wondering if there is a way to add external potential at atomic sites in siesta calculations.
I have calculated the potential distribution (using the Comsol electrostatics module) for a structure under consideration that produces a nonuniform potential distribution. I would like to add this potential to siesta calculations. Is it possible?
Since, the potential is non-uniform, the external electric field approach does not work.

Any insights would be really helpful.


not sure but maybe something like “simulatedoping” might be the way