Energie des klassischen idealen Gases im großkanonischen Ensemble

Die kanonische Zustandssumme für ein ideales Gas ist

Z ( N , v , β ) = 1 N ! ( v λ 3 ) N
Wo λ = β H 2 2 π M die thermische De-Broglie-Wellenlänge ist. Es ist einfach zu bekommen
E = Protokoll Z β = 3 2 N k B T .

Aus Z die großkanonische Zustandssumme ist

Q ( μ , v , β ) = N = 0 1 N ! ( e β μ v λ 3 ) N = e e β μ v λ 3 .

Die durchschnittliche Teilchenzahl ist

N = Protokoll Q ( β μ ) = e β μ v λ 3 .

Um die durchschnittliche Energie zu erhalten, sollten wir ersetzen N ,

E = Protokoll Q β = 3 2 N k B T

aber das ist nur wahr, wenn wir das magisch ignorieren e β μ Faktor bei der Ableitung, sonst gibt es einen zusätzlichen (unsinnigen) Term. Ich habe ein paar Quellen überprüft und dies ist die akzeptierte Lösung (schließlich muss es diese sein, um mit dem kanonischen Ensemble-Ergebnis konsistent zu sein), obwohl sie das Problem auf mysteriöse Weise beschönigen, also fehlt mir etwas. Danke.

Antworten (5)

Hinweis: Sie haben einen Fehler in Ihren Berechnungen. Insbesondere im großkanonischen Ensemble,

E Protokoll Q β .

Außerdem habe ich gerade die ganze Berechnung durchgeführt, nachdem ich diesen Fehler auf die entsprechende Weise korrigiert habe, und es hat so funktioniert, wie es sollte.

Nachtrag, 2019-02-02. Details über den Hinweis hinaus

Schritt 1. Erinnern Sie sich an die folgenden Definitionen der großkanonischen Zustandssumme Q , die Ensemble-Durchschnittsenergie E im großkanonischen Ensemble und die durchschnittliche Teilchenzahl des Ensembles N im großkanonischen Ensemble. Alle Summen sind über Staaten ich vom System:

Q ich e β ( E ich μ N ich ) , E ich e β ( E ich μ N ich ) Q E ich , N ich e β ( E ich μ N ich ) Q N ich

Schritt 2. Zeigen Sie, dass die folgende Identität aus den Definitionen in Schritt 1 folgt:

E = ln Q β + μ N .

Schritt 3. Zeigen Sie das, wenn wir nehmen

Q = v e β μ λ 3 ,
Dann

ln Q β = 3 2 N β μ N .

Schritt 4. Kombinieren Sie die Schritte 2 und 3, um das gewünschte Ergebnis zu erzielen.

Entschuldigung, aber ich denke nicht, dass ein Hinweis als Antwort gepostet werden sollte. Ich denke, Sie sollten mehr Inhalte hinzufügen.
@FGSUZ Fair genug. Ich habe einige Erläuterungen beigefügt.
Aussagen zu partiellen Ableitungen, ohne explizit anzugeben, welche Variablen festgehalten werden, lassen in diesem Zusammenhang alles auf ihr Gegenteil schließen.

Ja, es gibt ein subtiles Problem. Erinnern Sie sich an das allgemeine Problem, aus dem diese Formeln abgeleitet werden.

Ein Beispiel zuerst. Das sagt die Thermodynamik

U = Φ + T S + μ N

( Φ ist das großkanonische Potential)

von denen Sie erhalten könnten

S = ( Φ T ) v , μ

Und überprüfen Sie, dass wir in der Thermodynamik immer einen Dubindex in den Ableitungen angeben. Dieser Subindex ( v , μ ) , erinnert uns daran, dass diese Variablen konstant gehalten werden müssen.

Wenn Sie darüber nachdenken, sollten sie nicht notwendig sein, da eine partielle Ableitung bereits implizit davon ausgeht, dass "alle anderen Variablen" konstant gehalten werden. Aber die Thermodynamik kann ziemlich chaotisch werden, daher ist es wirklich wichtig, die konstant gehaltenen Variablen im Auge zu behalten.

