تمرین حل‌شدهٔ ترمودینامیک

غشای نیمه‌تراوا: کدام کمیت‌ها برابر می‌شوند؟

تمرین 25 · درس 6 — قانون دوم ترمودینامیک

  • تعادل، غشای نیمه‌تراوا، پتانسیل شیمیایی، فشار جزئی، دیوارهٔ ثابت، دیوارهٔ متحرک

صورت مسئله

ظرفی صلب و منزوی با حجم 2V02V_0 شامل دو مخلوط از گازهای ایده‌آل تک‌اتمی A وB است. غشایی گرماگذر آن‌ها را جدا می‌کند. A از آن می‌گذرد ولی B نمی‌گذرد. در آغاز غشا ثابت و نفوذناپذیر است و ظرف را به حجم‌های برابر V0V_0 در دمای یکسان T0T_0 تقسیم می‌کند. سمت چپ 2N02N_0 ذرهٔ A و N0N_0 ذرهٔ B و سمت راست N0N_0 ذرهٔ A و 2N02N_0 ذرهٔ B دارد. سپس تراوایی انتخابی فعال می‌شود. انرژی غشا و برهم‌کنش ذرات را نادیده بگیرید.

برای این تمرین، تعمیم همانی ترمودینامیکی به تعداد ذرات متغیر را بپذیرید:

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 پتانسیل شیمیایی هر ذره از گونهٔ ss است. برای هر مخلوط ایده‌آل نیز بپذیرید

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}.

پتانسیل شیمیایی گونهٔ s∈{A,B}s\in\{A,B\} برابر است با

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

که p∘p^\circ فشار مرجع است.

  1. ظرف و اجزای آن را رسم کنید. فشارهای کل آغازین و سپس فشارهای جزئی A را مقایسه کنید. با پتانسیل‌های شیمیایی جهت آغازین انتقال A را پیش‌بینی کنید.
  2. غشا ثابت می‌ماند. dStot\mathrm dS_{\mathrm{tot}} را برحسب بازتوزیع‌های مجاز dU1\mathrm dU_1 و dNA,1\mathrm dN_{A,1} بنویسید و شرط‌های ایستایی را نتیجه بگیرید. آیا فشارهای کل یا μB\mu_B باید برابر باشند؟
  3. دمای نهایی، تعداد ذرات A و فشارهای کل را بیابید. در دمای T0T_0 بررسی کنید که توزیع به‌دست‌آمدهٔ A آنتروپی را نسبت به انتقال A بیشینه می‌کند.
  4. از این حالت غشا را آزاد می‌کنیم. بدون اصطکاک یا نیروی اضافی حرکت می‌کند و گرماگذر و انتخابی می‌ماند. چه شرطی افزوده می‌شود؟ تعادل تازه را تعیین کنید: حجم‌ها، تعداد ذرات A، دما و فشار کل مشترک.
  5. کاربرد عددی: N0=NAn0N_0=\mathcal N_A n_0 با n0=1,0 moln_0=1{,}0\ \mathrm{mol}، V0=10 LV_0=10\ \mathrm L و T0=300 KT_0=300\ \mathrm K. رابطهٔ R=NAkBR=\mathcal N_A k_{\mathrm{B}} را یادآوری می‌کنیم و به کار می‌بریم 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}}.

    برای غشای ثابت و سپس متحرک، فشارهای نهایی و مقدار A در دو طرف (برحسب مول) و برای غشای متحرک حجم‌های نهایی را محاسبه کنید.

راهنمایی

راهنمایی
در حجم ثابت، اگر B از غشا نگذرد dS=dU/T−μA dNA/T\mathrm dS=\mathrm dU/T-\mu_A\,\mathrm dN_A/T. انتقال‌های دریافتی یک طرف از طرف دیگر کم می‌شوند. وقتی غشا متحرک است، جملهٔ (P1/T1−P2/T2) dV1(P_1/T_1-P_2/T_2)\,\mathrm dV_1 را بیفزایید.

پاسخ تفصیلی

پاسخ

Figure
شکل 1

پرسش 1. شکل دو گونه را متمایز می‌کند و نشان می‌دهد تنها A می‌تواند بگذرد. می‌گذاریم P∗=N0kBT0/V0P_*=N_0k_BT_0/V_0. هر محفظه در آغاز 3N03N_0 ذره دارد:

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_*.

