Debye-Huckel and Activity Coefficients

A pH meter reads 4.01 in a standard buffer. But the "concentration" of H⁺ implied by that number is not what a burette delivered. Real ions in solution feel each other. They experience electrostatic tugs from every neighboring charge. So the effective, thermodynamically active concentration of an ion is always somewhat less than its formal molarity. The gap between what you weighed out and what the ion actually "does" in an equilibrium expression is captured by a single multiplicative factor: the activity coefficient, γ. Debye-Hückel theory was born in 1923 from a deceptively simple electrostatic model. It gives us a working formula for γ. Any analytical chemist can compute this in seconds. It remains the backbone of every serious equilibrium calculation done today. This matters because analytical chemistry runs on equilibrium constants — $K_a$ , $K_{sp}$ , $K_f$ , $\beta_n$ . Every one of these is, strictly, a ratio of activities, not concentrations. Ignore activity coefficients, and several calculations can go badly wrong. A solubility calculation in seawater can be off by tens of percent. So can a pH calculation in a buffered serum matrix. So can a titration endpoint in the presence of a swamping electrolyte. The chemist who assumes ideality in a solution of ionic strength 0.1 M or higher is not doing approximate chemistry. They are doing wrong chemistry, silently. The big picture is this: every ion in solution sits inside a statistically averaged "ionic atmosphere" of opposite charge. This atmosphere is built up by thermal motion competing against Coulomb attraction. Debye and Hückel linearized the Poisson-Boltzmann equation describing that atmosphere. They extracted the extra electrostatic free energy per ion. Exponentiating that free energy gives the activity coefficient. The resulting formula depends on only three things. These are the ion's charge, its effective hydrated size, and the ionic strength of the whole solution. The extended (Debye-Hückel-Güntelberg-style) form pushes the useful range further. It extends from dilute solutions (< 0.001 M) out to the 0.1 M regime that dominates real laboratory work. This article builds the extended Debye-Hückel equation from the physical picture of the ionic atmosphere. It derives both the limiting law and the extended form with explicit algebra. It works three full numerical examples. These are a single-ion activity coefficient at low ionic strength, a realistic buffer/serum calculation, and a competition-style problem requiring iterative ionic-strength determination. It closes with the traps that catch students. These include charge-squared dependence, the difference between $\gamma_{\pm}$ and single-ion $\gamma_i$ , and the upper validity limit. It previews how these coefficients feed directly into Thermodynamic vs Concentration Equilibrium Constants and Ionic Strength Effects on Solubility and pH.

Key takeaways

Activity $a_i = \gamma_i [i]$ ; ions are never perfectly ideal once ionic strength exceeds ~10⁻⁴ M.
The extended Debye-Hückel equation gives $\log\gamma_i$ from charge $z_i$ , ionic strength $\mu$ , and hydrated-ion size $\alpha_i$ (in picometers via the parameter table).
$\gamma_i$ depends on $z_i^2$ , so multiply-charged ions (Fe³⁺, PO₄³⁻) deviate from ideality far faster than singly charged ions.
The equation is reliable up to $\mu \approx 0.1$ M; beyond that, specific-ion (Pitzer/Davies-type) corrections are needed.
Neutral molecules have $\gamma \approx 1$ at all but very high ionic strength — activity corrections apply to charged species.
Activity coefficients are the bridge between the thermodynamic constant $K$ (true constant) and the concentration-based $K_c$ used at the bench.

Key idea: every ion's activity coefficient falls below 1 as ionic strength rises, following $-\log\gamma_i \propto z_i^2\sqrt{\mu}/(1+\text{size term})$ , because each ion is electrostatically screened by an oppositely charged "ionic atmosphere" that lowers its effective chemical potential.

Where this lives today

