LSUBROT: Difference between revisions

From VASP Wiki
No edit summary
 
(7 intermediate revisions by 4 users not shown)
Line 1: Line 1:
{TAGDEF|LSUBROT|.FALSE. {{!}} .TRUE. }
{{TAGDEF|LSUBROT|.FALSE. {{!}} .TRUE. }}


{{DEF|SUBROT|.FALSE.|}}
{{DEF|SUBROT|.FALSE.|}}


Description: {{TAG|LSUBROT}} determines whether an optimal rotation matrix between the occupied and unoccupied block is sought or not.
Description: This flag can be set for [[:Category:Hybrid_functionals|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"/>. For {{TAG|LSUBROT}} =.TRUE. the rotation matrix is instead optimized by performing 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).
 
== References ==
<references>
<ref name="kresse:prb:96">[http://link.aps.org/doi/10.1103/PhysRevB.54.11169 G. Kresse and J. Furthmüller, Phys. Rev. B 54, 11169 (1996).]</ref>
</references>
----
 
[[Category:INCAR tag]][[Category:Electronic minimization]]

Latest revision as of 14:43, 8 April 2022

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]. For LSUBROT =.TRUE. the rotation matrix is instead optimized by performing 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).

References