Brownsche Bewegung: Erwarteter Wert gleichmäßiger Potenzen der Geschwindigkeitskorrelationsfunktion

Bei einem Brownschen Teilchen, das einer Reibung ausgesetzt ist γ γ und ein weißes Gaußsches Rauschen ξ ( t ) ξ ( t ) mit Varianz σ σ sind die Langevin-Gleichungen

d v ( t ) d t = - γ m v ( t ) + 1 m ξ ( t ) d v ( t ) d t = v ( t ) { d v ( t ) d t = - - γ m v ( t ) + 1 m ξ ( t ) d v ( t ) d t = v ( t )

Es ist möglich, diese Gleichungen zu lösen und vorausgesetzt v ( 0 ) = 0 v ( 0 ) = 0 , erhalten

⟨V 2 ( t ) ξ = g 2 m γ ( 1 - e - 2 γ m t ) , v 2 ( t ) ξ = G 2 m γ ( 1 - - e - - 2 γ m t ) ,

wo G G kommt aus dem Delta-korrelierten Rauschen,

⟨Ξ ( t 1 ) ξ ( t 2 ) ξ = g δ ( t 2 - t 1 ) . ξ ( t 1 ) ξ ( t 2 ) ξ = G δ ( t 2 - - t 1 ) .

Frage: Was ich jetzt wissen möchte, ist, ob es möglich ist zu berechnen ⟨V 4 ( t ) ξ v 4 ( t ) ξ , ⟨V 6 ( t ) ξ v 6 ( t ) ξ , ⟨V 8 ( t ) ξ v 8 ( t ) ξ , etc ..., von ⟨V 2 ( t ) ξ v 2 ( t ) ξ . Ist es?

Ich weiß, dass wir die daraus resultierenden Integrale direkt mit dem Wickschen Theorem berechnen können, z. B. haben wir die folgende Identität, die zur Berechnung verwendet werden kann ⟨V 4 ( t ) ξ v 4 ( t ) ξ

⟨Ξ ( t 1 ) ξ ( t 2 ) ξ ( t 3 ) ξ ( t 4 ) ξ = ⟨Ξ ( t 1 ) ξ ( t 2 ) ξ ⟨Ξ ( t 4 ) ξ ( t 4 ) ξ + ⟨Ξ ( t 1 ) ξ ( t 3 ) ξ ⟨Ξ ( t 2 ) ξ ( t 4 ) ξ + ⟨Ξ ( t 1 ) ξ ( t 4 ) ξ ⟨Ξ ( t 2 ) ξ ( t 3 ) ξ ξ ( t 1 ) ξ ( t 2 ) ξ ( t 3 ) ξ ( t 4 ) ξ = ξ ( t 1 ) ξ ( t 2 ) ξ ξ ( t 4 ) ξ ( t 4 ) ξ + ξ ( t 1 ) ξ ( t 3 ) ξ ξ ( t 2 ) ξ ( t 4 ) ξ + ξ ( t 1 ) ξ ( t 4 ) ξ ξ ( t 2 ) ξ ( t 3 ) ξ

aber ich würde gerne wissen, ob wir den langen Weg der Berechnung all dieser Integrale vermeiden können, indem wir nur das Ergebnis für verwenden ⟨V 2 ( t ) ξ v 2 ( t ) ξ .

BEARBEITEN:

Ich habe die integrale "Brute Force" unter Verwendung des Wickschen Theorems berechnet ⟨V 4 ( t ) ξ v 4 ( t ) ξ und erhalten

⟨V 4 ( t ) ξ = 3 ( g 2 m γ ) 2 e - 4 γ m t ( e 2 γ m t - 1 ) 2 v 4 ( t ) ξ = 3 ( G 2 m γ ) 2 e - - 4 γ m t ( e 2 γ m t - - 1 ) 2

Meine Berechnung ist lang und deshalb veröffentliche ich nur das Endergebnis. Ich weiß nicht, ob dies korrekt ist, aber wenn Sie es zufällig wissen oder wissen, welche Referenz zu überprüfen ist, können Sie mir bitte signalisieren, ob dieses Ergebnis korrekt ist? Ich wäre dankbar!

Wick's Theorem ist der richtige Weg. Sie können wahrscheinlich ein rekursives Muster finden, um Zeit zu sparen.
Vielen Dank, dass Sie @DanielSank. Ich würde es zwar gerne finden, aber in der Zwischenzeit habe ich versucht, es für die vierte Potenz brutal zu machen. Ich weiß nicht, ob mein Ergebnis korrekt ist. Haben Sie eine Möglichkeit zu überprüfen? Es wäre sehr hilfreich. Ich werde die Frage bearbeiten und mein Ergebnis veröffentlichen. Bitte geben Sie mir ein Feedback! :) :)
Ich bin mir nicht sicher, warum Sie Integrale erstellen müssen. v ( t ) v ( t ) wird nur als Gaußscher Wert mit dem Mittelwert Null verteilt, sodass der vierte Moment dreimal so groß ist wie der zweite Moment im Quadrat usw., wie Ihr Brute-Force-Ergebnis zeigt.

