Incondensed matter physics, thedensity of states (DOS) of a system describes the number of allowed modes orstates per unit energy range. The density of states is defined as${\displaystyle D(E)=N(E)/V}$, where${\displaystyle N(E)\delta E}$ is the number of states in the system of volume${\displaystyle V}$ whose energies lie in the range from${\displaystyle E}$ to${\displaystyle E+\delta E}$. It is mathematically represented as a distribution by aprobability density function, and it is generally an average over the space and time domains of the various states occupied by the system. The density of states is directly related to thedispersion relations of the properties of the system. High DOS at a specific energy level means that many states are available for occupation.

Generally, the density of states of matter is continuous. Inisolated systems however, such as atoms or molecules in the gas phase, the density distribution isdiscrete, like aspectral density. Local variations, most often due to distortions of the original system, are often referred to aslocal densities of states (LDOSs).

In quantum mechanical systems, waves, or wave-like particles, can occupy modes or states with wavelengths and propagation directions dictated by the system. For example, in some systems, the interatomic spacing and the atomic charge of a material might allow only electrons of certain wavelengths to exist. In other systems, the crystalline structure of a material might allow waves to propagate in one direction, while suppressing wave propagation in another direction. Often, only specific states are permitted. Thus, it can happen that many states are available for occupation at a specific energy level, while no states are available at other energy levels.

Looking at the density of states of electrons at the band edge between thevalence and conduction bands in a semiconductor, for an electron in the conduction band, an increase of the electron energy makes more states available for occupation. Alternatively, the density of states is discontinuous for an interval of energy, which means that no states are available for electrons to occupy within the band gap of the material. This condition also means that an electron at the conduction band edge must lose at least the band gap energy of the material in order to transition to another state in the valence band.

This determines if the material is aninsulator or ametal in the dimension of the propagation. The result of the number of states in aband is also useful for predicting the conduction properties. For example, in a one dimensional crystalline structure an odd number ofelectrons per atom results in a half-filled top band; there are free electrons at theFermi level resulting in a metal. On the other hand, an even number of electrons exactly fills a whole number of bands, leaving the rest empty. If then the Fermi level lies in an occupied band gap between the highest occupied state and the lowest empty state, the material will be an insulator orsemiconductor.

Depending on the quantum mechanical system, the density of states can be calculated forelectrons,photons, orphonons, and can be given as a function of either energy or thewave vectork. To convert between the DOS as a function of the energy and the DOS as a function of the wave vector, the system-specific energy dispersion relation betweenE andk must be known.

In general, the topological properties of the system such as the band structure, have a major impact on the properties of the density of states. The most well-known systems, likeneutron matter inneutron stars andfree electron gases in metals (examples ofdegenerate matter and aFermi gas), have a 3-dimensionalEuclidean topology. Less familiar systems, liketwo-dimensional electron gases (2DEG) ingraphite layers and thequantum Hall effect system inMOSFET type devices, have a 2-dimensional Euclidean topology. Even less familiar arecarbon nanotubes, thequantum wire andLuttinger liquid with their 1-dimensional topologies. Systems with 1D and 2D topologies are likely to become more common, assuming developments innanotechnology andmaterials science proceed.

The density of states related to volumeV andN countable energy levels is defined as:$${\displaystyle D(E)=\frac{1}{V}\sum _{i=1}^N\delta (E-E({\mathbf{k}}_i)).}$$Because the smallest allowed change of momentum${\displaystyle k}$ for a particle in a box of dimension${\displaystyle d}$ and length${\displaystyle L}$ is${\displaystyle (\mathrm{\Delta }k{)}^d=(2\pi /L{)}^d}$, the volume-related density of states for continuous energy levels is obtained in the limit${\displaystyle L\to \mathrm{\infty }}$ as$${\displaystyle D(E):={\int }_{{\mathbb{R}}^d}\frac{{\mathrm{d}}^{d}k}{(2\pi {)}^d}\cdot \delta (E-E(\mathbf{k})),}$$Here,${\displaystyle d}$ is the spatial dimension of the considered system and${\displaystyle \mathbf{k}}$ the wave vector.