برای فشار جزئی A تنها ذرات A را می‌شماریم: 2N02N_0 در چپ و N0N_0 در راست. در آغاز P1,i=P2,i=3P∗P_{1,i}=P_{2,i}=3P_*، ولی pA,1,i=2P∗p_{A,1,i}=2P_* و pA,2,i=P∗p_{A,2,i}=P_*. پس μA,1−μA,2=kBT0ln⁡2>0\mu_{A,1}-\mu_{A,2}=k_BT_0\ln2>0: A ابتدا از چپ به راست حرکت می‌کند. برای dNA,1<0\mathrm dN_{A,1}<0 در دماهای برابر، تغییر آنتروپی کل (μA,2−μA,1) dNA,1/T0(\mu_{A,2}-\mu_{A,1})\,\mathrm dN_{A,1}/T_0 مثبت است. برابری فشارهای کل برای توقف پخش کافی نیست. در اختلاف پتانسیل‌های شیمیایی، دو جملهٔ μA∘(T0)\mu_A^\circ(T_0) حذف می‌شوند؛ ترکیب لگاریتم‌ها می‌دهد kBT0ln⁡(pA,1,i/pA,2,i)=kBT0ln⁡2k_BT_0\ln(p_{A,1,i}/p_{A,2,i})=k_BT_0\ln2.

پرسش 2. انرژی کل و تعداد کل ذرات A پایسته‌اند؛ تعداد ذرات B در هر طرف ثابت است. با dU2=−dU1\mathrm dU_2=-\mathrm dU_1، dNA,2=−dNA,1\mathrm dN_{A,2}=-\mathrm dN_{A,1} و dNB,1=dNB,2=0\mathrm dN_{B,1}=\mathrm dN_{B,2}=0 داریم

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).

جملهٔ حجم صفر می‌شود چون غشا ثابت است. جمع می‌دهد

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}.

با جمع‌کردن ضرایب داریم

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}.

رسانش گرما و انتقال A مجازند: تغییراتشان را می‌توان مستقل آزمود. پس ایستایی ایجاب می‌کند 1/T1−1/T2=01/T_1-1/T_2=0، یعنی T1=T2T_1=T_2، سپس (μA,2−μA,1)/T=0(\mu_{A,2}-\mu_{A,1})/T=0، یعنی μA,1=μA,2\mu_{A,1}=\mu_{A,2}. نه P1=P2P_1=P_2 لازم است (غشا ثابت است) و نه μB,1=μB,2\mu_{B,1}=\mu_{B,2} (B نمی‌تواند بگذرد).

پرسش 3. انرژی کل آغازین 9N0kBT09N_0k_BT_0 است و تعداد کل ذرات 6N06N_0 می‌ماند. در واقع Ui=(3/2)(3N0+3N0)kBT0=9N0kBT0U_i=(3/2)(3N_0+3N_0)k_BT_0=9N_0k_BT_0. در دمای مشترک نهایی TT، داریم Uf=(3/2)(6N0)kBT=9N0kBTU_f=(3/2)(6N_0)k_BT=9N_0k_BT. پس برابری Uf=UiU_f=U_i می‌دهد T=T0T=T_0. برابری پتانسیل‌های شیمیایی A تبدیل می‌شود به

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}.

با تقسیم بر kBT0k_BT_0 و گرفتن نمایی، pA,1=pA,2p_{A,1}=p_{A,2} به دست می‌آید. چون حجم‌ها برابرند، عوامل kBT0/V0k_BT_0/V_0 حذف می‌شوند: NA,1=NA,2N_{A,1}=N_{A,2}. چون NA,1+NA,2=3N0N_{A,1}+N_{A,2}=3N_0، هر یک 3N0/23N_0/2 است. با افزودن ذرات B داریم

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}.

بنابراین،

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

تکیه‌گاه غشا این اختلاف فشار را موازنه می‌کند. با a=NA,1a=N_{A,1} و NA,2=3N0−aN_{A,2}=3N_0-a در دما و حجم‌های ثابت، پرسش 2 می‌دهد

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}.

جملهٔ انرژی حذف می‌شود چون دماها برابرند و dU1+dU2=0\mathrm dU_1+\mathrm dU_2=0. برای مشتق‌گیری دوباره می‌نویسیم ln⁡(3N0−a)−ln⁡a\ln(3N_0-a)-\ln a: مشتقش −1/(3N0−a)−1/a-1/(3N_0-a)-1/a است. پس،

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.

توزیع برابر یگانه بیشینه روی 0<a<3N00<a<3N_0 است.