Nun, wenn Sie das verstanden haben, kommen wir zum Punkt. Bei Ensembles verwenden wir die entropische Darstellung. Damit arbeitet man S als grundlegende Beziehung, statt U . Folglich müssen Sie den Ausdruck neu anordnen:

S = Φ T + U T + μ N ¯ T

Und jetzt werden wie immer Variablen erhalten:

U = Φ ( λ U )

Wo λ U ist die intensive Variable, die mit verbunden ist U , welches ist ( 1 / T ) . Also im Grunde muss man rechnen:

U = Φ ( 1 / T )

Aber diese partielle Ableitung muss durchgeführt werden, indem die anderen Variablen konstant gehalten werden . Dazu gehört die Konstanthaltung ( μ / T )

Ja, es ist seltsam: 1 / T variiert, aber μ / T darf nicht variieren. Das ist es.

... und da Φ = k B T ln ( Q )

U = [ k B ln ( Q ) ] ( 1 / T )

Und alles durch dividieren k B , du erhältst

U = [ ln ( Q ) ] [ 1 / ( k B T ) ] ) = ln ( Q ) β

Aber diese Ableitung wird beibehalten ( μ / T ) konstant. Wenn das konstant ist, dividieren durch k B wird konstant bleiben, und deshalb schreiben wir

U = ( ln ( Q ) β ) ( β μ )

Das heißt, die Ableitung muss beibehalten werden β μ Konstante. Das ist die Erklärung deines Kalküls.

Es gibt interessante Artikel über die mangelnde Strenge dieser Notation, haha.

Der Ausdruck, für den Sie sich entschieden haben E ist nicht mit einem temperaturunabhängigen chemischen Potential vereinbar!

Um seine Abhängigkeit zu finden, erinnern Sie sich, dass die Konstanten wie β ergeben sich aus der Lösung der stationären Punkte des Entropiefunktionals:

S = D Γ ρ ( Q , P ) ln ( ρ ( Q , P ) )
unterliegen verschiedenen Bedingungen, die zu Lagrange-Multiplikatoren führen, die später mit verschiedenen thermodynamischen Eigenschaften identifiziert werden:

  • λ = ln ( Z G ) 1 aus Normalisierungsbeschränkung;

  • β := T 1 von konstant E ;

  • γ := μ β von konstant N ;

Wo

Z G = N = 0 D Γ N e H β γ N

Es ist jetzt klar, dass die Flüchtigkeit e μ β = e γ ist wirklich unabhängig von β .

Nun ist der Erwartungswert der Energie gegeben durch

( ln Z ) β
mit dem Verständnis, dass γ ist konstant. Dies ist dasselbe wie die Formel für das kanonische Ensemble.

Die übliche Behandlung ist jedoch zu behandeln γ als abhängige Variable und μ als unabhängig, da μ hat eine offensichtlichere physikalische Interpretation.

Daher ist die durchschnittliche Energie

E   =   1 Z G N = 0 D Γ N H e H β + β μ N   =   1 Z G N = 0 D Γ N ( β + μ N ) e H β + β μ   =   μ N β ( ln Z G )
wo die partielle Ableitung gilt μ Konstante.

Also je nachdem, ob Sie sich für das Halten entscheiden γ oder chemisches Potential als Ihre unabhängige Konstante erhalten Sie eine andere Formel für E . Das Endergebnis Ihrer Berechnung sollte unabhängig von dieser Wahl sein, solange Sie mit allen Ihren Ableitungen konsistent sind.

Dies erklärt die "Magie" beim Ignorieren der e β μ Faktor. Das Konstanthalten des Arguments stimmt mit Ihrem Ausdruck for überein E . Der "korrekte" (konventionelle) großkanonische Ensemble-Ausdruck enthält einen Begriff, der den "zusätzlichen (unsinnigen) Begriff" genau aufhebt.

Kommentar zur FGSUZ-Antwort, das etwas seltsame Ergebnis, dass

E = ( Protokoll Q β ) β μ
ist kein Zufall, sondern hat eine ziemlich tiefe Interpretation in Bezug auf die Wahrscheinlichkeitstheorie.

