DFT+U: formalism: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
 
(10 intermediate revisions by 2 users not shown)
Line 20: Line 20:
\delta_{s_1 s_2} \delta_{s_3 s_4}
\delta_{s_1 s_2} \delta_{s_3 s_4}
</math>
</math>
:where <math>|m\rangle</math> represents a real spherical harmonics of angular momentum <math>\ell</math>={{TAG|LDAUL}}.
:where <math>|m\rangle</math> represents a real spherical harmonics of angular momentum <math>l</math>={{TAG|LDAUL}}.


:The unscreened electron-electron interaction <math>U_{\gamma_{1}\gamma_{3}\gamma_{2}\gamma_{4}}</math> can be written in terms of the Slater integrals <math>F^0</math>, <math>F^2</math>, <math>F^4</math>, and <math>F^6</math> (<math>f</math> electrons). Using values for the Slater integrals calculated from atomic orbitals, however, would lead to a large overestimation of the true electron-electron interaction, since in solids the Coulomb interaction is screened (especially <math>F^0</math>).
:The unscreened electron-electron interaction <math>U_{\gamma_{1}\gamma_{3}\gamma_{2}\gamma_{4}}</math> can be written in terms of the Slater integrals <math>F^0</math>, <math>F^2</math>, <math>F^4</math>, and <math>F^6</math> (<math>f</math> electrons). Using values for the Slater integrals calculated from atomic orbitals, however, would lead to a large overestimation of the true electron-electron interaction, since in solids the Coulomb interaction is screened (especially <math>F^0</math>).
Line 60: Line 60:
</math>
</math>


:This can be understood as adding a penalty functional to the semilocal total energy expression that forces the [[#occmat|on-site occupancy matrix]] in the direction of idempotency (''i.e.'' forcing orbitals to be either fully empty or fully occupied),
:This can be understood as adding a penalty functional to the semilocal total energy expression that forces the [[#occmat|on-site occupancy matrix]] in the direction of idempotency,
::<math>\hat n^{\sigma} = \hat n^{\sigma} \hat n^{\sigma}</math>.
::<math>\hat n^{\sigma} = \hat n^{\sigma} \hat n^{\sigma}</math>.


Line 66: Line 66:


:'''Note''': in Dudarev's approach the parameters <math>U</math> and <math>J</math> do not enter seperately, only the difference <math>U-J</math> is meaningful.
:'''Note''': in Dudarev's approach the parameters <math>U</math> and <math>J</math> do not enter seperately, only the difference <math>U-J</math> is meaningful.
*{{TAG|LDAUTYPE}}=3: This option is for the calculation of the parameter <math>U</math> using the linear response approach from Ref. {{cite|cococcioni:2005}}. The steps to use this method are shown for [[Calculate U for LSDA+U|the example of NiO]].


*{{TAG|LDAUTYPE}}=4: same as {{TAG|LDAUTYPE}}=1, but without exchange splitting (i.e., the total spin-up plus spin-down occupancy matrix is used). The double-counting term is given by
*{{TAG|LDAUTYPE}}=4: same as {{TAG|LDAUTYPE}}=1, but without exchange splitting (i.e., the total spin-up plus spin-down occupancy matrix is used). The double-counting term is given by
Line 73: Line 75:
\frac{J}{2} \sum_\sigma {\hat n}^\sigma_{\mathrm{tot}}({\hat n}^\sigma_{\mathrm{tot}}-1).
\frac{J}{2} \sum_\sigma {\hat n}^\sigma_{\mathrm{tot}}({\hat n}^\sigma_{\mathrm{tot}}-1).
</math>
</math>
== Related Tags and Sections ==
{{TAG|LDAU}},
{{TAG|LDAUTYPE}},
{{TAG|LDAUL}},
{{TAG|LDAUU}},
{{TAG|LDAUJ}},
{{TAG|LDAUPRINT}},
{{TAG|LMAXMIX}},


== References ==
== References ==
Line 78: Line 89:
----
----


[[Category:INCAR tag]][[Category:Exchange-correlation functionals]][[Category:DFT+U]]
[[Category:Exchange-correlation functionals]][[Category:DFT+U]]

Latest revision as of 15:35, 9 May 2023

DFT+U is a method that was proposed to improve the description of systems with strongly correlated or electrons, like antiferromagnetic NiO for instance, that are usually inaccurately described with the standard LDA and GGA functionals[1]. Several variants of the DFT+U method exist (see Refs. [2][3] for reviews) that differ for instance in the way the double counting term is calculated. Three variants of them are implemented in VASP, whose formalism is briefly summarized below.

  • LDAUTYPE=1: The rotationally invariant DFT+U introduced by Liechtenstein et al.[4]
This particular flavour of DFT+U is of the form
and is determined by the PAW on-site occupancies
and the (unscreened) on-site electron-electron interaction
where represents a real spherical harmonics of angular momentum =LDAUL.
The unscreened electron-electron interaction can be written in terms of the Slater integrals , , , and ( electrons). Using values for the Slater integrals calculated from atomic orbitals, however, would lead to a large overestimation of the true electron-electron interaction, since in solids the Coulomb interaction is screened (especially ).
In practice these integrals are often treated as parameters, i.e., adjusted to reach agreement with experiment for a property like for instance the equilibrium volume, the magnetic moment or the band gap. They are normally specified in terms of the effective on-site Coulomb- and exchange parameters, and (LDAUU and LDAUJ, respectively). and can also be extracted from constrained-DFT calculations[5][6].
These translate into values for the Slater integrals in the following way (as implemented in VASP at the moment):
- -
-
The essence of the DFT+U method consists of the assumption that one may now write the total energy as:
where the Hartree-Fock-like interaction replaces the semilocal on-site due to the fact that one subtracts a double-counting energy , which supposedly equals the on-site semilocal contribution to the total energy,
  • LDAUTYPE=2: The simplified (rotationally invariant) approach to the DFT+U, introduced by Dudarev et al.[7]
This flavour of DFT+U is of the following form:
This can be understood as adding a penalty functional to the semilocal total energy expression that forces the on-site occupancy matrix in the direction of idempotency,
.
Real matrices are only idempotent when their eigenvalues are either 1 or 0, which for an occupancy matrix translates to either fully occupied or fully unoccupied levels.
Note: in Dudarev's approach the parameters and do not enter seperately, only the difference is meaningful.
  • LDAUTYPE=3: This option is for the calculation of the parameter using the linear response approach from Ref. [8]. The steps to use this method are shown for the example of NiO.
  • LDAUTYPE=4: same as LDAUTYPE=1, but without exchange splitting (i.e., the total spin-up plus spin-down occupancy matrix is used). The double-counting term is given by

Related Tags and Sections

LDAU, LDAUTYPE, LDAUL, LDAUU, LDAUJ, LDAUPRINT, LMAXMIX,

References