Born effective charges: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
 
(4 intermediate revisions by the same user not shown)
Line 1: Line 1:
Born effective charges  
The change in polarization from the displacement of an atom is not uniquely defined in periodic systems, where atoms are repeated in different cells and the charge can be generalized.{{cite|ghosez:michenaud:gonze:1998}} Born effective charges are one way of defining this dynamical charge.


# https://journals.aps.org/prb/abstract/10.1103/PhysRevB.58.6224
= Introduction =
The dynamical charge is defined as the cell volume <span>'''&Omega;'''</span><sub>0</sub>, multiplied by the partial derivative of the macroscopic polarization '''P''' in the direction ''i'' with respect to a rigid displacement of the sublattice of atoms ''κ'' in the direction ''j''.


Equation:
However, the polarization is not uniquely defined in periodic systems and depends on the macroscopic electric field <math>\mathcal{E}_i</math> fixed by the periodic boundary conditions. The Born effective charge '''Z*''' is the partial derivative of the polarization with respect to position ''u'' for zero macroscopic electric field.{{cite|gonze:prb:1997}} As polarization is the first derivative of the total energy with respect to the macroscopic electric field, '''Z*''' may be rearranged in terms of the partial derivative of the force '''F''' in direction ''j'' on atom ''κ'' with respect to <math>\mathcal{E}_i</math>:


::<math>
::<math>
Z_{\kappa,\beta\alpha}^*
Z_{\kappa,ij}^*
=\Omega_0\frac{\partial \mathcal{P}_{\mathrm{mac},\beta}}{\partial \tau_{\kappa\alpha}(\textbf{q=0})}
=\frac{\Omega_0}{e} \frac{\partial \mathcal{P}_i} {\partial u_{\kappa,j}(\textbf{q=0})}
=\Omega_0\frac{\partial F_{\kappa\alpha}}{\partial \mathcal{E}_\beta}
=\frac{1}{e} \frac{\partial F_{\kappa,j}}{\partial \mathcal{E}_i}
\qquad
\qquad
{\alpha\beta=x,y,z,}.
{i,j=x,y,z}
</math>
</math>


# https://journals.aps.org/prb/abstract/10.1103/PhysRevB.55.10355
{{NB|mind|
*The '''*''' does not denote a complex conjugate, '''Z*''' is always a real quantity.
*'''Z*''' is given in units of <math>\vert e \vert</math> in VASP.
*VASP outputs <span><b>Z</b></span><sub><i>ij</i></sub><sup><span>&#8902;</span></sup> with ''i'' for the macroscopic electric field, and ''j'' for the direction of the force. In literature, <span><b>Z</b></span><sub><i>ji</i></sub><sup><span>&#8902;</span></sup> is commonly seen, with the force direction ''j'' followed by the electric field direction ''i''. Note, py4vasp follows the latter notation <span><b>Z</b></span><sub><i>ji</i></sub><sup><span>&#8902;</span></sup> for historic reasons.
}}
 
= How to calculate =
There are two ways of computing Born effective charges in VASP. The first is done using {{TAG|LCALCEPS}}, where a finite electric field is applied along the three cartesian directions and the resultant forces on the atoms are calculated:
 
{{TAG|LCALCEPS}} = .TRUE.
 
The other approach is done using  {{TAG|LEPSILON}}, where the derivative of the wavefunction with respect to an electric field is calculated using density functional perturbation theory (DFPT):
 
{{TAG|LEPSILON}} = .TRUE.
 
These may be used in combination with {{TAG|IBRION}} to obtain additional dielectric properties:


{{TAG|IBRION}} = 5 or 6 ! Calculated using finite differences.
{{TAG|IBRION}} = 7 or 8 ! Calculated using DFPT


= Introduction =  
For more details, see the pages for each tag. The Born effective charges including local field effects will be given in the {{TAG|OUTCAR}} file:
Insert something from paper/ thesis then edit
 
BORN EFFECTIVE CHARGES (including local field effects) (in |e|, cummulative output)
 
