Deriving the Debye–Hückel Limiting Law

Introduction

The Debye–Hückel limiting law is more than an empirical straight line. It follows from a model in which mobile ions respond to the potential of a central ion while the solvent acts as a dielectric continuum. The derivation explains the characteristic square root of ionic strength and also reveals why the law is restricted to sufficiently dilute solutions.

Core explanation

Choose an ion of charge zi e at the origin. A mean electrostatic potential ψ(r) changes the local number density of each other ionic species j. A Boltzmann factor gives nj(r) = nj,bulk exp[−zj e ψ(r)/(kBT)] in a simple mean-field picture. The exponential expresses the competition between electrostatic energy and thermal motion. It assumes point-like ions in a uniform dielectric and ignores specific short-range association.

Poisson's equation relates potential curvature to charge density: ∇²ψ = −ρ/ε. Substitute the Boltzmann distributions into ρ = Σj zj e nj(r), together with the central charge. Bulk electroneutrality ensures Σj zj nj,bulk = 0 away from the central ion. If zj e ψ is much smaller than kBT in the region where the approximation is used, expand exp(−x) ≈ 1 − x. The constant terms cancel by bulk electroneutrality, leaving a charge response proportional to −ψ Σj nj,bulk zj².

That linearized equation contains κ² = e²Σj nj,bulk zj²/(εkBT). Converting number densities to molar concentrations gives κ² = F²Σj cj zj²/(εRT), or 2F²I/(εRT) with a concentration-scale ionic strength in SI units. The potential outside the central ion then has a screened form proportional to exp(−κr)/r. Its characteristic decay length is κ−1, the Debye length.

To obtain an activity coefficient, calculate the reversible electrical work associated with building up the central ion's charge in the potential created by the surrounding atmosphere. In the point-ion, infinite-dilution limit, the excess chemical potential is proportional to −zi²κ. Since κ is proportional to √I, ln γi is proportional to −zi²√I. Expressing the logarithm in base ten gives log10 γi = −A zi²√I when the units and solvent-dependent constant A are specified. For a neutral salt, the stoichiometric mean becomes log10 γ± = −A z+z− √I in the simple limiting law.

The negative sign reflects stabilization by the countercharge cloud in the dilute model. The square of charge comes from stronger electrostatic coupling for higher valence; the square root comes from the screening parameter. The derivation does not predict detailed ion-specific behavior at moderate concentration, because its point-ion, continuum-solvent and linearized assumptions begin to fail.

Step-by-step reasoning

Start from a central charge and the Boltzmann response of mobile ions. Insert the distribution into Poisson's equation, use bulk electroneutrality and linearize for small dimensionless potential. Identify κ² and connect it to ionic strength. Compute the atmosphere's charging contribution to chemical potential and divide by RT to obtain ln γ. State the limiting result and every approximation before using it numerically.

Visual explanation

Draw a flow chart: central-ion potential → Boltzmann redistribution → charge density → Poisson equation → linearized screened potential → atmosphere charging work → log γ ∝ −z²√I. Show a small graph of log γ versus √I that begins at zero and slopes downward near infinite dilution.

Real-world analogy

A disturbance in a crowd causes nearby people to redistribute, and their redistribution partly counteracts the original disturbance. A linear model works when the disturbance is small. The ionic theory is more precise because it couples a Boltzmann population response to an electrostatic field equation; at large disturbances the linear approximation fails.

Real-world example

In very dilute aqueous NaCl, adding a small amount of salt raises ionic strength and the limiting law predicts a mean coefficient below one. For equally dilute CaCl2, the product of ionic charge magnitudes is two rather than one, so the leading log-coefficient correction is larger at the same ionic strength. This comparison isolates charge dependence, not differences in specific hydration.

Why?

The atmosphere contains excess opposite charge, lowering the reversible work of inserting the central ion relative to an ideal noninteracting solution. The strength of this effect scales with central-ion charge squared, while the density of available screening charge sets κ ∝ √I. These two relationships produce the distinctive limiting-law form.

Common misconception

The limiting law is not an exact formula for all electrolyte concentrations. Its derivation assumes diffuse point ions, a dielectric continuum and linear response. Another mistake is to use total analytical salt concentration in place of ionic strength for a multivalent electrolyte.

Worked example

Question: Without calculating a numerical coefficient, compare the leading limiting-law correction for a monovalent ion and a divalent ion at the same very low ionic strength.

Reasoning: The individual-ion limiting expression has log10 γi proportional to −zi²√I. At fixed I, changing z from 1 to 2 multiplies z² from 1 to 4. Therefore the divalent ion's leading logarithmic deviation from ideality has four times the magnitude in this model. The actual coefficients may diverge from this relation outside the dilute limit.

Answer: The divalent ion has four times the leading magnitude of log10 γ correction.

Quick check

1. Which mathematical approximation turns the Poisson–Boltzmann equation into the linear Debye–Hückel form? Answer: Expand exp(−x) ≈ 1 − x when the ion's electrostatic energy in the local potential is small compared with kBT.

Exam focus

Present the derivation as linked assumptions and results rather than memorizing only the final equation. Show how bulk electroneutrality removes constant terms, how κ² contains Σcj zj², and why the charging work yields a negative z²√I correction. State the ionic-strength scale used.

Advanced insight

The limiting law becomes exact only as ionic strength approaches zero within its continuum framework. Single-ion coefficients still rely on a convention; neutral-salt mean coefficients connect more directly to thermodynamic measurements. A finite-size model changes higher-order behavior while retaining the leading limiting slope.

Summary

Debye–Hückel theory linearizes the Poisson–Boltzmann response of a dilute ionic atmosphere. The resulting screened potential has κ proportional to √I, and charging an ion in its atmosphere gives log γ proportional to −z²√I. The law's assumptions explain both its usefulness near infinite dilution and its failure at higher concentration.

Practice questions

1. What does bulk electroneutrality remove in the linearized charge-density expansion? Answer: The zeroth-order constant term in the sum of ionic charges. 2. What sets the screening parameter κ in the dilute model? Answer: Solvent permittivity, temperature and the squared-charge-weighted ion concentrations. 3. Why is the correction larger for higher-valence ions? Answer: The leading excess chemical potential scales with charge squared. 4. Does the limiting law include specific ion pairing? Answer: No. It models long-range diffuse electrostatics of effectively point-like ions.