Umrechnung zwischen stetigen und diskreten Zustandssummen

Ich bin etwas verwirrt über die beiden Notationen für Partitionsfunktionen, dh die stetige Form

Z = 1 H 3 D P D R e β H ( P , R )

und die diskrete Form

Z = ich e β H ich .

Zunächst einmal, wie kommt man von einem dieser Ausdrücke zum jeweils anderen?

Angenommen, wir haben ein einzelnes Teilchen auf einem diskreten Gitter (jede Gitterzelle mit einer Ortslänge von A ) ohne kinetische Energie. Wir können dann den 2. Ausdruck verwenden und erhalten die Partitionssumme als

Z = ich e β v ich .

mit v ich die potentielle Energie ( T ich = 0 angenommen) vor Ort ich und die Summe über alle Gitterplätze. Wie würde der kontinuierliche Ausdruck aussehen?

Z = 1 H 3 ( D P ) ( D R e β v ( R ) ) = 1 H 3 ( D P ) ( ich e β v ich A 3 )

und können wir diesen Ausdruck wieder in eine Summe umwandeln, um auf den diskreten Ausdruck zurückzukommen?

Antworten (5)

Für einen Physiker ist ein Integral oft nichts anderes als eine „kontinuierliche Summe“. Die Anleitung X F ( X ) bedeutet oft eine tatsächliche Summe, wenn X aus einer abzählbaren Menge gezogen wird, aber ein ganzzahliges When X wird aus einer unzählbaren Menge gezogen.

Daher gibt es keine "Umwandlung" zwischen Ihren Ausdrücken - sie alle sind gültige Möglichkeiten, die Partitionsfunktion für verschiedene Arten von Systemen aufzuschreiben:

  • Klassisches statistisches System, endlich viele Mikrozustände ich { 1 , , N } mit Energie E ich : Z = ich e β E ich

  • Klassisches statistisches System, unendlich viele kontinuierliche Mikrozustände, gegeben durch Punkte ( Q , P ) R 6 N im Phasenraum für N Partikel: Z = H 3 e β H ( Q , P ) D 3 N X D 3 N P .

Das Integral ist einfach die einzige Möglichkeit, etwas zu bekommen, das einer "Summe über alles" nahe kommt ( X , P ) " das mathematisch wohldefiniert ist und sich wie gewünscht verhält.

Danke, aber sollte es nicht eine Möglichkeit geben, die Ausdrücke umzuwandeln, wenn sie konsistent sind? Angenommen, wir verringern die Gitterzellenlänge in dem Beispiel, das ich oben erklärt habe, dh nehmen die Grenze lim A Z D ich S C R e T e sollte dies nicht wieder den Continuos-Ausdruck ergeben?
@Guiste Sie können sie nicht "konvertieren", da sie für verschiedene Systeme gelten . Natürlich sollte an der Grenze, wo das diskrete System stetig wird, seine Zustandssumme zur Zustandssumme des kontinuierlichen Systems werden. Leider sind Details von Kontinuumsgrenzen schwierig und ich vermute, dass man dies tatsächlich nicht für alle Systeme garantieren kann. Warum wollen Sie diese Ausdrücke "konvertieren"?
Ich habe Probleme mit dem Beispiel, das ich unten in der Frage erklärt habe. Ich habe eine Referenz, die die diskrete Summe verwendet, um die Partitionsfunktion abzuleiten. Allerdings möchte ich aus summentechnischen Gründen von der kontinuierlichen Summe ausgehen. Da sich das System auf einem Gitter befindet, sollte es möglich sein, das Integral über den Raum als Summe über Gitterplätze mit dem jeweiligen Gewicht auszuführen (ich habe es in die Frage eingefügt). Dies sollte dann vergleichbar mit der Lösung der Referenz sein, die direkt die diskrete Summe nimmt ...

Gute Frage. Der Punkt ist das Rezept von Gibbs:

ρ e β H ,
um die Wahrscheinlichkeitsverteilung für ein klassisches System im thermischen Gleichgewicht zu erhalten, muss durch eine Vorschrift ergänzt werden, die uns sagt, wie wir Zustände zählen sollen, sodass die genaue Version des Rezepts lauten sollte:
D P ( S ) = e β H ( S ) D N ( S ) .
Mit anderen Worten, wir sollten immer ein Maß angeben D N ( S ) auf dem Raum der zulässigen Zustände S des Systems, für das der Boltzmann-Faktor e β H ( S ) wirkt als Dichte.

