Density of States
1.1 Introduction
The density of states (DOS) is essentially the number of different
states at a particular energy level that electrons are allowed to
occupy, i.e. the number of electron states per unit volume per unit
energy. Bulk properties such as specific heat,
paramagnetic susceptibility, and other transport phenomena of
conductive solids depend on this function. DOS calculations allow one to
determine the general distribution of states as a function of energy
and can also determine the spacing between energy bands in
semi-conductors [1].
1.2 Density of States for Waves
Before we get involved in the derivation of the DOS of electrons in a
material, it may be easier to first consider just an elastic wave
propagating through a solid. Elastic waves are in reference to the
lattice vibrations of a solid comprised of discrete atoms. Though, when
the wavelength is very long, the atomic nature of the solid can be
ignored and we can treat the material as a continuous medium [2].
We begin with the 1-D wave equation:
With which we then have a solution for a propagating plane wave:
(1)
= wave number
, =amplitude, = frequency, = velocity of sound
Dispersion Relation:
In equation(1), the temporal factor, 
can be omitted because it is not relevant to the derivation of the DOS[2]. So now we will use the solution:
(2)
To begin, we must apply some type of boundary conditions to the system. The easiest way to do this is to consider a periodic boundary condition. With a periodic boundary condition we can imagine our system having two ends, one being the origin, 0, and the other, . We now say that the origin end is constrained in a way that it is always at the same state of oscillation as end L[2].
This boundary condition is represented as:
Now we apply the boundary condition to equaiton(2) to get:
Now, using Euler’s identity; 
we can see that there are certain values of which satisfy the above equation. Those values are 
for any integer, . Leaving the relation:
If you choose integer values for and plot them along an axis you get a 1-D line of points, known as modes, with a spacing of 
between each mode.
We now have that the number of modes in an interval in -space equals:
Using the dispersion relation we can find the number of modes within a frequency range that lies within
. This number of modes in that range is represented by , where is defined as the density of states.
So now we see that =
which we turn into:
We do so in order to use the relation:
and obtain:
we multiply by a factor of two be cause there are modes in positive and negative -space, and we get the density of states for a phonon in 1-D:
2-D:
We can now derive the density of states for two dimensions. Equation(2) becomes:
now apply the same boundary conditions as in the 1-D case:
We now consider an area for each point in -space = and find the number of modes that lie within a flat ring with thickness , a radius and area:
Number of modes inside interval:
Now account for transverse and longitudinal modes (multiply by a factor of 2) and set equal to We get 
and apply dispersion
relation to get 
which simplifies to the 2-D result:
3-D:
We can now derive the density of states for three dimensions. Equation(2) becomes:
now apply the same boundary conditions as in the 1-D case to get:
We now consider a volume for each point in -space = and find the number of modes that lie within a spherical shell, thickness , with a radius and volume:
Number of modes inside shell:
Assuming a common velocity for transverse and longitudinal waves we can
account for one longitudinal and two transverse modes for each value of (multiply by a factor of 3) and set equal to :
Apply dispersion relation and let to get 
and simplify for the 3-D result:
1.3 Density of States for Electrons
Now that we have seen the distribution of modes for waves in a
continuous medium, we move to electrons. The calculation of some
electronic processes like absorption, emission, and the general
distribution of electrons in a material require us to know the number of
available states per unit volume per unit energy. The density of states
is once again represented by a function which this time is a function of energy and has the relation = the number of states per unit volume in the energy range: .
We begin by observing our system as a free electron gas confined to points
contained within the surface. We do this so that the electrons in our
system are free to travel around the crystal without being influenced by
the potential of atomic nuclei[3]. Using the Schrödinger wave equation we can determine that the solution of electrons confined in a box with rigid walls, i.e. the Particle in a box problem, gives rise to standing waves for which the allowed values of are expressible in terms of three nonzero integers, [1].
Taking a step back, we look at the free electron, which has a momentum, and velocity, , related by . Since the energy of a free electron is entirely kinetic we can disregard the potential energy term and state that the energy, =
Using De-Broglie’s particle-wave duality theory we can assume that
the electron has wave-like properties and assign the electron a wave
number :
h is the reduced Planck’s constant:
Now we substitute terms:
Substitute term into the equation for energy: 
We are now left with the dispersion relation for electron energy: 
Where is the effective mass of an electron.
As for the case of a phonon which we discussed earlier, the equation for allowed values of is found by solving the Schrödinger wave equation with the same boundary conditions that we used earlier.
We are left with the solution:
And after applying the same boundary conditions used earlier:
We can consider each position in -space being filled with a cubic unit cell volume of:
= making the number of allowed values per unit volume of -space :
We can picture the allowed values from as a sphere near the origin with a radius and thickness . The allowed states are now found within the volume contained between and , see Fig.1.
Fig.1 [3]. Spherical shell showing values of as points. The points contained within the shell and are the allowed values.
The volume of the shell with radius and thickness can be calculated by simply multiplying the surface area of the sphere, 
, by the thickness, :
Now we can form an expression for the number of states in the shell by combining the number of allowed states per unit volume of -space with the volume of the spherical shell seen in Fig.1.
Number of states:
Substitute in the dispersion relation for electron energy:
which leads to 
now substitute the expressions obtained for and in terms of back into the expression for the number of states:
We have now represented the electrons in a 3 dimensional -space, similar to our representation of the elastic waves in -space, except this time the shell in -space has its surfaces defined by the energy contours and , thus the number of allowed values within this shell gives the number of available states and when divided by the shell thickness, , we obtain the function [2].
we must now account for the fact that any state can contain two electrons, spin-up and spin-down, so we multiply by a factor of two to get:
The above expression for the DOS is valid only for the region in -space where the dispersion relation applies. As the energy increases the contours described by become non-spherical, and when the energies are large enough the shell will intersect the boundaries of the first Brillouin zone, causing the shell volume to decrease which leads to a decrease in the number of states. If the volume continues to decrease, goes to zero and the shell no longer lies within the zone. The energy at which becomes zero is the location of the top of the valance band and the range from where remains zero is the band gap[2].
Fig.2 [1]. The left hand side shows a two-band diagram and a DOS vs. plot for no band overlap. The right hand side shows a two-band diagram and a DOS vs. plot for the case when there is a band overlap.
In simple metals the DOS can be calculated for most of the energy band, using:
however when we reach energies near the top of the band we must use a
slightly different equation. To derive this equation we can consider
that the next band is ev below the minimum of the first band[1]. This is illustrated in the upper left plot in Fig.2. The energy of this second band is:
. Now we can derive the density of states in this region in the same
way that we did for the rest of the band and get the result:
Solving for the DOS in the other dimensions will be similar to what we did for the waves. i.e. for 2-D we would consider an area element in -space , and for 1-D a line element in -space . Fig.3 lists the equations for the density of states in 4 dimensions, (a quantum dot would be considered 0-D), along with corresponding plots of DOS vs. energy. Fig.4 plots DOS vs. energy over a range of values for each dimension and super-imposes the curves over each other to further visualize the different behavior between dimensions.
Fig.3 [4]
Fig.4 [3]
References
[1] Sachs, M., Solid State Theory, (New York, McGraw-Hill Book Company, 1963),pp159-160;238-242.
[2] Omar, Ali M., Elementary Solid State Physics, (Pearson Education, 1999), pp68- 75;213-215.
[3] Semiconductor Physics (online) http://britneyspears.ac/physics/dos/dos.htm
[4] Density of States (online) www.ecse.rpi.edu/~schubert/Course-ECSE-6968%20Quantum%20mechanics/Ch12%20Density%20of%20states.pdf
Good article keep it up
ReplyDeleteYou can also reach out to :- https://www.mastersindia.co/blog/gst-portal-login/