Antworten (1)

Ich hoffe, Sie entschuldigen mich, wenn Sie meinen Kommentar als Antwort kopieren. Die Langevin-Gleichung ist linear, daher impliziert die Gaußsche Natur der Zufallskraft dies v ( t ) v ( t ) wird als Gaußscher verteilt. Daher ist bei einer Anfangsgeschwindigkeit von Null das vierte Moment dreimal so groß wie das Quadrat des zweiten Moments usw. Dies haben Sie bei der Berechnung Ihrer Brute Force festgestellt.

[Nach OP-Kommentar bearbeiten]

Mehr Details. Für eine Gauß-verteilte Variable werden alle Momente durch die ersten beiden bestimmt. Eine Tabelle finden Sie beispielsweise auf der Wikipedia-Seite im Abschnitt "Momente" (wo sie auch eine allgemeine Formel angeben). Für den Fall, dass die Verteilung den Mittelwert Null hat ⟨V ( t ) = 0 v ( t ) = 0 Es gelten die Formeln in der rechten Spalte dieser Tabelle (zentrale Momente): if ⟨V ( t ) 2 = Σ 2 v ( t ) 2 = σ 2 , dann ⟨V ( t ) 4 = 3 σ 4 v ( t ) 4 = 3 σ 4 , ⟨V ( t ) 6 = 15 σ 6 v ( t ) 6 = fünfzehn σ 6 usw. Wenn die Anfangsgeschwindigkeit nicht Null ist, bedeutet dies, dass ⟨V ( t ) = μ 0 v ( t ) = μ 0 gelten die Formeln in der mittleren Spalte: etwas komplizierter, aber bei weitem nicht so schlimm wie das tatsächliche Aufschreiben von Integralen über die Zufallskraft.

Eine andere Möglichkeit, dies auszudrücken, besteht darin, zu sagen, dass die Kumulanten höherer Ordnung der Gaußschen Verteilung Null sind (dh eine höhere Ordnung als der zweite Kumulant).

Danke! Ich fürchte, ich verstehe Ihre Argumentation nicht. Warum ist der vierte Moment dreimal so groß wie das Quadrat des zweiten Augenblicks? Und wie verallgemeinert sich dies für höhere Mächte?
Ich habe meine Antwort entsprechend bearbeitet!
Vielen Dank. Dann nur zur Verdeutlichung, da Sie das gleiche Symbol verwendet haben: Sie definieren ⟨V 2 ( t ) ξ = σ 2 v 2 ( t ) ξ = σ 2 , aber das ist nicht dasselbe σ σ für die Varianz der zufälligen Kraft, die ich oben definiert habe, richtig?
Tut mir leid, dass ich das nicht klargestellt habe. Ich habe nur benutzt σ 2 σ 2 (und μ μ ), um mit der Tabelle auf der Wikipedia-Seite übereinzustimmen, und natürlich meinte ich die Varianz (und den Mittelwert) der Variablen v ( t ) v ( t ) . Ich hätte sagen sollen, dass dies nicht der Parameter ist, der mit der Verteilung der zufälligen Kraft verbunden ist.