Dieser Punkt ist nicht ganz trivial. Zum Beispiel das kanonische Maß D P D Q für ein hamiltonisches System mit Koordinaten P Und Q hat den Vorteil, dass es in Bezug auf die kanonische Transformation forminvariant ist, und ist der kinematischen Beschreibung des Systems in gewisser Weise eigen (ich denke, das ist der Grund, es als "kanonisch" zu bezeichnen). Daher ist es ein gutes, vielleicht das einzige, Kandidatenmaß für das Zählen von Zuständen.

Allerdings kanonisch ( P , Q ) -System sind nicht die einzigen klassischen Systeme, die es gibt. Vergleichen Sie zum Beispiel, was Feynman hier tut, um die Polarisation für ein Gas von Molekülen mit einem intrinsischen Dipol in einem externen elektrischen Feld zu berechnen. Vor Gl. ( 1.15 ) , Er behauptet, dass

[...] die relative Anzahl der Moleküle mit der potentiellen Energie U ist proportional zu"

e U k T .

Anschließend spezialisiert er sich auf das jeweilige Problem:

[...] die Anzahl der Moleküle bei θ pro Raumwinkeleinheit ist proportional zu e U k T .

Hier die Staaten S sind die möglichen Orientierungen eines Einheitsvektors im dreidimensionalen Raum, dh der Einheitskugel S 2 . Die natürliche Wahl für das Maß ist eindeutig der Raumwinkel D Ω , die die Hypothese widerspiegelt, dass jede Orientierung im Raum a priori äquivalent sein sollte (oder anders ausgedrückt, die Existenz des elektrischen Feldes, das zu ungleichen Orientierungen im Raum führt, muss vollständig durch den Boltzmann-Faktor berücksichtigt werden E ).

Auch beim Zählen von Zuständen von kanonischen (d. h. P , Q )-Systemen gibt es eine schon Gibbs bekannte Mehrdeutigkeit, die sich aus der Frage ergibt, wie man Zustände identischer Teilchen zählen kann, und die zu den 1 N ! Faktoren in der Zustandssumme.

Jedenfalls sind diese Probleme meines Wissens allesamt etwas akademisch, da in der statistischen Quantenmechanik die Angabe des Dichteoperators völlig eindeutig ist:

ρ ^ = 1 Z e β H ^ .
Daher hat man zumindest im Prinzip einen soliden Weg, um das "klassische" Maß zu erhalten D N ohne auf all diese (zugegebenermaßen schwammigen) argumente einzugehen: man muss das maß so wählen, dass es die quantenzustandssumme im grenzbereich hoher temperaturen wiedergibt.


In Anbetracht dessen denke ich, dass das, was Sie in Ihrem speziellen Beispiel versuchen, nicht viel Sinn macht, weil Ihr System kein natürliches hat ( P , Q ) Koordinaten. Vielleicht könnten Sie Ihr Problem als Annäherung bei niedriger Temperatur für ein Teilchen auf einem Gitter mit einem Nahbereichsbindungspotential betrachten v ich an jedem Ort, so dass wegen der niedrigen Temperatur die einzigen zugänglichen Zustände die Grundzustände jedes Potentials sind.

Das wird so aussehen, als käme es aus dem linken Feld, aber ich verspreche, es ist relevant!

Wahrscheinlichkeit und diskrete/kontinuierliche Zufallsvariablen.

Wahrscheinlich ist es sehr üblich, mit kontinuierlichen Zufallsvariablen zu arbeiten. Die Definition davon ist, dass sie eher eine Wahrscheinlichkeitsdichtefunktion als diskrete Wahrscheinlichkeiten für jedes gegebene Ergebnis haben. Die Definition ist also so F ist die Wahrscheinlichkeitsdichtefunktion für eine Zufallsvariable X wenn für ein kleines Δ X , die Wahrscheinlichkeit, dass X < X < X + Δ X wird von gegeben F ( X )   Δ X . Das ist die Definition.

