Wie wird die Dissipationsfunktion des Fluktuationssatzes zur Entropie?

In "The Fluctuation Theorem" von Evans and Searles leiten sie den Satz der transienten Fluktuation aus dem Satz von Liouville ab (S. 1541). Nach ihrer Notation Γ = ( q , p ) Γ = ( q , p ) Verwenden sie Liouvilles Theorem und die Kontinuitätsgleichung, um daraus zu schließen

d f d t 0 d f d t = f t + d f d Γ ⋅ ⋅ ˙ = f t + d d Γ ( f Γ ˙ ) = - f d d Γ ⋅ ⋅ ˙ . d f d t = f t + d f d Γ Γ ˙ 0 = f t + d d Γ ( f Γ ˙ ) d f d t = - f d d Γ Γ ˙ .

Sie definieren dann Λ ( Γ ) = d d Γ ⋅ ⋅ ˙ Λ ( Γ ) = d d Γ Γ ˙ und leiten Sie den Fluktuationssatz her, der die Dissipationsfunktion definiert als

Ω t t = t 0 Ω ( Γ ( s ) )   d s ln [ f ( Γ ( 0 ) , 0 ) f ( Γ ( t ) , 0 ) ] - t 0 Λ ( Γ ( s ) )   d s . Ω t t = 0 t Ω ( Γ ( s ) ) d s ln [ f ( Γ ( 0 ) , 0 ) f ( Γ ( t ) , 0 ) ] - 0 t Λ ( Γ ( s ) ) d s .

Wikipedia definiert,

Ω t ( Γ ) = t 0 d s Ω ( Γ ; s ) ln [ f ( Γ ( 0 ) , 0 ) f ( Γ ( t ) , 0 ) ] + ΔQ ( Γ ; t ) k T . Ω t ( Γ ) = 0 t d s Ω ( Γ ; s ) ln [ f ( Γ ( 0 ) , 0 ) f ( Γ ( t ) , 0 ) ] + Δ Q. ( Γ ; t ) k T .

Welche Magie erlaubt diesen Sprung? Warum fehlt Wikipedia der Faktor t t vor?

Eine Sache, die mir aufgefallen ist, ist die für nicht konservative Systeme

q ˙ ich = H p ich + F q   ;   p ˙ ich = - H q ich + F p q ˙ ich = H p ich + F q ; p ˙ ich = - H q ich + F p

Λ ( Γ ( s ) ) Δ Q ( Γ ; t ) k T = d d Γ ⋅ ⋅ ˙ = ich ( 2 H q ich p ich - 2 H p ich q ich + ( F q q ich + F p p ich ) ) = ich ( F q q ich + F p p ich ) = - t 0 Λ ( Γ ( s ) )   d s = - ich t 0 ( F q q ich + F p p ich ) . Λ ( Γ ( s ) ) = d d Γ Γ ˙ = ich ( 2 H q ich p ich - 2 H p ich q ich + ( F q q ich + F p p ich ) ) = ich ( F q q ich + F p p ich ) Δ Q. ( Γ ; t ) k T = - 0 t Λ ( Γ ( s ) ) d s = - ich 0 t ( F q q ich + F p p ich ) .

Also für rein konservative Systeme F q = F p = 0 F q = F p = 0 ,

P r ( Ω t = A ) P r ( Ω t = - A ) = 1 . P r ( Ω t = EIN ) P r ( Ω t = - EIN ) = 1 .

Die Entropieerzeugung erfolgt also durch Forcierung oder Dissipation.

