Dear Support,
I am resubmitting a post that I submitted previously (https://vasp.at/forum/viewtopic.php?p=18825#p18825). An updated version of it is given below:
I have a couple of questions about correcting the energy of a surface or slab model with a net dipole moment along the surface normal direction.
When performing slab calculations, a post hoc correction to the energy of a slab model with a net dipole moment along the surface normal can be made by adding a corrective term to the computed energy. In VASP, this correction is done using the VASP keyword IDIPOL = 1, 2, or 3. For example, if the surface normal is parallel to the c lattice vector, then the correction is done using IDIPOL = 3. A self-consistent correction to the potential of the system can also be made by setting the VASP keywords IDIPOL, LDIPOL, and DIPOL. However, I would like to ask about the correction to the energy.
My understanding is that the dipole moment of a periodic system is dependent on the choice of origin. If I am correct, then the size of the correction to the energy is also dependent on the choice of origin. My questions are:
1) What is the theoretical justification or basis for this correction if the correction is dependent on the choice of origin?
2) As discussed in the VASP manual under the DIPOL keyword (https://www.vasp.at/wiki/index.php/DIPOL#), to make the correction to the energy, the origin is chosen to include the plane of maximum charge density parallel to the surface, unless it is explicitly set using the DIPOL keyword. Why is this choice of origin made? Is it based on some convention?
Thanks,
Yves