Vacuum level / band alignment for isolated polar molecules from LOCPOT
Dear VASP users and developers,
I would like to ask for advice on determining the vacuum level for isolated molecule calculations in VASP, especially for polar molecules.
Our goal is to perform vacuum-level band alignment between several isolated molecules and a substrate. In practice, we want to compare molecular frontier levels, such as HOMO/LUMO, with substrate band edges. For the substrate/slab, the usual procedure is clear: compute the planar-averaged electrostatic potential from `LOCPOT`, identify the vacuum plateau, and align the band edges to that vacuum level.
However, we are less certain about the correct protocol for isolated molecules. Some of the molecules are polar, for example H₂O. We tested an isolated molecule in a relatively large box (approximately 20 × 21 × 22 ų).
We also tested calculations both without and with dipole correction:
LDIPOL = .TRUE.
IDIPOL = 4
DIPOL = 0.5 0.5 0.5
The molecule is placed near the center of the box. Setting LVHAR = True, we output `LOCPOT` and then computed the planar-averaged electrostatic potential along one direction using vaspkit or in-house script (we also triedvtotav) / planar averaging.
The issue is that even with LDIPOL = .TRUE. and IDIPOL = 4, we do not obtain a clear flat vacuum plateau in the planar-averaged potential. Instead, the planar average still shows a slope/curvature in the vacuum region. This happens both without and with the dipole correction, although the shape changes.
I understand that for an isolated polar molecule the electrostatic potential is directional, and the far-field potential should behave approximately like a dipole field. Therefore, compressing the full 3D potential into a 1D planar average may not produce a unique plateau, and the result may depend on whether the average is taken along x, y, or z.
My questions are:
1. For isolated polar molecules, is it expected that the planar-averaged potential does not show a clear vacuum plateau, even when using `LDIPOL = .TRUE.` and `IDIPOL = 4`? If so, how can we obtain such plateau?
2. Is there a recommended VASP protocol for obtaining a vacuum reference for molecular orbital energies of isolated molecules? Some guidelines, in particular, for band aligment?
I will upload representative input/output files and plots of the planar-averaged potential for reference.
Thank you very much for any advice or clarification.