For isotropic one-dimensional systems with parabolic energy dispersion, the density of states is$D_{1D}(E)=\frac{1}{2\pi \hslash }(\frac{2m}{E}{)}^{1/2}$. In two dimensions the density of states is a constant${\displaystyle D_{2D}=\frac{m}{2\pi {\hslash }^2}}$, while in three dimensions it becomes${\displaystyle D_{3D}(E)=\frac{m}{2{\pi }^2{\hslash }^3}(2mE{)}^{1/2}}$.

Equivalently, the density of states can also be understood as the derivative of the microcanonical partition function${\displaystyle Z_m(E)}$ (that is, the total number of states with energy less than${\displaystyle E}$) with respect to the energy:$${\displaystyle D(E)=\frac{1}{V}\cdot \frac{\mathrm{d}{Z}_m(E)}{\mathrm{d}E}.}$$

The number of states with energy${\displaystyle E^{\prime }}$ (degree of degeneracy) is given by:$${\displaystyle g\left(E^{\prime }\right)=\underset{\mathrm{\Delta }E\to 0}{lim}{\int }_{E^{\prime }}^{E^{\prime }+\mathrm{\Delta }E}D(E)\mathrm{d}E=\underset{\mathrm{\Delta }E\to 0}{lim}D\left(E^{\prime }\right)\mathrm{\Delta }E,}$$where the last equality only applies when the mean value theorem for integrals is valid.

There is a large variety of systems and types of states for which DOS calculations can be done.

Some condensed matter systems possess astructuralsymmetry on the microscopic scale which can be exploited to simplify calculation of their densities of states. In spherically symmetric systems, the integrals of functions are one-dimensional because all variables in the calculation depend only on the radial parameter of the dispersion relation.Fluids,glasses andamorphous solids are examples of a symmetric system whosedispersion relations have a rotational symmetry.

Measurements on powders or polycrystalline samples require evaluation and calculation functions and integrals over the wholedomain, most often aBrillouin zone, of the dispersion relations of the system of interest. Sometimes the symmetry of the system is high, which causes the shape of the functions describing the dispersion relations of the system to appear many times over the whole domain of the dispersion relation. In such cases the effort to calculate the DOS can be reduced by a great amount when the calculation is limited to a reduced zone orfundamental domain.^[1] The Brillouin zone of theface-centered cubic lattice (FCC) in the figure on the right has the 48-fold symmetry of thepoint groupO_h with fulloctahedral symmetry. This configuration means that the integration over the whole domain of the Brillouin zone can be reduced to a 48-th part of the whole Brillouin zone. As acrystal structure periodic table shows, there are many elements with a FCC crystal structure, likediamond,silicon andplatinum and their Brillouin zones and dispersion relations have this 48-fold symmetry. Two other familiar crystal structures are the body-centered cubic lattice (BCC) and hexagonal closed packed structures (HCP) with cubic and hexagonal lattices, respectively. The BCC structure has the 24-foldpyritohedral symmetry of the point groupT_h. The HCP structure has the 12-foldprismatic dihedral symmetry of the point groupD_3h. A complete list of symmetry properties of a point group can be found inpoint group character tables.

In general it is easier to calculate a DOS when the symmetry of the system is higher and the number of topological dimensions of the dispersion relation is lower. The DOS of dispersion relations with rotational symmetry can often be calculated analytically. This result is fortunate, since many materials of practical interest, such as steel and silicon, have high symmetry.

Inanisotropic condensed matter systems such as asingle crystal of a compound, the density of states could be different in one crystallographic direction than in another. These causes the anisotropic density of states to be more difficult to visualize, and might require methods such as calculating the DOS for particular points or directions only, or calculating the projected density of states (PDOS) to a particular crystal orientation.

