LCALCEPS: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
(14 intermediate revisions by 4 users not shown)
Line 31: Line 31:
from the [[Berry_phases_and_finite_electric_fields#Self-consistent_response_to_finite_electric_fields|self-consistent response to a finite electric field]] <span style="font-size:16pt">''&epsilon;''</span>.
from the [[Berry_phases_and_finite_electric_fields#Self-consistent_response_to_finite_electric_fields|self-consistent response to a finite electric field]] <span style="font-size:16pt">''&epsilon;''</span>.
In this case, the "response" of the system is the change in the polarization '''P''', the Hellmann-Feynman forces '''F''', and the stress tensor &sigma;.
In this case, the "response" of the system is the change in the polarization '''P''', the Hellmann-Feynman forces '''F''', and the stress tensor &sigma;.
If this is combined with {{TAG|IBRION}}=6, the contribution from the ionic relaxations to the piezoelectric and dielectric tensors are calculated as well.
To this end VASP will perform essentially three successive calculations, with:
{{TAG|EFIELD_PEAD}}= <span style="font-size:16pt">''&epsilon;''</span><sub>x</sub> 0 0
{{TAG|EFIELD_PEAD}}= 0 <span style="font-size:16pt">''&epsilon;''</span><sub>y</sub> 0
{{TAG|EFIELD_PEAD}}= 0 0 <span style="font-size:16pt">''&epsilon;''</span><sub>z</sub>
where, by default, VASP chooses <span style="font-size:16pt">''&epsilon;''</span><sub>x</sub>=<span style="font-size:16pt">''&epsilon;''</span><sub>y</sub>=<span style="font-size:16pt">''&epsilon;''</span><sub>z</sub>=0.01 eV/&Aring;.
This default may be overwritten by specifying
{{TAG|EFIELD_PEAD}}= <span style="font-size:16pt">''&epsilon;''</span><sub>x</sub> <span style="font-size:16pt">''&epsilon;''</span><sub>y</sub> <span style="font-size:16pt">''&epsilon;''</span><sub>z</sub>
in the {{FILE|INCAR}} file.
The relevant output is found in the {{FILE|OUTCAR}} file, immediately following the lines (see the description of {{TAG|LEPSILON}}=.TRUE. as well):
MACROSCOPIC STATIC DIELECTRIC TENSOR (including local field effects)
BORN EFFECTIVE CHARGES (including local field effects)
PIEZOELECTRIC TENSOR (including local field effects)
In the above, "including local field effects" pertains to the fact that changes in the orbitals due to the electric field induce changes in the Hartree- and exchange-correlation potential. One may choose to limit this to changes in the Hartree potential alone, by specifying:
{{TAG|LRPA}}=.TRUE.
This is commonly referred to as the response within the "Random Phase Approximation" (RPA), or the "neglect of local field effects". The OUTCAR file will now contain additional sections, headed by the lines:
MACROSCOPIC STATIC DIELECTRIC TENSOR (excluding local field effects)
BORN EFFECTIVE CHARGES (excluding local field effects)
PIEZOELECTRIC TENSOR (excluding local field effects)
----
'''Note''': For standard DFT functionals, &epsilon;<sub>&infin;</sub>, ''Z''<sup>*</sup>, and ''e''<sup>(0)</sup> may be more easily calculated from density functional perturbation theory (see {{TAG|LEPSILON}}=.TRUE.). For functionals that depend not only on the density but also explicitly on the orbitals, like hybrid functionals, density functional perturbation theory is presently not implemented and {{TAG|LEPSILON}}=.TRUE. is not applicable.
'''Note''': The piezoelectric tensor has the wrong sign in Vasp 5.4.4 and older. The bug is fixed with [http://cms.mpi.univie.ac.at/patches/patch.5.4.4.16052018.gz patch.5.4.4.16052018.gz].


== Related Tags and Sections ==
== Related Tags and Sections ==
Line 42: Line 85:
{{TAG|NPPSTR}},
{{TAG|NPPSTR}},
{{TAG|DIPOL}},
{{TAG|DIPOL}},
{{TAG|IBRION}},
[[Berry_phases_and_finite_electric_fields|Berry phases and finite electric fields]]
[[Berry_phases_and_finite_electric_fields|Berry phases and finite electric fields]]
{{sc|LCALCEPS|Examples|Examples that use this tag}}
----
----
[[The_VASP_Manual|Contents]]
[[The_VASP_Manual|Contents]]


[[Category:INCAR]][[Category:Berry phases]]
[[Category:INCAR]][[Category:Dielectric Properties]][[Category:Berry phases]][[Category:Howto]]

Revision as of 13:50, 2 July 2019

LCALCEPS = .TRUE. | .FALSE.
Default: LCALCEPS = .FALSE. 

Description: for LCALCEPS=.TRUE. the macroscopic ion-clamped static dielectric tensor, Born effective charge tensors, and the ion-clamped piezoelectric tensor of the system are determined from the response to finite electric fields.


For LCALCEPS=.TRUE., VASP calculates the ion-clamped static dielectric tensor

the Born effective charge tensors

and the ion-clamped piezoelectric tensor of the system

from the self-consistent response to a finite electric field ε. In this case, the "response" of the system is the change in the polarization P, the Hellmann-Feynman forces F, and the stress tensor σ.

If this is combined with IBRION=6, the contribution from the ionic relaxations to the piezoelectric and dielectric tensors are calculated as well.

To this end VASP will perform essentially three successive calculations, with:

EFIELD_PEAD= εx 0 0
EFIELD_PEAD= 0 εy 0
EFIELD_PEAD= 0 0 εz 

where, by default, VASP chooses εx=εy=εz=0.01 eV/Å.

This default may be overwritten by specifying

EFIELD_PEAD= εx εy εz

in the INCAR file.

The relevant output is found in the OUTCAR file, immediately following the lines (see the description of LEPSILON=.TRUE. as well):

MACROSCOPIC STATIC DIELECTRIC TENSOR (including local field effects) 
BORN EFFECTIVE CHARGES (including local field effects) 
PIEZOELECTRIC TENSOR (including local field effects) 

In the above, "including local field effects" pertains to the fact that changes in the orbitals due to the electric field induce changes in the Hartree- and exchange-correlation potential. One may choose to limit this to changes in the Hartree potential alone, by specifying:

LRPA=.TRUE.

This is commonly referred to as the response within the "Random Phase Approximation" (RPA), or the "neglect of local field effects". The OUTCAR file will now contain additional sections, headed by the lines:

MACROSCOPIC STATIC DIELECTRIC TENSOR (excluding local field effects) 
BORN EFFECTIVE CHARGES (excluding local field effects) 
PIEZOELECTRIC TENSOR (excluding local field effects)

Note: For standard DFT functionals, ε, Z*, and e(0) may be more easily calculated from density functional perturbation theory (see LEPSILON=.TRUE.). For functionals that depend not only on the density but also explicitly on the orbitals, like hybrid functionals, density functional perturbation theory is presently not implemented and LEPSILON=.TRUE. is not applicable.

Note: The piezoelectric tensor has the wrong sign in Vasp 5.4.4 and older. The bug is fixed with patch.5.4.4.16052018.gz.

Related Tags and Sections

LEPSILON, LCALCPOL, EFIELD_PEAD, LPEAD, IPEAD, LBERRY, IGPAR, NPPSTR, DIPOL, IBRION, Berry phases and finite electric fields

Examples that use this tag


Contents