Langevin-Gleichung

Stochastische Differentialgleichung
(Weitergeleitet von Langevingleichung)

Eine Langevin-Gleichung (nach Paul Langevin) ist eine stochastische Differentialgleichung, welche die Dynamik einer Teilmenge der Freiheitsgrade eines physikalischen Systems beschreibt. Dabei handelt es sich typischerweise um „langsame“ (makroskopische) Freiheitsgrade, die „schnellen“ (mikroskopischen) Freiheitsgrade sind verantwortlich für die stochastische Natur der Differentialgleichung.

Brownsche Bewegung als Prototyp

Bearbeiten

Die ursprüngliche Langevingleichung beschreibt die Bewegung eines in einer Flüssigkeit schwebenden Teilchens infolge von Zusammenstößen mit Flüssigkeitsmolekülen,

 

Die Gleichung wurde aufgrund heuristischer Überlegungen von Paul Langevin aufgestellt. Die langsame Variable ist hier die Geschwindigkeit   des Teilchens,   ist die Teilchenmasse, die Konstante   ist der Reibungskoeffizient. Die Größe   ist die sogenannte fluktuierende Kraft, in diesem Fall ein Gauß'sches weißes Rauschen mit Korrelationsfunktion

 

Hierbei ist   die Boltzmannkonstante,   die Temperatur und   die i-te Komponente des Vektors  . Die  -Distribution für die Korrelation in der Zeit bedeutet, dass die Kraft zur Zeit   völlig unkorreliert ist mit der Kraft zu einer anderen Zeit. Das ist natürlich eine Näherung. Die tatsächliche fluktuierende Kraft ist zumindest über ein Zeitintervall das der Stoßdauer entspricht korreliert. Jedoch wird die Langevingleichung zur Beschreibung der Bewegung eines „makroskopischen“ Teilchens über viel größere Zeitskalen verwendet, und in diesem Grenzfall ergeben die  -Korrelation und die Langevingleichung den korrekten stochastischen Prozess.

Eine weitere prototypische Eigenschaft der Langevingleichung ist das Auftreten des Dämpfungskoeffizienten   in der Korrelationsfunktion der fluktuierenden Kraft. Der technische Terminus dafür ist Einstein-Smoluchowski-Beziehung.

Mathematische Aspekte

Bearbeiten

Eine in der Zeit exakt  -korrelierte fluktuierende Kraft   ist keine mathematische Funktion im üblichen Sinn und auch die Ableitung   ist in diesem Fall nicht definiert. Dieses Problem verschwindet wenn die Langevingleichung in integraler Form

 

geschrieben wird, und eine Langevingleichung sollte im einfachsten Fall immer als Abkürzung für ihr Zeitintegral interpretiert werden. Die mathematische Bezeichnung für Gleichungen dieser Art ist stochastische Differentialgleichung.

Ein anderes mathematisches Problem tritt auf für (ziemlich spezielle[1]) Langevingleichungen mit einer multiplikativen fluktuierenden Kraft, d. h. Termen wie   auf der rechten Seite. Solche Gleichungen können nach Ito- oder Stratonovich-Schema interpretiert werden (nach Itō Kiyoshi, Ruslan Stratonovich), und wenn die Herleitung der Langevingleichung hierzu keine Information liefert, ist die Herleitung sowieso fragwürdig.

Generische Langevingleichung

Bearbeiten

Eine generische Langevingleichung lässt sich mit Projektionsoperator-Methoden aus der klassischen Mechanik herleiten.[2][3] Diese generische Langevingleichung spielt eine zentrale Rolle in der Theorie der kritischen Dynamik[4] und in anderen Bereichen der Nicht-Gleichgewichts-Thermodynamik. Die Gleichung für Brownsche Bewegung weiter oben ist ein Spezialfall.

Eine wesentliche Voraussetzung der Herleitung ist ein Kriterium zur Aufteilung der Freiheitsgrade in die Kategorien langsam und schnell. Z. B. wird in einer Flüssigkeit lokales thermisches Gleichgewicht innerhalb weniger Stoßzeiten erreicht. Es dauert dagegen viel länger, bis Dichten von Erhaltungsgrößen wie Masse oder Energie zum Gleichgewichtswert relaxieren. Dichten von Erhaltungsgrößen und insbesondere deren Anteile mit großen Wellenlängen sind also Kandidaten für langsame Variablen. Technisch wird die Unterteilung mit dem Zwanzig-Projektionsoperator[5] realisiert, dem wesentlichen Werkzeug der Herleitung. Die Herleitung verwendet wenige (plausible) Annahmen welche in ähnlicher Weise auch anderswo in der statistischen Physik benötigt werden und ist daher nicht streng mathematisch.