Betrachten Sie die großkanonische Verteilung

w k = 1 Q e β ( μ N k E k ) Wo  Q = k e β ( μ N k E k )
gibt die Wahrscheinlichkeit des Zustands an k mit Energie E k Und N k Teilchen bei der gegebenen β , μ .

In der Wahrscheinlichkeitstheorie bei gegebener Zufallsvariable E , definiert man seine momenterzeugende Funktion als Erwartungswert von e S E , M ( S ) = E [ e S E ] . Indem wir die Exponentialfunktion erweitern, sehen wir warum:

M ( S ) = E [ e S E ] = E [ N = 0 + ( S E ) N N ! ] = N = 0 + E [ E N ] S N N !
Also die n-te Ableitung von M bewertet bei S = 0 gibt den n-ten Moment M N = E [ E N ] von E .

Bei der großkanonischen Verteilung M kann als Zustandssummenfunktion ausgedrückt werden Q ( β , μ ) allein¹:

M ( S ) = k w k e S E k = 1 Q k e β μ N k ( β S ) E k = 1 Q k e ( β S ) / ( β S ) β μ N k ( β S ) E k = Q ( β S , β μ / ( β S ) ) / Q ( β , μ )
Das Ableiten davon ist ziemlich ärgerlich, weil β Und μ im zweiten Argument gemischt erscheinen. Es ist vorteilhaft, Variablen zu ändern, um sie zu entkoppeln: define λ = e β μ , die sogenannte Fugazität ². Die Partitionsfunktion wird dann
Q ( β , λ ) = k e β μ N k e β E k = k λ N k e β E k
Und M ( S ) = Q ( β S , λ ) / Q ( β , λ ) , was viel einfacher zu handhaben ist.

Wir können jetzt leicht das erste Moment, also die mittlere Energie berechnen:

E = M ' ( 0 ) = S [ Q ( β S , λ ) Q ( β , λ ) ] S = 0 = 1 Q Q β ( β , λ ) = ( Protokoll Q β ) λ
Zum Vergleich, nach langwieriger Berechnung Q ( β , μ ) hätte die Formel gegeben:
E = ( Protokoll Q β ) μ + μ β ( Protokoll Q μ ) β
was dem oben genannten vollständig entspricht, aber hässlicher und weniger nützlich ist.


Bonus : Höhere Momente können nicht so einfach als Ableitungen bzgl. ausgedrückt werden β , aber wir können stattdessen die Kumulanten berechnen . Diese k N sind Kombinationen der üblichen Momente M N und sind gleichermaßen gut darin, eine Wahrscheinlichkeitsverteilung zu beschreiben. Sie werden von der Funktion generiert Protokoll M ( S ) , wie vorher:

K ( S ) = Protokoll M ( S ) = N = 1 + k N S N N ! k N = K ( N ) ( S )

Diese Wahl ist sehr bequem, weil wir (bis auf die belanglose Konstante Protokoll Q ( β , λ ) )

K ( S ) = Protokoll Q ( β S , λ )
Daraus ergeben sich alle Kumulanten (die ersten drei sind der Durchschnitt E , Schwankungen Δ E 2 Und Δ E 3 ) einfach folgen:
k N = ( N Protokoll Q β N ) λ


¹ Dies ist nur eine Folge von w k exponentiell sein.

² Das Exponential ist nicht wichtig, aber es ist konventionell so definiert und ein bisschen netter.

Sie können eine Ableitung nehmen, indem Sie eine konstante Flüchtigkeit annehmen, dh e β μ . Denn die Hauptformel zur Berechnung der Zustandssumme im großkanonischen Ensemble ist Q = e β ( E μ N ) . Wenn Sie also die mittlere Energie berechnen möchten oder E = P ich E ich = e β ( E ich μ N ich ) Q E ich . Wir wissen das P ich = e β ( E ich μ N ich ) Q . Sie sollten also eine Ableitung von nehmen ln ( Q ) in einer Weise, dass für jeden Begriff nur E ich komm daneben P ich es bedeutet, dass Sie Ableitungen in konstanter Flüchtigkeit oder Konstante nehmen sollten e β μ