All requests for technical support from the VASP group must be addressed to: vasp.materialphysik@univie.ac.at
Difference between revisions of "LSUBROT"
Line 3: | Line 3: | ||
{{DEF|SUBROT|.FALSE.|}} | {{DEF|SUBROT|.FALSE.|}} | ||
− | Description: This flag can be set for hybrid functionals (HF-type calculations). {{TAG|LSUBROT}} determines whether an optimal rotation matrix between the occupied and unoccupied block is sought, when a direct optimization of the energy functional is performed (i.e. {{TAG| ALGO}}=All {{!}} Damped). The corresponding algorithm is unpublished. {{TAG|LSUBROT}} =.FALSE. is the standard algorithm, in which the rotation matrix between occupied and unoccupied orbitals is determined essentially using Loewdin perturbation theory, as for instance explained in Ref. <ref name="kresse:prb:96"/>. If {{TAG|LSUBROT}} =.TRUE., the rotation matrix is instead optimized. To do this, VASP performs a few standard SCF steps, in which the orbitals are kept fixed, but rotations between the occupied and unoccupied manifold are allowed. Once satisfactory convergence has been reached, the optimized density matrix (rotation matrix between occupied and unoccupied block) is passed back to the direct optimization routine and a rotation along the suggested direction is performed alongside an update of the orbitals. This generally speeds up calculations for small gap systems as well as metals. However, in rare cases, we have observed instabilities, so be careful when selecting {TAG|LSUBROT}} =.TRUE. | + | Description: This flag can be set for hybrid functionals (HF-type calculations). {{TAG|LSUBROT}} determines whether an optimal rotation matrix between the occupied and unoccupied block is sought, when a direct optimization of the energy functional is performed (i.e. {{TAG| ALGO}}=All {{!}} Damped). The corresponding algorithm is unpublished. {{TAG|LSUBROT}} =.FALSE. is the standard algorithm, in which the rotation matrix between occupied and unoccupied orbitals is determined essentially using Loewdin perturbation theory, as for instance explained in Ref. <ref name="kresse:prb:96"/>. If {{TAG|LSUBROT}} =.TRUE., the rotation matrix is instead optimized. To do this, VASP performs a few standard SCF steps, in which the orbitals are kept fixed, but rotations between the occupied and unoccupied manifold are allowed. Once satisfactory convergence has been reached, the optimized density matrix (rotation matrix between occupied and unoccupied block) is passed back to the direct optimization routine and a rotation along the suggested direction is performed alongside an update of the orbitals. This generally speeds up calculations for small gap systems as well as metals. However, in rare cases, we have observed instabilities, so be careful when selecting {{TAG|LSUBROT}} =.TRUE. |
Although the flag can be set for standard functionals, it is only efficient for hybrid functionals (HF-type calculations). | Although the flag can be set for standard functionals, it is only efficient for hybrid functionals (HF-type calculations). |
Revision as of 14:47, 12 August 2019
LSUBROT = .FALSE. | .TRUE.
Default: SUBROT | = .FALSE. |
Description: This flag can be set for hybrid functionals (HF-type calculations). LSUBROT determines whether an optimal rotation matrix between the occupied and unoccupied block is sought, when a direct optimization of the energy functional is performed (i.e. ALGO=All | Damped). The corresponding algorithm is unpublished. LSUBROT =.FALSE. is the standard algorithm, in which the rotation matrix between occupied and unoccupied orbitals is determined essentially using Loewdin perturbation theory, as for instance explained in Ref. ^{[1]}. If LSUBROT =.TRUE., the rotation matrix is instead optimized. To do this, VASP performs a few standard SCF steps, in which the orbitals are kept fixed, but rotations between the occupied and unoccupied manifold are allowed. Once satisfactory convergence has been reached, the optimized density matrix (rotation matrix between occupied and unoccupied block) is passed back to the direct optimization routine and a rotation along the suggested direction is performed alongside an update of the orbitals. This generally speeds up calculations for small gap systems as well as metals. However, in rare cases, we have observed instabilities, so be careful when selecting LSUBROT =.TRUE.
Although the flag can be set for standard functionals, it is only efficient for hybrid functionals (HF-type calculations).