پرسش 4. ایستایی نسبت به حجم P1=P2P_1=P_2 را می‌افزاید: جملهٔ اضافی در dStot\mathrm dS_{\mathrm{tot}} برابر (P1/T1−P2/T2) dV1(P_1/T_1-P_2/T_2)\,\mathrm dV_1 است. ضریب آن باید صفر شود؛ با T1=T2T_1=T_2 این امر برابری فشارهای کل را ایجاب می‌کند. دمای نهایی دوباره T0T_0 است، با همان موازنهٔ انرژی پرسش 3. چون pA,1=pA,2p_{A,1}=p_{A,2}، برابری فشارهای کل نتیجه می‌دهد pB,1=pB,2p_{B,1}=p_{B,2}. B نمی‌تواند بگذرد:

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

پس V2=2V1V_2=2V_1 و سپس 3V1=2V03V_1=2V_0: V1=2V0/3V_1=2V_0/3 و V2=4V0/3V_2=4V_0/3. برای A برابری فشارهای جزئی می‌دهد NA,1/V1=NA,2/V2N_{A,1}/V_1=N_{A,2}/V_2، پس NA,2=2NA,1N_{A,2}=2N_{A,1}. چون مجموع 3N03N_0 است، داریم NA,1=N0N_{A,1}=N_0 و NA,2=2N0N_{A,2}=2N_0. سرانجام Pf=(N0+N0)kBT0/(2V0/3)=3P∗P_f=(N_0+N_0)k_BT_0/(2V_0/3)=3P_*. پس نتایج عبارت‌اند از

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_*.

برابری فشارهای جزئی B اینجا از شرط‌های دیگر و مدل مخلوط ایده‌آل نتیجه می‌شود، نه از تبادل مجاز B. مقدارهای A دیگر برابر نیستند: غلظت‌ها و بنابراین فشارهای جزئی برابر می‌شوند.

پرسش 5. برای محاسبهٔ فشارها حجم‌ها را به یکاهای SI تبدیل می‌کنیم: V0=10 L=10−2 m3V_0=10\ \mathrm L=10^{-2}\ \mathrm{m^3}. ثابت NA\mathcal N_A مقدار ماده را به تعداد ذرات تبدیل می‌کند: N0=NAn0N_0=\mathcal N_A n_0 با n0=1,0 moln_0=1{,}0\ \mathrm{mol}. ابتدا محاسبه می‌کنیم

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}}.

با N0kB=n0RN_0k_B=n_0R، مقیاس فشار تعریف‌شده در پرسش 1 برابر است با

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},

زیرا 1 bar=105 Pa1\ \mathrm{bar}=10^5\ \mathrm{Pa}. این فشار یک مقیاس محاسباتی است؛ خودش فشار کل هیچ محفظه‌ای نیست.

غشای ثابت. هر حجم 10 L10\ \mathrm L می‌ماند. پرسش 3 می‌دهد NA,1=NA,2=3N0/2N_{A,1}=N_{A,2}=3N_0/2. برای تبدیل ذرات به مول بر NA\mathcal N_A تقسیم می‌کنیم:

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 نمی‌تواند از غشا بگذرد: 1 mol1\ \mathrm{mol} از B در چپ و 2 mol2\ \mathrm{mol} در راست می‌ماند. پس مقدارهای کل 2,5 mol2{,}5\ \mathrm{mol} و 3,5 mol3{,}5\ \mathrm{mol} هستند. در T0=300 KT_0=300\ \mathrm K، معادلهٔ P=nRT/VP=nRT/V جداگانه می‌دهد

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}.

فشارهای کل متفاوت‌اند؛ ثابت‌کردن غشا این حالت را نگه می‌دارد. فشارهای جزئی A برابرند، چون A با مقدارها و دماهای برابر حجم‌های یکسان را اشغال می‌کند.

غشای متحرک. از پرسش 4،

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.

مجموع همچنان 20 L20\ \mathrm L، حجم ظرف، است. جمعیت‌های A اکنون NA,1=N0N_{A,1}=N_0 و 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}.

مجموع واقعاً 3 mol3\ \mathrm{mol} از A است. با افزودن B، دو محفظه به‌ترتیب 22 و 4 mol4\ \mathrm{mol} گاز دارند. پیش از گردکردن، حجم‌های دقیق را به کار می‌بریم:

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}.

هر دو محاسبه فشار کل یکسان می‌دهند، مطابق الزام تعادل مکانیکی غشای آزادشده.