DIPOL

DIPOL = [real array]
Default: DIPOL = Template:FIXME

For VASP.4.X behavior please refer to Monopole, Dipole and Quadrupole corrections.

It is possible to calculate the total dipole-moment in the cell, using the option

DIPOL=R_x R_y R_z

where R_x, R_y and R_z specify the center of the cell in direct lattice coordinates.

Mind: the calculation of the dipole requires a definition of the center of the cell, and results might differ for different positions. You should use this option only for surfaces and isolated molecules. In this case use the center of mass for the position (for surface only the component normal to the surface is meaningful).

The main problem is that the definition of the dipole 'destroys' the translational symmetry, i.e. the dipole is defined as

{\displaystyle \int ({\mathbf r}-{\mathbf R}_{\rm center}) \rho_{\rm ions+valence}{\mathbf r} d^3 {\mathbf r}. }

Now, this makes only sense if $\rho_{\rm ions+valence}$ drops to zero at some distance from $\mathbf R_{\rm center}$ . If this is not the case, the values are extremely sensible with respect to changes in $\mathbf R_{\rm center}$ .

Related Tags and Sections

Contents