Skip to main content
Chemistry LibreTexts

12.5 Solid-Liquid Equilibria

  • Page ID
    20643
  • \( \newcommand{\tx}[1]{\text{#1}}      % text in math mode\)
     \( \newcommand{\subs}[1]{_{\text{#1}}} % subscript text\)
     \( \newcommand{\sups}[1]{^{\text{#1}}} % superscript text\)
     \( \newcommand{\st}{^\circ}            % standard state symbol\)
     \( \newcommand{\id}{^{\text{id}}}      % ideal\)
     \( \newcommand{\rf}{^{\text{ref}}}     % reference state\)
     \( \newcommand{\units}[1]{\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}[1]{\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}[1]{\Delta_{\text{#1}}}\)
     \( \newcommand{\pd}[3]{(\partial #1 / \partial #2 )_{#3}} % \pd{}{}{} - partial derivative, one line\)
     \( \newcommand{\Pd}[3]{\left( \dfrac {\partial #1} {\partial #2}\right)_{#3}} % Pd{}{}{} - Partial derivative, built-up\)
     \( \newcommand{\bpd}[3]{[ \partial #1 / \partial #2 ]_{#3}}\)
     \( \newcommand{\bPd}[3]{\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}[1]{\\[-2.5pt]{}\tag*{#1}}\)
     \( \newcommand{\nextcond}[1]{\\[-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} \) 

    A freezing-point curve (freezing point as a function of liquid composition) and a solubility curve (composition of a solution in equilibrium with a pure solid as a function of temperature) are different ways of describing the same physical situation. Thus, strange as it may sound, the composition \(x\A\) of an aqueous solution at the freezing point is the mole fraction solubility of ice in the solution.

    12.5.1 Freezing points of ideal binary liquid mixtures

    Section 12.2.1 described the use of freezing-point measurements to determine the solvent chemical potential in a solution of arbitrary composition relative to the chemical potential of the pure solvent. The way in which freezing point varies with solution composition in the limit of infinite dilution was derived in Sec. 12.4.1. Now let us consider the freezing behavior over the entire composition range of an ideal liquid mixture.

    Figure 12.4 Dependence on composition of the freezing point of binary liquid mixtures with benzene as component A (experimental data from R. Negishi, Rev. Phys. Chem. Japan, 15, 98–116, 1941). Solid curve: calculated for an ideal liquid mixture (Eq. 12.5.2), taking the temperature variation of \(\Delsub{sol,A}H\) into account. Open circles: B = toluene. Open triangles: B = cyclohexane.

    The general relation between temperature and the composition of a binary liquid mixture, when the mixture is in equilibrium with pure solid A, is given by Eq. 12.3.6: \begin{equation} \Pd{T}{x\A}{\!p} = \frac{T^2}{\Delsub{sol,A}H} \bPd{(\mu\A/T)}{x\A}{T,p} \tag{12.5.1} \end{equation} We can replace \(T\) by \(T\fA\) to indicate this is the temperature at which the mixture freezes to form solid A. From the expression for the chemical potential of component A in an ideal liquid mixture, \(\mu\A = \mu\A^* + RT\ln x\A\), we have \(\bpd{(\mu\A/T)}{x\A}{T,p}=R/x\A\). With these substitutions, Eq. 12.5.1 becomes \begin{gather} \s{ \Pd{T\fA}{x\A}{\!p} = \frac{RT\fA^2}{x\A\Delsub{sol,A}H} } \tag{12.5.2} \cond{(ideal liquid mixture)} \end{gather} Figure 12.4 compares the freezing behavior of benzene predicted by this equation with experimental freezing-point data for mixtures of benzene–toluene and benzene–cyclohexane. Any constituent that forms an ideal liquid mixture with benzene should give freezing points for the formation of solid benzene that fall on the curve in this figure. The agreement is good over a wide range of compositions for benzene–toluene mixtures (open circles), which are known to closely approximate ideal liquid mixtures. The agreement for benzene–cyclohexane mixtures (open triangles), which are not ideal liquid mixtures, is confined to the ideal-dilute region.

    Figure 12.5 Freezing-point curves of ideal binary liquid mixtures. The solid is component A. Each curve is calculated from Eq. 12.5.4 and is labeled with the value of \(\Delsub{fus,A}H/RT\fA^*\).

    If we make the approximation that \(\Delsub{sol,A}H\) is constant over the entire range of mixture composition, we can replace it by \(\Delsub{fus,A}H\), the molar enthalpy of fusion of pure solid A at its melting point. This approximation allows us to separate the variables in Eq. 12.5.2 and integrate as follows from an arbitrary mixture composition \(x'\A\) at the freezing point \(T'\fA\) to pure liquid A at its freezing point \(T\fA^*\): \begin{equation} \int_{T'\fA}^{T\fA^*} \frac{\dif T}{T^2} = \frac{R}{\Delsub{fus,A}H}\int_{x'\A}^{1}\frac{\dx\A}{x\A} \tag{12.5.3} \end{equation} The result, after some rearrangement, is \begin{gather} \s{ \ln x\A = \frac{\Delsub{fus,A}H}{R} \left( \frac{1}{T\fA^*} - \frac{1}{T\fA} \right) } \tag{12.5.4} \cond{(ideal liquid mixture,} \nextcond{\(\Delsub{sol,A}H{=}\Delsub{fus,A}H\))} \end{gather} This equation was used to generate the curves shown in Fig. 12.5. Although the shape of the freezing-point curve (\(T\fA\) versus \(x\B\)) shown in Fig. 12.4 is concave downward, Fig. 12.5 shows this is not always the case. When \(\Delsub{fus,A}H/RT\fA^*\) is less than \(2\), the freezing-point curve at low \(x\B\) is concave upward.

    12.5.2 Solubility of a solid nonelectrolyte

    Suppose we find that a solution containing solute B at a particular combination of temperature, pressure, and composition can exist in transfer equilibrium with pure solid B at the same temperature and pressure. This solution is said to be saturated with respect to the solid. We can express the solubility of the solid in the solvent by the value of the mole fraction, concentration, or molality of B in the saturated solution. We can also define solubility as the maximum value of the solute mole fraction, concentration, or molality that can exist in the solution without the possibility of spontaneous precipitation.

    This section considers the solubility of a solid nonelectrolyte. For the solution process B(s) \(\ra\) B(sln), the general expression for the thermodynamic equilibrium constant is \(K=a\B\sln/a\B\solid\). (In this and other expressions for equilibrium constants in this chapter, activities will be assumed to be for equilibrium states, although not indicated by the “eq” subscripts used in Chap. 11.) The activity of the pure solid is \(a\B\solid =\G\B\solid\). Let us use a solute standard state based on mole fraction; then the solute activity is \(a\B\sln=\G\xbB \g\xbB x\B\). From these relations, the solubility expressed as a mole fraction is \begin{equation} x\B = \frac{\G\B\solid K}{\G\xbB \g\xbB} \tag{12.5.5} \end{equation}

    If we measure the solubility at the standard pressure, the pressure factors \(\G\B\solid\) and \(\G\xbB\) are unity and the solubility is given by \begin{gather} \s{ x\B = \frac{K}{\g\xbB} } \tag{12.5.6} \cond{(solubility of solid B, \(p{=}p\st\))} \end{gather} If the pressure is not exactly equal to \(p\st\), but is not very much greater, the values of the pressure factors are close to unity and Eq. 12.5.6 is a good approximation.

    We can find the standard molar enthalpy of solution of B from the temperature dependence of the solubility. Combining Eqs. 12.1.12 and 12.5.6, we obtain \begin{gather} \s{ \Delsub{sol,B}H\st = RT^2\frac{\dif\ln(\g\xbB x\B)}{\dif T} } \tag{12.5.7} \cond{(\(p{=}p\st\))} \end{gather} The solubility may be small enough for us to be able to set the solute activity coefficient equal to \(1\), in which case Eq. 12.5.7 becomes \begin{gather} \Delsub{sol,B} \tag{12.5.8} H\st = RT^2 \frac{\dif\ln x\B}{\dif T} \cond{(\(p{=}p\st\), \(\g\xbB{=}1\))} \end{gather}

    If the solubility \(x\B\) increases with increasing temperature, \(\Delsub{sol,B}H\st\) must be positive and the solution process is endothermic. A decrease of solubility with increasing temperature implies an exothermic solution process. These statements refer to a solid of low solubility; see Sec. 11.9 for a discussion of the general relation between the temperature dependence of solubility and the sign of the molar differential enthalpy of solution at saturation.

    For a solute standard state based on molality, we can derive equations like Eqs. 12.5.7 and 12.5.8 with \(\g\xbB\) replaced by \(\g\mbB\) and \(x\B\) replaced by \(m\B/m\st\). If we use a solute standard state based on concentration, the expressions become slightly more complicated. The solubility in this case is given by \begin{equation} c\B =\frac{\G\B\solid Kc\st}{\G\cbB \g\cbB} \tag{12.5.9} \end{equation} From Eq. 12.1.11, we obtain, for a nonelectrolyte solid of low solubility, the relation \begin{gather} \Delsub{sol,B} \tag{12.5.10} H\st = RT^2 \left( \frac{\dif\ln(c\B/c\st)}{\dif T} + \alpha\A^* \right) \cond{(\(p{=}p\st\), \(\g\cbB{=}1\))} \end{gather}

    12.5.3 Ideal solubility of a solid

    The ideal solubility of a solid at a given temperature and pressure is the solubility calculated on the assumptions that (1) the liquid is an ideal liquid mixture, and (2) the molar differential enthalpy of solution equals the molar enthalpy of fusion of the solid (\(\Delsub{sol,B}H{=}\Delsub{fus,B}H\)). These were the assumptions used to derive Eq. 12.5.4 for the freezing-point curve of an ideal liquid mixture. In Eq. 12.5.4, we exchange the constituent labels A and B so that the solid phase is now component B: \begin{gather} \s{ \ln x\B = \frac{\Delsub{fus,B}H}{R} \left( \frac{1}{T\fB^*}-\frac{1}{T} \right) } \tag{12.5.11} \cond{(ideal solubility of solid B)} \end{gather} Here \(T\fB^*\) is the melting point of solid B.

    Figure 12.6 Ideal solubility of solid B as a function of \(T\). The curves are calculated for two solids having the same molar enthalpy of fusion (\(\Delsub{fus,B}H=20\units{kJ mol\(^{-1}\)}\)) and the values of \(T\fB^*\) indicated.

    According to Eq. 12.5.11, the ideal solubility of a solid is independent of the kind of solvent and increases with increasing temperature. For solids with similar molar enthalpies of fusion, the ideal solubility is less at a given temperature the higher is the melting point. This behavior is shown in Fig. 12.6. In order for the experimental solubility of a solid to agree even approximately with the ideal value, the solvent and solute must be chemically similar, and the temperature must be close to the melting point of the solid so that \(\Delsub{sol,B}H\) is close in value to \(\Delsub{fus,B}H\).

    From the freezing behavior of benzene–toluene mixtures shown by the open circles in Fig. 12.4, we can see that solid benzene has close to ideal solubility in liquid toluene at temperatures not lower than about \(20\K\) below the melting point of benzene.

    12.5.4 Solid compound of mixture components

    Binary liquid mixtures are known in which the solid that appears when the mixture is cooled is a compound containing both components in a fixed proportion. This kind of solid is called a solid compound, or stoichiometric addition compound. Examples are salt hydrates (salts with fixed numbers of waters of hydration in the formula unit) and certain metal alloys.

    The composition of the liquid mixture in this kind of system is variable, whereas the composition of the solid compound is fixed. Suppose the components are A and B, present in the liquid mixture at mole fractions \(x\A\) and \(x\B\), and the solid compound has the formula A\(_a\)B\(_b\). We assume that in the liquid phase the compound is completely dissociated with respect to the components; that is, that no molecules of formula A\(_a\)B\(_b\) exist in the liquid. The reaction equation for the freezing process is \[ a\tx{A(mixt)} + b\tx{B(mixt)} \arrow \tx{A\(_a\)B\(_b\)(s)} \] When equilibrium exists between the liquid and solid phases, the temperature is the freezing point \(T\subs{f}\) of the liquid. At equilibrium, the molar reaction Gibbs energy defined by \(\Delsub{r}G = \sum_i\!\nu_i \mu_i\) is zero: \begin{equation} -a\mu\A-b\mu\B+\mu\solid =0 \tag{12.5.12} \end{equation} Here \(\mu\A\) and \(\mu\B\) refer to chemical potentials in the liquid mixture, and \(\mu\solid\) refers to the solid compound.

    How does the freezing point of the liquid mixture vary with composition? We divide both sides of Eq. 12.5.12 by \(T\) and take differentials: \begin{gather} \s{ -a\dif(\mu\A/T)-b\dif(\mu\B/T)+\dif[\mu\solid /T] = 0 } \tag{12.5.13} \cond{(phase equilibrium)} \end{gather} The pressure is constant. Then \(\mu\A/T\) and \(\mu\B/T\) are functions of \(T\) and \(x\A\), and \(\mu\solid /T\) is a function only of \(T\). We find expressions for the total differentials of these quantities at constant \(p\) with the help of Eq. 12.1.3: \begin{equation} \dif(\mu\A/T) = -\frac{H\A}{T^2}\dif T + \frac{1}{T}\Pd{\mu\A}{x\A}{T,p} \dx\A \tag{12.5.14} \end{equation} \begin{equation} \dif(\mu\B/T) = -\frac{H\B}{T^2}\dif T + \frac{1}{T}\Pd{\mu\B}{x\A}{T,p}\dx\A \tag{12.5.15} \end{equation} \begin{equation} \dif[\mu\solid /T] = -\frac{H\m\solid }{T^2}\dif T \tag{12.5.16} \end{equation} When we substitute these expressions in Eq. 12.5.13 and solve for \(\dif T/\dx\A\), setting \(T\) equal to \(T\subs{f}\), we obtain \begin{equation} \frac{\dif T\subs{f}}{\dx\A} = \frac{T\subs{f}}{aH\A+bH\B-H\m\solid } \left[ a\Pd{\mu\A}{x\A}{T,p}+b\Pd{\mu\B}{x\A}{T,p} \right] \tag{12.5.17} \end{equation}

    The quantity \(aH\A+bH\B-H\m\solid\) in the denominator on the right side of Eq. 12.5.17 is \(\Delsub{sol}H\), the molar differential enthalpy of solution of the solid compound in the liquid mixture. The two partial derivatives on the right side are related through the Gibbs–Duhem equation \(x\A\dif\mu\A+x\B\dif\mu\B=0\) (Eq. 9.2.27), which applies to changes at constant \(T\) and \(p\). We rearrange the Gibbs–Duhem equation to \(\dif\mu\B=-(x\A/x\B)\dif\mu\A\) and divide by \(\dx\A\): \begin{equation} \Pd{\mu\B}{x\A}{T,p} = -\frac{x\A}{x\B} \Pd{\mu\A}{x\A}{T,p} \tag{12.5.18} \end{equation} Making this substitution in Eq. 12.5.17, we obtain the equation \begin{equation} \frac{\dif T\subs{f}}{\dx\A} = \frac{x\A T\subs{f}}{\Delsub{sol}H} \left( \frac{a}{x\A}-\frac{b}{x\B} \right)\Pd{\mu\A}{x\A}{T,p} \tag{12.5.19} \end{equation} which can also be written in the slightly rearranged form \begin{equation} \frac{\dif T\subs{f}}{\dx\A} = \frac{b T\subs{f}}{\Delsub{sol}H} \left( \frac{a}{b}-\frac{x\A}{1-x\A} \right)\Pd{\mu\A}{x\A}{T,p} \tag{12.5.20} \end{equation}

    Suppose we heat a sample of the solid compound to its melting point to form a liquid mixture of the same composition as the solid. The molar enthalpy change of the fusion process is the molar enthalpy of fusion of the solid compound, \(\Delsub{fus}H\), a positive quantity. When the liquid has the same composition as the solid, the dissolution and fusion processes are identical; under these conditions, \(\Delsub{sol}H\) is equal to \(\Delsub{fus}H\) and is positive.

    Equation 12.5.20 shows that the slope of the freezing-point curve, \(T\subs{f}\) versus \(x\A\), is zero when \(x\A/(1-x\A)\) is equal to \(a/b\), or \(x\A=a/(a+b)\); that is, when the liquid and solid have the same composition. Because \(\pd{\mu\A}{x\A}{T,p}\) is positive, and \(\Delsub{sol}H\) at this composition is also positive, we see from the equation that the slope decreases as \(x\A\) increases. Thus, the freezing-point curve has a maximum at the mixture composition that is the same as the composition of the solid compound. This conclusion applies when both components of the liquid mixture are nonelectrolytes, and also when one component is an electrolyte that dissociates into ions.

    Now let us assume the liquid mixture is an ideal liquid mixture of nonelectrolytes in which \(\mu\A\) obeys Raoult’s law for fugacity, \(\mu\A=\mu\A^*+RT\ln x\A\). The partial derivative \(\pd{\mu\A}{x\A}{T,p}\) then equals \(RT/x\A\), and Eq. 12.5.19 becomes \begin{equation} \frac{\dif T\subs{f}}{\dx\A} = \frac{RT\subs{f}^2}{\Delsub{sol}H} \left( \frac{a}{x\A}-\frac{b}{x\B} \right) \tag{12.5.21} \end{equation} By making the approximations that \(\Delsub{sol}H\) is independent of \(T\) and \(x\A\), and is equal to \(\Delsub{fus}H\), we can separate the variables and integrate as follows: \begin{equation} \int_{T'\subs{f}}^{T''\subs{f}}\frac{\dif T\subs{f}}{T\subs{f}^2} = \frac{R}{\Delsub{fus}H}\left( \int_{x'\A}^{x''\A}\frac{a}{x\A}\dx\A +\int_{x'\B}^{x''\B}\frac{b}{x\B}\dx\B \right) \tag{12.5.22} \end{equation} (The second integral on the right side comes from changing \(\dx\A\) to \(-\dx\B\).) The result of the integration is \begin{gather} \s{ \frac{1}{T'\subs{f}} = \frac{1}{T''\subs{f}} + \frac{R}{\Delsub{fus}H} \left( a\ln\frac{x''\A}{x'\A}+b\ln\frac{x''\B}{x'\B} \right) } \tag{12.5.23} \cond{(ideal liquid mixture in} \nextcond{equilibrium with solid} \nextcond{compound, \(\Delsub{sol}H{=}\Delsub{fus}H\))} \end{gather}

    Figure 12.7 Solid curve: freezing-point curve of a liquid melt of Zn and Mg that solidifies to the solid compound Zn\(_2\)Mg (Rodney P. Elliot, Constitution of Binary Alloys, First Supplement, McGraw-Hill, New York, 1965, p. 603). The curve maximum (open circle) is at the compound composition \(x''\subs{Zn}=2/3\) and the solid compound melting point \(T''\subs{f}=861\K\). Dashed curve: calculated using Eq. 12.5.23 with \(\Delsub{fus}H = 15.8\units{kJ mol\(^{-1}\)}\).

    Let \(T'\subs{f}\) be the freezing point of a liquid mixture of composition \(x'\A\) and \(x'\B=1-x'\A\), and let \(T''\subs{f}\) be the melting point of the solid compound of composition \(x''\A=a/(a+b)\) and \(x''\B=b/(a+b)\). Figure 12.7 shows an example of a molten metal mixture that solidifies to an alloy of fixed composition. The freezing-point curve of this system is closely approximated by Eq. 12.5.23.

    12.5.5 Solubility of a solid electrolyte

    Consider an equilibrium between a crystalline salt (or other kind of ionic solid) and a solution containing the solvated ions: \[ \tx{M\(_{\nu_+}\)X\(_{\nu_-}\)(s)} \arrows \nu_+\tx{M\(^{z_+}\)(aq)} + \nu_-\tx{X\(^{z_-}\)(aq)} \] Here \(\nu_+\) and \(\nu_-\) are the numbers of cations and anions in the formula unit of the salt, and \(z_+\) and \(z_-\) are the charge numbers of these ions. The solution in equilibrium with the solid salt is a saturated solution. The thermodynamic equilibrium constant for this kind of equilibrium is called a solubility product, \(K\subs{s}\).

    We can readily derive a relation between \(K\subs{s}\) and the molalities of the ions in the saturated solution by treating the dissolved salt as a single solute substance, B. We write the equilibrium in the form B\(^*\)(s)\(\arrows\)B(sln), and write the expression for the solubility product as a proper quotient of activities: \begin{equation} K\subs{s} = \frac{a\mbB}{a\B^*} \tag{12.5.24} \end{equation} From Eq. 10.3.16, we have \(a\mbB= \G\mbB \g_{\pm}^\nu(m_+/m\st)^{\nu_+}(m_-/m\st)^{\nu_-}\). This expression is valid whether or not the ions M\(^{z_+}\) and X\(^{z_-}\) are present in solution in the same ratio as in the solid salt. When we replace \(a\mbB\) with this expression, and replace \(a\B^*\) with \(\G\B^*\) (Table 9.5), we obtain \begin{equation} K\subs{s} = \left(\frac{\G\mbB}{\G\B^*}\right) \g_{\pm}^\nu \left( \frac{m_+}{m\st} \right)^{\nu_+} \left( \frac{m_-}{m\st} \right)^{\nu_-} \tag{12.5.25} \end{equation} where \(\nu=\nu_+ + \nu_-\) is the total number of ions per formula unit. \(\g_{\pm}\) is the mean ionic activity coefficient of the dissolved salt in the saturated solution, and the molalities \(m_+\) and \(m_-\) refer to the ions M\(^{z_+}\) and X\(^{z_-}\) in this solution.

    The first factor on the right side of Eq. 12.5.25, the proper quotient of pressure factors for the reaction B\(^*\)(s)\(\ra\)B(sln), will be denoted \(\G\subs{r}\) (the subscript “r” stands for reaction). The value of \(\G\subs{r}\) is exactly \(1\) if the system is at the standard pressure, and is otherwise approximately \(1\) unless the pressure is very high.

    If the aqueous solution is produced by allowing the salt to dissolve in pure water, or in a solution of a second solute containing no ions in common with the salt, then the ion molalities in the saturated solution are \(m_+=\nu_+m\B\) and \(m_-=\nu_-m\B\) where \(m\B\) is the solubility of the salt expressed as a molality. Under these conditions, Eq. 12.5.25 becomes \begin{gather} \s{ K\subs{s} = \G\subs{r} \g_{\pm}^{\nu} \left(\nu_+^{\nu_+}\nu_-^{\nu_-}\right) \left( \frac{m\B}{m\st} \right)^{\nu} } \tag{12.5.26} \cond{(no common ion)} \end{gather} We could also have obtained this equation by using the expression of Eq. 10.3.10 for \(a\mbB\).

    If the ionic strength of the saturated salt solution is sufficiently low (i.e., the solubility is sufficiently low), it may be practical to evaluate the solubility product with Eq. 12.5.26 and an estimate of \(\g_{\pm}\) from the Debye–Hückel limiting law (see Prob. 12.19). The most accurate method of measuring a solubility product, however, is through the standard cell potential of an appropriate galvanic cell (Sec. 14.3.3).

    Since \(K\subs{s}\) is a thermodynamic equilibrium constant that depends only on \(T\), and \(\G\subs{r}\) depends only on \(T\) and \(p\), Eq. 12.5.26 shows that any change in the solution composition at constant \(T\) and \(p\) that decreases \(\g_{\pm}\) must increase the solubility. For example, the solubility of a sparingly-soluble salt increases when a second salt, lacking a common ion, is dissolved in the solution; this is a salting-in effect.

    Equation 12.5.25 is a general equation that applies even if the solution saturated with one salt contains a second salt with a common ion. For instance, consider the sparingly-soluble salt M\(_{\nu_+}\)X\(_{\nu_-}\) in transfer equilibrium with a solution containing the more soluble salt M\(_{\nu'_+}\)Y\(_{\nu'_-}\) at molality \(m\C\). The common ion in this example is the cation M\(^{z_+}\). The expression for the solubility product is now \begin{gather} \s{ K\subs{s} = \G\subs{r} \g_{\pm}^{\nu} (\nu_+m\B+\nu'_+m\C)^{\nu_+}(\nu_-m\B)^{\nu_-}/(m\st)^{\nu} } \tag{12.5.27} \cond{(common cation)} \end{gather} where \(m\B\) again is the solubility of the sparingly-soluble salt, and \(m\C\) is the molality of the second salt. \(K\subs{s}\) and \(\G\subs{r}\) are constant if \(T\) and \(p\) do not change, so any increase in \(m\C\) at constant \(T\) and \(p\) must cause a decrease in the solubility \(m\B\). This is called the common ion effect.

    From the measured solubility of a salt in pure solvent, or in an electrolyte solution with a common cation, and a known value of \(K\subs{s}\), we can evaluate the mean ionic activity coefficient \(\g_{\pm}\) through Eq. 12.5.26 or 12.5.27. This procedure has the disadvantage of being limited to the value of \(m\B\) existing in the saturated solution.

    We find the temperature dependence of \(K\subs{s}\) by applying Eq. 12.1.12: \begin{equation} \frac{\dif\ln K\subs{s}}{\dif T} = \frac{\Delsub{sol,B}H\st}{RT^2} \tag{12.5.28} \end{equation} At the standard pressure, \(\Delsub{sol,B}H\st\) is the same as the molar enthalpy of solution at infinite dilution, \(\Delsub{sol,B}H^{\infty}\).

    Contributors

    • Was this article helpful?