Es bezeichne   die langsamen Variablen. Die generische Langevingleichung lautet dann

 

Die fluktuierende Kraft   gehorcht einer Normalverteilung mit Korrelationsfunktion

 

Dies impliziert die Onsagersche Reziprozitätsbeziehungen   für die Dämpfungskoeffizienten  . Die Abhängigkeit   von   von   ist in den meisten Fällen vernachlässigbar. Das Symbol   bezeichnet die Hamiltonian des Systems,   ist die Gleichgewichts-Wahrscheinlichkeitsverteilung der Variablen  . Die Klammer   schließlich ist die Projektion der Poissonklammer der langsamen Variablen   und   in den Raum der langsamen Variablen.

Im Fall der Brownschen Bewegung hätte man  ,   oder   und  . Die Bewegungsgleichung   für   ist exakt, es gibt keine fluktuierende Kraft   und keinen Dämpfungskoeffizienten  .

Beispiele

Bearbeiten
 
Ein elektrischer Stromkreis mit einem Widerstand und einem Kondensator.

Widerstandsrauschen

Bearbeiten

Die Abbildung rechts zeigt einen elektrischen Stromkreis bestehend aus einem Widerstand mit elektrischem Widerstand   und einem Kondensator mit Kapazität  . Die langsame Variable ist die Spannung   am Kondensator oder Widerstand. Die Hamiltonian lautet  , und die Langevingleichung wird

 

Die daraus folgende Korrelationsfunktion

 

wird zu einem weißen Rauschen (Nyquist-Rauschen, Johnson-Rauschen) wenn die Kapazität klein wird.

Überdämpfte Brownsche Bewegung

Bearbeiten

Ein Grenzfall der klassischen Langevingleichung ist der Fall starker Reibung (oder kleiner Masse  ). In diesem Fall ist der Trägheitsterm   vernachlässigbar und man spricht von Überdämpfung. Die Langevingleichung erhält die Form

 

Integration und Mittelung des Quadrats über die fluktuierende Kraft liefert   wo   die Raumdimension ist. Mit Hilfe dieser Gleichung kann bei bekanntem Reibungskoeffizienten aus experimentell gemessener Diffusion die Boltzmannkonstante   bestimmt werden. Die nicht überdämpfte Langevingleichung liefert übrigens dasselbe Diffusionsgesetz.

Mikrodrehwaage (harmonischer Oszillator)

Bearbeiten

Hierbei handelt es sich um einen an einem Quarzfaden hängenden Spiegel. Ein am Spiegel reflektierter Lichtstrahl ermöglicht genaue Messungen. Die langsamen Variablen sind Drehimpuls   und Auslenkungswinkel   mit Hamiltonian

 

Die Dämpfungskonstante ist über den Luftdruck adjustierbar. Mit Hilfe der entsprechenden Langevingleichung und Experimenten erhält man einen genauen Wert für die Boltzmannkonstante. Umgekehrt beschreibt die Langevingleichung, wie thermische Fluktuationen der Empfindlichkeit von Zeigerinstrumenten eine Grenze setzen.

Kritische Dynamik

Bearbeiten

Der Ordnungsparameter eines kontinuierlichen Phasenübergangs ändert sich in der Nähe des kritischen Punktes nur mehr langsam und seine Dynamik ist mit einer i. A. nichtlinearen Langevingleichung beschreibbar.[4] Der einfachste Fall ist die Universalitätsklasse „Modell A“ mit einem nicht erhaltenen skalaren Ordnungsparameter  , realisiert z. B. in einem axialen Ferromagneten. Die Langevingleichung lautet

 

Andere Universalitätsklassen (die Nomenklatur ist „Modell A“, …, „Modell J“)[4] haben einen diffundierenden Ordnungsparameter, Ordnungsparameter mit mehreren Komponenten, zusätzliche andere langsame Variable und meistens auch Beiträge von Poissonklammern.

Anderson-Lokalisierung in einer Dimension

Bearbeiten

Langevingleichungen treten manchmal in unerwartetem anderem Zusammenhang auf. Ein Beispiel ist die Anderson-Lokalisierung. Die zeitunabhängige Schrödingergleichung für ein Teilchen der Masse   in einem Potential   in einer Dimension lässt sich schreiben

 

Hierbei ist   die Wellenfunktion,   das Wirkungsquantum und   ein Energie-Eigenwert. Mit   wird die Schrödingergleichung zu[6]

 

