Gelöste Thermodynamikaufgabe

Semipermeable Membran: Welche Größen gleichen sich an?

Übung 25 · Lektion 6 — Der zweite Hauptsatz der Thermodynamik

  • Gleichgewicht
  • semipermeable Membran
  • chemisches Potential
  • Partialdruck
  • feste Trennwand
  • bewegliche Trennwand

Aufgabenstellung

Ein starrer, isolierter Behälter mit dem Volumen 2V02V_0 enthält zwei Mischungen einatomiger idealer Gase A und B. Sie sind durch eine diathermische Membran getrennt, die A durchlässt, B jedoch nicht. Anfangs ist sie fixiert und undurchlässig gemacht und teilt den Behälter in gleiche Volumina V0V_0 bei derselben Temperatur T0T_0. Links befinden sich 2N02N_0 Teilchen von A und N0N_0 von B, rechts N0N_0 von A und 2N02N_0 von B. Dann wird die selektive Durchlässigkeit aktiviert. Membranenergie und Wechselwirkungen zwischen den Teilchen werden vernachlässigt.

Für diese Aufgabe wird die Erweiterung der thermodynamischen Identität auf veränderliche Teilchenzahlen als gegeben angenommen:

dU=T dS−P dV+μA dNA+μB dNB.\boxed{\mathrm dU=T\,\mathrm dS-P\,\mathrm dV +\mu_A\,\mathrm dN_A+\mu_B\,\mathrm dN_B.}

μs\mu_s ist das chemische Potential pro Teilchen der Spezies ss. Für jede ideale Mischung gelte außerdem

U=32(NA+NB)kBT,P=pA+pB,ps=NskBTV.U=\frac32(N_A+N_B)k_BT,\quad P=p_A+p_B,\quad p_s=\frac{N_sk_BT}{V}.

Das chemische Potential der Spezies s∈{A,B}s\in\{A,B\} ist

μs(T,ps)=μs∘(T)+kBTln⁡psp∘,\mu_s(T,p_s)=\mu_s^\circ(T)+k_BT\ln\frac{p_s}{p^\circ},

wobei p∘p^\circ ein Referenzdruck ist.

  1. Skizzieren Sie den Behälter und seine Bestandteile. Vergleichen Sie zuerst die anfänglichen Gesamtdrücke, dann die Partialdrücke von A. Bestimmen Sie aus den chemischen Potentialen die anfängliche Transportrichtung von A.
  2. Die Membran bleibt fest. Drücken Sie dStot\mathrm dS_{\mathrm{tot}} durch die erlaubten Umverteilungen dU1\mathrm dU_1 und dNA,1\mathrm dN_{A,1} aus und bestimmen Sie die Stationaritätsbedingungen. Müssen die Gesamtdrücke oder die μB\mu_B gleich sein?
  3. Bestimmen Sie Endtemperatur, A-Teilchenzahlen und Gesamtdrücke. Überprüfen Sie bei der Temperatur T0T_0, dass die gefundene A-Verteilung die Entropie bezüglich des A-Transports maximiert.
  4. Ausgehend von diesem Zustand wird die Membran freigegeben. Sie bewegt sich reibungsfrei und ohne zusätzliche Kraft und bleibt diathermisch und selektiv durchlässig. Welche zusätzliche Bedingung gilt? Bestimmen Sie das neue Gleichgewicht: Volumina, A-Teilchenzahlen, Temperatur und gemeinsamen Gesamtdruck.
  5. Zahlenbeispiel: N0=NAn0N_0=\mathcal N_A n_0, mit n0=1,0 moln_0=1{,}0\ \mathrm{mol}, V0=10 LV_0=10\ \mathrm L und T0=300 KT_0=300\ \mathrm K. Verwenden Sie die Beziehung R=NAkBR=\mathcal N_A k_{\mathrm{B}} sowie NA=6,022×1023 mol−1,kB=1,380649×10−23 J K−1.\mathcal N_A=6{,}022\times10^{23}\ \mathrm{mol^{-1}},\qquad k_{\mathrm{B}}=1{,}380649\times10^{-23}\ \mathrm{J\,K^{-1}}.

    Berechnen Sie für die feste und anschließend für die bewegliche Membran die Enddrücke und die Stoffmengen von A auf beiden Seiten (in Mol) sowie die Endvolumina bei beweglicher Membran.

