How do I plot a total or partial density of states (DOS)?
Set tasks to 10 and run the code. This will produce a total and partial spin-resolved density of states. Note that to improve the quality of the DOS, you can first run one self-consistent loop (set maxscl to 1 and tasks to 1) with a larger k-point set. This will not update the density in STATE.OUT but will produce a new set of eigenvalues and Fermi energy. Compute the DOS as usual afterwards. See also the dos block in the Users’ Guide.