14.7: At Constant P and T, ∆rµ is the Change in Gibbs Free Energy
- Page ID
- 152394
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)At constant temperature and pressure, the chemical potential of component \(A\) is the contribution that one mole of \(A\) makes to the Gibbs free energy of the system. The Gibbs free energy of the system is the sum of the contributions made by all of its components. The chemical potential change that occurs during a reaction,
\[\mathrm{\ }{\Delta }_r\mu = \sum^{\omega }_{j=1}{{\mu }_j{\nu }_j} \nonumber \]
is the same thing as the Gibbs free energy change. To establish these points, we introduce Euler’s theorem on homogenous functions. For simplicity, we consider systems in which only pressure–volume work occurs.
A function \(f\left(x,y,z\right)\) is said to be homogeneous of order \(n\) if
\[{\lambda }^nf\left(x,y,z\right)=f\left(\lambda x,\lambda y,\lambda z\right) \nonumber \] Then,
\[\frac{d}{d\lambda }\left[{\lambda }^nf\left(x,y,z\right)\right]=n{\lambda }^{n-1}f\left(x,y,z\right)={\left(\frac{\partial f}{\partial \left(\lambda x\right)}\right)}_{yz}\frac{d\left(\lambda x\right)}{d\lambda }+{\left(\frac{\partial f}{\partial \left(\lambda y\right)}\right)}_{xz}\frac{d\left(\lambda y\right)}{d\lambda }+{\left(\frac{\partial f}{\partial \left(\lambda z\right)}\right)}_{xy}\frac{d\left(\lambda z\right)}{d\lambda } \nonumber \] \[=x{\left(\frac{\partial f}{\partial \left(\lambda x\right)}\right)}_{yz}+y{\left(\frac{\partial f}{\partial \left(\lambda y\right)}\right)}_{xz}+z{\left(\frac{\partial f}{\partial \left(\lambda z\right)}\right)}_{xy} \nonumber \]
Since this must be true for any \(\lambda\), it must be true for \(\lambda =1\). Making this substitution, we have Euler’s theorem for order \(n\):
\[nf\left(x,y,z\right)=x{\left(\frac{\partial f}{\partial x}\right)}_{yz}+y{\left(\frac{\partial f}{\partial y}\right)}_{xz}+z{\left(\frac{\partial f}{\partial z}\right)}_{xy} \nonumber \]
An extensive state function is homogeneous of order one in its extensive variables. In particular, we have
\[\lambda G\left(P,T,n_A,n_B,n_C,n_D\right)=G\left(P,T,{\lambda n}_A,{\lambda n}_B,{\lambda n}_C,{\lambda n}_D\right) \nonumber \]
for any \(\lambda\). (For, say, \(\lambda ={1}/{2}\), this says only that, if we divide a homogeneous equilibrium system into two equal portions, all else remaining constant, the Gibbs free energy of each portion will be half of that of the original system. The pressure and temperature are independent of \(\lambda\).) Taking the derivative with respect to \(\lambda\), we find
\[\frac{d}{d\lambda }\left[\lambda G\left(P,T,n_A,n_B,n_C,n_D\right)\right]=G\left(P,T,n_A,n_B,n_C,n_D\right)=n_A{\left(\frac{\partial G}{\partial \left(\lambda n_A\right)}\right)}_{PT}+n_B{\left(\frac{\partial G}{\partial \left(\lambda n_B\right)}\right)}_{PT}+n_C{\left(\frac{\partial G}{\partial \left(\lambda n_c\right)}\right)}_{PT}+n_D{\left(\frac{\partial G}{\partial \left(\lambda n_D\right)}\right)}_{PT} \nonumber \]
Since this must be true for any \(\lambda\), it must be true for \(\lambda =1\). We find
\[G\left(P,T,n_A,n_B,n_C,n_D\right)={\mu }_An_A+{\mu }_Bn_B+{\mu }_Cn_C+{\mu }_Dn_D \nonumber \]
The same equation follows from applying Euler’s theorem to other state functions. For example, viewing the internal energy as a function of entropy, volume, and composition, we have
\[\lambda E\left(S,V,n_A,n_B,n_C,n_D\right)=E\left(\lambda S,\lambda V,{\lambda n}_A,{\lambda n}_B,{\lambda n}_C,{\lambda n}_D\right) \nonumber \]
so that
\[\frac{d}{d\lambda }\left[\lambda E\left(S,V,n_A,n_B,n_C,n_D\right)\right]=E\left(S,V,n_A,n_B,n_C,n_D\right)=S{\left(\frac{\partial E}{\partial \left(\lambda S\right)}\right)}_V+V{\left(\frac{\partial E}{\partial \left(\lambda V\right)}\right)}_S+n_A{\left(\frac{\partial E}{\partial \left(\lambda n_A\right)}\right)}_{SV}+n_B{\left(\frac{\partial E}{\partial \left(\lambda n_B\right)}\right)}_{SV}+n_C{\left(\frac{\partial E}{\partial \left(\lambda n_c\right)}\right)}_{SV}+n_D{\left(\frac{\partial E}{\partial \left(\lambda n_D\right)}\right)}_{SV} \nonumber \]
Setting \(\lambda =1\), we have
\[E\left(S,V,n_A,n_B,n_C,n_D\right) \nonumber \] \[=TS-PV+{\mu }_An_A+{\mu }_Bn_B+{\mu }_Cn_C+{\mu }_Dn_D \nonumber \]
which we can rearrange to \[G=E-TS+PV={\mu }_An_A+{\mu }_Bn_B+{\mu }_Cn_C+{\mu }_Dn_D \nonumber \]
This equation describes any system whose thermodynamic functions are continuous functions of one another. Evidently, we can model the thermodynamic functions of any such system by modeling the chemical potentials of its components. In the remainder of this chapter, we develop this idea.
By definition, \({\Delta }_rG\) is the Gibbs free energy change for converting \(a\) moles of \(A\) and \(b\) moles of \(B\) to \(c\) moles of \(C\) and \(d\) moles of \(D\) at constant pressure and temperature and while the composition remains constant at \(n_A\), \(n_B\), \(n_C\), and \(n_D\). That is,
\[{\Delta }_rG=G\left(n_A,n_B,n_C+c,n_D+d\right)-G\left(n_A+a,n_B+b,n_C,n_D\right) \nonumber \]
where \(n_A\gg a\), \(n_B\gg b\), \(n_C\gg c\), and \(n_D\gg d\). Using the relationship we have just found between \(G\) and the \(n_i\), we have
\[{\Delta }_rG=\left[{\mu }_An_A+{\mu }_Bn_B+{\mu }_C\left(n_C+c\right)+{\mu }_D\left(n_D+d\right)\right] \nonumber \] \[-\left[{\mu }_A\left(n_A+a\right)+{\mu }_B\left(n_B+b\right)+{\mu }_Cn_C+{\mu }_Dn_D\right] \nonumber \] \[=c{\mu }_C+d{\mu }_D-a{\mu }_A-d{\mu }_d \nonumber \] \[={\Delta }_r\mu \nonumber \]
Since any other extensive thermodynamic function is also homogeneous of order one in the composition variables, a similar relationship will exist between the change in the function itself and its partial molar derivatives. For example, expressing the system volume as a function of pressure, temperature, and composition, we have
\[\lambda V\left(P,T,n_A,n_B,n_C,n_D\right)=V\left(P,T,{\lambda n}_A,{\lambda n}_B,{\lambda n}_C,{\lambda n}_D\right) \nonumber \]
and a completely parallel derivation shows that the volume of the system is related to the partial molar volumes of the components by
\[V\left(P,T,n_A,n_B,n_C,n_D\right)={\overline{V}}_An_A+{\overline{V}}_Bn_B+{\overline{V}}_Cn_C+{\overline{V}}_Dn_D \nonumber \]
The volume change for the reaction is given by
\[{\Delta }_rV=c{\overline{V}}_C+d{\overline{V}}_D-a{\overline{V}}_A-b{\overline{V}}_B \nonumber \]