NATURALO: Difference between revisions

Latest revision as of 15:12, 29 November 2024

NATURALO = [integer]
Default: NATURALO = 0; for low scaling GW type calculations the default is 4

Description: calculate RPA natural orbitals.

This flag should be used in combination with ALGO = G0W0R or ALGO = scGW0R. The VASP code diagonalizes the RPA density matrix and writes the final natural orbitals to the WAVECAR file. The one-electron occupancies on the WAVECAR file can also be updated to the eigenvalues of the RPA density matrix. For ALGO = G0W0R, the interacting Green's function is approximated as

${\displaystyle G = G_0 + G_0 \Sigma G_0 }$

whereas for ALGO = scGW0R the Dyson equation is solved

${\displaystyle G = G_0 + G \Sigma G_0. }$

In both cases, the RPA density matrix is determined as $\gamma= \lim_{\tau \to 0^-} G(\tau)$ . More details on the use of RPA natural orbitals can be found in Ref. ^[1].

The following settings are currently supported

NATURALO=0 calculate the density matrix, diagonalize the matrix and write the natural orbitals and eigenvalues of the density matrix to WAVECAR. The eigenvalues of the density matrix are stored in the occupancy entries of the WAVECAR file, whereas the one-electron DFT eigenvalues remain untouched. This setting is usually not particularly useful for practical calculations and should be used only by experts. Furthermore, LFINITE_TEMPERATURE should not be combined with this setting.

NATURALO=1 calculate the density matrix, diagonalize the matrix only in the sub-block of unoccupied states, and write the occupied Kohn Sham orbitals, as well as the natural orbitals corresponding to unoccupied states to the file WAVECAR. The unoccupied orbitals are ordered according to their occupancies in the RPA density matrix. The one-electron occupancies and KS-DFT eigenvalues are not updated from their KS values (the occupancies will remain 1 for occupied Kohn-Sham orbitals and 0 for natural orbitals representing the virtual manifold). This setting has been used in Ref. ^[1]. See also Ref. ^[2] for further information. Note that all orbitals- even those with a tiny fractional occupancy -are treated as occupied orbitals and not updated: the algorithm should hence even work for metallic systems.

NATURALO<0. Similar to NATURALO=1 but additionally conserves ABS|NATURALO| unoccupied Kohn-Sham states. This is expedient, for subsequent GW and BSE calculations to conserve few unoccupied orbitals to their Kohn-Sham states.

If 10 is added (e.g. NATURALO=10, NATURALO=11) the density matrix is diagonalizes using a perturbative Loewdin algorithm that attempts to keep the orbital order strictly conserved: E.g. the natural orbital matching closest to each Kohn-Sham orbital will be determined and stored. Use this tag for metals.

NATURALO=2 (or 12) is similar to 0, but the one-electron occupancies are not updated. In rare cases this might lead to inconsistencies, if the orbital order changes between DFT and the RPA density matrix (i.e. a previously occupied DFT orbitals posses a smaller occupation in the RPA density matrix than some unoccupied Kohn-Sham orbitals and are moved into the unoccupied block). This problem can be reduced using NATURALO=12, as described above. This flag, in combination with ALGO = scGW0R, can be used to evaluate the GW-singles contribution to the correlation energy.^[3] One can deduct the HF singles and the GW singles energies from the energies after

    Energies after diagonalization of HF Hamiltonian (single shot)
    Hartree-Fock free energy of the ion-electron system (eV)

    Energies after update of density matrix
    Hartree-Fock free energy of the ion-electron system (eV)

Experience has shown that there is very little difference between the natural orbitals obtained using ALGO = G0W0R and ALGO = scGW0R. We strongly recommend to use the more efficient and better tested algorithm ALGO = G0W0R (with the exception of GW-singles) to determine natural orbitals. Furthermore, perform careful tests for NOMEGA: the RPA total energy converges much faster then the natural orbitals. Using a too small NOMEGA can yield natural orbitals that are non-optimal, leading to very slow convergence of correlated calculations with respect to the number of natural orbitals. A crucial test is that the following line

correlated contrib. to density matrix         0.0000004037        0.0000000000

in the stdout and OUTCAR file shows values close to zero (for ALGO = G0W0R). The above value is perfectly acceptable and the value decreases as NOMEGA increases.

NATURALO=4 preserves original (DFT) orbitals but updates the eigenvalues in the WAVECAR file to the QP energies. This mode is useful if only quasi-particle energies are corrected with the GW method, for instance when selecting ALGO=EVGW0R.

Mind: available as of VASP.6.4:

References

[ramberger2019rpa-1] B. Ramberger, Z. Sukurma, T. Schäfer, G. Kresse, J. Chem. Phys. 151, 214106 (2019).

[GruneisNO-2] A. Grüneis et al. J. Chem. Theory Comput. 7, 2780 (2011).