Hinweis

Hinweis
Bei festem Volumen gilt dS=dU/T−μA dNA/T\mathrm dS=\mathrm dU/T-\mu_A\,\mathrm dN_A/T, wenn B die Membran nicht durchquert. Die von einer Seite aufgenommenen Größen gehen der anderen verloren. Bei beweglicher Membran kommt der Term (P1/T1−P2/T2) dV1(P_1/T_1-P_2/T_2)\,\mathrm dV_1 hinzu.

Ausführliche Lösung

Lösung

Figure
Abbildung 1

Frage 1. Die Skizze unterscheidet beide Spezies und zeigt, dass nur A hindurchtreten kann. Wir setzen P∗=N0kBT0/V0P_*=N_0k_BT_0/V_0. Jeder Teilraum enthält anfangs 3N03N_0 Teilchen:

P1,i=(2N0+N0)kBT0V0=3P∗,P2,i=(N0+2N0)kBT0V0=3P∗.P_{1,i}=\frac{(2N_0+N_0)k_BT_0}{V_0}=3P_*,\qquad P_{2,i}=\frac{(N_0+2N_0)k_BT_0}{V_0}=3P_*.

Für den Partialdruck von A werden nur A-Teilchen gezählt: 2N02N_0 links und N0N_0 rechts. Anfangs gilt P1,i=P2,i=3P∗P_{1,i}=P_{2,i}=3P_*, aber pA,1,i=2P∗p_{A,1,i}=2P_* und pA,2,i=P∗p_{A,2,i}=P_*. Daher μA,1−μA,2=kBT0ln⁡2>0\mu_{A,1}-\mu_{A,2}=k_BT_0\ln2>0: A bewegt sich zunächst von links nach rechts. Für dNA,1<0\mathrm dN_{A,1}<0 bei gleichen Temperaturen ist die gesamte Entropieänderung (μA,2−μA,1) dNA,1/T0(\mu_{A,2}-\mu_{A,1})\,\mathrm dN_{A,1}/T_0 positiv. Gleiche Gesamtdrücke reichen nicht aus, um Diffusion zu unterbinden. In der Differenz der chemischen Potentiale heben sich beide Terme μA∘(T0)\mu_A^\circ(T_0) auf; Zusammenfassen der Logarithmen ergibt kBT0ln⁡(pA,1,i/pA,2,i)=kBT0ln⁡2k_BT_0\ln(p_{A,1,i}/p_{A,2,i})=k_BT_0\ln2.

Frage 2. Gesamtenergie und gesamte A-Teilchenzahl bleiben erhalten; die B-Teilchenzahlen auf beiden Seiten sind konstant. Mit dU2=−dU1\mathrm dU_2=-\mathrm dU_1, dNA,2=−dNA,1\mathrm dN_{A,2}=-\mathrm dN_{A,1} und dNB,1=dNB,2=0\mathrm dN_{B,1}=\mathrm dN_{B,2}=0 folgt

dSj=dUjTj−μA,jTj dNA,j(j=1,2).\mathrm dS_j=\frac{\mathrm dU_j}{T_j} -\frac{\mu_{A,j}}{T_j}\,\mathrm dN_{A,j} \quad (j=1,2).

Der Volumenterm verschwindet, weil die Membran fest ist. Die Addition ergibt

dStot=dU1T1−dU1T2−μA,1T1dNA,1+μA,2T2dNA,1.\mathrm dS_{\mathrm{tot}} =\frac{\mathrm dU_1}{T_1}-\frac{\mathrm dU_1}{T_2} -\frac{\mu_{A,1}}{T_1}\mathrm dN_{A,1} +\frac{\mu_{A,2}}{T_2}\mathrm dN_{A,1}.

Zusammenfassen der Koeffizienten liefert

