Why VASP gives no straight vacuum line for the surface local potential?
Posted: Thu Apr 08, 2021 2:46 pm
Dear colleagues,
I'm calculating local potential (LOCPOT) data in VASP for the semiconductor/metal interface in order to extract work functions (WF), VBM, CBM, etc. It's known that in the case of two different WFs a step will appear when dipole correction is included (e.g. see https://www.researchgate.net/post/Who-c ... k-function)
Indeed, I easily got such a step for the case of Ag/h-BN interface, for example (see Fig. 1) when activate LDIPOL = .TRUE., IDIPOL = 3 (i.e. Z direction). However, for some other semiconductor/metal interfaces when using the same dipole correction I have very strange curve at the vacuum positions that looks like apple cores or oscillations.
Nevertheless, in Fig. 2 one can still see a step despite strange oscillations. However, in Fig. 3 no step is presented at all. So, my two questions are:
1) Why do these strange curves may appear at the vacuum level? (instead of a simple step as in Fig.1)
2) What should I change/control to fix it and find a step
P.S.
I have already played with vacuum size. In Fig.2 and 3 it's twice higher than the structure thickness. No evidences that further increase of vacuum size can help me.
Also, I was tested a bit DIPOL flag (center mass position for dipole correction) and found that VASP can set it automatically without any problem when the vacuum is high enough. But maybe I'm wrong and I have to set up manually DIPOL flag in some trouble cases.
I will highly appreciate any ideas and advices! Thank you.
I'm calculating local potential (LOCPOT) data in VASP for the semiconductor/metal interface in order to extract work functions (WF), VBM, CBM, etc. It's known that in the case of two different WFs a step will appear when dipole correction is included (e.g. see https://www.researchgate.net/post/Who-c ... k-function)
Indeed, I easily got such a step for the case of Ag/h-BN interface, for example (see Fig. 1) when activate LDIPOL = .TRUE., IDIPOL = 3 (i.e. Z direction). However, for some other semiconductor/metal interfaces when using the same dipole correction I have very strange curve at the vacuum positions that looks like apple cores or oscillations.
Nevertheless, in Fig. 2 one can still see a step despite strange oscillations. However, in Fig. 3 no step is presented at all. So, my two questions are:
1) Why do these strange curves may appear at the vacuum level? (instead of a simple step as in Fig.1)
2) What should I change/control to fix it and find a step
P.S.
I have already played with vacuum size. In Fig.2 and 3 it's twice higher than the structure thickness. No evidences that further increase of vacuum size can help me.
Also, I was tested a bit DIPOL flag (center mass position for dipole correction) and found that VASP can set it automatically without any problem when the vacuum is high enough. But maybe I'm wrong and I have to set up manually DIPOL flag in some trouble cases.
I will highly appreciate any ideas and advices! Thank you.