Debye-Modell

Das Debye-Modell beschreibt eine Methode zur Berechnung des Beitrags zur Wärmekapazität eines kristallinen Festkörpers, den die quantisierten Schwingungen in Kristallgittern, die sog. Phononen, ergeben; es stellt sich u. a. heraus, dass dies in der Regel der wesentliche Beitrag ist. Die von Peter Debye 1911 und 1912 entwickelte „Theorie der spezifischen Wärme von Kristallen“ gilt überdies als eine der ersten theoretischen Bestätigungen der 1900 von Max Planck erstmals vorgestellten Quantenthese^[1].

Grundlagen des Modells

Dispersionsrelation im Vergleich mit dem Ergebnis von einfachen harmonischen Oszillatoren

Gegenüber dem Einstein-Modell (von 1906), welches $N$ unabhängige Oszillatoren mit identischer Frequenz annimmt, geht das Debye-Modell (1912) von einer Vielzahl möglicher Frequenzen und von Null verschiedener Ausbreitungsgeschwindigkeit aller Wellen bzw. Phononen aus. Jedoch wird durchgehend die Langwellennäherung vorausgesetzt, d.h. es wird der Einfachheit halber angenommen, dass bis zu einer Grenzfrequenz, der sog. Debyefrequenz ${\omega }_D$, immer strenge Proportionalität zwischen Frequenz und Wellenvektor (also eine lineare Dispersionsrelation) gilt, wobei ein longitudinaler und zwei transversale Schallwellen-Freiheitsgrade vorausgesetzt werden.

Bemerkenswert an diesem Ansatz ist, dass er (abgesehen von der Nichtexistenz longitudinaler Lichtwellen) mit den Annahmen Plancks zur Berechnung der Hohlraumstrahlung identisch ist, wenn man die Schallgeschwindigkeit durch die Lichtgeschwindigkeit ersetzt. Somit ergeben sich für einen strahlenden Hohlraum (→Plancksches Strahlungsgesetz) Formeln mit dem selben Aufbau, wie für einen erwärmten Festkörper, bei dem Teilchen in gitterförmiger Anordnung schwingen. In beiden Fällen folgen charakteristische „T³-Gesetze“.^[2] (s. u.).

Phononen existieren aber nur bis zu einer Maximalfrequenz (im Debye-Modell also bis zu ${\omega }_D$). Diese Maximalfrequenz ergibt sich aus der Summe aller möglichen Schwingungsmoden, da deren Gesamtzahl höchstens gleich dem Dreifachen der Anzahl der schwingenden Gitterteilchen (Atome) sein kann. Daraus folgt auch, dass ${\omega }_D$ grundsätzlich etwas niedriger ist (${\omega }_D⪅{\omega }_{\mathrm{m}\mathrm{a}\mathrm{x}}$, mit der unten benutzten Bezeichnung ${\omega }_{\mathrm{m}\mathrm{a}\mathrm{x}}$) als die eines entsprechenden harmonischen Oszillators (siehe Bild) ohne Frequenzbegrenzung.

Ergebnisse des Modells

Temperaturbereiche

Das Modell macht korrekte Voraussagen über die Temperaturabhängigkeit der Wärmekapazität, und zwar sowohl im Niedrigtemperaturlimes als auch für den Hochtemperaturgrenzwert.

Im Niedrigtemperaturbereich, d. h. für $T\ll {\mathrm{\Theta }}_{\mathrm{D}}$ (${\mathrm{\Theta }}_{\mathrm{D}}$ ist die sog. Debye-Temperatur) gilt für den Phononen-Anteil der spezifischen Wärmekapazität:

$C_{\mathrm{v}}=\frac{12{\pi }^{4}N{k}_{\mathrm{B}}}{5}\cdot \frac{T^3}{{\mathrm{\Theta }}_{\mathrm{D}}^3},$

wobei ${\mathrm{\Theta }}_D=\frac{\hslash {\omega }_D}{k_B}$ ist, mit dem reduzierten planckschen Wirkungsquantum $\hslash$ und der Boltzmann-Konstante $k_{\mathrm{B}}$.

Die Debye-Temperatur ist dabei proportional zu einer effektiven Schallgeschwindigkeit, zu der die transversalen Schallwellen zu 2/3 und die longitudinalen Schallwellen zu 1/3 beitragen; und zwar gilt präzise

$\frac{1}{{\mathrm{\Theta }}_{\mathrm{D}}^3}\propto \frac{1}{c_{\mathrm{e}\mathrm{f}\mathrm{f}}^3}:=\frac{1}{3}(\frac{2}{c_t^3}+\frac{1}{c_l^3}).$

Im Hochtemperaturbereich, d. h. für $T\gg {\mathrm{\Theta }}_D$, gilt für die Innere Energie die Beziehung $U=3N{k}_{\mathrm{B}}T$ und somit für die Wärmekapazität

$C_{\mathrm{v}}=3N{k}_{\mathrm{B}}$.

In diesem Limes ergibt sich also, wie schon beim Einstein-Modell, das Gesetz von Dulong-Petit.