Dann können wir alle Wahrscheinlichkeiten zusammenzählen, um eine zu bekommen, D X   F ( X ) = 1 , und wir können auch damit beginnen, Dinge wie die Erwartungsspanne zu definieren, Q ( X ) = D X   F ( X )   Q ( X ) . Für etwas Übung mit diesen Ideen, bevor wir weitermachen:

  1. Beweisen Sie die übliche Abkürzungsformel für die Varianz: dass wenn X ¯ = X und wir definieren Var [ X ] = ( X X ¯ ) 2 , Dann Var [ X ] = X 2 X 2 .
  2. Gegeben zwei Zufallsvariablen, X , Y , wir haben einen "gemeinsamen Vertrieb" F ( X , j ) so dass X < X < X + Δ X , j < Y < j + Δ j wird von gegeben F ( X , j )   Δ X   Δ j . Nun könnte dies als faktorisierbar sein G ( X ) H ( j ) bei "unabhängigen Zufallsvariablen". Beweisen Sie das dann X Y = X Y .
  3. Um zwei Zufallsvariablen zu addieren, Z = X + Y , können Sie die gemeinsame Verteilung falten; H ( z ) = D X   F ( X , z X ) . Beweisen Sie, dass der Erwartungswert von a Z wird von gegeben X + Y = X + Y und das für unabhängige Zufallsvariablen, Var [ X + Y ] = Var [ X ] + Var [ Y ] .

Okay, jetzt wo Sie wissen, was diese Definition ist, haben wir Physiker ein grobes Werkzeug, um diese kontinuierlichen Wahrscheinlichkeiten zu verwenden, um diskrete Wahrscheinlichkeiten zu modellieren. Es ist eine Familie glatter Funktionen, die technisch nicht existieren, genannt "Dirac δ -Funktion." In Worten: δ ( X ) ist eine unendlich hohe, unendlich dünne Spitze, die bei Null zentriert ist, so dass das Integral darüber eins ist. In Mathematik,