Every certified pH buffer (NIST/NBS traceable standards) is manufactured and assigned its pH value using Debye-Hückel-based (specifically, Bates-Guggenheim convention) single-ion activity coefficients for chloride, since pH is operationally defined via a cell containing Cl⁻ — the entire international pH scale rests on this formalism. Clinical chemistry laboratories correct ion-selective electrode readings for sodium, potassium, and calcium in blood serum using activity coefficients computed from extended or Davies-type equations, because serum's ionic strength (~0.15 M, dominated by NaCl) is far from the ideal-dilute limit assumed by naive Nernst-equation calculations. Environmental and geochemical modeling software (PHREEQC, MINTEQA2, Visual MINTEQ) computes speciation and mineral saturation indices in natural waters — freshwater, seawater, acid mine drainage — using Debye-Hückel or Davies-equation activity corrections as the default engine for ionic strengths up to about 0.5–1 M, switching to Pitzer parameters for brines and seawater beyond that. In pharmaceutical formulation, the solubility and stability of ionic active pharmaceutical ingredients in buffered, salt-containing injectable solutions are predicted using exactly the kind of ionic-strength-dependent $K_{sp}'$ calculation performed in Example 3, since drug salts must remain soluble (or controllably precipitate) across physiological ionic strengths. Analytical method development for potentiometric titrations and ion-selective electrode calibration curves routinely uses total ionic strength adjustment buffers (TISAB) specifically to fix μ at a known value, so that a single activity coefficient correction (rather than a shifting one) applies throughout a calibration — a direct practical descendant of the extended equation's μ-dependence. Finally, every serious equilibrium calculation distinguishing thermodynamic $K$ from concentration-based $K_c$ — the subject taken up next in Thermodynamic vs Concentration Equilibrium Constants — uses exactly the machinery built here, and the same activity corrections drive the shifts in solubility and pH catalogued in Ionic Strength Effects on Solubility and pH.

Ask the tutor: How do I decide, for a specific real solution, whether to use the simple limiting law, the extended Debye-Hückel equation, or the Davies equation?

The physical picture

Picture a Cl⁻ ion sitting in a beaker of 0.05 M NaCl. It is not alone: within a few tenths of a nanometer, the time-averaged charge density around it is slightly enriched in Na⁺ and slightly depleted in Cl⁻, simply because opposite charges attract and like charges repel, moderated by thermal (Brownian) randomization. This time-averaged, statistically smeared cloud of net opposite charge is the ionic atmosphere. It has a characteristic radius — the Debye length, $1/\kappa$ — inversely related to $\sqrt{\mu}$ : the more concentrated the electrolyte, the tighter and more effective the screening. The consequence for thermodynamics is that the central ion's electrostatic potential energy is lowered relative to an isolated ion in vacuum (or relative to the hypothetical ideal solution used as reference state). Lowered energy means lowered chemical potential, which by definition means an activity coefficient less than unity. Debye and Hückel calculated this energy by solving a linearized Poisson-Boltzmann equation for the potential around a central ion, then used the Debye charging process (a Gedanken experiment of slowly building up the ion's charge from zero while the atmosphere continuously re-equilibrates) to convert that potential into a free energy, and hence into $\ln\gamma_i$ . Two physical facts fall directly out of the model. First, the effect must scale as $z_i^2$ : reversing the sign of the charge doesn't change whether it's stabilized by its atmosphere (both cations and anions are stabilized, since each attracts an atmosphere of opposite sign), and doubling the charge quadruples the Coulomb interaction energy. Second, the ionic atmosphere cannot pack in arbitrarily close — real ions have finite size, and their innermost coordinated water molecules keep other ions at a minimum distance of closest approach, $\alpha$ (in picometers, tabulated per ion). This finite distance is what the "extended" term in the extended Debye-Hückel equation accounts for, and it's why extending the theory beyond the ultra-dilute limiting law requires an empirical per-ion size parameter.

Figure

The ionic atmosphere around a central Cl⁻ ion in NaCl solution

The size of the finite ion also matters mechanically: it sets the distance of closest approach between the central ion and the nearest atmosphere charge, which is why the extended equation includes $\alpha_i$ (often quoted in Ångströms or picometers) in its denominator. Small, highly charged, poorly hydrated ions (Al³⁺) behave very differently from large, singly charged, weakly hydrated ions (I⁻, ClO₄⁻) at the same ionic strength — a fact table values of $\alpha_i$ directly encode.

Figure

log γ vs √μ: limiting law (dashed) vs extended equation (solid), z = ±1

The core ideas

The Debye-Hückel limiting law, valid strictly only for $\mu \lesssim 0.001$ – $0.005$ M, is

\log \gamma_i = -0.51\, z_i^2 \sqrt{\mu}

where:

$\gamma_i$ — single-ion activity coefficient (dimensionless), defined by $a_i = \gamma_i [i]$
$z_i$ — signed charge number of ion $i$ (the sign drops out because it is squared)
$\mu$ — ionic strength, $\mu = \tfrac{1}{2}\sum_j c_j z_j^2$ , in mol/L
$0.51$ — the Debye-Hückel constant $A$ for water at 25 °C (units mol⁻¹/² L¹/²; it depends on temperature and solvent dielectric constant)

Reading: the log of the activity coefficient falls linearly with $\sqrt{\mu}$ , faster for higher-charge ions because of the $z^2$ dependence — a triply charged ion deviates from ideality nine times faster (in log space) than a singly charged one at the same $\mu$ . The extended Debye-Hückel equation, the workhorse for $0.001 \lesssim \mu \lesssim 0.1$ M, adds the finite-ion-size correction:

\log \gamma_i = \frac{-0.51\, z_i^2 \sqrt{\mu}}{1 + \dfrac{\alpha_i \sqrt{\mu}}{305}}

where:

$\alpha_i$ — the effective hydrated-ion-size parameter, in picometers (pm); tabulated per ion (e.g., H⁺: 900 pm, Mg²⁺: 800 pm, Cl⁻: 300 pm, K⁺: 300 pm)
$305$ — a collected numerical constant equal to $1/(B \times 10^{-8}\,\text{cm/pm conversion factor})$ where $B \approx 0.33 \times 10^8 \,\text{cm}^{-1}\text{mol}^{-1/2}\text{L}^{1/2}$ at 25 °C — it converts $\alpha_i$ in pm consistently with $\sqrt{\mu}$ in mol $^{1/2}$ L $^{-1/2}$
all other symbols as above

Reading: as $\mu \to 0$ , the denominator $\to 1$ and the extended equation collapses exactly onto the limiting law — the size correction only matters once the ionic atmosphere shrinks to a radius comparable to $\alpha_i$ itself, which happens as $\mu$ grows.

Derivation AFrom the ionic atmosphere to log γᵢ

Goal. Show how linearizing the Poisson-Boltzmann equation for the potential around a central ion leads to the limiting-law form $\log\gamma_i = -A z_i^2\sqrt{\mu}$ , then explain where the size parameter $\alpha_i$ enters to give the extended form.

Step 1 — Electrostatic potential from Poisson's equation. Around a central ion of charge $z_i e$ , the mean electrostatic potential $\psi(r)$ obeys the Poisson equation relating it to the local charge density $\rho(r)$ from all other ions:

\nabla^2 \psi = -\frac{\rho(r)}{\varepsilon_0 \varepsilon_r}

The charge density itself is built from a Boltzmann distribution of every ion type $j$ around the central ion, weighted by its own electrostatic energy $z_j e \psi(r)$ :

\rho(r) = \sum_j n_j^0 z_j e \, \exp\!\left(\frac{-z_j e \psi(r)}{k_BT}\right)

This is the full nonlinear Poisson-Boltzmann equation — exact in principle, intractable in closed form.

Step 2 — Linearization (the "Debye-Hückel approximation"). For dilute solutions, $z_j e \psi(r) \ll k_BT$ almost everywhere (the electrostatic energy is small compared to thermal energy except very close to the ion). Expanding the exponential to first order and using electroneutrality ( $\sum_j n_j^0 z_j e = 0$ , which kills the zeroth-order term) gives a linear equation:

\nabla^2 \psi = \kappa^2 \psi, \qquad \kappa^2 = \frac{e^2}{\varepsilon_0\varepsilon_r k_BT}\sum_j n_j^0 z_j^2

The quantity $\kappa$ has units of inverse length; $1/\kappa$ is the Debye screening length, and $\kappa^2$ is directly proportional to ionic strength $\mu$ (since $\sum_j n_j^0 z_j^2 \propto \mu$ ). Solving this linear (Helmholtz-type) equation with the boundary condition that $\psi \to$ bare Coulomb potential as $r\to 0$ gives the screened Coulomb form:

\psi(r) = \frac{z_i e}{4\pi\varepsilon_0\varepsilon_r}\cdot\frac{e^{-\kappa r}}{r}

This is the potential of the central ion plus the polarization response of its atmosphere — an exponentially screened Coulomb potential rather than the bare $1/r$ .

Step 3 — Extract the atmosphere's contribution and charge it up. Subtracting the bare-ion potential from the total gives the potential due to the atmosphere alone, evaluated at the ion's own surface. Using the Debye charging process — imagining the central ion's charge built up continuously from 0 to $z_ie$ while the atmosphere stays equilibrated at every step — and integrating the resulting work done against the atmosphere's field gives the extra (non-ideal) electrostatic free energy per ion:

\Delta G_i^{el} = -\frac{z_i^2 e^2 \kappa}{8\pi\varepsilon_0\varepsilon_r}

This is negative (stabilizing), confirming both cations and anions are lowered in energy by their atmospheres, and scales as $z_i^2$ , confirming the charge-squared dependence anticipated physically.

