Mikrokanonischer Zustand

Der mikrokanonische Zustand ist ein Gleichgewichtszustand der statistischen Physik mit der exakten Nebenbedingung einer konstanten Gesamtenergie. Die einzige Information über ein Quantensystem sei, dass die Gesamtenergie $E=\langle \hat H\rangle$ gleich E₀ bzw. aus dem Intervall [E₀,E₀ + ΔE] mit $\Delta E\ll E$ ist, wobei die Zustände mit von außen vorgegebenen Parametern, wie Volumen und Teilchenzahl, verträglich sein müssen. Dies entspricht einem System innerhalb eines total abgeschlossenen Kastens (kein Energie- oder Teilchenaustausch mit der Umgebung, keine äußeren Felder). Das Potential außerhalb des Kastens ist unendlich, wodurch der Hamilton-Operator nur diskrete Energie-Eigenwerte besitzt und die Zustände abzählbar sind.

Quantenmechanisch

Im Gleichgewicht ändert sich der Dichteoperator des Systems nicht $\tfrac{\partial\hat{\rho}}{\partial t}=0$. Nach der von-Neumann-Gleichung $\tfrac{\partial\hat{\rho}}{\partial t}=-\tfrac{i}{\hbar}[\hat{H},\hat{\rho}]$ vertauscht im Gleichgewicht der Dichteoperator mit dem Hamiltonoperator (Kommutator gleich Null). Daher lassen sich gemeinsame Eigenzustände wählen, d.h. die Energieeigenzustände von $\hat{H}$ sind auch Eigenzustände von $\hat{\rho}$.

Man schränkt den Hilbertraum $\mathcal{H}$ auf einen Teilraum ein, der von den Zustandsvektoren mit Eigenwert E₀ aufgespannt wird (Eigenraum). Sei $\left\{ |E,n\rangle\right\}$ ein vollständiges Orthonormalsystem (VONS), nämlich die Eigenzustände des Hamiltonoperators $\hat{H}|E,n\rangle=E|E,n\rangle$, so wird der Unterraum $\mathcal{H}_{E_{0}}$ von den Basisvektoren $|E,n\rangle$ aufgespannt, für die E = E₀. Die Energie E₀ ist im Allgemeinen entartet (deswegen ist ein Zustand nicht eindeutig durch Angabe von E bestimmt, sondern eine weitere Quantenzahl n ist nötig); der Entartungsgrad g entspricht der Dimension des Unterraums $\textrm{dim}\mathcal{H}_{E_{0}}$.

Der Hamiltonoperator $\hat{H}$ kann nicht zwischen den g Basiszuständen zu E₀ unterscheiden (Entartung); da $\hat{H}$ keinen Basiszustand bevorzugt, wird somit jedem Basiszustand a priori dieselbe Wahrscheinlichkeit 1 / g zugeordnet: Nach der Maximum-Entropie-Methode ist das System durch den Zustand zu beschreiben, welcher die Entropie maximiert. Die Entropie wird genau dann maximal, wenn jeder Basisvektor die gleiche Wahrscheinlichkeit $1/\textrm{dim}\mathcal{H}_{E_{0}}$ hat.

Daher ergibt sich der Dichteoperator des mikrokanonischen Zustands zu

$\hat{\rho}_{E_{0}}=\frac{1}{Z_m(E_{0})}\sum_{n=1}^{\textrm{dim}\mathcal{H}_{E_{0}}}|E_{0},n\rangle\langle E_{0},n|$

mit der mikrokanonischen Zustandssumme (oft auch mit Ω bezeichnet)

$Z_m(E_{0})=\textrm{dim}\mathcal{H}_{E_{0}}=\textrm{Tr}\left(\sum_{n=1}^{\textrm{dim}\mathcal{H}_{E_{0}}}|E_{0},n\rangle\langle E_{0},n|\right),$

wobei die Spur eines Operators folgendermaßen definiert ist: $\textrm{Tr}(O)=\sum_{m}\langle m|O|m\rangle$ für beliebiges VONS $\left\{ |m\rangle\right\}$ von $\mathcal{H}$

Dass jeder Energieeigenzustand mit Energie E₀ dieselbe Wahrscheinlichkeit besitzt, ist die grundlegende Annahme der Gleichgewichts-Statistik. Aus ihr können alle Gleichgewichtseigenschaften für abgeschlossene oder offene Systeme hergeleitet werden (z. B. kanonischer oder großkanonischer Zustand).

