This page explains how to calculate electronic density of states [1] using density functional theory. We study silicon in the standard diamond-centered cubic structure and use Quantum ESPRESSO [2] as our simulation engine during this tutorial.

Accuracy of the results

Please note that this calculation is performed using Density Functional Theory and generalized gradient approximation [3] which is known to under-estimate the energy of unoccupied electronic states.

# Create job¶

Si is the default material, so if you choose "Create a Job" from the sidebar on the home page Si will automatically be loaded.

# Choose workflow¶

Density of states is calculated in conjunction with band structure calculation so under workflow chose the "Bandstructure + DOS" for either Quantum Espresso.

Click on the animation below to view:

## Adjust kpoints¶

It is critical to have a high sampling in reciprocal cell (high k-point density) to give enough detail to calculate an accurate density of states.

In QuantumEspresso, the band structure + DOS workflow has 5 units. The first unit specifies the settings for the self-consistent calculation of the eigenvalues and wave functions. The second unit calculation is a non self-consitent calculation using the wave functions and charge density of the previous calculation. Subsequent calculations calculate the density of states and also the projection of those states for partial density of states analysis

We set the kpoint density to 11x11x11 in the first workflow unit to provide sufficient density for the second non-consistent calculation step of the band structure.

# Submit job¶

Before submitting the calculation, click on the "Go to Compute" button and examine the compute parameters. This is a small structure so 1 core and 5 minutes are sufficient. Click "No" when it asks if you want to save a duplicate material.

# Monitor status¶

As each unit in the workflow is executing, you can monitor its progress live by viewing both the output of the executable as well as a graphical representation of the total energy convergence on the Status tab under each execution's sub-tab.

# Examine results¶

When all 5 units are complete, switching to the Results tab and the sub-tab for the final execution unit will show the density of states as well as the partial density of states due to each atom in the structure as well as their s and p electron-like character. By moving your cursor along each data series it will highlight which element electronic character that data series corresponds to.

Partial contributions

The numbers represent the order of the current orbital as included inside the pseudopotential, and __not__ the principal quantum number.