Solved thermodynamics exercise

Semipermeable membrane: which quantities become equal?

Exercise 25 · Lesson 6 — The Second Law of Thermodynamics

  • equilibrium
  • semipermeable membrane
  • chemical potential
  • partial pressure
  • fixed partition
  • movable partition

Statement

A rigid, isolated vessel of volume 2V02V_0 contains two mixtures of monatomic ideal gases A and B. A diathermal membrane separates them. It allows A through, but not B. Initially held fixed and made impermeable, it divides the vessel into equal volumes V0V_0 at the same temperature T0T_0. On the left are 2N02N_0 particles of A and N0N_0 of B; on the right, N0N_0 of A and 2N02N_0 of B. Its selective permeability is then activated. Neglect the membrane's energy and interactions between particles.

For this exercise, accept the extension of the thermodynamic identity to variable particle numbers:

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 is the chemical potential per particle of species ss. Also accept, for each ideal mixture,

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

The chemical potential of species s∈{A,B}s\in\{A,B\} is

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

where p∘p^\circ is a reference pressure.

  1. Sketch the vessel and its constituents. Compare the initial total pressures, then the partial pressures of A. Use the chemical potentials to predict the initial direction of A transfer.
  2. The membrane remains fixed. Express dStot\mathrm dS_{\mathrm{tot}} in terms of the allowed redistributions dU1\mathrm dU_1 and dNA,1\mathrm dN_{A,1}, then deduce the stationarity conditions. Must the total pressures or the μB\mu_B be equal?
  3. Determine the final temperature, the numbers of A particles and the total pressures. At temperature T0T_0, verify that the distribution of A found maximizes entropy with respect to A transfer.
  4. Starting from this state, release the membrane. It moves without friction or any additional force and remains diathermal and selective. What additional condition applies? Determine the new equilibrium: volumes, numbers of A particles, temperature and common total pressure.
  5. Numerical application: N0=NAn0N_0=\mathcal N_A n_0, with n0=1,0 moln_0=1{,}0\ \mathrm{mol}, V0=10 LV_0=10\ \mathrm L and T0=300 KT_0=300\ \mathrm K. Recall the relation R=NAkBR=\mathcal N_A k_{\mathrm{B}} and use 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}}.

    For the fixed and then the movable membrane, calculate the final pressures and the amounts of A on each side (in moles), together with the final volumes when the membrane is movable.

Hint

Hint
At fixed volume, dS=dU/T−μA dNA/T\mathrm dS=\mathrm dU/T-\mu_A\,\mathrm dN_A/T if B does not cross the membrane. Transfers received by one side are lost by the other. When the membrane is movable, add the term (P1/T1−P2/T2) dV1(P_1/T_1-P_2/T_2)\,\mathrm dV_1.

Detailed solution

Solution

Figure
Figure 1

Question 1. The diagram distinguishes the two species and shows that only A can cross. Set P∗=N0kBT0/V0P_*=N_0k_BT_0/V_0. Each compartment initially contains 3N03N_0 particles:

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

For A's partial pressure, count only A particles: 2N02N_0 on the left and N0N_0 on the right. Initially P1,i=P2,i=3P∗P_{1,i}=P_{2,i}=3P_*, but pA,1,i=2P∗p_{A,1,i}=2P_* and pA,2,i=P∗p_{A,2,i}=P_*. Thus μA,1−μA,2=kBT0ln⁡2>0\mu_{A,1}-\mu_{A,2}=k_BT_0\ln2>0: A initially moves from left to right. For dNA,1<0\mathrm dN_{A,1}<0 at equal temperatures, the total entropy change (μA,2−μA,1) dNA,1/T0(\mu_{A,2}-\mu_{A,1})\,\mathrm dN_{A,1}/T_0 is positive. Equal total pressures do not suffice to stop diffusion. In the chemical-potential difference, the two μA∘(T0)\mu_A^\circ(T_0) terms cancel; combining the logarithms gives kBT0ln⁡(pA,1,i/pA,2,i)=kBT0ln⁡2k_BT_0\ln(p_{A,1,i}/p_{A,2,i})=k_BT_0\ln2.

Question 2. Total energy and the total number of A particles are conserved; the numbers of B particles on each side are constant. Using dU2=−dU1\mathrm dU_2=-\mathrm dU_1, dNA,2=−dNA,1\mathrm dN_{A,2}=-\mathrm dN_{A,1} and dNB,1=dNB,2=0\mathrm dN_{B,1}=\mathrm dN_{B,2}=0 gives

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

The volume term vanishes because the membrane is fixed. Adding gives

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

Collecting the coefficients gives

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

Heat conduction and A transfer are allowed: their variations can be tested independently. Stationarity therefore requires 1/T1−1/T2=01/T_1-1/T_2=0, or T1=T2T_1=T_2, followed by (μA,2−μA,1)/T=0(\mu_{A,2}-\mu_{A,1})/T=0, or μA,1=μA,2\mu_{A,1}=\mu_{A,2}. It requires neither P1=P2P_1=P_2 (the membrane is fixed) nor μB,1=μB,2\mu_{B,1}=\mu_{B,2} (B cannot cross).