Klassisch

Analog ergibt sich der klassische mikrokanonische Zustand für N Teilchen (Phasenraumdichte)

$\rho_{E_{0}}(\vec{x}_1,\vec{p}_1,\ldots,\vec{x}_N,\vec{p}_N)=\frac{1}{Z_m(E_{0})} \delta (E_0 - H(\vec{x}_1,\vec{p}_1,\ldots,\vec{x}_N,\vec{p}_N))$

mit der klassischen mikrokanonischen Zustandssumme (Gesamtzahl der zugänglichen Mikrozustände, die dieselbe Gesamtenergie E₀ besitzen)

$Z_m(E_{0})=\int_{\mathbb R^{6N}}d\tau \;\delta (E_0 - H(\vec{x}_1,\vec{p}_1,\ldots,\vec{x}_N,\vec{p}_N))$

mit $d\tau = \frac{1}{\xi h^{3N}}\;d^3x_1d^3p_1\ldots d^3x_Nd^3p_N$

wobei für N identische Teilchen der Faktor ξ = N! die Mehrfachzählung ununterscheidbarer Teilchen verhindert,

und für n verschiedene Teilchensorten mit $N_1,\ldots,N_n$ Teilchenzahlen und $\sum^{n}_{i=1}N_i=N$ der Faktor $\xi=N_1!\ldots N_n!$.

Entropie

Die Entropie kann man durch die Zustandssumme ausdrücken:

$S(\hat{\rho}_{E_{0}})=k_{{\rm B}}\ln\left[Z_{m}(E_{0})\right]$

Dies lässt sich aus der Definition der Entropie herleiten, wobei die Zustandssumme gleich $Z_m(E_{0})=\textrm{dim}\mathcal{H}_{E_{0}}=g$ ist.

