# Ehrenfest-Modell

Das Ehrenfest-Modell (auch als Ehrenfest-Kette bekannt) ist ein stochastisches Modell, das den Stoffaustausch zwischen zwei durch eine Membran getrennte Behältnisse beschreibt. Das Modell wurde zuerst durch den österreichischen Physiker Paul Ehrenfest (1880–1933) vorgeschlagen und ist einer von vielen Beiträgen der Physik zur Entwicklung der mathematischen Theorie der stochastischen Prozesse.

## Das Modell

Bei verschiedenen Substanzen wurde beobachtet, dass die Verteilung der Substanz in einem solchen Experiment im Laufe der Zeit zwar einem Gleichgewichtszustand entgegenstrebt, aber dennoch auch nach Erreichen desselben stets unkontrollierbaren, scheinbar zufälligen Schwankungen ausgesetzt bleibt.
Diesen Umstand versuchte das folgende Modell zu erklären:
Zu Beginn befinden sich in beiden Behältern zusammen eine endliche Anzahl von $N$ Partikeln; etwa die einzelnen Moleküle des Stoffes, wovon sich anfangs $l_{0}\leq N$ im linken und analog $r_{0}=N-l_{0}\leq N$ im rechten Behälter aufhalten. In jedem Zeitschritt wird nun genau eines dieser $N$ Teilchen gleichverteilt ausgewählt, das den Behälter wechselt, sodass $l$ und $r$ in jedem Schritt genau um eins ansteigen oder fallen.
Mathematisch gesehen handelt es sich bei diesem zufälligen Vorgang um eine Markow-Kette $(l_{n})_{n\in \mathbb {N} _{0}}$ mit Zustandsraum $\{0,1,2,\ldots N\}$ und einer Übergangsmatrix $\Pi =\left(\pi _{a,b}\right)$, gegeben durch
    $\pi _{a,b}={\begin{cases}{\frac {a}{N}},&{\text{falls }}b=a-1\,,\\{\frac {N-a}{N}},&{\text{falls }}b=a+1\,,\\0&{\text{sonst}}.\end{cases}}$

## Mathematische Eigenschaften

* Die oben definierte Ehrenfest-Kette besitzt eine eindeutig bestimmte stationäre Verteilung: Ist die Anzahl $l_{n}$ der Teilchen im linken (oder rechten) Behälter binomialverteilt mit Parameter ${\frac {1}{2}}$, ist also $P(l_{n}=k)={\binom {N}{k}}{\frac {1}{2^{N}}}$ für $k=0,1,\ldots N$, so hat $l_{n+1}$ dieselbe Verteilung.
* Die Konvergenz der Kette gegen diese Verteilung ist allerdings nicht gegeben, da die Kette periodisch ist (das erkennt man daran, dass $l_{n}$ stets zwischen geraden und ungeraden Zahlen wechselt und somit $P(l_{n}=k)$ jedes zweite Mal gleich null ist). Dies kann man umgehen, indem man zur aperiodischen Version der Kette übergeht und die Übergangsmatrix $\Pi$ für einen festen Parameter $p\in0,1$ durch die Matrix ${\hat {\Pi }}:=pI_{N+1}+(1-p)\Pi$ ersetzt (dabei ist $I_{N+1}$ die Einheitsmatrix). Interpretation: mit Wahrscheinlichkeit $p$ bleibt die Anzahl der Teilchen in den Behältern unverändert, mit Wahrscheinlichkeit $1-p$ ändert sie sich nach dem oben beschriebenen Verfahren. Dadurch wird die Kette aperiodisch und konvergiert für $n\to \infty$ gegen die stationäre Verteilung, die sich durch diese Modifikation nicht ändert.

## Beispiel

Gegeben seien zwei Behälter, die durch eine Membran voneinander getrennt sind. In dem linken Behälter befinden sich zu Beginn des Experiments $10$ Moleküle und der rechte Behälter ist noch leer. Durch die Membran kann genau ein Molekül pro Zeiteinheit den Behälter wechseln.
Da der rechte Behälter zu Beginn noch leer ist, wird in der ersten Sekunde ein Molekül aus dem linken in den rechten Behälter fliegen. Anschließend befinden sich nur noch $9$ Moleküle in dem linken Behälter. Nun gibt es zwei Möglichkeiten: Entweder eines der verbleibenden $9$ Moleküle des linken Behälters fliegt in den rechten Bereich, oder das Molekül rechts fliegt wieder zurück in den linken Bereich. Jedes Molekül soll hierbei die gleiche Chance haben, den Behälter zu wechseln. Demnach beträgt die Wahrscheinlichkeit $90$ %, dass ein weiteres Molekül von links nach rechts fliegt. Bei $8$ Molekülen links beträgt diese Wahrscheinlichkeit nur noch $80$ % und so weiter.
Der Übergangsgraph enthält die Zustände $0$ bis $10$, welche die Anzahl der Moleküle im linken Behälter repräsentieren. Die Markow-Kette startet im Zustand $10$. Vervollständigt man den Übergangsgraphen und erstellt eine dazu passende Übergangsmatrix, kann man die Wahrscheinlichkeitsverteilungen für die Anzahl Moleküle im linken Behälter für jeden Zeitpunkt bestimmen. Nach $5$ Zeiteinheiten besteht mit einer Wahrscheinlichkeit von $30{,}2$ % zum ersten Mal die Möglichkeit zum physikalischen Gleichgewicht.
Die stationäre Verteilung lässt sich mit Hilfe der oben formulierten Formel
    $P(l_{n}=k)=P(l_{n+1}=k)={\binom {10}{k}}{\frac {1}{2^{10}}}$ für $k=0,1,\ldots ,10$
ermitteln. Dadurch ergibt sich die Wahrscheinlichkeitsverteilung
    $(0.10,0.98,4.39,11.72,20.51,24.61,20.51,11.72,4.39,0.98,0.10)$ %.