Äquipartitionstheorem

Das Äquipartitionstheorem (auch Gleichverteilungssatz genannt) besagt, dass im thermischen Gleichgewicht bei der Temperatur T im Mittel jeder Freiheitsgrad die gleiche Energie $⟨E⟩$ besitzt:

$⟨E⟩=\frac{1}{2}{k}_{B}T$

Dabei ist $k_B$ die Boltzmann-Konstante. Also gilt für Teilchen mit $f$ Freiheitsgraden:

$⟨E⟩=\frac{f}{2}{k}_{B}T$

Der Gleichverteilungssatz gilt nur für Freiheitsgrade, deren Variable im Ausdruck für die Energie, das heißt in der Hamilton-Funktion, als Quadrat vorkommen. Des Weiteren dürfen diese Freiheitsgrade nicht "eingefroren" sein, das heißt, dieser Freiheitsgrad muss tatsächlich angeregt werden. Beispielsweise werden Molekülschwingungen „kleiner Moleküle“ wie H₂ oder O₂ bei Raumtemperatur nicht angeregt, weil die für den Übergang auf den niedrigsten angeregten Zustand nötige Energie nicht erreicht wird.

Freiheitsgrade, deren Variable nicht in der Hamilton-Funktion vorkommen, führen natürlich zu keinem Beitrag zur Energie; für Freiheitsgrade, die anders als in rein quadratischer Form vorkommen, lässt sich die mittlere Energie nicht so einfach berechnen.

Beispiele

Spezifische Wärme von Gasen

Aus dem Äquipartitionstheorem lässt sich beispielsweise die Wärmekapazität (spezifische Wärme) $C_V$ eines idealen Gases bei konstantem Volumen berechnen. Wir betrachten zuerst ein einatomiges Gas (Edelgas):

Die Energie des Gases ist durch seine kinetische Energie der Atome gegeben; für jedes Atom gilt

$E=\frac{1}{2}m{v}^2=\frac{1}{2}m(v_x^2+v_y^2+v_z^2)$,

wobei $m$ die Masse des Atoms ist, und $v_i$ die Komponenten des Geschwindigkeits-Vektors. Es kommen also drei Freiheitsgrade je Atom als Quadrat vor, daher ist die mittlere Energie je Atom

$⟨E⟩=\frac{3}{2}{k}_{B}T$

Daraus ergibt sich durch Differenzieren nach der Temperatur $T$ die Wärmekapazität $C_V$ von $(3/2)k_B$ je Atom, also

$C_V=\frac{3}{2}N{k}_B$

für ein einatomiges Gas mit N Atomen.

Bei zweiatomigen Gasmolekülen, wie sie beispielsweise Wasserstoff, Sauerstoff und Stickstoff bilden, sind zusätzlich zwei Rotationsfreiheitsgrade zu berücksichtigen (die Rotation um die Molekülachse, also um die dritte Raumrichtung, ist nicht relevant, und die Molekülschwingungen sind „eingefroren“). Daraus ergeben sich 5 Freiheitsgrade je Molekül, und somit

$C_V=\frac{5}{2}N{k}_B$

für ein Gas mit N Molekülen. Schwingungsfreiheitsgrade ergeben sich erst bei Temperaturen über 1000 K.

Wärmekapazität von Festkörpern

Bei Festkörpern kann die Schwingung der Atome um ihre Ruheposition durch das Potential eines harmonischen Oszillators angenähert werden. Je Raumrichtung $i$ ist die dazugehörige Energie durch

$E_i=\frac{1}{2}m{v}_i^2+\frac{1}{2}m{x}_i^2{\omega }_0^2$

gegeben, wobei ${\omega }_0$ die Kreisfrequenz des Oszillator ist, und $x_i$ die Auslenkung des Atoms aus seiner Ruhelage in Richtung $i$ bedeutet. Der erste Summand ist die kinetische Energie , der zweite die potentielle Energie. Es kommen also zwei Freiheitsgrade pro Atom und Raumdimension als Quadrat vor, in drei Dimensionen also sechs Freiheitsgrade je Atom. Daher ist die mittlere Energie je Atom

$⟨E⟩=\frac{6}{2}{k}_{B}T=3k_{B}T$

Bei $N$ Atomen sind also $6N$ Freiheitsgrade zu berücksichtigen. Daraus ergibt sich unmittelbar die Wärmekapazität von

$C=3N{k}_B$.

Diese Gleichung ist als Dulong-Petit-Gesetz bekannt. Auch hier gilt, dass die Freiheitsgrade nicht „eingefroren“ sein dürfen (die Temperatur muss deutlich über der Debye-Temperatur liegen); ansonsten kann die Wärmekapazität nur mit dem Debye-Modell berechnet werden.

Herleitung