Die Kompressibilität des Phasenraums (dh d f / d t 0 d f / d t 0 ) impliziert selbst Irreversibilität, die oft als Entropieproduktion interpretiert wird (obwohl ich beachten sollte, dass Entropie und Irreversibilität möglicherweise nicht synonym sind).
@honeste_vivere Danke! Wie würden Sie in diesem Fall die Entropie berechnen? Die normale Boltzmann / Shannon-Entropie - f Log f - f Log f scheint nicht so, als könnte es jemals damit zusammenhängen Λ ( Γ ( s ) ) d s Λ ( Γ ( s ) ) d s da Λ = d d Γ ⋅ ⋅ ˙ Λ = d d Γ Γ ˙ hat keine Faktoren von f f
Ich frage, ob Sie es mit einer Art Summierung berechnen müssen ich G ( q ich , p ich ) ich G ( q ich , p ich ) über die Teilchen (dh Teilchenentropie, nicht über ähnlich präparierte Ensembles verteilt wie f ( Γ ) f ( Γ ) ), dann könnte ich mir möglicherweise vorstellen, eine Verbindung herzustellen.
Sie können etwas verwenden, das der Störungstheorie und / oder der mittleren Feldtheorie ähnelt, in der Sie konstruieren / annehmen f ⟨f + Δ f f f + δ f und dann die normale berechnen f   Log f f Log f zur ersten Bestellung. Technisch gesehen haben Sie durch Auferlegen / Unterstellen des Fluktuationsdissipationssatzes Irreversibilität in die Gleichungen "eingefügt". Dies ist eine der Konsequenzen der Zufallsphasenapproximation.
Wie hilft das, wenn es keine gibt? f f im Γ Γ ? Oder doch? f f irgendwie aufheben, um das zu bekommen Γ Γ nur Ausdruck?

Antworten (1)

Boltzmann, Wlasow und Entropie

Die Vlasov-Gleichung ist die kollisionsfreie Form der Boltzmann-Gleichung . Die Vlasov-Gleichung kann wie folgt geschrieben werden:

f s t + v f s + F m s ⋅ ⋅ v f s = 0 (1) (1) f s t + v f s + F m s v f s = 0
woher f s = f s ( x , v , t ) f s = f s ( x , v , t ) ist die Teilchengeschwindigkeitsverteilungsfunktion der Spezies s s (zB Maxwellian ), F F ist eine äußere Kraft und die k- te Komponente von v v ist gegeben durch k ^   / v k k ^ / v k . Ändern wir Gleichung 1 mit der Lorentzkraft für F F e m F F e m und linearisieren , dh annehmen Q ⟨Q⟩ + δ Q. Q. Q. + δ Q. , woher   ist ein Ensemble Durchschnitt und ⟨Δ Q⟩ = 0 δ Q. = 0 . Daraus können wir eine durchschnittliche und schwankende Form von Gleichung 1 konstruieren. Ich werde den Index fallen lassen s s aus Faulheit für den Rest der Ableitung.

Die durchschnittliche Form ergibt sich aus:

⟨F t + v ⟨f + ⟨F e m m ⋅ ⋅ v ⟨F = C (2a) (2a) f t + v f + F e m m v f = C
woher C C ist ein Proxy für einen Kollisionsbegriff von:
C = - ⟨δ   F e m m ⋅ ⋅ v δ f (2b) (2b) C = - δ F e m m v δ f
Die schwankende Form ist gegeben durch:
t δ f + v δ f + ⟨F e m ⟩⟩ v δ f + δ F e m ⋅ ⋅ v ⟨F = - δ F e m ⋅ ⋅ v δ f - C (3) (3) t δ f + v δ f + F e m v δ f + δ F e m v f = - δ F e m v δ f - C

Es gibt einen wichtigen Unterschied zwischen dem Kollisionsbegriff, C C und der klassische binäre Kollisionsterm, der in der Boltzmannschen Gleichung gefunden wird. C C bewahrt nicht den lokalen Impuls und die Energiedichte der Partikel. Wichtiger, C C bewirkt nicht, dass sich das Plasma zu einem Maxwellschen entspannt [siehe z. B. Tidman und Krall , 1971].

Gibbs 'Entropie kann wie folgt geschrieben werden:

S = - k B   a l l   d x   d v d x   d v   ⟨F   ln | ⟨F | (4) (4) S = - k B ein l l d x d v d x d v f ln | f |
Gibbs erkannte, dass Gleichung 4 zur negativen Unendlichkeit divergieren kann, wenn man nicht die Einteilchenwahrscheinlichkeitsverteilung verwendet , die ein Durchschnitt über Ensemblezustände ist (also die   s). Eine interessante Beobachtung ist, dass die Liouville-Gleichung auch vorhersagt, dass Gleichung 4 ohne die Verwendung von Ensemble-Durchschnitten zur negativen Unendlichkeit abweicht [z. B. Evans und Morriss , 1990]. Der Grund für die Divergenz liegt darin, dass Störungen im Phasenraum immer kleiner werden, was immer höher dimensionale Formen erfordert f ( x , v , t ) f ( x , v , t ) . Eine der Konsequenzen (möglicherweise indirekt?) Dieses intuitiven Sprunges war die Entwicklung von Dingen wie der Mittelfeldtheorie .