The density of states is dependent upon the dimensional limits of the object itself. In a system described by three orthogonal parameters (3 Dimension), the units of DOS is[Energy]⁻¹[Volume]⁻¹, in a two dimensional system, the units of DOS is[Energy]⁻¹[Area]⁻¹, in a one dimensional system, the units of DOS is[Energy]⁻¹[Length]⁻¹. The referenced volume is the volume ofk-space; the space enclosed by theconstant energy surface of the system derived through adispersion relation that relatesE tok. An example of a 3-dimensionalk-space is given in Fig. 1. It can be seen that the dimensionality of the system confines the momentum of particles inside the system.

The calculation for DOS starts by counting theN allowed states at a certaink that are contained within[k,k + dk] inside the volume of the system. This procedure is done by differentiating the whole k-space volume${\displaystyle {\mathrm{\Omega }}_{n,k}}$ in n-dimensions at an arbitraryk, with respect tok. The volume, area or length in 3, 2 or 1-dimensional sphericalk-spaces are expressed by$${\displaystyle {\mathrm{\Omega }}_n(k)=c_n{k}^n}$$

for an-dimensionalk-space with the topologically determined constants$${\displaystyle c_1=2,c_2=\pi ,c_3=\frac{4\pi }{3}}$$for linear, disk and spherical symmetrical shaped functions in 1, 2 and 3-dimensional Euclideank-spaces respectively.

According to this scheme, the density of wave vector statesN is, through differentiating${\displaystyle {\mathrm{\Omega }}_{n,k}}$ with respect tok, expressed by$${\displaystyle N_n(k)=\frac{\mathrm{d}{\mathrm{\Omega }}_n(k)}{\mathrm{d}k}=n{c}_n{k}^{n-1}}$$

The 1, 2 and 3-dimensional density of wave vector states for a line, disk, or sphere are explicitly written as

One state is large enough to contain particles having wavelength λ. The wavelength is related tok through the relationship.$${\displaystyle k=\frac{2\pi }{\lambda }}$$

In a quantum system the length of λ will depend on a characteristic spacing of the system L that is confining the particles. Finally the density of statesN is multiplied by a factor${\displaystyle s/V_k}$, wheres is a constant degeneracy factor that accounts for internal degrees of freedom due to such physical phenomena as spin or polarization. If no such phenomenon is present then${\displaystyle s=1}$.V_k is the volume in k-space whose wavevectors are smaller than the smallest possible wavevectors decided by the characteristic spacing of the system.

To finish the calculation for DOS find the number of states per unit sample volume at an energy${\displaystyle E}$ inside an interval${\displaystyle [E,E+\mathrm{d}E]}$. The general form of DOS of a system is given as$${\displaystyle D_n\left(E\right)=\frac{\mathrm{d}{\mathrm{\Omega }}_n(E)}{\mathrm{d}E}}$$The scheme sketched so faronly applies tomonotonically rising andspherically symmetric dispersion relations. In general the dispersion relation${\displaystyle E(k)}$ is not spherically symmetric and in many cases it isn't continuously rising either. To expressD as a function ofE theinverse of the dispersion relation${\displaystyle E(k)}$ has to be substituted into the expression of${\displaystyle {\mathrm{\Omega }}_n(k)}$ as a function ofk to get the expression of${\displaystyle {\mathrm{\Omega }}_n(E)}$ as a function of the energy. If the dispersion relation is not spherically symmetric or continuously rising and can't be inverted easily then in most cases the DOS has to be calculated numerically. More detailed derivations are available.^[2][3]

The dispersion relation for electrons in a solid is given by theelectronic band structure.

Thekinetic energy of a particle depends on the magnitude and direction of thewave vectork, the properties of the particle and the environment in which the particle is moving. For example, the kinetic energy of anelectron in aFermi gas is given by$${\displaystyle E=E_0+\frac{{\left(\hslash k\right)}^2}{2m},}$$

wherem is theelectron mass. The dispersion relation is a spherically symmetric parabola and it is continuously rising so the DOS can be calculated easily.