dStot=(1T1−1T2)dU1+(μA,2T2−μA,1T1)dNA,1.\mathrm dS_{\mathrm{tot}}= \left(\frac1{T_1}-\frac1{T_2}\right)\mathrm dU_1 +\left(\frac{\mu_{A,2}}{T_2}-\frac{\mu_{A,1}}{T_1}\right) \mathrm dN_{A,1}.

Wärmeleitung und A-Transport sind erlaubt: Ihre Variationen können unabhängig geprüft werden. Stationarität verlangt daher 1/T1−1/T2=01/T_1-1/T_2=0, also T1=T2T_1=T_2, und anschließend (μA,2−μA,1)/T=0(\mu_{A,2}-\mu_{A,1})/T=0, also μA,1=μA,2\mu_{A,1}=\mu_{A,2}. Sie verlangt weder P1=P2P_1=P_2 (die Membran ist fest) noch μB,1=μB,2\mu_{B,1}=\mu_{B,2} (B kann nicht hindurchtreten).

Frage 3. Anfangs beträgt die Gesamtenergie 9N0kBT09N_0k_BT_0, und die gesamte Teilchenzahl bleibt 6N06N_0. Denn Ui=(3/2)(3N0+3N0)kBT0=9N0kBT0U_i=(3/2)(3N_0+3N_0)k_BT_0=9N_0k_BT_0. Bei gemeinsamer Endtemperatur TT gilt Uf=(3/2)(6N0)kBT=9N0kBTU_f=(3/2)(6N_0)k_BT=9N_0k_BT. Die Gleichheit Uf=UiU_f=U_i liefert daher T=T0T=T_0. Die Gleichheit der chemischen Potentiale von A wird zu

kBT0ln⁡pA,1p∘=kBT0ln⁡pA,2p∘.k_BT_0\ln\frac{p_{A,1}}{p^\circ} =k_BT_0\ln\frac{p_{A,2}}{p^\circ}.

Division durch kBT0k_BT_0 und Exponentieren ergibt pA,1=pA,2p_{A,1}=p_{A,2}. Da die Volumina gleich sind, heben sich die Faktoren kBT0/V0k_BT_0/V_0 auf: NA,1=NA,2N_{A,1}=N_{A,2}. Wegen NA,1+NA,2=3N0N_{A,1}+N_{A,2}=3N_0 beträgt jede Zahl 3N0/23N_0/2. Nach Addition der B-Teilchen folgt

P1,f=(3N0/2+N0)kBT0V0,P2,f=(3N0/2+2N0)kBT0V0.P_{1,f}=\frac{(3N_0/2+N_0)k_BT_0}{V_0},\qquad P_{2,f}=\frac{(3N_0/2+2N_0)k_BT_0}{V_0}.

Somit

P1,f=52P∗,P2,f=72P∗.P_{1,f}=\frac52P_*,\qquad P_{2,f}=\frac72P_*.

Die Membranhalterung gleicht diesen Druckunterschied aus. Mit a=NA,1a=N_{A,1} und NA,2=3N0−aN_{A,2}=3N_0-a liefert Frage 2 bei fester Temperatur und festen Volumina

dStotda=μA,2−μA,1T0=kBln⁡pA,2pA,1=kBln⁡3N0−aa.\frac{\mathrm dS_{\mathrm{tot}}}{\mathrm da} =\frac{\mu_{A,2}-\mu_{A,1}}{T_0} =k_{\mathrm{B}}\ln\frac{p_{A,2}}{p_{A,1}} =k_{\mathrm{B}}\ln\frac{3N_0-a}{a}.

Der Energieterm verschwindet wegen der gleichen Temperaturen und dU1+dU2=0\mathrm dU_1+\mathrm dU_2=0. Zum erneuten Differenzieren schreiben wir ln⁡(3N0−a)−ln⁡a\ln(3N_0-a)-\ln a; die Ableitung lautet −1/(3N0−a)−1/a-1/(3N_0-a)-1/a. Somit