Liouvilles Gleichung

Der erste Teil des Folgenden stammt aus meiner Antwort unter https://physics.stackexchange.com/a/177972/59023 .

Wir wissen das ⟨F f erfüllt Liouvilles Gleichung , oder besser gesagt, ⟨F f / t = 0 t = 0 . Im Allgemeinen lautet die Bewegungsgleichung:

f t = f [ ( q d q d t ) + ( p d p d t ) ] + [ d q d t ⋅ ⋅ f q + d p d t ⋅ ⋅ f p ] (5) (5) f t = f [ ( q d q d t ) + ( p d p d t ) ] + [ d q d t f q + d p d t f p ]
wo ich den kanonischen Phasenraum von definiert habe ( q , p ) ( q , p ) . Wenn ich die Begriffe vereinfache d Q / d t d Q. / d t zu Q. ˙ Q. ˙ und lass Γ = ( q , p ) Γ = ( q , p ) dann finde ich:
f t = - f Γ ⋅ ⋅ ˙ - Γ ˙ ⋅ ⋅ f Γ = - Γ ( Γ ˙   f ) (6a) (6b) (6a) f t = - f Γ Γ ˙ - Γ ˙ f Γ (6b) = - Γ ( Γ ˙ f )
wo man sieht, dass die letzte Form wie die Kontinuitätsgleichung aussieht. Wenn ich die Gesamtzeitableitung wie folgt definiere:
d d t = t + Γ ˙ ⋅ ⋅ Γ (7) (7) d d t = t + Γ ˙ Γ
dann kann ich zeigen, dass die zeitliche Änderungsrate der Verteilungsfunktion gegeben ist durch:
d f d t = f t + Γ ˙ ⋅ ⋅ f Γ = - [ f   Γ ⋅ ⋅ ˙ + Γ ˙ ⋅ ⋅ f Γ ] + Γ ˙ ⋅ ⋅ f Γ = - f   Γ ⋅ ⋅ ˙ - f   Λ ( Γ ) (8a) (8b) (8c) (8d) (8a) d f d t = f t + Γ ˙ f Γ (8b) = - [ f Γ Γ ˙ + Γ ˙ f Γ ] + Γ ˙ f Γ (8c) = - f Γ Γ ˙ (8d) - f Λ ( Γ )
woher Λ ( Γ ) Λ ( Γ ) wird der Phasenraumkomprimierungsfaktor genannt [z. B. Evans und Morriss , 1990]. Es ist zu beachten, dass die Gleichungen 8a bis 8d verschiedene Formen der Liouville-Gleichung sind, die ohne Bezugnahme auf die Bewegungsgleichungen erhalten wurden und die Existenz eines Hamilton-Operators nicht erfordern. Ich kann Gleichung 8d in der folgenden Form umschreiben:
d d t ln | f | = - Λ ( Γ ) (9) (9) d d t ln | f | = - Λ ( Γ )

Wenn die Bewegungsgleichungen aus einem Hamilton-Operator generiert werden können, dann Λ ( Γ ) = 0 Λ ( Γ ) = 0 auch in Gegenwart von äußeren Feldern, die das System aus dem Gleichgewicht bringen. Beachten Sie, dass die Existenz eines Hamiltonianers eine ausreichende, aber nicht notwendige Bedingung für ist Λ ( Γ ) = 0 Λ ( Γ ) = 0 .

Entropie-Produktion

Denken Sie daran, dass wir definiert haben ⟨Q⟩ / t = 0 Q. / t = 0 und wir wissen aus Gleichung 9, dass d / d t ln ⟨F = 0 d / d t ln f = 0 für konservative Systeme (dh solche mit inkompressiblem Phasenraum). Wenn wir definieren F ⟨f   ln | ⟨F | F f ln | f | Dann können wir zeigen, dass:

F t ≡ ≡ t ( ⟨F Ln | ⟨F | ) F t + ( v F ) 0 + ( v F ) = ( 1 + ln | ⟨F | ) ⟨F t = 0 = 0 = d F d t + F ( v ) = ⟨F Ln | ⟨F | ( v ) + ( 1 + ln | ⟨F | ) [ v ⟨f ] (10 A) (10b) (10c) (10d) (10 A) F t t ( f ln | f | ) = ( 1 + ln | f | ) f t = 0 (10b) F t + ( v F ) = 0 (10c) = d F d t + F ( v ) (10d) 0 + ( v F ) = f ln | f | ( v ) + ( 1 + ln | f | ) [ v f ]
Mit diesen Beziehungen kann man zeigen, dass:
( v ) F ( v ) ( ⟨f Ln | ⟨F | ) = - v F = - v ( ⟨f Ln | ⟨F | ) = - ⟨f V ⋅ ⋅ ln | ⟨F | - In | ⟨F | v ⟨f = - v ⟨f - In | ⟨F | v ⟨f = - ( 1 + ln | ⟨F | ) v ⟨f (11a) (11b) (11c) (11d) (11e) (11a) ( v ) F = - v F (11b) ( v ) ( f ln | f | ) = - v ( f ln | f | ) (11c) = - f v ln | f | - ln | f | v f (11d) = - v f - ln | f | v f (11e) = - ( 1 + ln | f | ) v f

Wir wissen auch, dass in der Grenze von x ± x ± , der Begriff ⟨F e m 0 F e m 0 weil wir annehmen, dass alle Steigungen in dieser Grenze asymptotisch gegen Null gehen und C C -term geht ebenfalls auf Null. Daher finden wir, dass Gleichung 2a mit der Menge arbeitet:

  d v   ( 1 + ln | ⟨F | ) (12) (12) d v ( 1 + ln | f | )
führt zu folgenden Ergebnissen:
[ t + v + ⟨F e m ⟩⟩ v ] d v   ( 1 + ln | ⟨F | ) [ v ] d v   ( 1 + ln | ⟨F | ) ∇ ∇ d v   v ⟨f Ln | ⟨F | = d v   C ( 1 + ln | ⟨F | ) = d v   C ( 1 + ln | ⟨F | ) = d v   C ( 1 + ln | ⟨F | ) (13a) (13b) (13c) (13a) [ t + v + F e m v ] d v ( 1 + ln | f | ) = d v C ( 1 + ln | f | ) (13b) [ v ] d v ( 1 + ln | f | ) = d v C ( 1 + ln | f | ) (13c) d v v f ln | f | = d v C ( 1 + ln | f | )
Nun fügen wir das Formular für ein C C aus Gleichung 2b zu finden:
∇ ∇   d v   v   ⟨F Ln | ⟨F | = e m d v   ⟨Δ F e m   δ f ⟩⟩ v ln | ⟨F | (14) (14) d v v f ln | f | = e m d v δ F e m δ f v ln | f |
wobei der Ausdruck auf der linken Seite die Divergenz des Entropieflusses ist .

Gleichung 14 ist ein Beispiel dafür, wie Entropie in einem kollisionsfreien Medium mit inkompressiblem Phasenraum erzeugt werden kann, der sich aus der Zeitabhängigkeit in ergibt ⟨F f eingeführt durch die Schwankungen ungleich Null, die durch das gegeben sind ⟨Δ F e m   δ f δ F e m δ f -Begriff. Auch wenn ⟨F f Entlang von Phasenverläufen aus der inkompressiblen Liouville-Gleichung erhalten, entwickelt es immer feinere Merkmale, die der Phasenmischung entsprechen. Die Einführung der Irreversibilität (dh hier auch der Entropie) in dieses System ist größtenteils eine Folge des Ansatzes, der dem Grobkörnungsverfahren entspricht, das zur Herleitung der Boltzmann-Gleichung aus der reversiblen Liouville-Gleichung verwendet wird [z. B. Evans und Morriss , 1990] ; Tidman und Krall , 1971].

Ein ähnlicher Ansatz kann verwendet werden, wenn wir nicht annehmen d / d t ln ⟨F = 0 d / d t ln f = 0 eine Form für die Entropie ableiten.

Verweise

  • Evans, DJ und G. Morriss Statistical Mechanics of Nonquilibrium Liquids, 1. Auflage , Academic Press, London, 1990.
  • Tidman, DA und NA Krall Stoßwellen in kollisionsfreien Plasmen , Wiley-Reihe in Plasmaphysik, New York: Wiley-Interscience, 1971.
Geniale Antwort! Es wird eine Weile dauern, bis ich mich eingehend damit befasst habe. Danke auch für die Hinweise