For longitudinalphonons in a string of atoms the dispersion relation of the kinetic energy in a 1-dimensionalk-space, as shown in Figure 2, is given by$${\displaystyle E=2\hslash {\omega }_0\left|\sin \left(\frac{ka}{2}\right)\right|}$$where${\omega }_0=\sqrt{k_{\text{F}}/m}$ is the oscillator frequency,${\displaystyle m}$ the mass of the atoms,${\displaystyle k_{\text{F}}}$ the inter-atomic force constant and${\displaystyle a}$ inter-atomic spacing. For small values of${\displaystyle k\ll \pi /a}$ the dispersion relation is linear:$${\displaystyle E=\hslash {\omega }_{0}ka}$$

When${\displaystyle k\approx \pi /a}$ the energy is$${\displaystyle E=2\hslash {\omega }_0\left|\cos \left(\frac{\pi -ka}{2}\right)\right|}$$

With the transformation${\displaystyle q=k-\pi /a}$ and small${\displaystyle q}$ this relation can be transformed to$${\displaystyle E\approx 2\hslash {\omega }_0\left[1-{\left(\frac{qa}{2}\right)}^2\right]}$$

The two examples mentioned here can be expressed like$${\displaystyle E=E_0+c_k{k}^p}$$

This expression is a kind ofdispersion relation because it interrelates two wave properties and it isisotropic because only the length and not the direction of the wave vector appears in the expression. The magnitude of the wave vector is related to the energy as:$${\displaystyle k={\left(\frac{E-E_0}{c_k}\right)}^{1/p},}$$

Accordingly, the volume of n-dimensionalk-space containing wave vectors smaller thank is:$${\displaystyle {\mathrm{\Omega }}_n(k)=c_n{k}^n}$$

Substitution of the isotropic energy relation gives the volume of occupied states$${\displaystyle {\mathrm{\Omega }}_n(E)=\frac{c_n}{{c_k}^{\frac{n}{p}}}{\left(E-E_0\right)}^{n/p},}$$

Differentiating this volume with respect to the energy gives an expression for the DOS of the isotropic dispersion relation$${\displaystyle D_n\left(E\right)=\frac{\mathrm{d}}{\mathrm{d}E}{\mathrm{\Omega }}_n(E)=\frac{n{c}_n}{p{c_k}^{\frac{n}{p}}}{\left(E-E_0\right)}^{\frac{n}{p}-1}}$$

In the case of a parabolic dispersion relation (p = 2), such as applies to free electrons in a Fermi gas, the resulting density of states,${\displaystyle D_n\left(E\right)}$, for electrons in a n-dimensional systems is

for${\displaystyle E>E_0}$, with${\displaystyle D(E)=0}$ for.

In 1-dimensional systems the DOS diverges at the bottom of the band as${\displaystyle E}$ drops to${\displaystyle E_0}$. In 2-dimensional systems the DOS turns out to be independent of${\displaystyle E}$. Finally for 3-dimensional systems the DOS rises as the square root of the energy.^[4]

Including the prefactor${\displaystyle s/V_k}$, the expression for the 3D DOS is$${\displaystyle N(E)=\frac{V}{2{\pi }^2}{\left(\frac{2m}{{\hslash }^2}\right)}^{\frac{3}{2}}\sqrt{E-E_0},}$$

where${\displaystyle V}$ is the total volume, and${\displaystyle N(E-E_0)}$ includes the 2-fold spin degeneracy.

In the case of a linear relation (p = 1), such as applies tophotons,acoustic phonons, or to some special kinds of electronic bands in a solid, the DOS in 1, 2 and 3 dimensional systems is related to the energy as:

The density of states plays an important role in thekinetic theory of solids. The product of the density of states and theprobability distribution function is the number of occupied states per unit volume at a given energy for a system in thermal equilibrium. This value is widely used to investigate various physical properties of matter. The following are examples, using two common distribution functions, of how applying a distribution function to the density of states can give rise to physical properties.

Fermi–Dirac statistics: The Fermi–Dirac probability distribution function, Fig. 4, is used to find the probability that a fermion occupies a specific quantum state in a system at thermal equilibrium.Fermions are particles which obey thePauli exclusion principle (e.g. electrons, protons, neutrons). The distribution function can be written as$${\displaystyle f_{\mathrm{F}\mathrm{D}}(E)=\frac{1}{\exp \left(\frac{E-\mu }{k_{\mathrm{B}}T}\right)+1}.}$$

