Self-consistency cycle: Difference between revisions

Latest revision as of 21:27, 13 November 2023

The term self-consistency cycle (SCC) denotes a category of algorithms that determine the electronic ground state by a combination of iterative matrix diagonalization and density mixing.

Figure 1. shows a procedural flowchart of the self-consistency cycle:

When starting from scratch, the SCC starts with an initial guess for the electronic density of the system under consideration: VASP uses the approximation of overlapping atomic charge densities. The orbitals are initialized with random numbers. Alternatively, the SCC may (re-)start from the orbitals and/or electronic density obtained in a previous calculation.
The density defines the Hamiltonian.
By means of iterative matrix diagonalization techniques, one obtains the NBANDS lowest lying eigenstates of the Hamiltonian. The iterative matrix diagonalization algorithms implemented in VASP are the blocked-Davidson algorithm and the residual-minimization method with direct inversion in the iterative subspace (RMM-DIIS). Per default VASP uses the blocked-Davidson algorithm (ALGO = Normal). This step is often referred to as iterative optimization/refinement of the orbitals.
After the eigenstates and eigenvalues of the Hamiltonian, i.e., orbitals and one-electron energies, have been determined with sufficient accuracy, the corresponding partial occupancies of the orbitals are calculated.
From the one-electron energies and partial occupancies, the free energy of the system is computed.
From the orbitals and partial occupancies, a new electronic density is constructed.
In principle, the new density could be directly used to define a new Hamiltonian. In most cases, however, this does not lead to a stable algorithm (on account of, e.g., charge sloshing). Instead, the new density is not used directly but is mixed with the old density. By default VASP uses a Broyden mixer (IMIX). The resulting density then defines the new Hamiltonian for the next round of iterative matrix diagonalization.

Steps 2-7 are repeated until the change in the free energy from one cycle to the next drops below a specific threshold (EDIFF).

Note that when starting from scratch (ISTART = 0), the self-consistency cycle procedure of VASP always begins with several (NELMDL) cycles where the density is kept fixed at the initial approximation (overlapping atomic charge densities). This ensures that the wavefunctions that are initialized with random numbers have converged to something sensible before they are used to construct a new charge density.

Especially in case of the RMM-DIIS the initial set of orbitals plays a critical role. Therefore, either the number of non-selfconsistent cycles is chosen to be large (NELMDL = 12, for ALGO = VeryFast), or the non-selfconsistent cycles are done with the blocked-Davidson algorithm before switching over to the use of the RMM-DIIS (ALGO = Fast).

Furthermore, note that per default (LDIAG=.TRUE.) the iterative refinement of the orbitals is preceded (RMM-DIIS) or followed (blocked-Davidson) by a diagonalization of the subspace spanned by the current orbitals. In case of the RMM-DIIS, the optimization step is additionally followed by an orthogonalization of the refined orbitals.

With respect to the aforementioned it should be emphasized that, in principle, the RMM-DIIS method should also converge without any explicit subspace diagonalization and/or re-orthonormalization. However, in our experience their inclusion speeds up the convergence of the self-consistency cycle so substantially that it shortens the time-to-solution of most calculations, even though these operations scale as $O(N^3)$ .

For (a lot) more details on the self-consistency cycle and associated algorithms in VASP, we recommend the seminal papers by Kresse and Furthmüller.^[1]^[2]

References

[kresse:cms:1996-1] G. Kresse and J. Furthmüller, Comp. Mater. Sci. 6, 15 (1996)

[kresse:prb:96-2] G. Kresse and J. Furthmüller, Phys. Rev. B 54, 11169 (1996).

[1]

[2]