Die Zustandsdichte gemäß dem Debye-Modell ergibt sich aus folgender Rechnung: $g(\omega )d\omega =g(k)dk$. Oder umgeschrieben: $g(\omega )=g(k)\cdot \frac{d(k)}{d(\omega )}$. Nun gilt aber allgemein im $k$-Raum $g(k)=\frac{L}{\pi }$ und nach dem Debye-Modell ist $\omega =v_s\cdot k$, also $\frac{d(k)}{d(\omega )}=\frac{1}{v_s}$ und damit insgesamt: $g(\omega )=\frac{L}{\pi }\cdot \frac{1}{v_s}.$

In beiden Grenzfällen ist die Theorie völlig korrekt, während das intermediäre Verhalten durch die Debye-Theorie nur im Sinne einer „vernünftigen Interpolation“ beschrieben wird, die man gegebenenfalls verbessern kann (s. u.). Und zwar ist das Tieftemperaturverhalten deshalb korrekt, weil im Limes $\omega \ll {\omega }_D$ die Debye-Näherung mit dem exakten $g(\omega )$ übereinstimmt; das Hochtemperaturverhalten ist deshalb korrekt, weil die Debye-Näherung per constructionem auch die Summenregel

${\int }_0^{{\omega }_{\mathrm{m}\mathrm{a}\mathrm{x}}}g(\omega )\mathrm{d}\omega \equiv 3N$

erfüllt.

Zur Begründung der angegebenen Resultate

Das Debye-Modell nähert die Dispersionsrelation von Phononen in der angegebenen Weise linear an. Die Berechnung, die auch für den (realistischen!) Fall elementar durchgeführt werden kann, dass longitudinale und transversale Schallgeschwindigkeit sich erheblich unterscheiden, dauert lange, so dass Details hier nur aus Platzgründen unterbleiben^[3].

Da in einem Festkörper höchstens dreimal so viele Schwingungsmoden wie Atome vorhanden sein können, die Zustandsdichte für hohe $\omega$ jedoch divergiert, muss die Dichte bei einer bestimmten (materialabhängigen) Frequenz ${\omega }_{\mathrm{m}\mathrm{a}\mathrm{x}}$ (in der Debye-Näherung: bei ${\omega }_{\mathrm{m}\mathrm{a}\mathrm{x}}={\omega }_D$) abgeschnitten werden.

Ausgehend von der exakten Formel für die Schwingungsenergie U,

$U={\int }_0^{{\omega }_{\mathrm{m}\mathrm{a}\mathrm{x}}}\frac{g(\omega )\hslash \omega }{e^{\frac{\hslash \omega }{k_{B}T}}-1}\mathrm{d}\omega$

ergibt sich obige Wärmekapazität $C_v$ explizit durch Ausführung des Integrals und Differentiation nach der Temperatur T,

$C_v=\frac{\partial U}{\partial T}.$

Hierin ist $g(\omega )\mathrm{d}\omega$ die Zahl der Schwingungsmoden mit Kreisfrequenzen $\in [\omega ,\omega +\mathrm{d}\omega ].$

(Man beachte, dass oben statt der Debye-Näherung $g_D$ das exakte $g(\omega )$ steht, und statt ${\omega }_D$ die exakte Maximalfrequenz. In der Tieftemperaturnäherung benutzt man, dass man in dieser Näherung die obere Integrationsgrenze durch $\mathrm{\infty }$ ersetzen kann und dass die niedrigsten nicht-trivialen Terme der Taylorentwicklungen von g und $g_D$ bei $\omega \to 0$ übereinstimmen. Für das Hochtemperaturverhalten ersetzt man im Nenner den Term $e^x-1$ einfach durch x und berechnet das verbleibende Integral mit der Summenregel. Die Zustandsdichte $g(\omega )$ (die für die Tieftemperaturnäherung explizit benötigt wird) kann im Debye-Modell angegeben werden, wobei ${\omega }_{\mathrm{m}\mathrm{a}\mathrm{x}}$ zu ${\omega }_D$ wird.

Die konkrete, über die Debye-Näherung hinausgehende Berechnung der Zustandsdichte g, ist allerdings nicht allgemein analytisch, sondern nur numerisch oder für Teile der Temperaturskala genähert lösbar, wie oben für tiefe Temperaturen. Hier liegen auch die oben angedeuteten Verbesserungsmöglichkeiten für das intermediäre Verhalten.

Verallgemeinerung für andere Quasiteilchen

Das Debye'sche Verfahren kann für andere bosonische Quasiteilchen im Festkörper leicht in analoger Weise durchgeführt werden, z. B. in ferromagnetischen Systemen für sog. Magnonen anstelle der Phononen. Man hat jetzt andere Dispersionsrelationen für $\omega \to 0$, z. B. $\omega \propto k^2$ im genannten Fall, und andere Summenregeln, z. B. ${\int }_0^{{\omega }_{\mathrm{m}\mathrm{a}\mathrm{x}}}g(\omega )\mathrm{d}\omega =N.$ Auf diese Weise ergibt sich in Ferromagneten bei tiefen Temperaturen ein Magnonenbeitrag $\propto T^{3/2}$ zur Wärmekapazität, der gegenüber dem Phononenbeitrag , $\propto T^3,$ dominiert. In Metallen dagegen kommt der Hauptbeitrag , $\propto T,$ natürlich von den Elektronen. Er ist fermionisch und wird mit anderen Methoden berechnet, die auf Arnold Sommerfeld zurückgehen.

Einzelnachweise