$\begin{array}{rcl} S(\hat{\rho}_{E_{0}}) & = & -k_{{\rm B}}\left\langle \ln\hat{\rho}_{E_{0}}\right\rangle =-k_{{\rm B}}\textrm{Tr}\left[\hat{\rho}_{E_{0}}\ln\hat{\rho}_{E_{0}}\right]=-k_{{\rm B}}\sum_{E,m}\langle E,m|\hat{\rho}_{E_{0}}\ln\hat{\rho}_{E_{0}}|E,m\rangle\\ & = & -k_{{\rm B}}\sum_{E,m}\langle E,m|\sum_{n=1}^{g}|E_{0},n\rangle\frac{1}{g}\langle E_{0},n|\ln\left(\sum_{n'=1}^{g}|E_{0},n'\rangle\frac{1}{g}\langle E_{0},n'|\right)|E,m\rangle\\ & = & -k_{{\rm B}}\sum_{n=1}^{g}\frac{1}{g}\ln(\frac{1}{g})=-k_{{\rm B}}\ln(\frac{1}{g})=k_{{\rm B}}\ln\left(g\right)\end{array}$

Thermodynamik

Die Entropie ist von der Gesamtenergie E₀ (im folgenden E genannt), vom Volumen V und der Teilchenzahl N abhängig (weil der Hamiltonoperator von V und N abhängig ist): S(E,V,N). Die Gesamtenergie wird in der Thermodynamik meist innere Energie genannt. Nun werden die Ableitungen der Entropie untersucht.

Temperatur

Bringt man zwei abgeschlossene Systeme in schwachen thermischen Kontakt, dann kann pro Zeiteinheit nur wenig Energie ausgetauscht werden, so dass die Teilsysteme etwa im Gleichgewicht bleiben und die Entropie des Gesamtsystems sich additiv schreiben lässt: S_ges = S₁ + S₂. Die Gesamtenergie bleibt stets konstant E_ges = E₁ + E₂. Nimmt die Energie des einen Systems zu, so muss im gleichen Maß die des anderen abnehmen: dE₁ = − dE₂. Durch Energieaustausch wird das Gleichgewicht des Gesamtsystems erreicht, dort nimmt die Gesamtentropie ein Maximum ein:

$0=dS_{\text{ges}}=\frac{\partial S_{1}}{\partial E_{1}}dE_{1}+\frac{\partial S_{2}}{\partial E_{2}}dE_{2}=\left(\frac{\partial S_{1}}{\partial E_{1}}-\frac{\partial S_{2}}{\partial E_{2}}\right)dE_{1}\quad\Rightarrow\quad\frac{\partial S_{1}}{\partial E_{1}}=\frac{\partial S_{2}}{\partial E_{2}}$

Im Gleichgewicht werden die Ableitungen der Entropie nach der Energie für die Teilsysteme gleich groß. Man definiert hierüber die Temperatur:

$\frac{1}{T}:=\left(\frac{\partial S}{\partial E}\right)_{N,V}\quad\Rightarrow\quad T:=\left(\frac{\partial E}{\partial S}\right)_{N,V}$

Somit werden im Gleichgewicht die Temperaturen beider Teilsysteme gleich groß.

Druck

Der Druck wird definiert über

$p:=-\left\langle\frac{\partial\hat{H}}{\partial V}\right\rangle=-\left(\frac{\partial E}{\partial V}\right)_{S,N}$

also die isentrope (S=const) Energieänderung pro Volumen. Die Entropie S(E(V,N),V,N) wird nach dem Volumen differenziert:

$0=\frac{dS}{dV}=\left(\frac{\partial E}{\partial V}\right)_{S,N}\left(\frac{\partial S}{\partial E}\right)_{N,V}+\left(\frac{\partial S}{\partial V}\right)_{E,N}=-p\frac{1}{T}+\left(\frac{\partial S}{\partial V}\right)_{E,N}$

Somit erhält man

$\left(\frac{\partial S}{\partial V}\right)_{E,N}=\frac{p}{T}$

Chemisches Potential

Das chemische Potential wird definiert über

$\mu:=\left\langle\frac{\partial\hat{H}}{\partial N}\right\rangle=\left(\frac{\partial E}{\partial N}\right)_{S,V}$

also die isentrope Energieänderung pro Teilchen. Die Entropie S(E(V,N),V,N) wird nach der Teilchenzahl differenziert:

$0=\frac{dS}{dN}=\left(\frac{\partial E}{\partial N}\right)_{S,V}\left(\frac{\partial S}{\partial E}\right)_{N,V}+\left(\frac{\partial S}{\partial N}\right)_{E,V}=\mu\frac{1}{T}+\left(\frac{\partial S}{\partial N}\right)_{E,V}$

Daraus folgt

$\left(\frac{\partial S}{\partial N}\right)_{E,V}=-\frac{\mu}{T}$

Allgemein lässt sich festhalten: Ist der Hamiltonoperator von einem äußeren Parameter a abhängig (z. B. Volumen oder Teilchenzahl), so ist die Ableitung der Entropie nach a bei konstanter Energie gleich:

$\left(\frac{\partial S}{\partial a}\right)_{E}=-\frac{1}{T}\left\langle\frac{\partial\hat{H}}{\partial a}\right\rangle$

Thermodynamisches Potential

Zusammenfassend lassen sich die Ableitungen nach Energie, Volumen und Teilchenzahl darstellen:

$\frac{1}{T}\begin{pmatrix}1\\ p\\ -\mu\end{pmatrix}=\begin{pmatrix}\partial_{E}\\ \partial_{V}\\ \partial_{N}\end{pmatrix}S(E,V,N)$

Das totale Differential der Entropie lautet:

$\mathrm{d}S=\frac{1}{T}\,\mathrm{d}E+\frac{p}{T}\,\mathrm{d}V-\frac{\mu}{T}\,\mathrm{d}N$

Die Entropie S(E,V,N) lässt sich nach der Energie auflösen E(S,V,N). Die Energie ist das thermodynamische Potential der Mikrokanonik. Mit ihm lassen sich die obigen Ableitungen kompakt als Gradient des Potentials schreiben:

$\begin{pmatrix}T\\ -p\\ \mu\end{pmatrix}=\begin{pmatrix}\partial_{S}\\ \partial_{V}\\ \partial_{N}\end{pmatrix}E(S,V,N)$

Das totale Differential der Energie lautet somit:

$\mathrm{d}E=T\,\mathrm{d}S-p\,\mathrm{d}V+\mu\,\mathrm{d}N$

Dies ist die Fundamentalgleichung der Thermodynamik.

Beispiele

Ideales Gas

Ein Beispiel eines mikrokanonisch präparierten Systems, das mit den klassischen Methoden berechnet werden kann, ist das ideale Gas; siehe Herleitung unter Sackur-Tetrode-Gleichung.

Ungekoppelte Oszillatoren

Ein weiteres Beispiel ist ein System aus N nicht-wechselwirkenden gleichartigen harmonischen Oszillatoren. Deren Gesamt-Hamiltonoperator ist

$\hat{H}=\sum_{i=1}^{N}\hbar\omega\left(\hat{n}_{i}+\frac{1}{2}\right)$

wobei $\hat{n}_i$ der Besetzungszahloperator des i-ten Oszillators ist. Zu gegebener Gesamtenergie

$E=\sum_{i=1}^{N}\hbar\omega\left(n_{i}+\frac{1}{2}\right)=\hbar\omega\left(\sum_{i=1}^{N}n_{i}+\frac{N}{2}\right)\quad\Rightarrow\quad n:=\sum_{i=1}^{N}n_{i}=\frac{E}{\hbar\omega}-\frac{N}{2}$

soll nun die Zustandssumme berechnet werden. Diese ist gleich dem Entartungsgrad zur Energie E bzw. gleich der Anzahl der Möglichkeiten n ununterscheidbare Energiequanten auf N mehrfach besetzbare Oszillatoren zu verteilen bzw. der Anzahl der Möglichkeiten n ununterscheidbare Energiequanten auf N + n − 1 einfach besetzbare Oszillatoren zu verteilen:

$Z_{m}=\binom{N+n-1}{n}=\frac{(N+n-1)!}{n!\,(N-1)!}$

Hieraus lässt sich die Entropie berechnen:

$S=k_{{\rm B}}\ln(Z_{m})=k_{{\rm B}}\left\{ \ln\left((N+n-1)!\right)-\ln(n!)-\ln\left((N-1)!\right)\right\}$

Für große N und n kann man den Logarithmus der Fakultät mit der Stirling-Formel bis zur ersten Ordnung $\ln(N!)\approx N\ln(N)-N$ entwickeln und außerdem die 1 gegenüber der riesigen Zahl N vernachlässigen:

$\begin{align} S & =k_{{\rm B}}\left\{ (N+n-1)\,\ln(N+n-1)-n\,\ln(n)-(N-1)\,\ln(N-1)\right\} \\ & =k_{{\rm B}}\left\{ (N+n)\,\ln(N+n)-n\,\ln(n)-N\,\ln(N)\right\} \end{align}$

Umsortieren und Verwenden von $\epsilon:=\frac{E}{N\hbar\omega}=\frac{n}{N}+\frac{1}{2}$ liefert die Entropie, die wegen des Vorfaktors N extensiv ist (n / N und $\epsilon$ sind nämlich intensiv):

$\begin{align} S & =Nk_{{\rm B}}\left\{ \left(\frac{n}{N}+1\right)\,\ln\left(\frac{n}{N}+1\right)-\left(\frac{n}{N}\right)\,\ln\left(\frac{n}{N}\right)\right\} \\ & =Nk_{{\rm B}}\left\{ \left(\epsilon+\frac{1}{2}\right)\,\ln\left(\epsilon+\frac{1}{2}\right)-\left(\epsilon-\frac{1}{2}\right)\,\ln\left(\epsilon-\frac{1}{2}\right)\right\} \end{align}$

Nun kann die Temperatur berechnet werden:

$\frac{1}{T}=\frac{\partial S}{\partial E}=\frac{\partial\epsilon}{\partial E}\frac{\partial S}{\partial\epsilon}=\frac{k_{B}}{\hbar\omega}\ln\frac{\epsilon+\frac{1}{2}}{\epsilon-\frac{1}{2}}=\frac{k_{B}}{\hbar\omega}\ln\frac{E+\frac{N\hbar\omega}{2}}{E-\frac{N\hbar\omega}{2}}$

Schließlich kann man noch nach der Energie auflösen. Die Energie steigt monoton mit der Temperatur an:

$E=N\hbar\omega\left\{ \frac{1}{2}+\frac{1}{\exp(\frac{\hbar\omega}{k_{B}T})-1}\right\}$

Siehe auch

• Ensemble (Physik)
• Kanonischer Zustand
• Zustandssumme

Literatur

• Balian: From Microphysics to Macrophysics 1. 2. Auflage. Springer-Verlag, Berlin 2006, ISBN 3-540-45469-1