dStotda=kBln⁡3N0−aa,d2Stotda2=−kB(1a+13N0−a)<0.\frac{\mathrm dS_{\mathrm{tot}}}{\mathrm da} =k_{\mathrm{B}}\ln\frac{3N_0-a}{a},\qquad \frac{\mathrm d^2S_{\mathrm{tot}}}{\mathrm da^2} =-k_{\mathrm{B}}\left(\frac1a+\frac1{3N_0-a}\right)<0.

Die Gleichverteilung ist das einzige Maximum auf 0<a<3N00<a<3N_0.

Frage 4. Die Stationarität bezüglich des Volumens ergänzt P1=P2P_1=P_2: Der zusätzliche Term in dStot\mathrm dS_{\mathrm{tot}} ist (P1/T1−P2/T2) dV1(P_1/T_1-P_2/T_2)\,\mathrm dV_1. Sein Koeffizient muss verschwinden; mit T1=T2T_1=T_2 erzwingt dies gleiche Gesamtdrücke. Die Endtemperatur ist nach derselben Energiebilanz wie in Frage 3 wieder T0T_0. Wegen pA,1=pA,2p_{A,1}=p_{A,2} implizieren gleiche Gesamtdrücke pB,1=pB,2p_{B,1}=p_{B,2}. B kann nicht hindurchtreten:

N0V1=2N0V2,V1+V2=2V0,\frac{N_0}{V_1}=\frac{2N_0}{V_2},\qquad V_1+V_2=2V_0,

also V2=2V1V_2=2V_1 und dann 3V1=2V03V_1=2V_0: V1=2V0/3V_1=2V_0/3 und V2=4V0/3V_2=4V_0/3. Für A liefern gleiche Partialdrücke NA,1/V1=NA,2/V2N_{A,1}/V_1=N_{A,2}/V_2, somit NA,2=2NA,1N_{A,2}=2N_{A,1}. Da die Summe 3N03N_0 ist, folgen NA,1=N0N_{A,1}=N_0 und NA,2=2N0N_{A,2}=2N_0. Schließlich Pf=(N0+N0)kBT0/(2V0/3)=3P∗P_f=(N_0+N_0)k_BT_0/(2V_0/3)=3P_*. Die Ergebnisse lauten somit

V1=23V0,V2=43V0,NA,1=N0,NA,2=2N0,Pf=3P∗.V_1=\frac23V_0,\quad V_2=\frac43V_0,\quad N_{A,1}=N_0,\quad N_{A,2}=2N_0,\quad P_f=3P_*.

Die Gleichheit der B-Partialdrücke folgt hier aus den übrigen Bedingungen und dem Modell der idealen Mischung, nicht aus einem erlaubten B-Austausch. Die A-Stoffmengen sind nicht mehr gleich: Ihre Konzentrationen und damit ihre Partialdrücke gleichen sich an.

Frage 5. Zur Druckberechnung werden die Volumina in SI-Einheiten umgerechnet: V0=10 L=10−2 m3V_0=10\ \mathrm L=10^{-2}\ \mathrm{m^3}. Die Konstante NA\mathcal N_A überführt eine Stoffmenge in eine Teilchenzahl: N0=NAn0N_0=\mathcal N_A n_0 mit n0=1,0 moln_0=1{,}0\ \mathrm{mol}. Zunächst berechnen wir

R=NAkB=(6,022×1023)(1,380649×10−23)≃8,314 J mol−1 K−1.R=\mathcal N_A k_{\mathrm{B}} =(6{,}022\times10^{23})(1{,}380649\times10^{-23}) \simeq8{,}314\ \mathrm{J\,mol^{-1}\,K^{-1}}.

Mit N0kB=n0RN_0k_B=n_0R beträgt die in Frage 1 eingeführte Druckskala

P∗=N0kBT0V0=n0RT0V0≃1,0×8,314×30010−2=249420 Pa≃2,494 bar,P_* =\frac{N_0k_BT_0}{V_0} =\frac{n_0RT_0}{V_0} \simeq\frac{1{,}0\times8{,}314\times300}{10^{-2}} =249420\ \mathrm{Pa}\simeq2{,}494\ \mathrm{bar},

