Zustandsdichte

Die Zustandsdichte D(E) bzw. $D(\omega)$ (engl. density of states, abgekürzt DOS) ist eine physikalische Größe, die angibt, wie viele Zustände pro Energieintervall ${\mathrm {d}}E$ bzw. pro Frequenzintervall ${\mathrm {d}}\omega$ in einem physikalischen System existieren.

Im Allgemeinen wird die Zustandsdichte für Vielteilchensysteme im Rahmen eines Modells unabhängiger Teilchen betrachtet. Dann beziehen sich die Variablen $\omega$ bzw. $E=\hbar \omega$ auf die Energie der 1-Teilchenzustände. Häufig wird die Zustandsdichte dann auch in Abhängigkeit vom Impuls ${\vec {p}}=\hbar {\vec {k}}$ bzw. Wellenvektor ${\vec {k}}$ der 1-Teilchenzustände betrachtet und gibt deren Anzahl pro Volumenintervall des Impulsraums ( $\mathrm{d}^3\vec{p}$ ) bzw. des reziproken Raums ( $\mathrm{d}^3\vec{k}$ ) an. Die Zustandsdichte kann sich auf verschiedene Teilchensorten beziehen, z.B. auf Photonen, Phononen, Elektronen, Magnonen, Quasiteilchen, und wird pro Einheit des räumlichen Volumens angegeben. Für freie Teilchen ohne Spin lässt sich die Zustandsdichte daraus berechnen, dass im Phasenraum jeder quantenmechanische Zustand das Volumen $(2 \pi \hbar)^3$ einnimmt. Die Zustandsdichte (pro Volumen) $D(\vec{k})$ ist dann konstant

$D(\vec{k}) = \frac{1}{(2 \pi)^3} \ .$

Im Falle von Wechselwirkungen der Teilchen, sei es untereinander oder mit vorgegebenen Potentialen, kann die Zustandsdichte stark davon abweichen (siehe z.B. Bändermodell).

Definition

Allgemein ist die – auf das Volumen bezogene – Zustandsdichte für eine abzählbare Anzahl an Energieniveaus definiert durch:

$D(E) = \frac{1}{V} \cdot \sum_{i=1}^N \delta(E - E(\vec{k}_i))$

mit der Delta-Distribution $\delta$ .

Daraus erhält man durch Erweitern mit $(\Delta k)^{d}=\left({\frac {2\pi }{L}}\right)^{d}$ (der kleinsten erlaubten Änderung von k für ein Teilchen in einem Kasten der Dimension d und Länge L) und Übergang zu einem Riemann-Integral (Limes $L\to \infty$ ) die auf das Volumen bezogene Zustandsdichte für kontinuierliche Energiniveaus:

$D(E):=\int _{\mathbb {R} ^{d}}{\frac {\mathrm {d} ^{d}k}{(2\pi )^{d}}}\cdot \delta (E-E({\vec {k}}))\qquad (*)$

mit

der räumlichen Dimension des betrachteten Systems
dem Betrag des Wellenvektors.

Äquivalent kann die Zustandsdichte auch als Ableitung der mikrokanonischen Zustandssumme $Z_{m}(E)=N(E)$ nach der Energie aufgefasst werden:

$D(E)={\frac {1}{V}}\cdot {\frac {\mathrm {d} N(E)}{\mathrm {d} E}}$

Die Zahl der Zustände mit Energie (Entartungsgrad) ist gegeben durch:

$g(E')=\lim _{\Delta E\to 0}\int _{E'}^{E'+\Delta E}D(E)\mathrm {d} E=\lim _{\Delta E\to 0}D(E')\Delta E$ , wobei das letzte Gleichheitszeichen nur dann gilt, wenn der Mittelwertsatz der Integralrechnung für das Integral anwendbar ist.

Anschauung

Anschaulich zählt man die Mikrozustände für eine vorgegebene Energie : betrachtet man ein System mit N Mikrozuständen , so wird die Zustandsdichte beschrieben durch

$D(E)=\sum _{{i=1}}^{N}\delta (E-E_{i})$

da das Integral über die Zustandsdichte gerade die Gesamtzahl der Mikrozustände liefert:

$\int _{\mathbb {R} }\sum _{i=1}^{N}\delta (E-E_{i})\cdot \mathrm {d} E=\sum _{i=1}^{N}=N$

und außerdem liefert folgendes Integral die Zahl n(E) der Mikrozustände bei Energie E:

$\lim _{\Delta E\to 0}\int _{E_{i}-\Delta E}^{E_{i}+\Delta E}\sum _{i=1}^{N}\delta (E-E_{i})\cdot \mathrm {d} E=n(E)$

In obiger Formel (*) ist zumindest für die Anschauung die Eigenschaft $\delta (g(x))=\sum _{{i=1}}^{{n}}{\frac {\delta (x-x_{{i}})}{|g^{{\prime }}(x_{{i}})|}}$ der Deltadistribution wichtig, die jedoch nur für endlich viele und einfache Nullstellen $x_{i}$ von g(x) gilt.

