IMIX: Difference between revisions

From VASP Wiki
No edit summary
 
(2 intermediate revisions by 2 users not shown)
Line 1: Line 1:
{{TAGDEF|IMIX|0 {{!}} 1 {{!}} 2 {{!}} 4|4}}
{{TAGDEF|IMIX|0 {{!}} 1 {{!}} 2 {{!}} 4|4}}


Description: {{TAG|IMIX}} specifies the type of [[Density Mixing|density mixing]].
Description: {{TAG|IMIX}} specifies the type of [[:Category:Density mixing|density mixing]].
----
----
=={{TAG|IMIX}}=0: No mixing ==
=={{TAG|IMIX}}=0: No mixing ==
Line 33: Line 33:
:For {{TAG|WC}}=0, VASP uses Broyden's 2<sup>nd</sup> method,<ref name="bluegel:thesis:88"/><ref name="johnson:prb:88"/> and, for {{TAG|WC}}>0, VASP uses Pulay-mixing method<ref name="pulay:cpl:80"/>.
:For {{TAG|WC}}=0, VASP uses Broyden's 2<sup>nd</sup> method,<ref name="bluegel:thesis:88"/><ref name="johnson:prb:88"/> and, for {{TAG|WC}}>0, VASP uses Pulay-mixing method<ref name="pulay:cpl:80"/>.
:The default is a Pulay mixer with an initial approximation for the charge-dielectric function according to Kerker<ref name="kerker:prb:81"/>
:The default is a Pulay mixer with an initial approximation for the charge-dielectric function according to Kerker<ref name="kerker:prb:81"/>
::<math>A\times\min\left(\frac{G^2}{G^2+B^2},A_{\rm min}\right)</math>
::<math>A\times\max\left(\frac{G^2}{G^2+B^2},A_{\rm min}\right)</math>
:where <math>A</math>={{TAG|AMIX}}, <math>B</math>={{TAG|BMIX}}, and <math>A_{\rm min}</math>={{TAG|AMIN}}.
:where <math>A</math>={{TAG|AMIX}}, <math>B</math>={{TAG|BMIX}}, and <math>A_{\rm min}</math>={{TAG|AMIN}}.


Line 62: Line 62:
----
----


[[Category:INCAR tag]][[Category:Electronic minimization]][[Category:Density mixing]]
[[Category:INCAR tag]][[Category:Density mixing]]

Latest revision as of 11:35, 18 October 2024

IMIX = 0 | 1 | 2 | 4
Default: IMIX = 4 

Description: IMIX specifies the type of density mixing.


IMIX=0: No mixing

IMIX=1: Kerker mixing

For Kerker mixing[1], the mixed density is given by
with =AMIX and =BMIX. If BMIX is very small, e.g., BMIX=0.0001, a straight mixing is obtained.
Mind: BMIX=0 might cause floating-point exceptions on some platforms.

IMIX=2: Variant of Tchebycheff mixing

VASP uses a variant of the popular Tchebycheff-mixing scheme[2]. Here, the following second order equation of motion is used:
with =AMIX, =BMIX, and =AMIN. A velocity Verlet algorithm is used to integrate this equation. The discretized equation reads:
where
and
,
where the index N is the electronic iteration, and F is the force acting on the charge.
For BMIX≈0, no model for the dielectric matrix is used. For a simple straight mixing is obtained. Therefore, corresponds to maximal damping, while implies no damping. To determine the optimal parameters for and AMIX, first converge to the ground state with the Pulay mixer (IMIX=4). Then, search for the the eigenvalues of the charge-dielectric matrix in the OUTCAR file at the last occurrence of
eigenvalues of (default mixing * dielectric matrix)
The optimal parameters are then given by:
AMIX
AMIN

IMIX=4: Broyden's 2nd method and Pulay-mixing method (default)

For WC=0, VASP uses Broyden's 2nd method,[3][4] and, for WC>0, VASP uses Pulay-mixing method[5].
The default is a Pulay mixer with an initial approximation for the charge-dielectric function according to Kerker[1]
where =AMIX, =BMIX, and =AMIN.
AMIN=0.4 usually yields good convergence. AMIX strongly depends on the system, for instance, it should be small, e.g., AMIX= 0.02, for metals.
In the Broyden scheme, the functional form of the initial mixing matrix is determined by AMIX and BMIX or the INIMIX tag. The metric used in the Broyden scheme is specified through MIXPRE.

Related tags and articles

INIMIX, MAXMIX, AMIX, BMIX, AMIX_MAG, BMIX_MAG, AMIN, MIXPRE, WC

Examples that use this tag

References