${\displaystyle \mu }$ is thechemical potential (also denoted as E_F and called theFermi level whenT=0),${\displaystyle k_{\mathrm{B}}}$ is the Boltzmann constant, and${\displaystyle T}$ is temperature. Fig. 4 illustrates how the product of the Fermi-Dirac distribution function and the three-dimensional density of states for a semiconductor can give insight to physical properties such as carrier concentration and Energy band gaps.

Bose–Einstein statistics: The Bose–Einstein probability distribution function is used to find the probability that a boson occupies a specific quantum state in a system at thermal equilibrium.Bosons are particles which do not obey the Pauli exclusion principle (e.g. phonons and photons). The distribution function can be written as$${\displaystyle f_{\mathrm{B}\mathrm{E}}(E)=\frac{1}{\exp \left(\frac{E-\mu }{k_{\text{B}}T}\right)-1}.}$$

From these two distributions it is possible to calculate properties such as theinternal energy per unit volume${\displaystyle u}$, the number of particles${\displaystyle N}$,specific heat capacity${\displaystyle c}$, andthermal conductivity${\displaystyle k}$. The relationships between these properties and the product of the density of states and the probability distribution, denoting the density of states by${\displaystyle g(E)}$ instead of${\displaystyle D(E)}$, are given by

${\displaystyle d}$ is dimensionality,${\displaystyle \nu }$ is sound velocity and${\displaystyle \mathrm{\Lambda }}$ ismean free path.

The density of states appears in many areas of physics, and helps to explain a number of quantum mechanical phenomena.

Calculating the density of states for small structures shows that the distribution of electrons changes as dimensionality is reduced. Forquantum wires, the DOS for certain energies actually becomes higher than the DOS for bulk semiconductors, and forquantum dots the electrons become quantized to certain energies.

The photon density of states can be manipulated by using periodic structures with length scales on the order of the wavelength of light. Some structures can completely inhibit the propagation of light of certain colors (energies), creating a photonic band gap: the DOS is zero for those photon energies. Other structures can inhibit the propagation of light only in certain directions to create mirrors, waveguides, and cavities. Such periodic structures are known asphotonic crystals.^[5][6][7][8] In nanostructured media the concept oflocal density of states (LDOS) is often more relevant than that of DOS, as the DOS varies considerably from point to point.

Interesting systems are in general complex, for instance compounds, biomolecules, polymers, etc. Because of the complexity of these systems the analytical calculation of the density of states is in most of the cases impossible. Computer simulations offer a set of algorithms to evaluate the density of states with a high accuracy. One of these algorithms is called theWang and Landau algorithm.^[9]

Within the Wang and Landau scheme any previous knowledge of the density of states is required. One proceeds as follows: the cost function (for example the energy) of the system is discretized. Each time the bini is reached one updates a histogram for the density of states,${\displaystyle g(i)}$, by$${\displaystyle g(i)\to g(i)+f}$$wheref is called the modification factor. As soon as each bin in the histogram is visited a certain number of times (10-15), the modification factor is reduced by some criterion, for instance,$${\displaystyle f_{n+1}\to \frac{1}{2}{f}_n}$$wheren denotes then-th update step. The simulation finishes when the modification factor is less than a certain threshold, for instance.

The Wang and Landau algorithm has some advantages over other common algorithms such asmulticanonical simulations andparallel tempering. For example, the density of states is obtained as the main product of the simulation. Additionally, Wang and Landau simulations are completely independent of the temperature. This feature allows to compute the density of states of systems with very rough energy landscape such as proteins.^[10]

Mathematically the density of states is formulated in terms of a tower of covering maps.^[11]

An important feature of the definition of the DOS is that it can be extended to any system. One of its properties are the translationally invariability which means that the density of the states ishomogeneous and it's the same at each point of the system. But this is just a particular case and the LDOS gives a wider description with aheterogeneous density of states through the system.

Local density of states (LDOS) describes a space-resolved density of states. In materials science, for example, this term is useful when interpreting the data from ascanning tunneling microscope (STM), since this method is capable of imaging electron densities of states with atomic resolution. According to crystal structure, this quantity can be predicted by computational methods, as for example withdensity functional theory.

