Density Of States

To calculate the density of States different switches and modes are present in the inp-file:

Note, that this is in Hartree units and not referred to some Fermi-level. The DOS is broadened by a Gaussian with the width sig_dos (also in htr.).

Note that you have to delete the fl7para and kpts files if you want to calculate the DOS.

Producing Raw-data For A DOS

  1. Setting dos=t and ndir=0 in the inp-file switches into the DOS-mode of FLEUR. In this mode a dosinp-file is created which contains the eigenvalues and thier relative weights for all kpts.
  2. Similarly, one obtains a vacDOS file containing the LDOS in the vacuum by setting vacdos=t.
  3. By setting dos=t and ndir>0 one obtains a dosinp-file with additional symmetry-information which is used for band-structure plotting

IntegratedDOSProgramm

total & partial DOS

By setting dos=t and ndir<0 the integrated DOS program is used. Three additional parameters in the inp file have to be specified in this case. The two energies emin_dos and emax_dos specify the energy-range in which the DOS is evaluated. The parameter sig_dos specifies the broadening, as shown below (where the last few lines of the input file are given.

vacdos=T,layers= 1,integ=F,star=T,nstars= 2
 50
 iplot=F,score=F,plpot=F
 0  0.000000  0.000000,nnne=  0,pallst=F
 xa=   2.00000,thetad= 300.00000,epsdisp=   0.00010,epsforce=   0.00010
 relax 000
 emin_dos=    -1.000,emax_dos=   0.00000,sig_dos=   0.00050

Three different negative values for ndir are implemented at the moment:

For information about the output see the documentation of the DOS.x file.

vacuum DOS

If one wishes to calculate the DOS in selected parts of the vacuum region, the parameter vacdos=T has to be set. In combination with dos=T and ndir=-1, this will generally produce a file VACDOS.[1,2] which has a format as describec in the documentation of the DOS.x file.

This is done by integrating over the 2-D unit cell and in z-direction one can either integrate (integ=T) up several layers (defined by layers= in the inp-file), or choose planes, where the DOS is evaluated.

i.e. in the sample inp the layer # 50 is at 11 a.u. or about 3.7 Angstroem above the topmost Cu atom. The format is (in FORTRAN) (20(i3,1x)), i.e. a maximum of 20 planes can be defined. By default, the largest number (N) that can be entered is 250.

Notice, that there is a parameter layerd in the fl7para file. Typically, this should be increased to the number of layers.

Orbital Decomposition Of The DOS

In cases it might be useful to obtain also the orbital decomposed density of states of a certain atom. To get a DOS.x file with this information, you have to specify

dos=T

and

ndir=-3

in the inp file. Additionally, a file named orbcomp has to be present in the working directory, where the number of an atom (counted in the order as it appears in the inp file) is specified. Then, the DOS-file contains information about the selected atom (not atom-type!) in the following form (the energy is in eV):

energy  total-DOS                     selected atom-type
                   s   p   p   p   d    d    d    d      d   f   f   f   f   ...
                       x   y   z   xy   yz   zx   x²-y²  z²  x³  y³  z³  x²y

(Note: you have to provide k-points in the whole Brillouin-zone to get a correct DOS here! Proper symmetrization will be implemented in future releases.) The other f-components can be found in the file ek_orco.[1,2], where also the orbital composition is given for each k-point. To get this file,

cdinf=T

should be set in the inp-file.

rotated coordinate frames

Since the decomposition in spherical harmonics is per default performed in the global coordinate frame, it can be sometimes useful to have a tool to perform the decomposition in a rotated coordinate frame. For this, another file orbcomprot has to be specified that contains three lines, specifying the Euler-angles \alpha, \beta and \gamma for the rotation.

layer resolved DOS

If a zero is specified in the orbcomp file, or no such file has been provided, then the DOS.x files will contain the information about the layer resolved DOS. This is especially meaningful for a film-calculation, where the film can be divided into layers; the boundaries of the layers are defined as (x,y)-planes that intersects halfway between two neighbouring atoms of different height in the film. The results are listed in the following form:

energy  total-DOS   layer            interstitial         muffin-tin 
                                contribution to layer
                     1, 2, ... n      1, 2, ... n          1, 2, ... n