Wenn nun das Potential von zufällig verteilten Verunreinigungen herrührt, so dass es mit einem Gaußschen weißen Rauschen mit   und   beschreibbar ist, dann ist die Gleichung für   eine Langevingleichung, wobei die Koordinate   die Rolle der Zeit spielt.

Äquivalente Techniken

Bearbeiten

In der Regel ist eine Lösung einer Langevingleichung zu einer bestimmten Realisierung der fluktuierenden Kraft uninteressant. Stattdessen interessiert man sich für Korrelationsfunktionen der langsamen Variablen nach Mittelung über die fluktuierende Kraft. Solche Mittelwerte lassen sich auch auf anderem Weg erhalten.

Fokker-Planck-Gleichung

Bearbeiten

Eine Fokker-Planck-Gleichung ist eine deterministische Gleichung für die zeitabhängige Wahrscheinlichkeitsverteilung   der langsamen Variablen  . Die der generischen Langevingleichung weiter oben entsprechende Fokker-Planck-Gleichung lässt sich mit Standard-Techniken erhalten (z. B. Ref.[7]),

 

Die Gleichgewichtsverteilung   ist eine stationäre Lösung.

Pfadintegral

Bearbeiten

Ein zu einer Langevingleichung äquivalentes Pfadintegral kann man aus der entsprechenden Fokker-Planck-Gleichung oder direkt aus der Langevingleichung[8] erhalten. Das der generischen Langevingleichung entsprechende Pfadintegral ist

 

wo   ein Normierungsfaktor ist und

 

Die Hilfsvariablen   heißen Antwort-Variablen. Die Pfadintegraldarstellung vereinfacht die Anwendung von Techniken der Quantenfeldtheorie, wie etwa der Störungsrechnung und der Renormierungsgruppe.

Simulationen

Bearbeiten

Simulationen der Langevin-Gleichung sind eine Art von Monte-Carlo-Simulation[9]

Literatur

Bearbeiten
  • Don S. Lemons, Anthony Gythiel: Paul Langevin’s 1908 paper „On the Theory of Brownian Motion“ [„Sur la théorie du mouvement brownien,“ C. R. Acad. Sci. (Paris) 146, 530–533 (1908)], Am. J. Phys. 65, 1079 (1997), DOI:10.1119/1.18725
  • N.G. Van Kampen: „Stochastic Processes in Physics and Chemistry.“. 3. Auflage. North Holland, 2007.
  • Schwabl, Franz: Statistische Mechanik. Springer, ISBN 3-540-31095-9
  • Huang, Kerson: Statistical Mechanics. Wiley, ISBN 978-81-265-1849-4
  • Huang, Kerson: Introduction to Statistical Physics. CRC Press, ISBN 0-7484-0942-4

Einzelnachweise

Bearbeiten
  1. N. G. Van Kampen: Stochastic Processes in Physics and Chemistry. Elsevier, 2007, ISBN 978-0-444-52965-7, doi:10.1016/b978-0-444-52965-7.x5000-4.
  2. K. Kawasaki: Simple derivations of generalized linear and nonlinear Langevin equations. In: J. Phys. A: Math. Nucl. Gen. 6. Jahrgang, 1973, S. 1289, doi:10.1088/0305-4470/6/9/004, bibcode:1973JPhA....6.1289K.
  3. arxiv:1506.02650v2
  4. a b c P. C. Hohenberg, B. I. Halperin: Theory of dynamic critical phenomena. In: Reviews of Modern Physics. 49. Jahrgang, Nr. 3, 1977, S. 435?479, doi:10.1103/RevModPhys.49.435, bibcode:1977RvMP...49..435H.
  5. R. Zwanzig: Memory effects in irreversible thermodynamics. In: Phys. Rev. 124. Jahrgang, Nr. 4, 1961, S. 983?992, doi:10.1103/PhysRev.124.983, bibcode:1961PhRv..124..983Z.
  6. B. I. Halperin: Green's Functions for a Particle in a One-Dimensional Random Potential. In: Phys. Rev. 139. Jahrgang, 1A, 1965, S. A104, doi:10.1103/PhysRev.139.A104.
  7. S. Ichimaru: Basic Principles of Plasma Physics. Benjamin, 1973, ISBN 3-540-61236-X.
  8. H. K. Janssen: Lagrangean for Classical Field Dynamics and Renormalization Group Calculations of Dynamical Critical Properties. In: Z. Phys. B. 23. Jahrgang, 1976, S. 377, doi:10.1007/BF01316547, bibcode:1976ZPhyB..23..377J.
  9. Mikio Namiki: Stochastic Quantization. Springer Science & Business Media, 2008, ISBN 978-3-540-47217-9, S. 176 (englisch, google.com).