Dear experts,
I am confused regarding the steps when we use ICHARG=11.
This is used to plot the band structure.
This means we are starting the calculation from a preconverged
CHGCAR file and keeping the density fixed throughout the calculation.
If we know the density, we can calculate the potential (at least upto
GGA level), and after solving the KS equation in one step, we can get
the eigenvalue at any kpoint we want. But, when running with ICHARG=11,
there are some scf like steps in the output. What are these steps?
With regards,
Bikash Patra
TIFR, Mumbai
Steps in case of ICHAR=11.
Moderators: Global Moderator, Moderator

 Global Moderator
 Posts: 341
 Joined: Fri Nov 08, 2019 7:18 am
Re: Steps in case of ICHAR=11.
This is a common misunderstanding about the way that VASP works. If you would use an exact diagonalization (ALGO = Exact), you could indeed stop after a single diagonalization and obtain converged eigenvalues. However, VASP obtains the orbitals with an iterative diagonalizer, i.e., there are two selfconsistencies, one with respect to the orbitals and one with respect to the density. You need multiple steps to get good enough orbitals so that your subspace diagonalization yields the correct eigenvalues. By the way, this is also the reason that for delaying the density mixing (NELMDL) in the beginning of every calculation.