Im Folgenden wird der Gleichverteilungssatz für klassische Systeme hergeleitet. Ausgangspunkt ist ein abgeschlossenes System, das energetisch an ein Wärmebad gekoppelt ist. Es bietet sich daher die Betrachtung des kanonischen Ensembles an.

Wir betrachten den Mittelwert der Größe $x_i\frac{\partial H}{\partial x_j}$, wobei $x_i$ für Orts- ($q_i$) oder Impulskoordinaten ($p_i$) stehen kann. $H$ bezeichnet die Hamiltonfunktion des Systems.

$⟨x_i\frac{\partial H}{\partial x_j}⟩=\frac{1}{Z}\int d\mathrm{\Gamma }{x}_i\frac{\partial H}{\partial x_j}{e}^{-\beta H}=-\frac{1}{\beta Z}\int d\mathrm{\Gamma }{x}_i\frac{\partial e^{-\beta H}}{\partial x_j}$

Die Integration erfolgt über den zugänglichen Phasenraum. $\beta =\frac{1}{kT}$ bezeichnet die inverse Temperatur, $Z=\int d\mathrm{\Gamma }{e}^{-\beta H}$ die kanonische Zustandssumme. Partielle Integration führt auf:

$⟨x_i\frac{\partial H}{\partial x_j}⟩=\frac{1}{\beta Z}\int d\mathrm{\Gamma }\frac{\partial x_i}{\partial x_j}{e}^{-\beta H}-\frac{1}{\beta Z}\int d{\mathrm{\Gamma }}_k{\left[x_i{e}^{-\beta H}\right]}_{x_j=a}^{x_j=b}$

wobei wir annehmen, dass $e^{-\beta H}$ hinreichend schnell für große $x_j$ abfällt, sodass die Randterme vernachlässigt werden können:

$⟨x_i\frac{\partial H}{\partial x_j}⟩=\frac{{\delta }_{ij}}{\beta Z}\int d\mathrm{\Gamma }{e}^{-\beta H}=\frac{{\delta }_{ij}}{\beta }$

Die allgemeinste Formulierung des Gleichverteilungssatzes für klassische Systeme im thermischen Gleichgewicht lautet

$⟨x_i\frac{\partial H}{\partial x_j}⟩=kT{\delta }_{ij}$

Die Herleitung wurde hier mit Hilfe des kanonischen Ensembles durchgeführt, sie ist auch mittels mikrokanonischen Ensembles möglich.

Anwendungen

Aus dem allgemeinen Gleichverteilungssatz lässt sich folgern: Jede Variable, die quadratisch in die Hamiltonfunktion eingeht, trägt mit $kT/2$ zur mittleren Energie bei:

$H=a{x}_i^2+f(x_j)\Rightarrow kT=⟨x_i\frac{\partial H}{\partial x_i}⟩=⟨2a{x}_i^2⟩\Rightarrow ⟨a{x}_i^2⟩=\frac{kT}{2}$

Dies gilt auch allgemeiner für die quadratische Form:

$H=\sum _{i,j}{x}_i{a}_{ij}{x}_j\Rightarrow \sum _{n}kT=⟨\sum _n{x}_n\frac{\partial H}{\partial x_n}⟩=⟨\sum _n{x}_n(\sum _j{a}_{nj}{x}_j+\sum _i{x}_i{a}_{in})⟩=2⟨\sum _{i,j}{x}_i{a}_{ij}{x}_j⟩\Rightarrow ⟨H⟩=N\frac{kT}{2}$

Einatomiges ideales Gas

Für $N$ nicht wechselwirkende Teilchen (einatomiges ideales Gas) in drei Raumdimensionen besteht die Hamiltonfunktion nur aus dem kinetischen Anteil:

$H=E_{\mathrm{kin}}=\sum _{i=1}^N\frac{{\mathbf{p}}_i^2}{2m}=\sum _{j=1}^{3N}\frac{p_j^2}{2m}$

Die Anwendung des obigen Ergebnisses

$kT=⟨p_j\frac{\partial H}{\partial p_j}⟩=2⟨\frac{p_j^2}{2m}⟩$

liefert

$⟨E_{\mathrm{kin}}⟩=\sum _{j=1}^{3N}⟨\frac{p_j^2}{2m}⟩=3N\frac{kT}{2}$

Das heißt, pro Translationsfreiheitsgrad (hier $3N$) ist die mittlere kinetische Energie $kT/2$.

Zweiatomiges ideales Gas

Für ein zweiatomiges ideales Gas, das heißt, die einzelnen Moleküle wechselwirken nicht miteinander, lautet die Hamiltonfunktion unter Vernachlässigung der Rotations-Vibrationskopplung (konstantes Trägheitsmoment)