u v D X   F ( X )   δ ( X X 0 ) = { F ( X 0 )  Wenn  u < X 0 < v , F ( X 0 )  Wenn  v < X 0 < u , 0  ansonsten.
(Obwohl Sie einen unbestimmten Fall erhalten, wenn u = X 0 oder v = X 0 . )

Einer davon kann beispielsweise ein Grenzwert eines Gaußschen Integrals mit Einheitsfläche sein,

δ σ ( X ) = 1 2 π σ 2   e X 2 / ( 2 σ 2 ) ,
und Sie können sich vorstellen, dass wir die gesamte mathematische Theorie mit diesem Ding als glatte Funktion ( σ > 0 ) aber am Ende nehmen wir immer eine Grenze σ 0 um die diskreten Antworten zu erhalten, die wir erwarten. Allerdings gibt es viele Funktionen, die sich auf den Dirac beschränken δ -Funktion, es ist nicht nur die Gaußsche.

Jetzt, da Sie diesen Trick kennen, muss ich Ihnen nur noch sagen, dass Ihr normaler sechsseitiger Würfel eine Wahrscheinlichkeitsdichtefunktion hat,

1 6 [ δ ( X 1 ) + δ ( X 2 ) + δ ( X 3 ) + δ ( X 4 ) + δ ( X 5 ) + δ ( X 6 ) ] .
Und dann, wenn Sie trainieren X Sie werden Ihren Vertrauten bekommen 1 / 6 + 2 / 6 + 3 / 6 + 4 / 6 + 5 / 6 + 6 / 6 = 3.5 Ausdruck, wie bei allen anderen Ausdrücken.

Ensembles und geheime Distributionen

Jetzt behaupte ich Ihnen, dass wir insgeheim eine Wahrscheinlichkeitsverteilung im Phasenraum gewählt haben, und das bedeutet, dass wir ... eine Wahrscheinlichkeitsdichtefunktion haben!

Wenn wir in der statistischen Mechanik von der Annahme der Entropiemaximierung ausgehen, müssen wir in einem ganz bestimmten Umstand beginnen, und diese Umstände werden "Ensembles" genannt. Dasjenige, von dem wir ausgehen müssen, ist das "mikrokanonische Ensemble", es sagt: "Es gibt eine sehr solide Box mit festem Volumen, keine Partikel gehen in die Box hinein oder aus ihr heraus, keine Energie geht in die Box hinein oder aus ihr heraus." Und dann sagen wir in diesem speziellen Fall, dass die Wahrscheinlichkeitsverteilung für ein Bündel von Teilchen ich ist flach über dem darunter liegenden Phasenraum, F ( X 1 , X 2 , , P 1 , P 2 , ) = C . Die winzige Konstante C muss unter Umständen am Ende sogar auf 0 begrenzt werden! Aber das ist in Ordnung, was wir wirklich meinen, ist, dass die relativen Wahrscheinlichkeiten, A oder B zu sehen, proportional zu ihren Phasenraumvolumina sind.

Die Tatsache, dass dies konstant ist , ist der Grund, warum es nicht in Ihrem ersten Ausdruck für die Partitionsfunktion erscheint, aber ich behaupte, dass es da ist und dass die generische Form vorhanden ist

Z = D D P   D D X   F ( X 1 , , X D , P 1 , , P D )   e β H ( X 1 , , P 1 , ) .
Sie sehen es nicht, weil es eine Konstante ist, die aus dem Integral herausgezogen und stillschweigend in den anderen Teil der Wahrscheinlichkeitsrechnung geschoben wurde (denken Sie daran, dass Wahrscheinlichkeiten immer sind e β   H / Z ...).

Wenn Sie eine Summe von gleichgewichteten Dirac schlagen δ -Funktionen in für F , entsprechend dem diskreten Wahrscheinlichkeitsraum von "dieses Ding kann in Konformation X oder Konformation Y sein, aber es gibt keinen kontinuierlichen Pfad zwischen den beiden), erhalten Sie stattdessen Ihren zweiten Ausdruck, wieder mit einer multiplikativen Konstante, die aus der Definition für gezogen wird der Bequemlichkeit halber.

Die folgende Annahme wird in einigen Büchern der statistischen Mechanik ganz am Anfang verwendet und hat nichts mit der kanonischen Gesamtheit zu tun:

"Die Phase kann in kleine Teile aufgeteilt werden C ich jeweils mit Volumen H 3 so dass H ist jeweils konstant C ich und nimmt Wert H ich drauf".

Mit dieser Annahme lässt sich die stetige Formel vereinfachen:

1 H 3 D P D R e β H = 1 H 3 ich C ich D P D R e β H = ich 1 H 3 C ich D P D R e β H ich = ich e β H ich 1 H 3 C ich D P D R

C ich hat Volumen H 3 ist äquivalent zu C ich D P D R = H 3

Endlich:

1 H 3 D P D R e β H = ich e β H ich

Es geht nur darum, eine stetige Funktion durch eine Stufenfunktion zu approximieren .

Zunächst einmal kann der Ausdruck für die Zustandssumme in der Quantenmechanik in Form eines Integrals ausgedrückt werden. Die Partitionsfunktion ist

Z = Tr   e β H .

Dieser Trace kann auf beliebiger Basis ausgewertet werden. In die Basis, wo der Hamiltonoperator diagonal ist, schreiben wir Z = ich e β E ich , Wo E ich sind die Eigenwerte der Energien und die Summe ist über alle Zustände. Allerdings auf anderer Basis, zum Beispiel der X -Basis werden wir schreiben,

Z = D X X | e β H | X .
Mein Punkt ist, dass wir selbst in der Quantenmechanik eine integrale Form erhalten, wenn wir die Berechnung auf kontinuierlicher Basis durchführen.

Nun, ich hoffe, Sie stimmen mir zu, dass die Frage nicht wirklich lautet, wie man von einer diskreten Form zu einer kontinuierlichen Form kommt, sondern eher, wie man von einer quantenmechanischen Form zu einer klassischen Form kommt. Ich denke, der einfachste Weg, diese Verbindung zu erkennen, besteht darin, einen Pfad zu formulieren, der integral zur Partitionsfunktion ist, und die klassische Grenze zu nehmen. Dies wurde in dem klassischen Buch von Feynman und Hibbs in Kapitel 9 über statistische Mechanik diskutiert.