LDAUTYPE: Difference between revisions

Latest revision as of 16:22, 27 February 2025

LDAUTYPE = 1 | 2 | 4
Default: LDAUTYPE = 2

Description: LDAUTYPE specifies the DFT+U variant that will be used.

The following variants of the DFT+U are available:

LDAUTYPE=1: The rotationally invariant DFT+U introduced by Liechtenstein et al.^[1]

LDAUTYPE=2: The simplified (rotationally invariant) approach to DFT+U, introduced by Dudarev et al.^[2]

LDAUTYPE=3: Linear response ansatz of Cococcioni et al. ^[3] to compute U. See how to calculate U.

Mind: For LDAUTYPE=3, the LDAUU and LDAUJ tags specify the strength of the spherical potential acting on the spin-up and spin-down manifolds, respectively.

LDAUTYPE=4: Same as LDAUTYPE=1, but without exchange splitting.

A method to estimate the parameters for DFT+U is the constrained-random-phase approximation. Another method is the linear response ansatz with LDAUTYPE=3, mentioned above. On the other hand, in many applications, the DFT+U parameters are used as tuning parameters to fit experimental data.

Tip: For band-structure calculations, increase LMAXMIX to 4 ([math]\displaystyle{ d }[/math] elements) or 6 ([math]\displaystyle{ f }[/math] elements).

This is because the CHGCAR file contains only information up to angular momentum quantum number set by LMAXMIX for the on-site PAW occupancy matrices. When the CHGCAR file is read and kept fixed in the course of the calculations (ICHARG=11), the results will necessarily not be identical to a self-consistent run. The deviations are often large for DFT+U calculations.

Warning: The total energy will depend on the parameters [math]\displaystyle{ U }[/math] (LDAUU) and [math]\displaystyle{ J }[/math] (LDAUJ). It is, therefore, not meaningful to compare the total energies resulting from calculations with different [math]\displaystyle{ U }[/math] and/or [math]\displaystyle{ J }[/math]; or [math]\displaystyle{ U-J }[/math] in the case of Dudarev's approach (LDAUTYPE=2).

It is possible to use LDAUTYPE=1, 2, and 3 for a non–spin-polarized calculation with ISPIN=1.

References

[liechtenstein:prb:95-1] A. I. Liechtenstein, V. I. Anisimov, and J. Zaanen, Phys. Rev. B 52, R5467 (1995).

[dudarev:prb:98-2] S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys, and A. P. Sutton, Phys. Rev. B 57, 1505 (1998).

[cococcioni:2005-3] M. Cococcioni and S. de Gironcoli, Phys. Rev. B 71, 035105 (2005).

[1]

[2]

[3]

@@ Line 1: / Line 1: @@
 {{TAGDEF|LDAUTYPE|1 {{!}} 2 {{!}} 4|2}}
-Description: {{TAG|LDAUTYPE}} specifies which type of L(S)DA+U approach will be used.
+Description: {{TAG|LDAUTYPE}} specifies the DFT+U variant that will be used.
 ----
-*{{TAG|LDAUTYPE}}=1: The rotationally invariant LSDA+U introduced by Liechtenstein ''et al.''<ref name="liechtenstein:prb:95"/>
+The following variants of the [[:Category:DFT+U|DFT+U]] are available:
-*{{TAG|LDAUTYPE}}=2: The simplified (rotationally invariant) approach to the LSDA+U, introduced by Dudarev ''et al.''<ref name="dudarev:prb:98"/>
+*{{TAG|LDAUTYPE}}=1: The rotationally invariant DFT+U introduced by Liechtenstein ''et al.''{{cite|liechtenstein:prb:95}}
-*{{TAG|LDAUTYPE}}=4: same as {{TAG|LDAUTYPE}}=1, but LDA+U instead of LSDA+U (i.e. no LSDA exchange splitting).
+*{{TAG|LDAUTYPE}}=2: The simplified (rotationally invariant) approach to DFT+U, introduced by Dudarev ''et al.''{{cite|dudarev:prb:98}}
-== Related Tags and Sections ==
+*{{TAG|LDAUTYPE}}=3: Linear response ansatz of Cococcioni et al. {{cite|cococcioni:2005}} to compute U. See [[Calculate U for LSDA+U|how to calculate U]].
+{{NB|mind|For {{TAG|LDAUTYPE}}{{=}}3, the {{TAG|LDAUU}} and {{TAG|LDAUJ}} tags specify the strength of the spherical potential acting on the spin-up and spin-down manifolds, respectively.|:}}
+*{{TAG|LDAUTYPE}}=4: Same as {{TAG|LDAUTYPE}}=1, but without exchange splitting.
+A method to estimate the parameters for DFT+U is the [[Constrained-random-phase approximation|constrained-random-phase approximation]]. Another method is the linear response ansatz with {{TAG|LDAUTYPE}}=3, mentioned above. On the other hand, in many applications, the DFT+U parameters are used as tuning parameters to fit experimental data.
+{{NB|tip|For band-structure calculations, increase {{TAG|LMAXMIX}} to 4 (<math>d</math> elements) or 6 (<math>f</math> elements).}}
+This is because the {{FILE|CHGCAR}} file contains only information up to angular momentum quantum number set by {{TAG|LMAXMIX}} for the [[LDAUTYPE#occmat|on-site PAW occupancy matrices]]. When the {{FILE|CHGCAR}} file is read and kept fixed in the course of the calculations ({{TAG|ICHARG}}=11), the results will necessarily not be identical to a self-consistent run. The deviations are often large for DFT+U calculations.
+{{NB|warning|The total energy will depend on the parameters <math>U</math> ({{TAG|LDAUU}}) and <math>J</math> ({{TAG|LDAUJ}}). It is, therefore, not meaningful to compare the total energies resulting from calculations with different <math>U</math> and/or <math>J</math>; or <math>U-J</math> in the case of Dudarev's approach ({{TAG|LDAUTYPE}}{{=}}2).}}
+It is possible to use {{TAG|LDAUTYPE}}=1, 2, and 3 for a non–spin-polarized calculation with {{TAG|ISPIN}}=1.
+== Related tags and articles ==
 {{TAG|LDAU}},
 {{TAG|LDAUL}},
 {{TAG|LDAUU}},
 {{TAG|LDAUJ}},
-{{TAG|LDAUPRINT}}
+{{TAG|LDAUPRINT}},
+{{TAG|LMAXMIX}}
+{{sc|LDAUTYPE|Examples|Examples that use this tag}}
 == References ==
-<references>
+<references/>
-<ref name="liechtenstein:prb:95">[[http://link.aps.org/doi/10.1103/PhysRevB.52.R5467 A. I. Liechtenstein, V. I. Anisimov and J. Zaane, Phys. Rev. B 52, R5467 (1995).]]</ref>
-<ref name="dudarev:prb:98">[[http://link.aps.org/doi/10.1103/PhysRevB.57.1505 S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys and A. P. Sutton, Phys. Rev. B 57, 1505 (1998).]]</ref>
-</references>
 ----
-[[The_VASP_Manual|Contents]]
-[[Category:INCAR]][[Category:LDA+U]]
+[[Category:INCAR tag]][[Category:Exchange-correlation functionals]][[Category:DFT+U]]