$H=\sum _{i=1}^{N}h\mathrm{mit}h=\underset{\mathrm{Translation}}{\underbrace{\frac{1}{2m_{\mathrm{ges}}}(p_x^2+p_y^2+p_z^2)}}+\underset{\mathrm{Rotation}}{\underbrace{\frac{1}{2\mathrm{\Theta }}(\frac{p_{\phi }^2}{{\sin }^2\vartheta }+p_{\vartheta }^2)}}+\underset{\mathrm{Vibration}}{\underbrace{(\frac{p_{\xi }^2}{2m_{\mathrm{r}}}+\frac{m_{\mathrm{r}}{\omega }^2{\xi }^2}{2})}}$

wobei $m_{\mathrm{ges}}$ die Gesamtmasse, $m_{\mathrm{r}}$ die reduzierte Masse und $\mathrm{\Theta }$ das Trägheitsmoment eines Moleküls ist. $\xi$ beschreibt die Auslenkung aus dem Gleichgewichtsabstand. Insgesamt gehen also sieben Größen quadratisch in die Hamiltonfunktion ein: $p_x,p_y,p_z,p_{\phi },p_{\vartheta },p_{\xi },\xi$. Daraus folgt:

$⟨h⟩=\frac{7}{2}kT\Rightarrow ⟨H⟩=\frac{7}{2}NkT$

Diese mittlere Energie ist so nur bei hohen Temperaturen gültig, wenn auch Rotationen und Vibrationen thermisch angeregt werden.

Thermische Zustandsgleichung

Betrachten wir ein reales Gas in einem Behälter, so lautet der Hamiltonian

$H=\sum _{i=1}^{3N}\frac{p_i^2}{2m}+V_{\mathrm{Wand}}(q_1,\dots ,q_{3N})+\mathrm{\Phi }(q_1,\dots ,q_{3N})$

wobei $V_{\mathrm{Wand}}$ das Potential zwischen Wand und Teilchen ist und $\mathrm{\Phi }$ das Potential zwischen den Teilchen. Für einen würfelförmigen Behälter mit Seitenlänge $L$ schreibt sich das Wandpotential z. B. wie folgt

$V_{\mathrm{Wand}}=A\sum _{i=1}^{3N}[\mathrm{\Theta }(q_i-L)+\mathrm{\Theta }(-q_i)]$

Dabei wurde die Heaviside-Funktion $\mathrm{\Theta }$ verwendet. Die Anwendung des Gleichverteilungssatzes liefert:

$3NkT=\sum _{i=1}^{3N}⟨q_i\frac{\partial H}{\partial q_i}⟩=\sum _{i=1}^{3N}⟨q_i\frac{\partial V_{\mathrm{Wand}}}{\partial q_i}⟩+\sum _{i=1}^{3N}⟨q_i\frac{\partial \mathrm{\Phi }}{\partial q_i}⟩$

Betrachte nun den ersten Term auf der rechten Seite

$\begin{array}{rl}\sum _{i=1}^{3N}{q}_i\frac{\partial V_{\mathrm{Wand}}}{\partial q_i}& =\sum _{i=1}^{3N}{q}_{i}A[\delta (q_i-L)-\delta (q_i)]=\sum _{i=1}^{3N}LA\delta (q_i-L)\\ & =-L\frac{\partial V_{\mathrm{Wand}}}{\partial L}=-L\frac{\partial H}{\partial L}=-L\frac{\partial V}{\partial L}\frac{\partial H}{\partial V}=-3V\frac{\partial H}{\partial V}\end{array}$

Es wurde ausgenutzt, dass die distributive Ableitung der Heaviside-Funktion die Delta-Distribution ist. Im vorletzten Schritt konnte das Volumen $V=L^3$ eingeführt werden. Nun führt man noch eine Ensemblemittelung durch und verwendet, dass der Druck definiert ist durch $p=-⟨\partial H/\partial V⟩$ (siehe z.B. hier: Kanonisches Ensemble).

$\sum _{i=1}^{3N}⟨q_i\frac{\partial V_{\mathrm{Wand}}}{\partial q_i}⟩=-3⟨V\frac{\partial H}{\partial V}⟩=-3V\underset{-p}{\underbrace{⟨\frac{\partial H}{\partial V}⟩}}=3pV$

Somit erhält man die thermische Zustandsgleichung:

$NkT=pV+\frac{1}{3}\underset{\mathrm{Virial}}{\underbrace{\sum _{i=1}^{3N}⟨q_i\frac{\partial \mathrm{\Phi }}{\partial q_i}⟩}}$

Diese entspricht der idealen Gasgleichung, die um einen Zusatzterm – dem Virial – erweitert ist. Das Virial kann in Potenzen der Teilchendichte $N/V$ entwickelt werden (siehe: Virialentwicklung).

Literatur

• Schwabl: Statistische Mechanik. Springer-Verlag, Berlin, 3. Auflage 2006, ISBN 978-3-540-31095-2