n-dimensionales Elektronengas

Die folgenden Erläuterungen beziehen sich vorrangig auf Anwendungen in der Festkörperphysik.

Zustandsdichte über der Energie abhängig von der Dimension (3D = gepunktet, 2D = rot, 1D = grün, 0D = blau). Die Sprünge in den Zustandsdichten für die Dimensionen D=0 bis D=2 sind darin begründet, dass in diesen Fällen die Zustandsdichten um verschiedene Energiezustände gezeichnet sind. Um diese Energiezustände herum hat die Zustandsdichte dann die berechnete und in der Tabelle dargestellte Form.

In einem n-dimensionalen Elektronengas können sich Ladungsträger in den Dimensionen 1, …, n frei bewegen. Der entsprechende Anteil der Energie ist kontinuierlich und kann unter Nutzung der parabolischen Näherung angegeben werden:

$E={\frac {\hbar ^{2}k^{2}}{2m^{*}}}$

Dabei ist

m* die effektive Masse des Ladungsträgers im Festkörper, genauer ist m* die effektive Zustandsdichtemasse
$\hbar$ das (durch $2\pi$ geteilte) Plancksche Wirkungsquantum.

Im Gegensatz dazu ist die Energiekomponente der anderen Dimensionen diskretisiert in den Werten $E_{l}$ . Die (auf das Volumen V bezogene) Zustandsdichte kann allgemein beschrieben werden:

$D(E)=2\cdot {\frac {{\mathrm {d}}N(E)}{{\mathrm {d}}E}}{\frac {1}{V}}.$

Darin entspricht

der Vorfaktor 2 den zwei möglichen Spinzuständen (oft wird er aber in berücksichtigt, hier wurde dies nicht so gemacht)
$V=L_{x}\cdot L_{y}\cdot L_{z}$ dem Volumen des Festkörpers
N(E) der Anzahl aller Zustände mit Energie kleiner gleich E (vgl.: Mikrokanonische Zustandssumme $Z_{m}$ ):

$N(E)={\begin{cases}{\frac {V_{k}}{\Omega _{k}}}&{\text{wenn}}\quad n=3\\\sum _{l}\Theta (E-E_{l}){\frac {V_{k}}{\Omega _{k}}}&{\text{wenn}}\quad n=1,2\\\sum _{l}\Theta (E-E_{l})&{\text{wenn}}\quad n=0\end{cases}}$

- $V_{k}$ beschreibt im n-dimensionalen k-Raum das Gesamtvolumen aller Zustände, die bei der verbleibenden Energie $E-E_{l}$ zugänglich sind
- $\Omega _{k}$ ist das Volumen eines solchen Zustandes.
- $\Theta$ beschreibt hier die Heaviside-Funktion.

Werte für verschieden-dimensionale Elektronengase
	Gesamtvolumen aller Zustände	Volumen eines Zustandes ${\Omega_k}$	(auf das Volumen bezogene) Zustandsdichte
	im k-Raum bei der verbleibenden Energie $E-E_{l}$		(auf das Volumen bezogene) Zustandsdichte
3D – Bulk	${\frac {4}{3}}\pi k^{3}$	$\frac{(2\pi)^3}{L_x L_y L_z}$	${\frac {(2m^{*})^{{\frac {3}{2}}}}{2\pi ^{2}\hbar ^{3}}}{\sqrt {E}}$
2D – Quantentopf/Quantenfilm	$\pi k^2$	$\frac{(2\pi)^2}{L_x L_y}$	$\frac{m^*} {\pi\hbar^2 L_z} \sum_l \Theta(E-E_l)$
1D – Quantendraht		$\frac{2\pi}{L_x}$	$\frac{\sqrt{2m^*}}{\pi\hbar L_y L_z}\sum_l \frac{1}{\sqrt{E-E_l}}$
0D – Quantenpunkt			$\frac{2}{ L_x L_y L_z}\sum_l \delta (E-E_l)$

Im Halbleiter

Zustandsdichten (farbig) in einem undotierten Halbleiter mit direktem Bandübergang. Zusätzlich ist die Fermi-Verteilung bei Raumtemperatur nach links aufgetragen, als Energieniveaus das Fermi-Niveau E_F und die Leitungsbandenergie E_C.