= Excluding local field effects =  


Previously, the local field effects have been included, that is changes in the orbitals due to the electric field induce changes in the Hartree- and exchange-correlation potentials. This may be limited to changes in the Hartree potential, by specifying:


MIND THAT VASP outputs it oddly?
{{TAG|LRPA}} = .TRUE.
{{TAG|LCALCEPS}} = .TRUE. ! N.B. {{TAG|LEPSILON}} does not output the final Born effective charges.


= Calculating =
This prints out the Born effective charges excluding local field effects:
Taken from Dielctric properties:
There are two approaches to compute Born effective charges implemented in VASP: one is done by applying finite electric fields along the three cartesian directions and computing the forces on the atoms which are activated using LCALCEPS or by computing the derivating of the wavefunction with respect to an electric field using density functional perturbation theory (DFPT) using LEPSILON.


Note: This is different in py4vasp.
BORN EFFECTIVE CHARGES (excluding local field effects) (in |e|, cummulative output)


These are calculated normally but remain hidden unless explicitly specified.


== Related tags and articles ==


<!--[[Category:Dielectric properties]][[LEPSILON]][[LCALCEPS]][[Category:Howto]]-->
{{TAG|LEPSILON}},
{{TAG|LCALCEPS}},
{{TAG|IBRION}},
{{TAG|LRPA}}


==References==
==References==
[[Category:Dielectric properties]][[Category:Howto]]

Latest revision as of 12:54, 13 August 2024

The change in polarization from the displacement of an atom is not uniquely defined in periodic systems, where atoms are repeated in different cells and the charge can be generalized.[1] Born effective charges are one way of defining this dynamical charge.

Introduction

The dynamical charge is defined as the cell volume Ω0, multiplied by the partial derivative of the macroscopic polarization P in the direction i with respect to a rigid displacement of the sublattice of atoms κ in the direction j.

However, the polarization is not uniquely defined in periodic systems and depends on the macroscopic electric field fixed by the periodic boundary conditions. The Born effective charge Z* is the partial derivative of the polarization with respect to position u for zero macroscopic electric field.[2] As polarization is the first derivative of the total energy with respect to the macroscopic electric field, Z* may be rearranged in terms of the partial derivative of the force F in direction j on atom κ with respect to :


Mind:
  • The * does not denote a complex conjugate, Z* is always a real quantity.
  • Z* is given in units of in VASP.
  • VASP outputs Zij with i for the macroscopic electric field, and j for the direction of the force. In literature, Zji is commonly seen, with the force direction j followed by the electric field direction i. Note, py4vasp follows the latter notation Zji for historic reasons.

How to calculate

There are two ways of computing Born effective charges in VASP. The first is done using LCALCEPS, where a finite electric field is applied along the three cartesian directions and the resultant forces on the atoms are calculated:

LCALCEPS = .TRUE.

The other approach is done using LEPSILON, where the derivative of the wavefunction with respect to an electric field is calculated using density functional perturbation theory (DFPT):

LEPSILON = .TRUE.

These may be used in combination with IBRION to obtain additional dielectric properties:

IBRION = 5 or 6 ! Calculated using finite differences.
IBRION = 7 or 8 ! Calculated using DFPT

For more details, see the pages for each tag. The Born effective charges including local field effects will be given in the OUTCAR file:

BORN EFFECTIVE CHARGES (including local field effects) (in |e|, cummulative output)

Excluding local field effects

Previously, the local field effects have been included, that is changes in the orbitals due to the electric field induce changes in the Hartree- and exchange-correlation potentials. This may be limited to changes in the Hartree potential, by specifying:

LRPA = .TRUE.
LCALCEPS = .TRUE. ! N.B. LEPSILON does not output the final Born effective charges.

This prints out the Born effective charges excluding local field effects:

BORN EFFECTIVE CHARGES (excluding local field effects) (in |e|, cummulative output)

These are calculated normally but remain hidden unless explicitly specified.

Related tags and articles

LEPSILON, LCALCEPS, IBRION, LRPA

References