@@ Line 1: / Line 1: @@
-The term "Self-Consistency Cycle" (SCC) denotes a category of algorithms that calculate the electronic ground state by a combination of iterative matrix diagonalisation and density mixing.
+[[File:SCC.png|350px|thumb|Fig. 1: the self-consistency cycle]]
-# The SCC starts with an initial guess for the electronic density of the system under consideration: {{VASP}} uses the approximation of overlapping atomic charge densities. This density defines the initial Hamiltonian. The orbitals are initialised with random numbers.
+The term '''self-consistency cycle''' (SCC) denotes a category of algorithms that determine the electronic ground state by a combination of iterative matrix diagonalization and [[Density mixing|density mixing]].
-# By means of iterative matrix diagonalisation techniques one obtains the ''N'' lowest lying eigenstates of the Hamiltonian (where ''N'' is of the order of the number of electrons in the unit cell). The iterative matrix diagonalisation algorithms implemented in {{VASP}} are the blocked-Davidson algorithm and the Residual Minimization Method with Direct Inversion in the Iterative Subspace (RMM-DIIS). Per default {{VASP}} uses the blocked-Davidson algorithm ({{TAG|ALGO}} = Normal).
-# After the eigenstates and eigenvalues have been the determined with sufficient accuracy, they are used to to compute the total energy of the system, and to construct a new electronic density.
-# In principle, this new density could be taken to define a new Hamiltonian. However, in order to arrive at a stable algorithm this new density is not used as is, but is mixed with the old density. By default {{VASP}} uses a Broyden mixer. The resulting density then defines the new Hamiltonian for the next round of iterative metrix diagonalisation (step 2).
-Steps 2-4 are repeated until the change in the total energy from one cycle to the next drops below a specific threshold ({{TAG|EDIFF}}).
+Figure 1. shows a procedural flowchart of the self-consistency cycle:
+# When starting from scratch, the SCC starts with an initial guess for the electronic density of the system under consideration: {{VASP}} uses the approximation of overlapping atomic charge densities. The orbitals are initialized with random numbers. Alternatively, the SCC may (re-)start from the orbitals and/or electronic density obtained in a previous calculation.
+# The density defines the Hamiltonian.
+# By means of iterative matrix diagonalization techniques, one obtains the {{TAG|NBANDS}} lowest lying eigenstates of the Hamiltonian. The iterative matrix diagonalization algorithms implemented in {{VASP}} are the [[Blocked-Davidson algorithm|blocked-Davidson algorithm]] and the [[RMM-DIIS|residual-minimization method with direct inversion in the iterative subspace (RMM-DIIS)]]. Per default {{VASP}} uses the [[Blocked-Davidson algorithm|blocked-Davidson algorithm]] ({{TAG|ALGO}} = Normal). This step is often referred to as ''iterative optimization/refinement of the orbitals''.
+# After the eigenstates and eigenvalues of the Hamiltonian, i.e., orbitals and one-electron energies, have been determined with sufficient accuracy, the corresponding partial occupancies of the orbitals are calculated.
+# From the one-electron energies and partial occupancies, the free energy of the system is computed.
+# From the orbitals and partial occupancies, a new electronic density is constructed.
+# In principle, the new density could be directly used to define a new Hamiltonian. In most cases, however, this does not lead to a stable algorithm (on account of, e.g., charge sloshing). Instead, the new density is not used directly but is mixed with the old density. By default {{VASP}} uses a Broyden mixer ({{TAG|IMIX}}). The resulting density then defines the new Hamiltonian for the next round of iterative matrix diagonalization.
+Steps 2-7 are repeated until the change in the free energy from one cycle to the next drops below a specific threshold ({{TAG|EDIFF}}).
 Note that when starting from scratch ({{TAG|ISTART}} = 0), the self-consistency cycle procedure of {{VASP}} always begins with several ({{TAG|NELMDL}}) cycles where the density is kept fixed at the initial approximation (overlapping atomic charge densities).
-This ensures that the wave functions that are initialised with random numbers have converged to a something sensible before they are used to construct a new charge density.
+This ensures that the wavefunctions that are initialized with random numbers have converged to something sensible before they are used to construct a new charge density.
+Especially in case of the [[RMM-DIIS]] the initial set of orbitals plays a critical role.
+Therefore, either the number of non-selfconsistent cycles is chosen to be large ({{TAG|NELMDL}} = 12, for {{TAG|ALGO}} = VeryFast), or the non-selfconsistent cycles are done with the [[Blocked-Davidson algorithm|blocked-Davidson algorithm]] before switching over to the use of the [[RMM-DIIS]] ({{TAG|ALGO}} = Fast).
+Furthermore, note that per default ({{TAG|LDIAG}}=.TRUE.) the iterative refinement of the orbitals is preceded ([[RMM-DIIS]]) or followed ([[Blocked-Davidson algorithm|blocked-Davidson]]) by a diagonalization of the subspace spanned by the current orbitals.
+In case of the [[RMM-DIIS]], the optimization step is additionally followed by an orthogonalization of the refined orbitals.
+With respect to the aforementioned it should be emphasized that, in principle, the [[RMM-DIIS]] method should also converge without any explicit subspace diagonalization and/or re-orthonormalization. However, in our experience their inclusion speeds up the convergence of the [[Self-consistency_cycle|self-consistency cycle]] so substantially that it shortens the time-to-solution of most calculations, even though these operations scale as <math>O(N^3)</math>.
+For (a lot) more details on the self-consistency cycle and associated algorithms in {{VASP}}, we recommend the seminal papers by Kresse and Furthmüller.{{cite|kresse:cms:1996}}{{cite|kresse:prb:96}}
+<!--
+----
 The following section discusses the minimization
@@ Line 73: / Line 92: @@
 This residual vector is added to the wavefunction <math>\phi_n</math>, the algorithms differ in the way this is exactly done.
+-->
 <!--
 Alternately, the ground-state energy can be evaluated with the [[RPA/ACFDT: Correlation energy in the Random Phase Approximation| random phase approximation]].