Step 4 — Convert free energy to an activity coefficient. By definition $\Delta G_i^{el} = RT\ln\gamma_i$ (the excess chemical potential over the ideal reference state). Substituting $\kappa \propto \sqrt{\mu}$ explicitly and converting to base-10 log, collecting all physical constants (elementary charge, permittivity, Avogadro's number, $R$ , $T$ ) into a single numeric constant $A = 0.51\,\text{mol}^{-1/2}\text{L}^{1/2}$ for water at 25 °C yields exactly the limiting law:

\log\gamma_i = -A z_i^2 \sqrt{\mu}

Step 5 — Reintroducing finite ion size. Step 2's linearization implicitly treated the central ion as a point charge, letting the atmosphere approach arbitrarily close ( $r\to 0$ ). Real ions (with their hydration shells) exclude other ions within a distance of closest approach $\alpha_i$ . Repeating the boundary-value problem in Step 2 with the boundary condition applied at $r = \alpha_i$ instead of $r = 0$ modifies the exponential screening factor, and the free energy in Step 3 picks up an extra factor $(1+\kappa\alpha_i)^{-1}$ . Carrying this correction through Step 4, with the numeric collection $B\alpha_i\sqrt{\mu} = \alpha_i\sqrt{\mu}/305$ (using $\alpha_i$ in pm), gives the extended equation:

\log\gamma_i = \frac{-A z_i^2\sqrt{\mu}}{1+\alpha_i\sqrt{\mu}/305}

As $\alpha_i \to 0$ (point-charge limit) this reduces back exactly to Step 4's limiting law — confirming internal consistency.

Go deeper · Derivation B · 1 minMean ionic activity coefficient from single-ion values

Derivation BMean ionic activity coefficient from single-ion values

Goal. Show why experiments can only measure the geometric mean activity coefficient $\gamma_\pm$ of a salt, and connect it back to the single-ion $\gamma_i$ values the extended equation predicts.

Step 1 — Why single ions are unmeasurable alone. Any real solution must be electrically neutral; you cannot prepare a beaker containing only Na⁺ with no counter-ion. Consequently no experiment (electrode potential, freezing-point depression, vapor pressure) can isolate the chemical potential of Na⁺ alone from that of its accompanying anion — every measurable quantity is necessarily a combination of both.

Step 2 — Define the mean ionic activity coefficient. For a salt $M_{p}X_{q}$ dissociating into $p$ cations of charge $z_+$ and $q$ anions of charge $z_-$ , define the mean activity coefficient as the geometric mean weighted by stoichiometry:

\gamma_\pm = \left(\gamma_+^{\,p}\,\gamma_-^{\,q}\right)^{1/(p+q)}

This single number is what conductivity, EMF, and colligative-property experiments actually determine.

Step 3 — Substitute the Debye-Hückel expressions. Writing $\log\gamma_+ = -Az_+^2\sqrt{\mu}/(1+\alpha_+\sqrt{\mu}/305)$ and similarly for $\gamma_-$ , and taking the logarithm of Step 2's definition:

\log\gamma_\pm = \frac{p\log\gamma_+ + q\log\gamma_-}{p+q}

Substituting and using electroneutrality of the salt ( $pz_+ = qz_-$ in magnitude, so $pz_+^2 + qz_-^2 = (p+q)z_+|z_-|$ after simplification) collapses the weighted average to a single simplified relation using $z_+ z_-$ as $|z_+ z_-|$ :

\log\gamma_\pm = \frac{-A\,|z_+z_-|\sqrt{\mu}}{1+\bar\alpha\sqrt{\mu}/305}

where $\bar\alpha$ is an appropriately averaged size parameter (often just taken as the tabulated value for the "salt" itself, e.g., $\alpha = 400$ pm for NaCl).

Step 4 — Reconciling with single-ion coefficients used at the bench. Even though $\gamma_\pm$ is the only rigorously measurable quantity, analytical practice (titration curves, $K_{sp}$ calculations, pH definitions) requires individual ion activities. The practical convention — used throughout quantitative analysis and endorsed operationally by IUPAC for pH — is to compute $\gamma_i$ for each ion separately from the extended equation using its own tabulated $\alpha_i$ , accepting that this individual value carries a non-thermodynamic (convention-dependent) assumption. This is precisely the approximation used in every $K_{sp}$ and $K_a$ calculation done with individual ion activity coefficients, and it is what feeds directly into Thermodynamic vs Concentration Equilibrium Constants.

Worked examples

Example 1 — Foundational: single-ion activity coefficient of Mg²⁺ in a dilute solution. A solution is 3.00 × 10⁻³ M MgCl₂ with no other electrolyte. Find $\gamma_{\mathrm{Mg^{2+}}}$ using the extended Debye-Hückel equation ( $\alpha_{\mathrm{Mg^{2+}}} = 800$ pm from standard tables). Step 1 — Ionic strength. MgCl₂ gives 3.00 × 10⁻³ M Mg²⁺ ( $z=+2$ ) and 6.00 × 10⁻³ M Cl⁻ ( $z=-1$ ):

\mu = \tfrac12\left[(3.00\times10^{-3})(2)^2 + (6.00\times10^{-3})(1)^2\right] = \tfrac12(1.20\times10^{-2}+6.00\times10^{-3}) = 9.00\times10^{-3}\ \text{M}

Step 2 — √μ.

\sqrt{\mu} = \sqrt{9.00\times10^{-3}} = 0.0949

Step 3 — Numerator.

-A z^2\sqrt{\mu} = -(0.51)(2)^2(0.0949) = -(0.51)(4)(0.0949) = -0.1936

Step 4 — Denominator.

1 + \frac{\alpha\sqrt{\mu}}{305} = 1 + \frac{(800)(0.0949)}{305} = 1 + \frac{75.9}{305} = 1 + 0.2489 = 1.2489

Step 5 — Combine.

\log\gamma_{\mathrm{Mg^{2+}}} = \frac{-0.1936}{1.2489} = -0.1550

\gamma_{\mathrm{Mg^{2+}}} = 10^{-0.1550} = 0.700

So even in a modestly dilute 3 mM MgCl₂ solution, Mg²⁺'s activity is only 70.0% of its formal concentration — already a substantial correction, driven entirely by the $z^2=4$ factor. As a limiting-case check: had this been a singly charged ion of the same $\alpha$ and $\mu$ , $-\log\gamma = 0.51(1)(0.0949)/1.2489 = 0.0388$ , giving $\gamma = 0.915$ — a far smaller deviation, exactly as the $z^2$ scaling predicts (ratio of exponents $4:1$ matches the ratio of $\log\gamma$ magnitudes $0.155:0.0388\approx 4.0$ ).

Example 2 — Applied: pH of an acetate buffer in the presence of background KCl. A buffer is prepared to be 0.100 M acetic acid / 0.100 M sodium acetate, then 0.150 M KCl is added as a "swamping" electrolyte (common in polarography and ionic-strength control). Given $K_a$ (thermodynamic) $= 1.75\times10^{-5}$ for acetic acid, $\alpha_{\mathrm{H^+}} = 900$ pm, $\alpha_{\mathrm{OAc^-}} = 450$ pm, find the actual pH, accounting for activity effects. Step 1 — Ionic strength. Contributions: Na⁺ (0.100 M, z=1), OAc⁻ (0.100 M, z=1), K⁺ (0.150 M, z=1), Cl⁻ (0.150 M, z=1). Acetic acid itself is neutral and contributes nothing.

\mu = \tfrac12\left[(0.100)(1)+(0.100)(1)+(0.150)(1)+(0.150)(1)\right] = \tfrac12(0.500) = 0.250\ \text{M}

This is above the nominal 0.1 M comfort zone of the extended equation, but it is still the standard tool used in practice (the alternative, the Davies equation, is reserved for $\mu > 0.5$ M); we proceed with the extended equation as the working estimate, flagging the limitation explicitly. Step 2 — √μ.

\sqrt{0.250} = 0.5000

Step 3 — γ for H⁺ (z=1, α=900 pm).

\log\gamma_{H^+} = \frac{-(0.51)(1)(0.5000)}{1+(900)(0.5000)/305} = \frac{-0.2550}{1+1.475} = \frac{-0.2550}{2.475} = -0.1030

\gamma_{H^+} = 10^{-0.1030} = 0.789

Step 4 — γ for acetate (z=1, α=450 pm).

\log\gamma_{OAc^-} = \frac{-(0.51)(1)(0.5000)}{1+(450)(0.5000)/305} = \frac{-0.2550}{1+0.7377} = \frac{-0.2550}{1.7377} = -0.1467

\gamma_{OAc^-} = 10^{-0.1467} = 0.713

Step 5 — Concentration-based apparent equilibrium constant. Since $K_a = \dfrac{a_{H^+}a_{OAc^-}}{a_{HOAc}} = \dfrac{\gamma_{H^+}[H^+]\,\gamma_{OAc^-}[OAc^-]}{\gamma_{HOAc}[HOAc]}$ and $\gamma_{HOAc}\approx1$ (neutral species):

K_a = K_a' \,\gamma_{H^+}\gamma_{OAc^-} \quad\Rightarrow\quad K_a' = \frac{K_a}{\gamma_{H^+}\gamma_{OAc^-}} = \frac{1.75\times10^{-5}}{(0.789)(0.713)} = \frac{1.75\times10^{-5}}{0.5626} = 3.11\times10^{-5}

Step 6 — Apply Henderson-Hasselbalch to concentrations using $K_a'$ .

[H^+] = K_a' \times \frac{[HOAc]}{[OAc^-]} = 3.11\times10^{-5} \times \frac{0.100}{0.100} = 3.11\times10^{-5}\ \text{M}

Step 7 — Convert to true pH (defined via activity, not concentration).

\mathrm{pH} = -\log a_{H^+} = -\log(\gamma_{H^+}[H^+]) = -\log[(0.789)(3.11\times10^{-5})] = -\log(2.454\times10^{-5}) = 4.610

Compare with the naive ideal calculation ignoring activity entirely: $\mathrm{pH} = pK_a = -\log(1.75\times10^{-5}) = 4.757$ . The real, swamped-electrolyte buffer reads about 0.15 pH units lower than the "ideal" prediction — a difference easily seen on a calibrated pH meter and a standard illustration of why ionic strength control matters in quantitative analysis.

Example 3 — Exam/competition level: iterative ionic strength for a sparingly soluble salt with a common-ion background. Calculate the molar solubility of PbSO₄ ( $K_{sp}$ (thermodynamic) $= 1.6\times10^{-8}$ ) in a solution that is 0.0400 M Na₂SO₄, accounting fully for activity effects. Use $\alpha_{\mathrm{Pb^{2+}}} = 450$ pm, $\alpha_{\mathrm{SO_4^{2-}}} = 400$ pm. Step 1 — Set up with an unknown solubility $s$ . Let $s$ = mol/L of PbSO₄ dissolved. Then $[\mathrm{Pb^{2+}}] = s$ and $[\mathrm{SO_4^{2-}}] = 0.0400+s \approx 0.0400$ M (since $K_{sp}$ is tiny, $s\ll 0.0400$ ; verified at the end). Step 2 — First-pass ionic strength (ignoring the small Pb²⁺/extra SO₄²⁻ contribution from dissolution). Na₂SO₄ alone: Na⁺ at 0.0800 M (z=1), SO₄²⁻ at 0.0400 M (z=2):

\mu \approx \tfrac12\left[(0.0800)(1)^2+(0.0400)(2)^2\right] = \tfrac12(0.0800+0.1600) = 0.1200\ \text{M}

Step 3 — Activity coefficients at μ = 0.1200 M. $\sqrt{\mu}=0.3464$ . For Pb²⁺ (z=2, α=450 pm):

\log\gamma_{Pb^{2+}} = \frac{-(0.51)(4)(0.3464)}{1+(450)(0.3464)/305} = \frac{-0.7067}{1+0.5113} = \frac{-0.7067}{1.5113} = -0.4677

\gamma_{Pb^{2+}} = 10^{-0.4677} = 0.3406

For SO₄²⁻ (z=2, α=400 pm):

\log\gamma_{SO_4^{2-}} = \frac{-(0.51)(4)(0.3464)}{1+(400)(0.3464)/305} = \frac{-0.7067}{1+0.4545} = \frac{-0.7067}{1.4545} = -0.4860

\gamma_{SO_4^{2-}} = 10^{-0.4860} = 0.3266

Step 4 — Concentration-based conditional solubility product.

K_{sp}' = \frac{K_{sp}}{\gamma_{Pb^{2+}}\gamma_{SO_4^{2-}}} = \frac{1.6\times10^{-8}}{(0.3406)(0.3266)} = \frac{1.6\times10^{-8}}{0.1113} = 1.438\times10^{-7}

Step 5 — Solve for s.

K_{sp}' = [\mathrm{Pb^{2+}}][\mathrm{SO_4^{2-}}] = s(0.0400+s) \approx s(0.0400)

s = \frac{1.438\times10^{-7}}{0.0400} = 3.60\times10^{-6}\ \text{M}

Step 6 — Check self-consistency of the ionic-strength approximation. $s = 3.60\times10^{-6}$ M is indeed $\ll 0.0400$ M (four orders of magnitude smaller), so the assumption in Step 1 that $[\mathrm{SO_4^{2-}}]\approx 0.0400$ M, and the neglect of Pb²⁺'s tiny contribution to $\mu$ in Step 2, are both justified — no further iteration is needed. Comparison with the naive (no-activity) calculation: $s_{naive} = \sqrt{K_{sp}/0.0400} = \sqrt{1.6\times10^{-8}/0.0400} = \sqrt{4.0\times10^{-7}} = 6.32\times10^{-4}$ M — wait, this naive treatment (ignoring the common-ion background entirely and treating it as simple $K_{sp}=s^2$ ) is not the correct point of comparison; the fair "ideal" comparison keeping the common-ion suppression but assuming $\gamma=1$ gives $s_{ideal} = K_{sp}/0.0400 = 1.6\times10^{-8}/0.0400 = 4.00\times10^{-7}$ M. The real, activity-corrected solubility ( $3.60\times10^{-6}$ M) is nine times larger than the ideal prediction — the ionic atmosphere from the swamping Na₂SO₄ stabilizes both Pb²⁺ and SO₄²⁻ ions enough to pull noticeably more PbSO₄ into solution than naive $K_{sp}$ arithmetic would suggest. This "salt effect" (increased solubility with increased ionic strength, for salts of multiply-charged ions) is a classic and easily mismeasured phenomenon in quantitative analysis.

Limiting cases, scope & common traps

Regime / case	Behavior of the extended equation	Practical implication
$\mu \to 0$	Denominator $\to 1$ ; equation $\to$ limiting law $\log\gamma_i = -Az_i^2\sqrt\mu$	Correct low-concentration check; matches Step 4 of the first derivation exactly
$\mu > 0.1$ –0.5 M	Extended equation increasingly inaccurate; γ may even be predicted to keep falling when real γ turns upward	Switch to the Davies equation or specific-ion (Pitzer) parameters
$z_i = 0$ (neutral species, e.g. HOAc, H₂O)	$\log\gamma_i \approx 0$ regardless of μ (to good approximation)	Neutral-species activities can be set to 1 in most bench calculations
Very high charge ( $z=3$ , e.g. Fe³⁺, PO₄³⁻, Al³⁺)	$z^2=9$ — deviations from ideality appear even at μ as low as 0.001 M	Never assume ideality for trivalent ions except in extremely dilute standards
$\alpha_i \to \infty$ (very large hydrated ion, hypothetical)	Denominator dominated by size term; γ approaches 1 even at moderate μ	Large, weakly hydrated ions (e.g., large organic anions) are "buffered" against ionic-strength effects
Mixed electrolytes, same μ	Extended equation depends only on μ and z, not on which specific ions supply that μ	This "ionic strength principle" is itself an approximation; real specific-ion effects (Hofmeister-type) are not captured

The chemist behind it

Peter Debye trained as an engineer and physicist in Aachen and Munich, working under Arnold Sommerfeld before building an extraordinarily broad career spanning X-ray diffraction, dipole moments, specific heats of solids, and — with his junior collaborator Erich Hückel — the statistical-mechanical theory of strong electrolytes. In 1923, while at the ETH Zürich, Debye and Hückel tackled a problem that had stumped physical chemists for decades: why do strong electrolytes (fully dissociated salts) show conductivities and colligative properties that deviate systematically from the ideal, dilute-solution laws, even at concentrations where "chemical" incomplete dissociation seemed implausible? Their answer — that the deviation is purely electrostatic, arising from the ionic atmosphere rather than from any incomplete ionization — was a conceptual triumph, converting a messy empirical curve-fitting problem into a clean, derivable physical law. Debye later emigrated to the United States, working at Cornell University, and received the 1936 Nobel Prize in Chemistry for his broader work on molecular structure via dipole moments and electron/X-ray diffraction — the Debye-Hückel theory itself, while arguably his most cited single contribution to analytical and physical chemistry, was folded into a career recognized more for its breadth. Erich Hückel, the junior partner on the 1923 papers, went on to develop the Hückel molecular orbital method for conjugated π systems, an entirely separate and equally foundational tool in structural organic chemistry — a reminder that the same scientist can leave two independent, unrelated legacies.

Origins · 1 minPeter Debye and how this was discovered

Peter Debye

1884-1966 · ETH Zürich and Cornell University, era of early statistical mechanics of electrolytes

Debye's 1923 papers with Hückel treated the electrolyte solution as a solvable statistical-mechanical problem for the first time, replacing ad hoc empirical activity corrections with a derivable formula rooted in Coulomb's law and Boltzmann statistics — the same style of reasoning he later applied to X-ray scattering from gases and the dipole moments that carry his name (the debye unit).

History

- 1923 — Debye and Hückel publish "Zur Theorie der Elektrolyte," deriving the limiting law from the linearized Poisson-Boltzmann equation.

- 1926 — Erich Hückel and later collaborators extend the treatment with finite-ion-size corrections, giving the extended equation used today.

- 1936 — Peter Debye awarded the Nobel Prize in Chemistry for work on dipole moments and X-ray/electron diffraction in gases.

- 1938 — Carl Davies proposes an empirical additional linear term in μ, extending practical validity toward μ ≈ 0.5 M (the Davies equation).

- 1973 — Kenneth Pitzer introduces a virial-expansion framework for activity coefficients that remains accurate to several molar ionic strength, used in modern geochemical and industrial process modeling.

Why does activity, not concentration, belong in equilibrium constant expressions?

Thermodynamic equilibrium constants are rigorously derived from the condition that the sum of chemical potentials of products equals that of reactants at equilibrium, and chemical potential in real (non-ideal) solutions is defined in terms of activity, not concentration — $\mu_i = \mu_i^\circ + RT\ln a_i$ . Concentration only equals activity in the idealized limit of infinite dilution, where every ion is so far from its neighbors that electrostatic interactions vanish. In any real solution with non-negligible ionic strength, using concentration directly in a $K$ expression introduces a systematic error that grows with both ionic strength and the charges of the species involved, which is exactly what the activity coefficient corrects for.

Why is there a size parameter α in the extended equation but not in the simple limiting law?

The Debye-Hückel limiting law is derived by treating the central ion as a mathematical point charge, valid only when the ionic atmosphere's characteristic radius (the Debye length) is much larger than the ions' actual physical size — true only at very low ionic strength. As ionic strength increases, the atmosphere compresses toward the ion, and its finite size (plus its hydration shell, since analytical calculations use hydrated rather than crystallographic radii) starts to matter because other ions physically cannot approach closer than a minimum distance. The extended equation reintroduces this distance of closest approach, α, as a denominator correction, extending accuracy from roughly 0.001 M up to about 0.1 M.

Can I use the extended Debye-Hückel equation for seawater (μ ≈ 0.7 M)?

Not reliably. The extended equation's own derivation assumes the linearized Poisson-Boltzmann approximation, which breaks down once ionic strength exceeds roughly 0.1 M — the predicted γ often keeps decreasing monotonically with μ, while real activity coefficients for many salts actually turn around and increase at high ionic strength (the "salt effect" reversal). For seawater and brines, oceanographers and geochemists use the Pitzer ion-interaction model, a virial-type expansion with empirically fitted binary and ternary interaction parameters, which remains accurate to several molar ionic strength.

Do activity coefficients matter for neutral molecules like acetic acid or glucose?

Generally not much, at least at the ionic strengths typical of aqueous analytical chemistry. Neutral species have z = 0, so the Debye-Hückel electrostatic mechanism (which depends on z²) predicts essentially no deviation from γ = 1 due to ionic strength; deviations that do occur for neutrals (the "salting out/salting in" effects) arise from different, non-electrostatic mechanisms such as changes in solvent structure and are treated with separate empirical parameters (Setschenow coefficients), not Debye-Hückel theory.

How were the α (ion-size) parameters actually determined?

The α values in standard tables (Kielland's 1937 compilation is the classic reference) were back-calculated by fitting the extended Debye-Hückel equation to experimentally measured mean ionic activity coefficients (from EMF cells, freezing-point depression, and vapor-pressure osmometry) across a range of ionic strengths for many salts, then apportioning the fitted "salt size" between the two constituent ions using consistency across multiple salts sharing a common ion. The resulting α values correlate reasonably with hydrated ionic radii but are best regarded as empirically optimized fitting parameters rather than directly measured physical radii.

What's the practical difference between the extended Debye-Hückel and the Davies equation?

The Davies equation modifies the extended form by dropping the individual α parameter (setting an implicit average size) and adding an extra empirical linear term in μ: $\log\gamma_i = -Az_i^2\left[\frac{\sqrt\mu}{1+\sqrt\mu} - 0.3\mu\right]$ . This extra term lets Davies' equation track the experimentally observed upturn in γ at higher ionic strength, extending usable accuracy to roughly μ ≈ 0.5 M, at the cost of losing ion-specific size information — useful when you want one formula for a whole ionic-strength-adjusted matrix rather than per-ion size lookups.

Why does the acetate-buffer example show the pH is lower than the ideal pKa, rather than higher?

Because both H⁺ and acetate ion activity coefficients drop below 1 with added ionic strength, the product $\gamma_{H^+}\gamma_{OAc^-}$ in the denominator of the apparent constant $K_a'=K_a/(\gamma_{H^+}\gamma_{OAc^-})$ is less than 1, which makes $K_a'$ larger than the thermodynamic $K_a$ — a numerically larger apparent acid-dissociation constant at fixed formal concentrations means more dissociation, hence a lower (more acidic) equilibrium pH than the ideal $pK_a$ value would suggest.