DIPOL: Difference between revisions

From VASP Wiki
mNo edit summary
 
(11 intermediate revisions by 4 users not shown)
Line 1: Line 1:
{{TAGDEF|DIPOL|[real array]|{{FIXME|default none or the origin 0 0 0}}|}}
{{TAGDEF|DIPOL|[real array]}}


Description: calculate the total dipole-moment in the cell, using the option
Description: specifies the center of the cell in direct lattice coordinates with respect to which the total dipole-moment in the cell is calculated.
{{TAG|DIPOL}}='''R'''<sub>x</sub> '''R'''<sub>y</sub> '''R'''<sub>z</sub>,
where '''R'''<sub>x</sub>, '''R'''<sub>y</sub> and '''R'''<sub>z</sub>, to 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.
The center of the cell w.r.t. which the total dipole-moment in the cell is calculated is specified as
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
{{TAG|DIPOL}}='''R'''<sub>x</sub> '''R'''<sub>y</sub> '''R'''<sub>z</sub>


where '''R'''<sub>x</sub>, '''R'''<sub>y</sub> and '''R'''<sub>z</sub> are given in direct lattice coordinates.
Calculations using the dipole correction, i.e. using tags {{TAG|IDIPOL}} or {{TAG|LDIPOL}}, require a definition of the center of the cell. Results of the computed dipole moment might differ for different positions. The reason for this difference is that the definition of the dipole 'destroys' the translational symmetry, i.e., the dipole is defined as
:<math>
:<math>
\int ({\mathbf r}-{\mathbf R}_{\rm center}) \rho_{\rm ions+valence}{\mathbf r} d^3 {\mathbf r}.
\int ({\mathbf r}-{\mathbf R}_{\rm center}) \rho_{\rm ions+valence}({\mathbf r}) d^3 {\mathbf r}.
</math>
</math>  
This measure will provide consistent values only if <math>\rho_{\rm ions+valence}</math> drops to zero at some distance from <math>\mathbf R_{\rm center}</math>. If this is
not the case, the values are extremely sensitive with respect to changes in <math>\mathbf R_{\rm center}</math>. In such cases, it might be beneficial to increase the size of the cell along the vacuum dimension (for surfaces) or for the entire cell (for isolated molecules).
{{NB|mind| If the flag is not set, VASP determines where the charge density averaged over one plane drops to a minimum and calculates the center of the charge distribution by adding half of the lattice vector perpendicular to the plane where the charge density has a minimum (this is a rather reliable approach for orthorhombic cells)}}
{{NB|tip| For calculations of isolated molecules and surfaces with the dipole correction, use {{TAG|DIPOL}} as the center of mass of the atoms in your cell. Additionally, note that for surfaces, only the component normal to the surface is meaningful.}}


Now, this makes only sense if <math>\rho_{\rm ions+valence}</math> drops to zero at some distance from <math>\mathbf R_{\rm center}</math>. If this is
== Related tags and articles ==
not the case, the values are extremely sensible with respect to changes in <math>\mathbf R_{\rm center}</math>.
 
== Related Tags and Sections ==
{{TAG|NELECT}},
{{TAG|NELECT}},
{{TAG|EPSILON}},
{{TAG|EPSILON}},
Line 26: Line 26:
{{TAG|EFIELD}},
{{TAG|EFIELD}},
[[Monopole_Dipole_and_Quadrupole_corrections|Monopole, Dipole and Quadrupole corrections]]
[[Monopole_Dipole_and_Quadrupole_corrections|Monopole, Dipole and Quadrupole corrections]]
{{sc|DIPOL|Examples|Examples that use this tag}}
----
----
[[The_VASP_Manual|Contents]]


[[Category:INCAR]]
[[Category:INCAR tag]][[Category:Molecules]][[Category:Monopole Dipole and Quadrupole Corrections]]

Latest revision as of 12:08, 18 October 2023

DIPOL = [real array] 

Description: specifies the center of the cell in direct lattice coordinates with respect to which the total dipole-moment in the cell is calculated.


The center of the cell w.r.t. which the total dipole-moment in the cell is calculated is specified as

DIPOL=Rx Ry Rz

where Rx, Ry and Rz are given in direct lattice coordinates.

Calculations using the dipole correction, i.e. using tags IDIPOL or LDIPOL, require a definition of the center of the cell. Results of the computed dipole moment might differ for different positions. The reason for this difference is that the definition of the dipole 'destroys' the translational symmetry, i.e., the dipole is defined as

This measure will provide consistent values only if drops to zero at some distance from . If this is not the case, the values are extremely sensitive with respect to changes in . In such cases, it might be beneficial to increase the size of the cell along the vacuum dimension (for surfaces) or for the entire cell (for isolated molecules).

Mind: If the flag is not set, VASP determines where the charge density averaged over one plane drops to a minimum and calculates the center of the charge distribution by adding half of the lattice vector perpendicular to the plane where the charge density has a minimum (this is a rather reliable approach for orthorhombic cells)
Tip: For calculations of isolated molecules and surfaces with the dipole correction, use DIPOL as the center of mass of the atoms in your cell. Additionally, note that for surfaces, only the component normal to the surface is meaningful.

Related tags and articles

NELECT, EPSILON, IDIPOL, LDIPOL, LMONO, EFIELD, Monopole, Dipole and Quadrupole corrections

Examples that use this tag