Question 3. Initially the total energy is 9N0kBT09N_0k_BT_0 and the total number of particles remains 6N06N_0. Indeed, Ui=(3/2)(3N0+3N0)kBT0=9N0kBT0U_i=(3/2)(3N_0+3N_0)k_BT_0=9N_0k_BT_0. At a common final temperature TT, Uf=(3/2)(6N0)kBT=9N0kBTU_f=(3/2)(6N_0)k_BT=9N_0k_BT. Equality Uf=UiU_f=U_i therefore gives T=T0T=T_0. Equality of A's chemical potentials becomes

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

Dividing by kBT0k_BT_0 and exponentiating gives pA,1=pA,2p_{A,1}=p_{A,2}. Since the volumes are equal, the factors kBT0/V0k_BT_0/V_0 cancel: NA,1=NA,2N_{A,1}=N_{A,2}. Since NA,1+NA,2=3N0N_{A,1}+N_{A,2}=3N_0, each equals 3N0/23N_0/2. Adding the B particles gives

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

Thus,

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

The membrane's support balances this pressure difference. Setting a=NA,1a=N_{A,1} and NA,2=3N0−aN_{A,2}=3N_0-a, at fixed temperature and volumes, question 2 gives

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

The energy term cancels because the temperatures are equal and dU1+dU2=0\mathrm dU_1+\mathrm dU_2=0. To differentiate again, write ln⁡(3N0−a)−ln⁡a\ln(3N_0-a)-\ln a: its derivative is −1/(3N0−a)−1/a-1/(3N_0-a)-1/a. Thus,

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.

The equal distribution is the unique maximum on 0<a<3N00<a<3N_0.

Question 4. Stationarity with respect to volume adds P1=P2P_1=P_2: the additional term in dStot\mathrm dS_{\mathrm{tot}} is (P1/T1−P2/T2) dV1(P_1/T_1-P_2/T_2)\,\mathrm dV_1. Its coefficient must vanish; with T1=T2T_1=T_2, this requires equal total pressures. The final temperature is again T0T_0, from the same energy balance as in question 3. Since pA,1=pA,2p_{A,1}=p_{A,2}, equal total pressures imply pB,1=pB,2p_{B,1}=p_{B,2}. B cannot cross:

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

so V2=2V1V_2=2V_1, then 3V1=2V03V_1=2V_0: V1=2V0/3V_1=2V_0/3 and V2=4V0/3V_2=4V_0/3. For A, equal partial pressures give NA,1/V1=NA,2/V2N_{A,1}/V_1=N_{A,2}/V_2, hence NA,2=2NA,1N_{A,2}=2N_{A,1}. Since the sum is 3N03N_0, we obtain NA,1=N0N_{A,1}=N_0 and NA,2=2N0N_{A,2}=2N_0. Finally, Pf=(N0+N0)kBT0/(2V0/3)=3P∗P_f=(N_0+N_0)k_BT_0/(2V_0/3)=3P_*. The results are therefore

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

Equal B partial pressures follow here from the other conditions and the ideal-mixture model, rather than from an allowed exchange of B. The amounts of A are no longer equal: their concentrations, and hence their partial pressures, become equal.

Question 5. Convert the volumes to SI units to calculate pressures: V0=10 L=10−2 m3V_0=10\ \mathrm L=10^{-2}\ \mathrm{m^3}. The constant NA\mathcal N_A converts an amount of substance to a particle number: N0=NAn0N_0=\mathcal N_A n_0 with n0=1,0 moln_0=1{,}0\ \mathrm{mol}. First calculate

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

With N0kB=n0RN_0k_B=n_0R, the pressure scale introduced in question 1 is

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

because 1 bar=105 Pa1\ \mathrm{bar}=10^5\ \mathrm{Pa}. This pressure is a calculation scale; it is not itself a compartment's total pressure.

Fixed membrane. Each volume remains 10 L10\ \mathrm L. Question 3 gives NA,1=NA,2=3N0/2N_{A,1}=N_{A,2}=3N_0/2. To convert particles to moles, divide by 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 cannot cross the membrane: 1 mol1\ \mathrm{mol} of B remains on the left and 2 mol2\ \mathrm{mol} on the right. The total amounts are therefore 2,5 mol2{,}5\ \mathrm{mol} and 3,5 mol3{,}5\ \mathrm{mol}. At T0=300 KT_0=300\ \mathrm K, the equation P=nRT/VP=nRT/V gives separately

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

The total pressures differ; fixing the membrane maintains this state. A's partial pressures are equal because A occupies equal volumes with equal amounts and temperatures.

Movable membrane. From question 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.

The sum remains 20 L20\ \mathrm L, the vessel's volume. The A populations are now NA,1=N0N_{A,1}=N_0 and 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}.

The total is indeed 3 mol3\ \mathrm{mol} of A. Adding B, the compartments contain 22 and 4 mol4\ \mathrm{mol} of gas respectively. Use the exact volumes before rounding:

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

Both calculations give the same total pressure, as required for mechanical equilibrium of the released membrane.