[Klimessingles-3] Jiří Klimeš et al., J. Chem. Phys. 143, 102816 (2015).

[1]

[2]

[3]

@@ Line 1: / Line 1: @@
-{{TAGDEF|NATURALO|[integer]|0}}
+{{TAGDEF|NATURALO|[integer]|0; for low scaling GW type calculations the default is 4}}
 Description: calculate RPA natural orbitals.
 ----
-If {{TAG|ALGO}} = G0W0R or {{TAG|ALGO}} = GWR, the VASP code diagonalizes the RPA density
+This flag should be used in combination with {{TAG|ALGO}} = G0W0R or {{TAG|ALGO}} = scGW0R.  The VASP code diagonalizes the RPA density
 matrix and writes the final natural orbitals to the {{TAG|WAVECAR}} file.
-The one-electron occupancies on the {{TAG|WAVECAR}} file are also updated to the eigevalues
+The one-electron occupancies on the {{TAG|WAVECAR}} file can also be updated to the eigenvalues of the RPA density matrix. For {{TAG|ALGO}} = G0W0R, the interacting Green's function is approximated as
-of the RPA density matrix.
+<math>
+G = G_0 + G_0 \Sigma G_0
+</math>
+whereas for  {{TAG|ALGO}} = scGW0R the Dyson equation is solved
+<math>
+G = G_0 + G \Sigma G_0.
+</math>
+In both cases, the RPA density matrix is determined as <math> \gamma= \lim_{\tau \to 0^-} G(\tau) </math>. More details on the use of RPA natural orbitals can be found in Ref. <ref name="ramberger2019rpa"/>.
 The following settings are currently supported
