The Debye-Huckel theory

Result

The Debye-Huckel theory

T-102Home CU-307Threads energy · equilibrium

Statement

Activity coefficients in dilute electrolyte solutions.

Why it matters

Every equilibrium constant and every Nernst-equation calculation this network has treated so far has quietly used concentration as a stand-in for the thermodynamically correct quantity, activity. That substitution is an excellent approximation for dilute, uncharged species, but it fails increasingly badly for ions as concentration rises, because a charged ion's long-range electrostatic interactions with the surrounding cloud of other ions cannot be ignored the way a neutral solute's can. The Debye–Huckel theory supplies the correction: a way to compute an ion's activity coefficient from first principles, restoring the thermodynamic rigor that electrical-double-layer, butler-volmer, and kohlrausch-conductivity all implicitly assume is available whenever they invoke an ion's true chemical potential rather than its bare concentration.

It is also the clearest worked example, anywhere in physical chemistry, of ionic strength as a genuinely distinct and useful quantity from raw concentration — a multiply-charged ion at low concentration can dominate a solution's ionic strength just as much as a singly-charged ion at much higher concentration, a fact with real consequences for solubility, reaction rates, and electrode behaviour alike.

Hypotheses

Ions are treated as point charges immersed in a continuous, structureless dielectric medium (the solvent), interacting purely via long-range Coulomb electrostatics.This ignores the solvent's actual molecular structure (hydrogen bonding, finite molecular size, local ordering around an ion) entirely, replacing it with a single bulk dielectric constant; the approximation is reasonable at the low concentrations the theory targets, where ion–ion separations are large compared to solvent molecular dimensions. The solution is dilute enough that ion–ion correlations can be treated by linearising the Poisson–Boltzmann equation around the bulk ion concentrations.This linearisation is what makes the mathematics of the limiting law tractable in closed form; it is valid only when the electrostatic potential energy of a typical ion–ion interaction is small compared to \(kT\), which restricts the original Debye–Huckel limiting law to genuinely dilute solutions (Result's Scope), with the extended law relaxing this only partially by adding an empirical ion-size correction. Each ion is surrounded, on time-average, by a diffuse spherical "ionic atmosphere" of net opposite charge, rather than any specific, fixed arrangement of neighbouring ions.This statistical, time-averaged picture (rather than a fixed ion-pair or lattice-like arrangement) is what allows the theory to compute a single, well-defined activity coefficient for an ionic species at a given concentration, instead of requiring detailed knowledge of instantaneous ion positions.

Proof

I = \tfrac12\sum_i c_i z_i^2

Define the ionic strength \(I\) as one-half the sum, over every ionic species present, of its molar concentration \(c_i\) weighted by the square of its charge \(z_i\); squaring the charge (rather than using it linearly) is essential because it is exactly what the underlying electrostatic (Coulombic) interaction energy between ions depends on, so a doubly charged ion contributes four times as strongly to \(I\) as a singly charged ion at the same concentration. B

\log_{10}\gamma_{\pm} = -A\,|z_+z_-|\,\sqrt{I}

Solving the linearised Poisson–Boltzmann equation (Hypotheses) for the electrostatic free energy of a central ion surrounded by its diffuse ionic atmosphere yields the Debye–Huckel limiting law: the mean ionic activity coefficient \(\gamma_\pm\) falls below unity (activity below concentration) as ionic strength rises, with \(A\) a temperature- and solvent-dependent constant (\(A\approx0.509\) for water at \(25^\circ\text{C}\)) and the characteristic \(\sqrt{I}\) dependence a direct signature of the underlying electrostatic screening physics. B

\log_{10}\gamma_{\pm} = \frac{-A\,|z_+z_-|\,\sqrt{I}}{1+Ba\sqrt{I}}

The extended Debye–Huckel law introduces an empirical ion-size parameter \(a\) (the distance of closest approach between ions) in the denominator, correcting for the finite size of real ions, which the point-charge limiting law of Step 2 ignores entirely; this extends the law's useful range to somewhat higher ionic strengths than the bare limiting law alone can reach. B

a_i = \gamma_i\,c_i \quad\text{(activity as concentration corrected by the computed activity coefficient)}