In a local density of states the contribution of each state is weighted by the density of its wave function at the point.${\displaystyle N(E)}$ becomes${\displaystyle n(E,x)}$

$${\displaystyle n(E,x)=\sum _j|{\varphi }_j(x){|}^2\delta (E-{\epsilon }_j)}$$

the factor of${\displaystyle |{\varphi }_j(x){|}^2}$ means that each state contributes more in the regions where the density is high. An average over${\displaystyle x}$ of this expression will restore the usual formula for a DOS. The LDOS is useful in inhomogeneous systems, where${\displaystyle n(E,x)}$ contains more information than${\displaystyle n(E)}$ alone.

For a one-dimensional system with a wall, the sine waves give

$${\displaystyle n_{1D}(E,x)=\frac{2}{\pi \hslash }\sqrt{\frac{2m}{E}}{\sin }^{2}kx}$$

where$k=\sqrt{2mE}/\hslash$.

In a three-dimensional system with${\displaystyle x>0}$ the expression is

$${\displaystyle n_{3D}(E,x)=\left(1-\frac{\sin 2kx}{2kx}\right)n_{3D}(E)}$$

In fact, we can generalise the local density of states further to

$${\displaystyle n(E,x,x^{\prime })=\sum _j{\varphi }_j(x){\varphi }_j^{\ast }(x^{\prime })\delta (E-{\epsilon }_j)}$$

this is called thespectral function and it's a function with each wave function separately in its own variable. In more advanced theory it is connected with the Green's functions and provides a compact representation of some results such asoptical absorption.

LDOS can be used to gain profit into a solid-state device. For example, the figure on the right illustrates LDOS of atransistor as it turns on and off in a ballistic simulation. The LDOS has clear boundary in the source and drain, that corresponds to the location of band edge. In the channel, the DOS is increasing as gate voltage increase and potential barrier goes down.

This sectionrelies excessively onreferences toprimary sources. Please improve this section by addingsecondary or tertiary sources. Find sources: "Density of states" – news ·newspapers ·books ·scholar ·JSTOR(March 2024) (Learn how and when to remove this message)

Inoptics andphotonics, the concept of local density of states refers to the states that can be occupied by a photon. For light it is usually measured by fluorescence methods, near-field scanning methods or by cathodoluminescence techniques. Different photonic structures have different LDOS behaviors with different consequences for spontaneous emission. In photonic crystals, near-zero LDOS are expected, inhibiting spontaneous emission.^[12]Similar LDOS enhancement is also expected in plasmonic cavity.^[13]However, in disordered photonic nanostructures, the LDOS behave differently. They fluctuate spatially with their statistics, and are proportional to the scattering strength of the structures.^[14]In addition, the relationship with themean free path of the scattering is trivial as the LDOS can be still strongly influenced by the short details of strong disorders in the form of a strong Purcell enhancement of the emission.^[15] and finally, for the plasmonic disorder, this effect is much stronger for LDOS fluctuations as it can be observed as a strong near-field localization.^[16]

• Chen, Gang. Nanoscale Energy Transport and Conversion. New York: Oxford, 2005
• Streetman, Ben G. and Sanjay Banerjee. Solid State Electronic Devices. Upper Saddle River, NJ: Prentice Hall, 2000.
• Muller, Richard S. and Theodore I. Kamins. Device Electronics for Integrated Circuits. New York: John Wiley and Sons, 2003.
• Kittel, Charles and Herbert Kroemer. Thermal Physics. New York: W.H. Freeman and Company, 1980
• Sze, Simon M. Physics of Semiconductor Devices. New York: John Wiley and Sons, 1981

• Online lecture:ECE 606 Lecture 8: Density of States by M. Alam
• Scientists shed light on glowing materials How to measure the Photonic LDOS

• Statistical mechanics
• Physical quantities
• Electronic band structures

• Articles with short description
• Short description is different from Wikidata
• Articles lacking reliable references from March 2024
• All articles lacking reliable references