Discrepancy between Fermi levels and electrostatic potential from LOCPOT for heterostructure?
I thought I would start over with a simpler question since I think I scared everyone away with my last Related to the same topic, I am confused about how to reconcile two difference outputs I have from the same calculation for an insulator-heavy metal heterostructure with vacuum. If I calculate the layer-projected density of states for the heterostructure, I see that in the center of the insulator (Cr2O3) region of the slab, the highest occupied states are below the highest occupied states in the center of the heavy metal (Pt), which sit at the heterostructure Fermi level, marked by the dashed black line. Then I would expect electrons to flow from Pt to Cr2O3 at the interface to try to equilibrate chemical potentials (and indeed, this is what looks to happen for the layers near the interface).
If I plot the electrostatic potential as a function of distance along the slab, with LVTOT=.TRUE. however, I see that the average electrostatic potential (marked by the red line) in the Cr2O3 part of the slab is higher/less negative than the average electrostatic potential in the Pt region of the slab. Since VASP plots the electrostatic potential with respect to electrons, this seems to imply that electrons should flow from Cr2O3 (higher potential) to Pt (lower potential) at the interface. So the exact opposite of what is implied by the density of states. It seems to me like if the highest occupied states in the Cr2O3 region of the slab are lower in energy than in the Pt region as indicated by the DOS, this has to mean that the electrostatic potential of the Cr2O3 region is lower in energy (more tightly bound).
If you can explain how I should understand this apparent discrepancy in the output, I would greatly appreciate it. Also, I'm curious as to what the reference level is set as for the LOCPOT, since it seems that some areas of the Pt potential are actually above zero; perhaps this is because I've included the exchange-correlation energy?
I attached the input files and the image showing the layer-projected density of states and the LOCPOT output which to me seem to contradict each other.