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 ω_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 demselben 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 ω_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 ω_D grundsätzlich etwas niedriger ist ($\omega_D \lessapprox \omega_{{\rm max}}$, mit der unten benutzten Bezeichnung ω_max) als die eines entsprechenden harmonischen Oszillators (siehe Bild) ohne Frequenzbegrenzung.

Ergebnisse des Modells

Debye-Temperaturen verschiedener Materialien

	ΘD in K
Blei	95
Natrium	160
Gold	165
Silber	215
Kupfer	345
Aluminium	428
α-Eisen	464
Chrom	610
Diamant	1850

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 \Theta_\mathrm{D}$ (Θ_D ist die sog. Debye-Temperatur) gilt für den Phononen-Anteil der spezifischen Wärmekapazität:

$C_\mathrm{v} = {12 \pi^4 N k_\mathrm{B} \over 5} \cdot {T^3 \over \Theta_\mathrm{D}^3 },$

wobei $\Theta_D = \frac{\hbar\omega_D}{k_B}$ ist, mit dem reduzierten planckschen Wirkungsquantum $\hbar$ und der Boltzmann-Konstante k_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}{\Theta_\mathrm{D}^3}\propto\frac{1}{c_{\rm eff}^3} :=\frac{1}{3}\left(\frac{2}{c_t^3}+\frac{1}{c_l^3}\right).$

Sie ist somit ein Maß für die Härte des Kristalls.

Im Hochtemperaturbereich, d. h. für $T\gg \Theta_D$, gilt für die Innere Energie die Beziehung U = 3Nk_BT und somit für die Wärmekapazität

$C_\mathrm{v}=3Nk_\mathrm{B}\,$.

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

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(ω) übereinstimmt; das Hochtemperaturverhalten ist deshalb korrekt, weil die Debye-Näherung per constructionem auch die Summenregel

$\int_0^{\omega_{\rm max}} g(\omega )\,{\rm 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 ω jedoch divergiert, muss die Dichte bei einer bestimmten (materialabhängigen) Frequenz ω_max (in der Debye-Näherung: bei ω_max = ω_D) abgeschnitten werden.

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

$U = \int_0^{\omega_{\rm max}} {g(\omega) \hbar \omega \over e^{\hbar \omega \over 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(ω)dω 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(ω) steht, und statt ω_D die exakte Maximalfrequenz. In der Tieftemperaturnäherung benutzt man, dass man in dieser Näherung die obere Integrationsgrenze durch $\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(ω) (die für die Tieftemperaturnäherung explizit benötigt wird) kann im Debye-Modell angegeben werden, wobei ω_max zu ω_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_{\rm max}}\,g(\omega ) {\rm 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

1. ↑ Peter Debye (1884-1966): Nobelpreisträger für Chemie
2. ↑ Für eine ausführliche klassische Herleitung siehe z.B.: Georg Joos, Lehrbuch der theoretischen Physik - 12. Aufl. - Frankfurt am Main : Akademische Verlagsgesellschaft, 1970 - einerseits "Die Debyesche Theorie der spezifischen Wärme fester Körper" S. 566 ff, bzw. anderseits "Das Plancksche Strahlungsgesetz" S. 580 ff.
3. ↑ Weitere Details findet man z. B bei Werner Döring, Einführung in die Theoretische Physik, Bd. 5, §14, Sammlung Göschen, De Gruyter, Berlin 1957