-* {{TAG|NATURALO}}=0 calculates the density matrix, diagonalizes the matrix and writes the natural orbitals and eigenvalues of the density matrix to the file {{TAG|WAVECAR}}.
+* {{TAG|NATURALO}}=0 calculate the density matrix, diagonalize the matrix and write the natural orbitals and eigenvalues of the density matrix to {{TAG|WAVECAR}}. The eigenvalues of the density matrix are stored in the occupancy entries of the {{TAG|WAVECAR}} file, whereas the one-electron DFT eigenvalues remain untouched. This setting is usually not particularly useful for practical calculations and should be used only by experts. Furthermore, {{TAG|LFINITE_TEMPERATURE}} should not be combined with this setting.
+* {{TAG|NATURALO}}=1 calculate the density matrix, diagonalize the matrix only in the sub-block of unoccupied states, and write the occupied Kohn Sham orbitals, as well as the natural orbitals corresponding to unoccupied states to the file {{TAG|WAVECAR}}. The unoccupied orbitals are ordered according to their occupancies in the RPA density matrix. The one-electron occupancies and KS-DFT eigenvalues are not updated from their KS values (the occupancies will remain 1 for occupied Kohn-Sham orbitals and 0 for natural orbitals representing the virtual manifold). This setting has been used in Ref. <ref name="ramberger2019rpa"/>. See also Ref. <ref name="GruneisNO"/> for further information. Note that all orbitals- even those with a tiny fractional occupancy -are treated as occupied orbitals and not updated: the algorithm should hence even work for metallic systems.
-* {{TAG|NATURALO}}=1 calculates the density matrix, diagonalizes the matrix only in the sub-block of unoccupied states, and writes the occupied Kohn Sham as well as natural orbitals corresponding to unoccupied states to the file {{TAG|WAVECAR}}. This should be only used for gapped systems (insulators and semiconductors). The one-electron occupancies are not updated (should remain 1 and 0).
+* {{TAG|NATURALO}}<0. Similar to  {{TAG|NATURALO}}=1 but additionally conserves ABS|{{TAG|NATURALO}}| unoccupied Kohn-Sham states. This is expedient, for subsequent GW and BSE calculations to conserve  few unoccupied orbitals to their Kohn-Sham states.
-* {{TAG|NATURALO}}=negative value: Similar to  {{TAG|NATURALO}}=1 but additionally conserve |{{TAG|NATURALO}}| unoccupied states. This is expedient, for subsequent GW calculations, to conserve few unoccupied orbitals to their Kohn-Sham state.
+* If 10 is added (e.g. {{TAG|NATURALO}}=10, {{TAG|NATURALO}}=11) the density matrix is diagonalizes using a perturbative Loewdin algorithm that attempts to keep the orbital order strictly conserved: E.g. the natural orbital matching closest to each Kohn-Sham orbital will be determined and stored. Use this tag for metals.
-* If 10 is added (e.g. {{TAG|NATURALO}}=10, {{TAG|NATURALO}}=12, {{TAG|NATURALO}}=11) the density matrix is diagonalizes using a perturbative Loewdin algorithm that attempts to keep the orbital order strictly conserved: E.g. orbital, the natural orbital matching closest the Nth Kohn-Sham orbital will be stored.
+* {{TAG|NATURALO}}=2 (or 12) is similar to 0, but the one-electron occupancies are not updated. In rare cases this might lead to inconsistencies, if the orbital order changes between DFT and the RPA density matrix (i.e. a previously occupied DFT orbitals posses a smaller occupation in the RPA density matrix than some unoccupied Kohn-Sham orbitals and are moved into the unoccupied block). This problem can be reduced using {{TAG|NATURALO}}=12, as described above. This flag, in combination with {{TAG|ALGO}} = scGW0R, can be used to evaluate the GW-singles contribution to the correlation energy.<ref name="Klimessingles"/> One can deduct the HF singles and the GW singles energies from the energies after
-* {{TAG|NATURALO}}=2 (or 12) is similar to 0, but the one-electron occupancies are not updated. In rare cases this might lead to inconsistencies, if the orbital order changes between DFT and the RPA density matrix (i.e. a previously occupied DFT orbitals posses a small occupation in the RPA density matrix and is moved into the unoccupied block). Use this with care and only for gapped systems. This flag in combination with {{TAG|ALGO}} = GWR, can be used to evaluate the  GW-singles contribution to the correlation energy. One can deduct the singles energies from the energies after
+     Energies after diagonalization of HF Hamiltonian (single shot)
+     Hartree-Fock free energy of the ion-electron system (eV)
- Energies after update of density matrix
+     Energies after update of density matrix
- Hartree-Fock free energy of the ion-electron system (eV)
+     Hartree-Fock free energy of the ion-electron system (eV)
-Experience has shown that there is very little difference between the orbitals obtained using {{TAG|ALGO}} = G0W0R and {{TAG|ALGO}} = GWR. We strongly recommend to use the more efficient and better tested algorithm {{TAG|ALGO}} = G0W0R. Furthermore, perform careful tests for {{TAG|NOMEGA}}: the RPA total energy converges much faster then the natural orbitals. Using a too small {{TAG|NOMEGA}} can yields natural orbitals that are non-optimal, leading to very slow convergence of correlated calculations with respect to the number of natural orbitals.
+Experience has shown that there is very little difference between the natural orbitals obtained using {{TAG|ALGO}} = G0W0R and {{TAG|ALGO}} = scGW0R. We strongly recommend to use the more efficient and better tested algorithm {{TAG|ALGO}} = G0W0R (with the exception of GW-singles) to determine natural orbitals. Furthermore, perform careful tests for {{TAG|NOMEGA}}: the RPA total energy converges much faster then the natural orbitals. Using a too small {{TAG|NOMEGA}} can yield natural orbitals that are non-optimal, leading to very slow convergence of correlated calculations with respect to the number of natural orbitals.
+A crucial test is that the following line
+ correlated contrib. to density matrix         0.0000004037        0.0000000000
+in the stdout and OUTCAR file shows values close to zero (for {{TAG|ALGO}} = G0W0R). The above value is perfectly acceptable and the value decreases as {{TAG|NOMEGA}} increases.
+* {{TAG|NATURALO}}=4 preserves original (DFT) orbitals but updates the eigenvalues in the {{TAG|WAVECAR}} file to the QP energies. This mode is useful if only quasi-particle energies are corrected with the GW method, for instance when selecting {{TAG|ALGO}}=[[GW calculations#gw0|EVGW0R]].
+{{NB|mind|available as of VASP.6.4:}}
-== Related Tags and Sections ==
+== Related tags and articles ==
 * {{TAG|ALGO}} for response functions and ''RPA'' calculations
-* for an overview [[ACFDT/RPA calculations]]
+* for an overview on total energies using the [[ACFDT/RPA calculations|ACFDT/RPA formalism]]
+* for a practical guide to [[GW calculations]]
+== References ==
+<references>
+<ref  name="ramberger2019rpa">[https://doi.org/10.1063/1.5128415   B. Ramberger, Z. Sukurma, T. Schäfer, G. Kresse,  J. Chem. Phys. 151, 214106 (2019).]</ref>
+<ref name="Klimessingles">[https://doi.org/10.1063/1.4929346 Jiří Klimeš et al., J. Chem. Phys.  143, 102816 (2015).]</ref>
+<ref name="GruneisNO"> [https://doi.org/10.1021/ct200263g A. Grüneis et al. J. Chem. Theory Comput. 7, 2780 (2011).] </ref>
+</references>
 ----
-[[Category:INCAR]][[Category:Many-Body Perturbation Theory]] [[Category:GW]] [[Category:ACFDT]][[Category:Performance]][[Category:Parallelization]][[Category:Low-scaling GW and RPA]][[Category:VASP6]]
+[[Category:INCAR tag]][[Category:Many-body perturbation theory]][[Category:GW]] [[Category:ACFDT]][[Category:Low-scaling GW and RPA]]