21.6: Finding ß and the Thermodynamic Functions for Distinguishable Molecules

All of a substance’s thermodynamic functions can be derived from the molecular partition function. We begin with the entropy. We consider closed (constant N) systems of independent, distinguishable molecules in which only pressure–volume work is possible. In Sections 20.10 and 20.14, we find that two different approaches give the entropy of this system,

\[S=-Nk\sum^{\infty }_{i=1}{P_i}{ \ln \rho \left(\epsilon_i\right)\ }. \nonumber \]

In Sections 20.1, 20.3, and 20.4, we find that three different approaches give the Boltzmann equation,

\[P_i=g_i\rho \left(\epsilon_i\right)=z^{-1}g_i\ \mathrm{exp}\left(-\beta \epsilon_i\right). \nonumber \]

We have

\[ \ln \rho \left(\epsilon_i\right)=-\ln z -\beta \epsilon_i \nonumber \]

Substituting, and recognizing that the energy of the N-molecule system is \(E=N\left\langle \epsilon \right\rangle\), we find that the entropy of the system is

\[S=kN\sum^{\infty }_{i=1}{P_i}\left[{ \ln z\ }+\beta \epsilon_i\right]=kN{ \ln z\ }\sum^{\infty }_{i=1}{P_i}+k\beta N\sum^{\infty }_{i=1}{P_i}\epsilon_i=kN{ \ln z\ }+k\beta E \nonumber \]

In Section 10.1, we find that the fundamental equation implies that

\[{\left(\frac{\partial E}{\partial S}\right)}_V=T \nonumber \]

Since the \(\epsilon_i\) are fixed when the volume and temperature of the system are fixed, \({ \ln z\ }\) is constant when the volume and temperature of the system are constant. Differentiating \(S=kN{ \ln z\ }+k\beta E\) with respect to \(S\) at constant \(V\), we find

\[1=k\beta {\left(\frac{\partial E}{\partial S}\right)}_V=k\beta T \nonumber \] so that \[\beta =\frac{1}{kT} \nonumber \]

This is an important result: Because we have now identified all of the parameters in our microscopic model, we can write the results we have found in forms that are more useful:

1. \[\underbrace{z=\sum^{\infty}_{i=1}{g_i} \mathrm{exp}\left(\frac{-\epsilon_i}{kT}\right)}_{\text{molecular partition function}} \nonumber \]
2. \[\underbrace{P_i=g_i{\rho }\left({\epsilon}_i\right) =z^{-1}g_i\mathrm{exp}\left(\frac{-\epsilon_i}{kT}\right)}_{\text{Boltzmann’s equation}} \nonumber \]
3. \[\underbrace{S=kN \ln z+\frac{E}{T}}_{\text{Entropy of an N-molecule system}} \nonumber \]

To express the system energy in terms of the molecular partition function, we first observe that

\[E=N\left\langle \epsilon \right\rangle =N\sum^{\infty }_{i=1}{P_i\epsilon_i}=Nz^{-1}\sum^{\infty }_{i=1}{g_i\epsilon_i}\mathrm{exp}\left(\frac{-\epsilon_i}{kT}\right) \nonumber \]

Then we observe that

\[\begin{align*} \left(\frac{\partial \ln z}{\partial T}\right)_V &= z^{-1} \sum^{\infty}_{i=1} g_i \left(\frac{\epsilon_i}{kT^2}\right) \mathrm{exp}\left(\frac{-\epsilon_i}{kT}\right) \\[4pt] &= \left(\frac{1}{NkT^2}\right) N z^{-1} \sum^{\infty}_{i = 1} g_i \epsilon_i \mathrm{exp} \left(\frac{-\epsilon_i}{kT}\right) \\[4pt] &=\frac{E}{NkT^2} \end{align*} \nonumber \]

The system energy becomes

\[\underbrace{E=NkT^2{\left(\frac{\partial { \ln z\ }}{\partial T}\right)}_V}_{\text{energy of an N-molecule system}} \nonumber \]

By definition, \(A=E-TS\). Rearranging our entropy result, \(S=kN{ \ln z\ }+{E}/{T}\), we have \(E-TS=-NkT{ \ln z\ }\). Thus, \[A=-NkT{ \ln z\ } \nonumber \] (Helmholtz free energy of an N-molecule system)

From \(dA=-SdT-PdV\), we have

\[{\left(\frac{\partial A}{\partial V}\right)}_T=-P \nonumber \]

(Here, of course, \(P\) is the pressure of the system, not a probability.) Differentiating \(A=-NkT{ \ln z\ }\) with respect to \(V\) at constant \(T\), we find

\[P=NkT{\left(\frac{\partial { \ln z\ }}{\partial V}\right)}_T \nonumber \] (pressure of an N-molecule system)

The pressure–volume product becomes

\[PV=NkTV{\left(\frac{\partial { \ln z\ }}{\partial V}\right)}_T \nonumber \]

Substituting into \(H=E+PV\), the enthalpy becomes

\[H=NkT\left[T{\left(\frac{\partial { \ln z\ }}{\partial T}\right)}_V+V{\left(\frac{\partial { \ln z\ }}{\partial V}\right)}_T\right] \nonumber \] (enthalpy of an N-molecule system)

The Gibbs free energy is given by \(G=A+PV\). Substituting, we find

\[G=-NkT{ \ln z\ }+NkTV{\left(\frac{\partial { \ln z\ }}{\partial V}\right)}_T \nonumber \] (Gibbs free energy of an N-molecule system)

The chemical potential can be found from

\[\mu ={\left(\frac{\partial A}{\partial n}\right)}_{V,T} \nonumber \]

At constant volume and temperature, \(kT{ \ln z\ }\) is constant.

Substituting \(N=n\overline{N}\) into \(A=-NkT{ \ln z\ }\) and taking the partial derivative, we find

\[\underbrace{\mu =-\overline{N}kT{ \ln z\ }=-RT \ln z}_{\text{chemical potential of distinguishable molecules}} \nonumber \]

In statistical thermodynamics we frequently express the chemical potential per molecule, rather than per mole; then, \[\mu ={\left(\frac{\partial A}{\partial N}\right)}_{V,T} \nonumber \] and

\[\mu =-kT{ \ln z\ } \nonumber \](chemical potential per molecule)