The bulk of the discussion in this chapter dealt with ideal solutions. However, real solutions will deviate from this kind of behavior. So much as in the case of gases, where fugacity was introduced to allow us to use the ideal models, activity is used to allow for the deviation of real solutes from limiting ideal behavior. The activity of a solute is related to its concentration by
The problem that then remains is the measurement of the activity coefficients themselves, which may depend on temperature, pressure, and even concentration.
Activity Coefficients for Ionic Solutes
For an ionic substance that dissociates upon dissolving
\[ MX(s) \rightarrow M^+(aq) + X^-(aq) \nonumber \]
the chemical potential of the cation can be denoted \(\mu_+\) and that of the anion as \(\mu_-\). For a solution, the total molar Gibbs function of the solutes is given by
\[G = \mu_+ + \mu_- \nonumber \]
where
\[ \mu = \mu^* + RT \ln a \nonumber \]
where \(\mu^*\) denotes the chemical potential of an ideal solution, and \(a\) is the activity of the solute. Substituting his into the above relationship yields
\[G = \mu^*_+ + RT \ln a_+ + \mu_-^* + RT \ln a_- \nonumber \]
Using a molal definition for the activity coefficient
\[a_i = \gamma_im_i \nonumber \]
The expression for the total molar Gibbs function of the solutes becomes
\[G = \mu_+^* + RT \ln \gamma_+ m_+ + \mu_-^* + RT \ln \gamma_- m_- \nonumber \]
This expression can be rearranged to yield
\[ G = \mu_+^* + \mu_-^* + RT \ln m_+m_- + RT \ln \gamma_+\gamma _- \nonumber \]
where all of the deviation from ideal behavior comes from the last term. Unfortunately, it impossible to experimentally deconvolute the term into the specific contributions of the two ions. So instead, we use a geometric average to define the mean activity coefficient, \(\gamma _\pm\).
\[\gamma_{\pm} = \sqrt{\gamma_+\gamma_-} \nonumber \]
For a substance that dissociates according to the general process
\[ M_xX_y(s) \rightarrow x M^{y+} (aq) + yX^{x-} (aq) \nonumber \]
the expression for the mean activity coefficient is given by
\[ \gamma _{\pm} = (\gamma_+^x \gamma_-^y)^{1/x+y} \nonumber \]
Debeye-Hückel Law
In 1923, Debeye and Hückel (Debye & Hückel, 1923) suggested a means of calculating the mean activity coefficients from experimental data. Briefly, they suggest that
\[ \log _{10} \gamma_{\pm} = \dfrac{1.824 \times 10^6}{(\epsilon T)^{3/2}} |z_++z_- | \sqrt{I} \nonumber \]
where \(\epsilon\) is the dielectric constant of the solvent, \(T\) is the temperature in K, \(z_+\) and \(z_-\) are the charges on the ions, and \(I\) is the ionic strength of the solution. \(I\) is given by
\[ I = \dfrac{1}{2} \dfrac{m_+ z_+^2 + m_-z_-^2}{m^o} \nonumber \]
For a solution in water at 25 oC,