denn 1 bar=105 Pa1\ \mathrm{bar}=10^5\ \mathrm{Pa}. Dieser Druck dient als Rechenskala; er ist selbst kein Gesamtdruck eines Teilraums.

Feste Membran. Jedes Volumen bleibt 10 L10\ \mathrm L. Frage 3 liefert NA,1=NA,2=3N0/2N_{A,1}=N_{A,2}=3N_0/2. Zur Umrechnung der Teilchenzahl in Mol wird durch NA\mathcal N_A dividiert:

nA,1=nA,2=3N02NA=32×1 mol=1,5 mol.n_{A,1}=n_{A,2}=\frac{3N_0}{2\mathcal N_A} =\frac32\times1\ \mathrm{mol}=1{,}5\ \mathrm{mol}.

B durchquert die Membran nicht: Links bleiben 1 mol1\ \mathrm{mol} B, rechts 2 mol2\ \mathrm{mol}. Die gesamten Stoffmengen betragen daher 2,5 mol2{,}5\ \mathrm{mol} und 3,5 mol3{,}5\ \mathrm{mol}. Bei T0=300 KT_0=300\ \mathrm K liefert die Gleichung P=nRT/VP=nRT/V getrennt

P1,f=(1,5+1)×8,314×30010−2=623550 Pa≃6,24 bar,P2,f=(1,5+2)×8,314×30010−2=872970 Pa≃8,73 bar.P_{1,f}=\frac{(1{,}5+1)\times8{,}314\times300}{10^{-2}} =623550\ \mathrm{Pa}\simeq6{,}24\ \mathrm{bar}, P_{2,f}=\frac{(1{,}5+2)\times8{,}314\times300}{10^{-2}} =872970\ \mathrm{Pa}\simeq8{,}73\ \mathrm{bar}.

Die Gesamtdrücke unterscheiden sich; die Fixierung der Membran erhält diesen Zustand. Die A-Partialdrücke sind gleich, da A gleiche Volumina mit gleichen Stoffmengen und Temperaturen einnimmt.

Bewegliche Membran. Aus Frage 4 folgt

V1,f=23V0=23×10≃6,67 L,V2,f=43V0=43×10≃13,33 L.V_{1,f}=\frac23V_0=\frac23\times10\simeq6{,}67\ \mathrm L, \qquad V_{2,f}=\frac43V_0=\frac43\times10\simeq13{,}33\ \mathrm L.

Die Summe bleibt 20 L20\ \mathrm L, das Behältervolumen. Die A-Teilchenzahlen betragen nun NA,1=N0N_{A,1}=N_0 und NA,2=2N0N_{A,2}=2N_0:

nA,1=N0NA=1 mol,nA,2=2N0NA=2 mol.n_{A,1}=\frac{N_0}{\mathcal N_A}=1\ \mathrm{mol},\qquad n_{A,2}=\frac{2N_0}{\mathcal N_A}=2\ \mathrm{mol}.

Die Summe beträgt tatsächlich 3 mol3\ \mathrm{mol} A. Einschließlich B enthalten die Teilräume 22 beziehungsweise 4 mol4\ \mathrm{mol} Gas. Wir verwenden vor dem Runden die exakten Volumina:

P1,f=(1+1)×8,314×300(2/3)×10−2=748260 Pa≃7,48 bar,P2,f=(2+2)×8,314×300(4/3)×10−2=748260 Pa≃7,48 bar.P_{1,f}=\frac{(1+1)\times8{,}314\times300}{(2/3)\times10^{-2}} =748260\ \mathrm{Pa}\simeq7{,}48\ \mathrm{bar}, P_{2,f}=\frac{(2+2)\times8{,}314\times300}{(4/3)\times10^{-2}} =748260\ \mathrm{Pa}\simeq7{,}48\ \mathrm{bar}.

Beide Rechnungen ergeben denselben Gesamtdruck, wie es das mechanische Gleichgewicht der freigegebenen Membran verlangt.