Once \(\gamma_\pm\) (or the single-ion \(\gamma_i\), not independently measurable but usable within the mean ionic convention) is computed from Step 2 or 3, it converts a raw concentration into the thermodynamically correct activity, which is the quantity that properly belongs in an equilibrium constant expression, a Nernst equation, or any other rigorous thermodynamic relation involving that ion. A

Result

\log_{10}\gamma_{\pm} = -A\,|z_+z_-|\,\sqrt{I}, \qquad I=\tfrac12\textstyle\sum_i c_i z_i^2

Reading. An ion's activity coefficient falls below one as ionic strength rises, roughly as \(\sqrt{I}\), because each ion is progressively more effectively electrostatically screened by the diffuse cloud of opposite-charge ions surrounding it as the solution becomes more concentrated in charged species.

Scope. The bare limiting law (Step 2) is reliable only for \(I\) below roughly \(0.01\,\text{mol/L}\); the extended law (Step 3) pushes this to perhaps \(0.1\,\text{mol/L}\); beyond that, more elaborate, largely empirical models (outside this result's scope) are needed, since the underlying dilute, linearised-electrostatics approximation of the Hypotheses breaks down.

Corollaries & converses

kohlrausch-conductivity's own concentration-dependence of molar conductivity shares the identical physical origin: the same ionic-atmosphere screening that lowers \(\gamma_\pm\) here also retards an ion's mobility under an applied field, which is why both effects scale characteristically with \(\sqrt{c}\) at low concentration.
electrical-double-layer's treatment of the diffuse layer near a charged electrode surface uses the same linearised Poisson–Boltzmann framework (Hypotheses) as this bulk-solution theory, differing mainly in geometry (a planar or spherical charged surface rather than a point ion).
butler-volmer's electrode-kinetics expressions are, strictly, written in terms of activities rather than bare concentrations; the Debye–Huckel correction of Step 4 is what justifies treating concentration and activity interchangeably in the dilute-solution limit those expressions are usually applied in.

Fails without

Apply the bare limiting law (Step 2) well outside its dilute-solution range: the underlying linearised Poisson–Boltzmann approximation (Hypotheses) breaks down as ionic strength rises, and the predicted \(\gamma_\pm\) diverges increasingly from the true, experimentally measured activity coefficient.
Ignore finite ion size and use the point-charge limiting law at moderate ionic strength, where the extended law's correction actually matters: the bare law systematically over-predicts how far \(\gamma_\pm\) falls below unity, since it takes no account of how closely real, finite-sized ions can actually approach one another.

Common errors

Computing ionic strength by simply summing concentrations, forgetting the \(z_i^2\) weighting of Step 1 — a solution of a 2:2 electrolyte has a substantially higher ionic strength than a 1:1 electrolyte at the identical molar concentration.
Applying the bare limiting law (Step 2) at concentrations well above its dilute-solution range of validity, where the extended law (Step 3), or an even more elaborate model, is required instead.
Treating \(\gamma_\pm\) as always representing a "correction toward higher" activity; in the dilute regime this theory covers, \(\gamma_\pm\) is always less than one, reflecting net stabilisation of an ion by its oppositely-charged atmosphere, which lowers its effective activity below its raw concentration.
Forgetting that \(\gamma_\pm\) is a measurable mean ionic quantity, while individual single-ion activity coefficients are not independently measurable at all (no experiment can isolate one ion's chemical potential from its necessarily accompanying counter-ion).

Discussion

Peter Debye and Erich Huckel published their theory in 1923, providing the first rigorous statistical-mechanical explanation for a body of empirical deviations from ideal solution behaviour in electrolytes that had puzzled chemists since the recognition, decades earlier, that strong electrolytes are fully dissociated in solution and yet still show markedly non-ideal colligative and conductivity behaviour even at quite low concentrations.

A genuinely subtle point is that the theory's \(\sqrt{I}\) dependence, rather than a simple linear dependence on concentration, is itself direct evidence for the long-range nature of Coulombic interactions: short-range interactions between solute particles typically produce activity coefficient corrections that are approximately linear in concentration at low concentration, and the characteristically different square-root form is a fingerprint of the electrostatic, ionic-atmosphere mechanism specifically.

Common misconception: that activity coefficients only matter for "concentrated" solutions and can be safely ignored for anything a student would call dilute. In fact, measurable deviations from \(\gamma_\pm=1\) appear even in millimolar electrolyte solutions for multiply charged ions, precisely because of the \(z^2\) weighting in ionic strength (Step 1); the theory's own bare limiting law is calibrated for exactly this low-concentration regime, not for concentrated solutions.

Worked examples

\text{Compute } I \text{ for } 0.010\,\text{mol/L MgCl}_2 (aq)

MgCl2 dissociates fully into \(\text{Mg}^{2+}\) (\(c=0.010\,\text{mol/L}\), \(z=+2\)) and \(2\times0.010=0.020\,\text{mol/L}\) \(\text{Cl}^-\) (\(z=-1\)). Applying Step 1: \(I=\tfrac12\big[(0.010)(2)^2+(0.020)(1)^2\big]=\tfrac12[0.040+0.020]=0.030\,\text{mol/L}\). A

\log_{10}\gamma_{\pm} = -(0.509)(2)(1)\sqrt{0.030} = -0.176 \ \Rightarrow\ \gamma_{\pm}\approx0.667

Using \(A=0.509\) for water at \(25^\circ\text{C}\) and \(|z_+z_-|=|(+2)(-1)|=2\) in Step 2's limiting law, the computed mean ionic activity coefficient of roughly \(0.67\) shows the ions' true thermodynamic activity is noticeably (about a third) below their raw molar concentration, even at this fairly modest ionic strength — illustrating precisely why the \(z^2\) weighting of ionic strength matters so much for a 2:1 electrolyte like MgCl2. A

I=0.030\,\text{mol/L}, \quad \gamma_{\pm}\approx0.67 \ \text{ for } 0.010\,\text{mol/L MgCl}_2

Reading. A modest concentration of a 2:1 electrolyte already produces a substantial ionic strength, owing to the \(z^2\) weighting, and a correspondingly significant departure of activity from raw concentration.

Scope. The identical two-step procedure (compute \(I\), then apply the limiting or extended law) applies to any electrolyte concentration within the theory's dilute-solution range of validity.

Problems

Compute the ionic strength of a solution that is \(0.020\,\text{mol/L}\) in NaCl and simultaneously \(0.010\,\text{mol/L}\) in Na2SO4.
Solution
Species and concentrations: \(\text{Na}^+\) from both salts, \(0.020+2(0.010)=0.040\,\text{mol/L}\) (\(z=+1\)); \(\text{Cl}^-\), \(0.020\,\text{mol/L}\) (\(z=-1\)); \(\text{SO}_4^{2-}\), \(0.010\,\text{mol/L}\) (\(z=-2\)). \(I=\tfrac12\big[(0.040)(1)+(0.020)(1)+(0.010)(4)\big]=\tfrac12[0.040+0.020+0.040]=0.050\,\text{mol/L}\).
Without recomputing numerically, state whether \(\gamma_{\pm}\) for a 1:1 electrolyte (e.g. NaCl) or a 2:2 electrolyte (e.g. MgSO4) deviates further from unity at the same molar concentration, and justify your answer using the Result.
Solution
The 2:2 electrolyte deviates further from unity. Ionic strength scales with \(z^2\) (Step 1), so MgSO4 at a given molar concentration has a substantially higher \(I\) than NaCl at the same concentration; and \(\log_{10}\gamma_\pm\) itself also scales directly with the product \(|z_+z_-|\) (Step 2), which is \(4\) for a 2:2 electrolyte versus \(1\) for a 1:1 electrolyte. Both factors compound in the same direction, so MgSO4's activity coefficient falls far below unity at concentrations where NaCl's remains close to it.
Explain, using the Hypotheses, why the Debye–Huckel limiting law is expected to become progressively less accurate as ionic strength increases, even before invoking the extended law's ion-size correction.
Solution
The limiting law's derivation rests on linearising the Poisson–Boltzmann equation, an approximation valid only when the electrostatic interaction energy between a given ion and its surrounding atmosphere is small compared to \(kT\) (Hypotheses). As ionic strength (and hence the typical electrostatic interaction energy between neighbouring ions) increases, this small-perturbation assumption becomes progressively less valid, so the linearised solution increasingly departs from the true, nonlinear behaviour — independent of, and prior to, the additional finite-ion-size correction the extended law separately introduces.