# 9.6: Evaluation of Activity Coefficients

• • Howard DeVoe
• University of Maryland
$$\newcommand{\vecs}{\overset { \rightharpoonup} {\mathbf{#1}} }$$ $$\newcommand{\vecd}{\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 \|}$$ $$\newcommand{\inner}{\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 \|}$$ $$\newcommand{\inner}{\langle #1, #2 \rangle}$$ $$\newcommand{\Span}{\mathrm{span}}$$

$$\newcommand{\tx}{\text{#1}} % text in math mode$$
$$\newcommand{\subs}{_{\text{#1}}} % subscript text$$
$$\newcommand{\sups}{^{\text{#1}}} % superscript text$$
$$\newcommand{\st}{^\circ} % standard state symbol$$
$$\newcommand{\id}{^{\text{id}}} % ideal$$
$$\newcommand{\rf}{^{\text{ref}}} % reference state$$
$$\newcommand{\units}{\mbox{\thinspace#1}}$$
$$\newcommand{\K}{\units{K}} % kelvins$$
$$\newcommand{\degC}{^\circ\text{C}} % degrees Celsius$$
$$\newcommand{\br}{\units{bar}} % bar (\bar is already defined)$$
$$\newcommand{\Pa}{\units{Pa}}$$
$$\newcommand{\mol}{\units{mol}} % mole$$
$$\newcommand{\V}{\units{V}} % volts$$
$$\newcommand{\timesten}{\mbox{\,\times\,10^{#1}}}$$
$$\newcommand{\per}{^{-1}} % minus one power$$
$$\newcommand{\m}{_{\text{m}}} % subscript m for molar quantity$$
$$\newcommand{\CVm}{C_{V,\text{m}}} % molar heat capacity at const.V$$
$$\newcommand{\Cpm}{C_{p,\text{m}}} % molar heat capacity at const.p$$
$$\newcommand{\kT}{\kappa_T} % isothermal compressibility$$
$$\newcommand{\A}{_{\text{A}}} % subscript A for solvent or state A$$
$$\newcommand{\B}{_{\text{B}}} % subscript B for solute or state B$$
$$\newcommand{\bd}{_{\text{b}}} % subscript b for boundary or boiling point$$
$$\newcommand{\C}{_{\text{C}}} % subscript C$$
$$\newcommand{\f}{_{\text{f}}} % subscript f for freezing point$$
$$\newcommand{\mA}{_{\text{m},\text{A}}} % subscript m,A (m=molar)$$
$$\newcommand{\mB}{_{\text{m},\text{B}}} % subscript m,B (m=molar)$$
$$\newcommand{\mi}{_{\text{m},i}} % subscript m,i (m=molar)$$
$$\newcommand{\fA}{_{\text{f},\text{A}}} % subscript f,A (for fr. pt.)$$
$$\newcommand{\fB}{_{\text{f},\text{B}}} % subscript f,B (for fr. pt.)$$
$$\newcommand{\xbB}{_{x,\text{B}}} % x basis, B$$
$$\newcommand{\xbC}{_{x,\text{C}}} % x basis, C$$
$$\newcommand{\cbB}{_{c,\text{B}}} % c basis, B$$
$$\newcommand{\mbB}{_{m,\text{B}}} % m basis, B$$
$$\newcommand{\kHi}{k_{\text{H},i}} % Henry's law constant, x basis, i$$
$$\newcommand{\kHB}{k_{\text{H,B}}} % Henry's law constant, x basis, B$$
$$\newcommand{\arrow}{\,\rightarrow\,} % right arrow with extra spaces$$
$$\newcommand{\arrows}{\,\rightleftharpoons\,} % double arrows with extra spaces$$
$$\newcommand{\ra}{\rightarrow} % right arrow (can be used in text mode)$$
$$\newcommand{\eq}{\subs{eq}} % equilibrium state$$
$$\newcommand{\onehalf}{\textstyle\frac{1}{2}\D} % small 1/2 for display equation$$
$$\newcommand{\sys}{\subs{sys}} % system property$$
$$\newcommand{\sur}{\sups{sur}} % surroundings$$
$$\renewcommand{\in}{\sups{int}} % internal$$
$$\newcommand{\lab}{\subs{lab}} % lab frame$$
$$\newcommand{\cm}{\subs{cm}} % center of mass$$
$$\newcommand{\rev}{\subs{rev}} % reversible$$
$$\newcommand{\irr}{\subs{irr}} % irreversible$$
$$\newcommand{\fric}{\subs{fric}} % friction$$
$$\newcommand{\diss}{\subs{diss}} % dissipation$$
$$\newcommand{\el}{\subs{el}} % electrical$$
$$\newcommand{\cell}{\subs{cell}} % cell$$
$$\newcommand{\As}{A\subs{s}} % surface area$$
$$\newcommand{\E}{^\mathsf{E}} % excess quantity (superscript)$$
$$\newcommand{\allni}{\{n_i \}} % set of all n_i$$
$$\newcommand{\sol}{\hspace{-.1em}\tx{(sol)}}$$
$$\newcommand{\solmB}{\tx{(sol,\,m\B)}}$$
$$\newcommand{\dil}{\tx{(dil)}}$$
$$\newcommand{\sln}{\tx{(sln)}}$$
$$\newcommand{\mix}{\tx{(mix)}}$$
$$\newcommand{\rxn}{\tx{(rxn)}}$$
$$\newcommand{\expt}{\tx{(expt)}}$$
$$\newcommand{\solid}{\tx{(s)}}$$
$$\newcommand{\liquid}{\tx{(l)}}$$
$$\newcommand{\gas}{\tx{(g)}}$$
$$\newcommand{\pha}{\alpha} % phase alpha$$
$$\newcommand{\phb}{\beta} % phase beta$$
$$\newcommand{\phg}{\gamma} % phase gamma$$
$$\newcommand{\aph}{^{\alpha}} % alpha phase superscript$$
$$\newcommand{\bph}{^{\beta}} % beta phase superscript$$
$$\newcommand{\gph}{^{\gamma}} % gamma phase superscript$$
$$\newcommand{\aphp}{^{\alpha'}} % alpha prime phase superscript$$
$$\newcommand{\bphp}{^{\beta'}} % beta prime phase superscript$$
$$\newcommand{\gphp}{^{\gamma'}} % gamma prime phase superscript$$
$$\newcommand{\apht}{\small\aph} % alpha phase tiny superscript$$
$$\newcommand{\bpht}{\small\bph} % beta phase tiny superscript$$
$$\newcommand{\gpht}{\small\gph} % gamma phase tiny superscript$$

$$\newcommand{\upOmega}{\Omega}$$

$$\newcommand{\dif}{\mathop{}\!\mathrm{d}} % roman d in math mode, preceded by space$$
$$\newcommand{\Dif}{\mathop{}\!\mathrm{D}} % roman D in math mode, preceded by space$$
$$\newcommand{\df}{\dif\hspace{0.05em} f} % df$$

$$\newcommand{\dBar}{\mathop{}\!\mathrm{d}\hspace-.3em\raise1.05ex{\Rule{.8ex}{.125ex}{0ex}}} % inexact differential$$
$$\newcommand{\dq}{\dBar q} % heat differential$$
$$\newcommand{\dw}{\dBar w} % work differential$$
$$\newcommand{\dQ}{\dBar Q} % infinitesimal charge$$
$$\newcommand{\dx}{\dif\hspace{0.05em} x} % dx$$
$$\newcommand{\dt}{\dif\hspace{0.05em} t} % dt$$
$$\newcommand{\difp}{\dif\hspace{0.05em} p} % dp$$
$$\newcommand{\Del}{\Delta}$$
$$\newcommand{\Delsub}{\Delta_{\text{#1}}}$$
$$\newcommand{\pd}{(\partial #1 / \partial #2 )_{#3}} % \pd{}{}{} - partial derivative, one line$$
$$\newcommand{\Pd}{\left( \dfrac {\partial #1} {\partial #2}\right)_{#3}} % Pd{}{}{} - Partial derivative, built-up$$
$$\newcommand{\bpd}{[ \partial #1 / \partial #2 ]_{#3}}$$
$$\newcommand{\bPd}{\left[ \dfrac {\partial #1} {\partial #2}\right]_{#3}}$$
$$\newcommand{\dotprod}{\small\bullet}$$
$$\newcommand{\fug}{f} % fugacity$$
$$\newcommand{\g}{\gamma} % solute activity coefficient, or gamma in general$$
$$\newcommand{\G}{\varGamma} % activity coefficient of a reference state (pressure factor)$$
$$\newcommand{\ecp}{\widetilde{\mu}} % electrochemical or total potential$$
$$\newcommand{\Eeq}{E\subs{cell, eq}} % equilibrium cell potential$$
$$\newcommand{\Ej}{E\subs{j}} % liquid junction potential$$
$$\newcommand{\mue}{\mu\subs{e}} % electron chemical potential$$
$$\newcommand{\defn}{\,\stackrel{\mathrm{def}}{=}\,} % "equal by definition" symbol$$

$$\newcommand{\D}{\displaystyle} % for a line in built-up$$
$$\newcommand{\s}{\smash[b]} % use in equations with conditions of validity$$
$$\newcommand{\cond}{\\[-2.5pt]{}\tag*{#1}}$$
$$\newcommand{\nextcond}{\\[-5pt]{}\tag*{#1}}$$
$$\newcommand{\R}{8.3145\units{J\,K\per\,mol\per}} % gas constant value$$
$$\newcommand{\Rsix}{8.31447\units{J\,K\per\,mol\per}} % gas constant value - 6 sig figs$$

$$\newcommand{\jn}{\hspace3pt\lower.3ex{\Rule{.6pt}{2ex}{0ex}}\hspace3pt}$$
$$\newcommand{\ljn}{\hspace3pt\lower.3ex{\Rule{.6pt}{.5ex}{0ex}}\hspace-.6pt\raise.45ex{\Rule{.6pt}{.5ex}{0ex}}\hspace-.6pt\raise1.2ex{\Rule{.6pt}{.5ex}{0ex}} \hspace3pt}$$
$$\newcommand{\lljn}{\hspace3pt\lower.3ex{\Rule{.6pt}{.5ex}{0ex}}\hspace-.6pt\raise.45ex{\Rule{.6pt}{.5ex}{0ex}}\hspace-.6pt\raise1.2ex{\Rule{.6pt}{.5ex}{0ex}}\hspace1.4pt\lower.3ex{\Rule{.6pt}{.5ex}{0ex}}\hspace-.6pt\raise.45ex{\Rule{.6pt}{.5ex}{0ex}}\hspace-.6pt\raise1.2ex{\Rule{.6pt}{.5ex}{0ex}}\hspace3pt}$$

This section describes several methods by which activity coefficients of nonelectrolyte substances may be evaluated. Section 9.6.3 describes an osmotic coefficient method that is also suitable for electrolyte solutes, as will be explained in Sec. 10.6.

## 9.6.1 Activity coefficients from gas fugacities

Suppose we equilibrate a liquid mixture with a gas phase. If component $$i$$ of the liquid mixture is a volatile nonelectrolyte, and we are able to evaluate its fugacity $$\fug_i$$ in the gas phase, we have a convenient way to evaluate the activity coefficient $$\g_i$$ in the liquid. The relation between $$\g_i$$ and $$\fug_i$$ will now be derived.

When component $$i$$ is in transfer equilibrium between two phases, its chemical potential is the same in both phases. Equating expressions for $$\mu_i$$ in the liquid mixture and the equilibrated gas phase (from Eqs. 9.5.14 and 9.5.11, respectively), and then solving for $$\g_i$$, we have \begin{equation} \mu_i^* + RT\ln\left(\g_i x_i\right) = \mu_i\rf\gas + RT\ln (\fug_i/p) \tag{9.6.1} \end{equation} \begin{equation} \g_i = \exp \left[ \frac{\mu_i\rf\gas -\mu_i^*}{RT} \right] \times \frac{\fug_i}{x_i p} \tag{9.6.2} \end{equation} On the right side of Eq. 9.6.2, only $$\fug_i$$ and $$x_i$$ depend on the liquid composition. We can therefore write \begin{equation} \g_i = C_i\frac{\fug_i}{x_i} \tag{9.6.3} \end{equation} where $$C_i$$ is a factor whose value depends on $$T$$ and $$p$$, but not on the liquid composition. Solving Eq. 9.6.3 for $$C_i$$ gives $$C_i=\g_i x_i/\fug_i$$.

Now consider Eq. 9.5.20. It says that as $$x_i$$ approaches 1 at constant $$T$$ and $$p$$, $$\g_i$$ also approaches 1. We can use this limit to evaluate $$C_i$$: \begin{equation} C_i = \lim_{x_i \ra 1}\frac{\g_i x_i}{\fug_i} = \frac{1}{\fug_i^*} \tag{9.6.4} \end{equation} Here $$\fug_i^*$$ is the fugacity of $$i$$ in a gas phase equilibrated with pure liquid $$i$$ at the temperature and pressure of the mixture. Then substitution of this value of $$C_i$$ (which is independent of $$x_i$$) in Eq. 9.6.3 gives us an expression for $$\g_i$$ at any liquid composition: \begin{equation} \g_i=\frac{\fug_i}{x_i\fug_i^*} \tag{9.6.5} \end{equation}

We can follow the same procedure for a solvent or solute of a liquid solution. We replace the left side of Eq. 9.6.1 with an expression from among Eqs. 9.5.15–9.5.18, then derive an expression analogous to Eq. 9.6.3 for the activity coefficient with a composition-independent factor, and finally apply the limiting conditions that cause the activity coefficient to approach unity (Eqs. 9.5.21–9.5.24) and allow us to evaluate the factor. When we take the limits that cause the solute activity coefficients to approach unity, the ratios $$\fug\B/x\B$$, $$\fug\B/c\B$$, and $$\fug\B/m\B$$ become Henry’s law constants (Eqs. 9.4.19–9.4.21). The resulting expressions for activity coefficients as functions of fugacity are listed in Table 9.4.

### Examples

Figure 9.11(a) shows the function $$(\phi_m - 1)/m\B$$ for aqueous sucrose solutions over a wide range of molality. The dependence of the solute activity coefficient on molality, generated from Eq. 9.6.20, is shown in Fig. 9.11(b). Figure 9.11(c) is a plot of the effective sucrose molality $$\g\mbB m\B$$ as a function of composition. Note how the activity coefficient becomes greater than unity beyond the ideal-dilute region, and how in consequence the effective molality $$\g\mbB m\B$$ becomes considerably greater than the actual molality $$m\B$$.

## 9.6.4 Fugacity measurements

Section 9.6.1 described the evaluation of the activity coefficient of a constituent of a liquid mixture from its fugacity in a gas phase equilibrated with the mixture. Section 9.6.3 mentioned the use of solvent fugacities in gas phases equilibrated with pure solvent and with a solution, in order to evaluate the osmotic coefficient of the solution.

Various experimental methods are available for measuring a partial pressure in a gas phase equilibrated with a liquid mixture. A correction for gas nonideality, such as that given by Eq. 9.3.16, can be used to convert the partial pressure to fugacity.

If the solute of a solution is nonvolatile, we may pump out the air above the solution and use a manometer to measure the pressure, which is the partial pressure of the solvent. Dynamic methods involve passing a stream of inert gas through a liquid mixture and analyzing the gas mixture to evaluate the partial pressures of volatile components. For instance, we could pass dry air successively through an aqueous solution and a desiccant and measure the weight gained by the desiccant.

The isopiestic vapor pressure technique is one of the most useful methods for determining the fugacity of H$$_2$$O in a gas phase equilibrated with an aqueous solution. This is a comparative method using a binary solution of the solute of interest, B, and a nonvolatile reference solute of known properties. Some commonly used reference solutes for which data are available are sucrose, NaCl, and CaCl$$_2$$.

In this method, solute B can be either a nonelectrolyte or electrolyte. Dishes, each containing water and an accurately weighed sample of one of the solutes, are placed in wells drilled in a block made of metal for good thermal equilibration. The assembly is placed in a gas-tight chamber, the air is evacuated, and the apparatus is gently rocked in a thermostat for a period of up to several days, or even weeks. During this period, H$$_2$$O is transferred among the dishes through the vapor space until the chemical potential of the water becomes the same in each solution. The solutions are then said to be isopiestic. Finally, the dishes are removed from the apparatus and weighed to establish the molality of each solution. The H$$_2$$O fugacity is known as a function of the molality of the reference solute, and is the same as the H$$_2$$O fugacity in equilibrium with the solution of solute B at its measured molality.

The isopiestic vapor pressure method can also be used for nonaqueous solutions.

This page titled 9.6: Evaluation of Activity Coefficients is shared under a CC BY 4.0 license and was authored, remixed, and/or curated by Howard DeVoe via source content that was edited to the style and standards of the LibreTexts platform; a detailed edit history is available upon request.