next up previous contents index
Next: PROCAR file Up: Files used by VASP Previous: EIGENVALUE file   Contents   Index

N.B. This document is no longer maintained, please visit our wiki.


The file DOSCAR contains the DOS and integrated DOS. The units are ``number of states/unit cell''. For dynamic simulations and relaxations, an averaged DOS and an averaged integrated DOS is written to the file. For a description of how the averaging is done see 6.21, 6.37). The first few lines of the DOSCAR file are made up by a header which is followed by NDOS lines holding three data

 energy     DOS     integrated DOS
The density of states (DOS) $ \bar n$, is actually determined as the difference of the integrated DOS between two pins, i.e.

$\displaystyle \bar n(\epsilon_i) = (N(\epsilon_i) - N(\epsilon_{i-1})) / \Delta \epsilon,

where $ \Delta \epsilon$ is the distance between two pins (energy difference between two grid point in the DOSCAR file), and $ N(\epsilon_i)$ is the integrated DOS

$\displaystyle N (\epsilon_{i}) = \int_{-\infty}^{\epsilon_i} n(\epsilon) d \epsilon.

This method conserves the total number of electrons exactly. For spin-polarized calculations each line holds five data
 energy     DOS(up) DOS(dwn)  integrated DOS(up) integrated DOS(dwn)
If RWIGS or LORBIT (Wigner Seitz radii, see section 6.336.34) is set in the INCAR file, a lm- and site-projected DOS is calculated and also written to the file DOSCAR. One set of data is written for each ion, each set of data holds NDOS lines with the following data
 energy s-DOS p-DOS d-DOS
 energy s-DOS(up) s-DOS(down) p-DOS(up) p-DOS(dwn) d-DOS(up) d-DOS(dwn)
for the non spin-polarized and spin polarized case respectively. As before the written densities are understood as the difference of the integrated DOS between two pins.

For non-collinear calculations, the total DOS has the following format:

 energy     DOS(total)   integrated-DOS(total)
Information on the individual spin components is available only for the site projected density of states, which has the format:
 energy s-DOS(total) s-DOS(mx) s-DOS(my) s-DOS(mz) p-DOS(total) p-DOS(mx) ...
In this case, the (site projected) total density of states (total) and the (site projected) energy resolved magnetization density in the $ x$ (mx), $ y$ (my) and $ z$ (mz) direction are available.

In all cases, the units of the l- and site projected DOS are states/atom/energy.

The site projected DOS is not evaluated in the parallel version for the following cases:

vasp.4.5, NPAR$ \ne$1 no site projected DOS
vasp.4.6, NPAR$ \ne$1, LORBIT=0-5 no site projected DOS
In vasp.4.6 the site projected DOS can be evaluated for LORBIT=10-12, even if NPAR is not equal 1 (contrary to previous releases).

Mind: For relaxations, the DOSCAR is usually useless. If you want to get an accurate DOS for the final configuration, first copy CONTCAR to POSCAR and continue with one static (ISTART=1; NSW=0) calculation.

N.B. Requests for support are to be addressed to: