Zum Inhalt springen

Statistische Mechanik/ Fokker-Planck-Gleichung

Aus Wikibooks

Die Fokker-Planck-Gleichung ist eine stochastische Differentialgleichung, die die zeitliche Evolution eines Systems von Teilchens mit umgebenden Wärmebad (z.B. ein flüssiges Lösungsmittel) beschreibt. Die Größe, die von Interesse ist, ist die Wahrscheinlichkeit, ein Teilchen an einer Position zu finden. Die Wahrscheinlichkeitsdichte abhängig von der Zeit ist also gesucht. Für diese lässt sich eine Differentialgleichung herleiten. Wir gehen zunächst von der Langevin-Gleichung für ein kugelförmiges Teilchen mit Radius in einem Wärmebad mit Viskosität (diese ist als wesentlich größer als die Masse des Teilchens angenommen, sodass die Massenträgheit des Teilchens vernachlässigt werden kann) aus; diese lautet:

.

Hierbei ist die Summe aller auf das Teilchen wirkende Kräfte, die nicht vom Wärmebad ausgehen, eine Konstante und das gaußverteilte weiße Rauschen. Das weiße Rauschen macht obige Bewegungsgleichung vom Zufall abhängig, sodass man auf eine Teilchen-Aufenthaltswahrscheinlichkeit angewiesen ist. Wir betrachten nun zu zwei verschiedenen Zeitpunkten, nämlich und mit einer möglichst klein gewählten Zeitspanne . Somit können wir den Quotienten

betrachten. Wir analysieren nun, wie die Wahrscheinlichkeitsverteilung zu einem späteren Zeitpunkt aussehen müsste, falls wir bereits kennen. Aus der Wahrscheinlichkeitstheorie wissen wir, wie man eine Wahrscheinlichkeit berechnet werden kann, falls die Wahrscheinlichkeitsverteilung für ein bestimmtes Ereignis bekannt ist; siehe auch: Bayes'scher Satz. Es wird die Wahrscheinlichkeit, von Position zum Zeitpunkt zur Position zum Zeitpunkt zu gelangen, bezogen auf die Zeiteinheit , benötigt; diese bezeichnen wir als . Dann gilt folgende Beziehung (Chapman-Kolmogorov-Gleichung):

.

Es wird über die Menge aller möglichen Positionen integriert. Wir erhalten nun:

.

Innerhalb der kleinen Zeitspanne wird sich das Teilchen in der Nähe von der Position befinden, also es kann angenommen, dass die vorherige Position , hier auch Integrationsvariable nur geringfügig von abweicht. Wir machen also eine Taylor-Entwicklung um :

.

Nun ist und dessen partielle x-Ableitungen nicht mehr von abhängig und können dies vor das Integral ziehen. Es gilt noch , denn es wird auf alle Fälle einer der Zustände im Ereignisraum erreicht (sicheres Ereignis) und wir erhalten die Fokker-Planck-Gleichung:

.

Hierbei ist

(Driftkoeffizient),

(Diffusionskoeffizient).

Aus der obigen Langevin-Gleichung kann der Drift- und Diffusionskoeffizient in Abhängigkeit von der Teilchenposition berechnet werden. Dazu schaut man sich an, wie sich die Teilchenposition innerhalb der sehr kurzen Zeitspanne ändert. Sei die Anfangsposition, so folgt durch Integration der Langevin-Gleichung:

.

Wir können die gewichtete Mittelung über die Operation durch eine statistische Mittelung über das statistische Rauschen dividiert durch den Zeitschritt ersetzen; ein deterministischer Term bleibt nach der Mittelung unverändert. Es gilt (Rauschen ist im Mittel Null) und folglich für den Driftkoeffizient:

Hierbei haben wir im letzten Schritt angenommen, dass sich der deterministische Anteil über die Zeitspanne nicht verändert. Der Diffusionskoeffizient kann analog berechnet werden, hierbei verwendet man die Identität und führt den Grenzübergang durch. Dies resultiert in

.

Analog kann die Fokker-Planck-Gleichung für einen Phasenraum mit n verschiedenen Variablen für die entsprechende zeitabhängige Wahrscheinlichkeitsdichte verallgemeinert werden. In diesem Fall lautet sie

mit Driftkoeffizienten und Korrelationskoeffizienten .

Fluktuation-Dissipation-Theorem

[Bearbeiten]

Es soll noch ein Zusammenhang zwischen der Konstante und thermodynamischen Größen wie der Temperatur des Wärmebades hergestellt werden. Dieser Zusammenhang wird auch Fluktuation-Dissipation-Theorem genannt, da dieses Die Reibungskraft und die stochastische Kraft miteinander verknüpft. Nehmen wir an, dass keine äußeren Kräfte außer Reibungs- und Zufallskraft auf das Teilchen wirken (auch wenn andere Kräfte angreifen, ändern sie nicht die Varianz der Geschwindigkeit, sondern nur den Erwartungswert; daher brauchen wir nur dieses Szenario zum Beweis des Fluktuation-Dissipations-Theorems zu betrachten). Dann folgt (wenn die Teilchenmasse und die Geschwindigkeit ist):

.

Die Lösung dieser Differentialgleichung ist

und nach Quadrieren ergibt sich:

Wenn wir nun den Erwartungswert bilden, die Zeitintegrale auswerten und den Limes , erhalten wir:

.

Aus der klassischen Gleichgewichtsstatistik wissen wir, dass gilt (Äquipartitionstheorem). Hierbei ist die Boltzmann-Konstante und die Standardabweichung in der Zufallsgröße "Teilchen-Geschwindigkeit". Diese Beziehung muss im Gleichgewicht, d.h. aufgrund des Strebens ins Gleichgewicht eines Nichtgleichgewichtssystems für unendlich große Beobachtungszeiten, erfüllt sein. Schließlich ergibt sich

(*)

und nach Auflösen nach :

.

Die Intensität der Zufallskraft ist demnach der Reibungskonstante proportional; also sind Reibungs-und Zufallskraft miteinander verknüpft. Die Gleichung (*) wird auch Einstein-Beziehung zwischen Diffusionskoeffizienten, Temperatur und Reibungskoeffizienten genannt.