In Halbleitermaterialien wird wegen der periodisch auftretenden Atomkerne ein ähnlicher Ansatz für das Leitungs- und Valenzband gemacht (siehe Bändermodell). Halbleiter zeichnet aus, dass ihre Dispersionkurven oder auch Bandstruktur ein Maximum (Valenzband) und Minimum (Leitungsband) besitzt, welche nicht überlappen, sondern durch die Bandlücke getrennt sind. Dabei wird bei einer Versetzung der Extrema im k-Raum (Impulsraum) von einem indirekten, bei gleichem Impulsunterschied von einem direkten Halbleiter gesprochen. Das funktionale Verhalten um solche Extremwerte lässt sich parabolisch (quadratisch) nähern. Die Krümmung dieser zur Näherung verwendeten Form muss allerdings nicht mit der Krümmung der Dispersionskurve der oben besprochenen freien Elektronen übereinstimmen. Stattdessen weist man den Ladungsträgern, also Elektronen und Löchern, in den beiden Bändern bei diesen Extrema effektive Massen m* zu so, dass nun die funktionale Beschreibung identisch zu der, der echten freien Elektronen, ist.

Zustandsdichten (farbig) in einem n-dotierten Halbleiter mit direktem Bandübergang. Energieniveau der Dotieratome E_D.

Die Energie der Leitungsband-Unterkante sei $E_{{\mathrm {C}}}$ , die der Valenzband-Oberkante $E_{{\mathrm {V}}}$ , die Differenz ist gleich der Bandlückenenergie $E_{{\mathrm {G}}}=E_{{\mathrm {C}}}-E_{{\mathrm {V}}}$ . Die Zustandsdichte im Leitungsband ist ( $m_{{\mathrm {e,d}}}^{*}$ ist die Zustandsdichtemasse des Elektrons im Leitungsband, also seine gemittelte effektive Masse):

$D_{{\mathrm {C}}}(E)={\frac {(2m_{{\mathrm {e,d}}}^{*})^{{\frac {3}{2}}}}{2\pi ^{2}\hbar ^{3}}}{\sqrt {E-E_{{\mathrm {C}}}}}$

Die Zustandsdichte im Valenzband ist ( $m_{{\mathrm {p,d}}}^{*}$ ist die Zustandsdichtemasse des Lochs im Valenzband):

$D_{{\mathrm {V}}}(E)={\frac {(2m_{{\mathrm {p,d}}}^{*})^{{\frac {3}{2}}}}{2\pi ^{2}\hbar ^{3}}}{\sqrt {E_{{\mathrm {V}}}-E}}$

Bei dotierten Halbleitern treten zu diesen möglichen Zuständen noch Zustände in der Bandlücke auf. Diese sind bei n-Dotierung nahe am Leitungsband und bei p-Dotierung nahe am Valenzband. Durch Zuführen von Energie kann die Aktivierungsenergie überwunden werden und es bilden sich vermehrt besetzte Zustände in Leitungs- bzw. Valenzband. Darüber hinaus ändert sich durch Dotierung die Lage des Fermi-Niveaus: es wird bei n-Dotierung angehoben, bzw. senkt sich bei p-Dotierung zum Valenzband hin ab. Bei einer n-Dotierung sind damit bereits bei Raumtemperatur wegen der thermischen Energie weit mehr Zustände im Leitungsband besetzt als bei einem undotierten Material. Die zusätzlichen freien Ladungsträger können damit den Stromtransport erhöhen.

Die thermische Besetzung der Zustände wird durch die Fermi-Verteilung bestimmt. Die Wahrscheinlichkeitsdichte, dass ein Zustand mit der Energie $[E, E+\mathrm{d}E]$ besetzt ist, schreibt sich

$W_{{\mathrm {e}}}(E)={\frac {1}{\exp {\left({\frac {E-\mu }{k_{{\mathrm {B}}}T}}\right)}+1}}$

Die Wahrscheinlichkeitsdichte, dass ein Zustand mit der Energie $[E, E+\mathrm{d}E]$ nicht besetzt oder äquivalent ausgedrückt mit einem Loch besetzt ist, schreibt sich

$W_{{\mathrm {h}}}(E)=1-W_{{\mathrm {e}}}(E)={\frac {1}{\exp {\left(-{\frac {E-\mu }{k_{{\mathrm {B}}}T}}\right)}+1}}$

Damit lassen sich die Ladungsträgerdichten, also Elektronendichte im Leitungsband und Löcherdichte im Valenzband, angeben:

$n=\int _{{E_{{\mathrm {C}}}}}^{{\infty }}W_{{\mathrm {e}}}(E)\,D_{{\mathrm {C}}}(E)\,{\mathrm {d}}E$

sowie

$p=\int _{{-\infty }}^{{E_{{\mathrm {V}}}}}W_{{\mathrm {h}}}(E)\,D_{{\mathrm {V}}}(E)\,{\mathrm {d}}E$

Eigentlich sollten die Integrationsgrenzen nicht bis unendlich ausgedehnt sein, sondern nur bis zum Ende des jeweiligen Bandes. Allerdings ist dort die Fermi-Verteilung schon näherungsweise Null – das chemische Potential liegt nämlich im Bereich der Bandlücke – sodass der Fehler vernachlässigbar ist. Zur Berechnung dieser Integrale siehe Fermi-Dirac-Integral.

Basierend auf einem Artikel in:

Wikipedia.de