25.6: The Debye-Hückel Theory
- Page ID
- 14518
\( \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}\)As ionic solutions are very common in chemistry, having to measure all activity coefficients (\(γ_{\pm}\)) for all possible solute-solvent combinations is a pretty daunting task, even though in times past extensive tabulation has taken place. We should be grateful for the rich legacy that our predecessors have left us in this respect (it would be hard to get any funding to do such tedious work today). Of course, it would be very desirable to be able to calculate \(γ_{\pm}\) values from first principles or if that fails by semi-empirical means. Fortunately, considerable progress has been made on this front as well. We can only scratch the surface of that topic in this course and will briefly discuss the simplest approach due to Debye and Hückel
Debye and Hückel came up with a theoretical expression that makes is possible to predict mean ionic activity coefficients as sufficiently dilute concentrations. The theory considers the vicinity of each ion as an atmosphere-like cloud of charges of opposite sign that cancels out the charge of the central ion (Figure \(\PageIndex{1}\)). From a distance the cloud looks neutral. The quantity \(1/κ\) is a measure for the size of this cloud and \(\kappa\) is the Debye-length. Its size depends on the concentration of all other ions.
Ionic Strength
To take the effect from all other ions into account, it is useful to define the ionic strength (\(I\)) as:
\[I =\dfrac{1}{2} \sum m_iz_i^2 \nonumber \]
where \(m_i\) is the molality of ion \(i\) and \(z_i\) is its charge coefficient. Note that highly charged ions (e.g. \(z=3+\)) contribute strongly (nine times more than +1 ions), but the formula is linear in the molality. Using the ionic strength the Debye-length becomes:
\[κ^2 = constant \, I \nonumber \]
The constant contains \(kT\) and \(ε_rε_o\) in the denominator and the number of Avogadro \(N_A\) and the square of the charge of the electron \(e\) in the numerator:
\[constant= 2000 \dfrac{e^2N_A}{ε_rε_okT} \nonumber \]
The Debye length and the logarithmic mean ionic activity coefficient are proportional:
\[\ln γ_{\pm} \propto κ \nonumber \]
Again there are a number of factors in the proportionality constant:
\[\ln γ_± = -|q_+q_-| \dfrac{κ}{8πε_rε_okT} \nonumber \]
The factors \(ε_r\) and \(ε_o\) are the relative permittivity of the medium and the permittivity of vacuum, respectively. Note that the factor \(8πε_rε_o\) is specific to the SI system of units. In cgs units the expression would look different, because the permittivities are defined differently in that system
If there is only one salt being dissolved, the ionic strength depends linearly on its concentration, the Debye length κ and \(\ln γ_±\), therefore, go as the square root of concentration (usually molality):
\[\ln \gamma_{\pm} \propto \sqrt{m} \nonumber \]
If there are other ions present the ionic strength involves all of them. This fact is sometimes used to keep ionic strength constant while changing the concentration of one particular ion. Say we wish to lower the concentration of Cu2+ in a redox reaction but we want to keep activity coefficients the same as much as possible. We could then replace it by an ion of the same charge say Zn2+ that does not partake in the reaction. A good way to do that is to dilute the copper solution with a zinc solution of the same concentration instead of with just solvent. The mean activity coefficient is given by the logarithm of this quantity as follows:
\[\log _{10}\gamma _{\pm }=-Az_{j}^{2}{\frac {\sqrt {I}}{1+Ba_{0}{\sqrt {I}}}} \label{DH}\]
with:
\[ A={\frac {e^{2}B}{2.303\times 8\pi \epsilon _{0}\epsilon _{r}k_{\rm {B}}T}}\]
\[B=\left({\frac {2e^{2}N}{\epsilon _{0}\epsilon _{r}k_{\rm {B}}T}}\right)^{1/2}\]
where \(I\) is the ionic strength and \(a_0\) is a parameter that represents the distance of closest approach of ions. For aqueous solutions at 25 °C \(A = 0.51\, mol^{−1}/2dm^{3/2}\) and \(B = 3.29 nm^{−1}mol^{−1}/2dm^{3/2}\).
Unfortunately this theory only works at very low concentrations and is therefore also known as the Debye limiting law (Figure \(\PageIndex{2}\)). There are a number of refinements that aim at extending the range of validity of the theory to be able to work at somewhat higher concentrations. These are discussed in the next section.
The most significant aspect of Equation \ref{DH} is the prediction that the mean activity coefficient is a function of ionic strength rather than the electrolyte concentration. For very low values of the ionic strength the value of the denominator in the expression above becomes nearly equal to one. In this situation the mean activity coefficient is proportional to the square root of the ionic strength.
When a solid is formed by a reaction from solution it is sometimes possible that it remains dispersed as very small particles in the solvent. The sizes typically range in the nanometers This is why it has become fashionable to call them nanoparticles, although they had been known as colloidal particles since the mid nineteenth century. They are smaller than the wavelength of the visible reason. This causes liquids that contain them to remain clear, although they can at times be beautifully colored. A good example is the reduction of AuCl4- with citrate to metallic gold. This produces clear wine red solutions, even at tiny gold concentrations.
\[\ce{2n AuCl4(aq)^{-} + n \,citrate^{3-}(aq) + 2n\, H_2O(l) \rightarrow} \\[4pt] \ce{2n\, Au(colloid) + 3n\, CH2O(aq) + 3n\, CO2(g) + 8n \,Cl^{-}(aq) + 3n\, H^{+}(aq)} \nonumber \]
The reason the gold does not precipitate completely is typically that the nanoparticle (AuNP) formed during the reaction are charged by the attachment of some of the ionic species in solution to its surface. This results in an charged particle with an atmosphere with a certain Debye length around it (Figure \(\PageIndex{3}\)). This charged cloud prevents the particle form coalescing with other particles by electrostatic repulsion.
Such a system is called a colloid. Of course these systems are metastable. Often they have a pretty small threshold to crashing to a real precipitate under influence of the strong van der Waals interactions that the particles experience once they manage to get in close contact. Under the right conditions colloids can survive for a long time. Some gold colloids prepared by Faraday in the 1850's are still stable today.
It will be clear from the above that addition of a salt -particularly containing highly charged ions like 3+ or 3-- may destabilize the colloid because the ionic strength will changed drastically and this will affect the Debye length.