# Debye–Hückel equation

The chemists Peter Debye and Erich Hückel noticed that solutions that contain ionic solutes do not behave ideally even at very low concentrations. So, while the concentration of the solutes is fundamental to the calculation of the dynamics of a solution, they theorized that an extra factor that they termed gamma is necessary to the calculation of the activity coefficients of the solution. Hence they developed the **Debye–Hückel equation** and **Debye–Hückel limiting law**. The activity is only proportional to the concentration and is altered by a factor known as the activity coefficient . This factor takes into account the interaction energy of ions in solution.

## Contents

## Debye–Hückel limiting lawEdit

In order to calculate the activity of an ion C in a solution, one must know the concentration and the activity coefficient:

where

- is the activity coefficient of C,
- is the concentration of the chosen
*standard state*, e.g. 1 mol/kg if molality is used, - is a measure of the concentration of C.

Dividing with gives a dimensionless quantity.

The **Debye–Hückel limiting law** enables one to determine the activity coefficient of an ion in a *dilute* solution of known ionic strength. The equation is^{[1]}^{:section 2.5.2}

where

- is the charge number of ion species
*i*, - is the elementary charge,
- is the inverse of the Debye screening length (defined below),
- is the relative permittivity of the solvent,
- is the permittivity of free space,
- is Boltzmann's constant,
- is the temperature of the solution,
- is Avogadro's number,
- is the ionic strength of the solution (defined below),
- is a constant that depends on temperature. If is expressed in terms of molality, instead of molarity (as in the equation above and in the rest of this article), then an experimental value for
*of water*is at 25 °C. It is common to use a base-10 logarithm, in which case we factor , so*A*is .

It is important to note that because the ions in the solution act together, the activity coefficient obtained from this equation is actually a mean activity coefficient.

The excess osmotic pressure obtained from Debye–Hückel theory is in cgs units:^{[1]}

Therefore, the total pressure is the sum of the excess osmotic pressure and the ideal pressure . The osmotic coefficient is then given by

## Summary of Debye and Hückel's first article on the theory of dilute electrolytesEdit

The English title of the article is "On the Theory of Electrolytes. I. Freezing Point Depression and Related Phenomena". It was originally published in 1923 in volume 24 of a German-language journal *Physikalische Zeitschrift*. An English translation^{[2]}^{:217–63} of the article is included in a book of collected papers presented to Debye by "his pupils, friends, and the publishers on the occasion of his seventieth birthday on March 24, 1954".^{[2]}^{:xv} The article deals with the calculation of properties of electrolyte solutions that are not under the influence of net electric fields, thus it deals with electrostatics.

In the same year they first published this article, Debye and Hückel, hereinafter D&H, also released an article that covered their initial characterization of solutions under the influence of electric fields called "On the Theory of Electrolytes. II. Limiting Law for Electric Conductivity", but that subsequent article is not (yet) covered here.

In the following summary (as yet incomplete and unchecked), modern notation and terminology are used, from both chemistry and mathematics, in order to prevent confusion. Also, with a few exceptions to improve clarity, the subsections in this summary are (very) condensed versions of the same subsections of the original article.

### IntroductionEdit

D&H note that the Guldberg–Waage formula for electrolyte species in chemical reaction equilibrium in classical form is^{[2]}^{:221}

where

- is a notation for multiplication,
- is a dummy variable indicating the species,
- is the number of species participating in the reaction,
- is the mole fraction of species ,
- is the stoichiometric coefficient of species ,
*K*is the equilibrium constant.

D&H say that, due to the "mutual electrostatic forces between the ions", it is necessary to modify the Guldberg–Waage equation by replacing with , where is an overall activity coefficient, not a "special" activity coefficient (a separate activity coefficient associated with each species)—which is what is used in modern chemistry as of 2007^{[update]}.

The relationship between and the special activity coefficients is^{[2]}^{:248}

### FundamentalsEdit

D&H use the Helmholtz and Gibbs free entropies and to express the effect of electrostatic forces in an electrolyte on its thermodynamic state. Specifically, they split most of the thermodynamic potentials into classical and electrostatic terms:

where

- is Helmholtz free entropy,
- is entropy,
- is internal energy,
- is temperature,
- is Helmholtz free energy.

D&H give the total differential of as^{[2]}^{:222}

where

By the definition of the total differential, this means that

which are useful further on.

As stated previously, the internal energy is divided into two parts:^{[2]}^{:222}

where

- indicates the classical part,
- indicates the electric part.

Similarly, the Helmholtz free entropy is also divided into two parts:

D&H state, without giving the logic, that^{[2]}^{:222}

It would seem that, without some justification,

Without mentioning it specifically, D&H later give what might be the required (above) justification while arguing that , an assumption that the solvent is incompressible.

The definition of the Gibbs free entropy is^{[2]}^{:222–3}

where is Gibbs free energy.

D&H give the total differential of as^{[2]}^{:222}

At this point D&H note that, for water containing 1 mole per liter of potassium chloride (nominal pressure and temperature aren't given), the electric pressure amounts to 20 atmospheres. Furthermore, they note that this level of pressure gives a relative volume change of 0.001. Therefore, they neglect change in volume of water due to electric pressure, writing^{[2]}^{:223}

and put

D&H say that, according to Planck, the classical part of the Gibbs free entropy is^{[2]}^{:223}

where

- is a species,
- is the number of different particle types in solution,
- is the number of particles of species
*i*, - is the particle specific Gibbs free entropy of species
*i*, - is Boltzmann's constant,
- is the mole fraction of species
*i*.

Species zero is the solvent. The definition of is as follows, where lower-case letters indicate the particle specific versions of the corresponding extensive properties:^{[2]}^{:223}

D&H don't say so, but the functional form for may be derived from the functional dependence of the chemical potential of a component of an ideal mixture upon its mole fraction.^{[3]}

D&H note that the internal energy of a solution is lowered by the electrical interaction of its ions, but that this effect can't be determined by using the crystallographic approximation for distances between dissimilar atoms (the cube root of the ratio of total volume to the number of particles in the volume). This is because there is more thermal motion in a liquid solution than in a crystal. The thermal motion tends to smear out the natural lattice that would otherwise be constructed by the ions. Instead, D&H introduce the concept of an ionic atmosphere or cloud. Like the crystal lattice, each ion still attempts to surround itself with oppositely charged ions, but in a more free-form manner; at small distances away from positive ions, one is more likely to find negative ions and vice versa.^{[2]}^{:225}

### The potential energy of an arbitrary ion solutionEdit

Electroneutrality of a solution requires that^{[2]}^{:233}

where

- is the total number of ions of species
*i*in the solution, - is the charge number of species
*i*.

To bring an ion of species *i*, initially far away, to a point within the ion cloud requires interaction energy in the amount of , where is the elementary charge, and is the value of the scalar electric potential field at . If electric forces were the only factor in play, the minimal-energy configuration of all the ions would be achieved in a close-packed lattice configuration. However, the ions are in thermal equilibrium with each other and are relatively free to move. Thus they obey Boltzmann statistics and form a Boltzmann distribution. All species' number densities are altered from their bulk (overall average) values by the corresponding Boltzmann factor , where is the Boltzmann constant, and is the temperature.^{[4]} Thus at every point in the cloud^{[2]}^{:233}

Note that in the infinite temperature limit, all ions are distributed uniformly, with no regard for their electrostatic interactions.^{[2]}^{:227}

The charge density is related to the number density:^{[2]}^{:233}

When combining this result for the charge density with the Poisson equation from electrostatics, a form of the Poisson–Boltzmann equation results:^{[2]}^{:233}

This equation is difficult to solve and does not follow the principle of linear superposition for the relationship between the number of charges and the strength of the potential field. It has been solved by the Swedish mathematician Thomas Hakon Gronwall and his collaborators physicical chemists V. K. La Mer and Karl Sandved in a 1928 article from Physikalische Zeitschrift dealing with extensions to Debye–Huckel theory, which resorted to Taylor series expansion.

However, for sufficiently low concentrations of ions, a first-order Taylor series expansion approximation for the exponential function may be used ( for ) to create a linear differential equation (Hamann, Hamnett, and Vielstich. Electrochemistry. Wiley-VCH. section 2.4.2). D&H say that this approximation holds at large distances between ions,^{[2]}^{:227} which is the same as saying that the concentration is low. Lastly, they claim without proof that the addition of more terms in the expansion has little effect on the final solution.^{[2]}^{:227} Thus

The Poisson–Boltzmann equation is transformed to^{[2]}^{:233}

because the first summation is zero due to electroneutrality.^{[2]}^{:234}

Factor out the scalar potential and assign the leftovers, which are constant, to . Also, let be the ionic strength of the solution:^{[2]}^{:234}

So, the fundamental equation is reduced to a form of the Helmholtz equation:^{[5]}

Today, is called the Debye screening length. D&H recognize the importance of the parameter in their article and characterize it as a measure of the thickness of the ion atmosphere, which is an electrical double layer of the Guoy–Chapman type.^{[2]}^{:229}

The equation may be expressed in spherical coordinates by taking at some arbitrary ion:^{[6]}^{[2]}^{:229}

The equation has the following general solution (keep in mind that is a positive constant):^{[2]}^{:229}

where , , and are undetermined constants

The electric potential is zero at infinity by definition, so must be zero.^{[2]}^{:229}

In the next step, D&H assume that there is a certain radius , beyond which no ions in the atmosphere may approach the (charge) center of the singled out ion. This radius may be due to the physical size of the ion itself, the sizes of the ions in the cloud, and any water molecules that surround the ions. Mathematically, they treat the singled out ion as a point charge to which one may not approach within the radius .^{[2]}^{:231}

The potential of a point charge by itself is

D&H say that the total potential inside the sphere is^{[2]}^{:232}

where is a constant that represents the potential added by the ionic atmosphere. No justification for being a constant is given. However, one can see that this is the case by considering that any spherical static charge distribution is subject to the mathematics of the shell theorem. The shell theorem says that no force is exerted on charged particles inside a sphere (of arbitrary charge).^{[7]} Since the ion atmosphere is assumed to be (time-averaged) spherically symmetric, with charge varying as a function of radius , it may be represented as an infinite series of concentric charge shells. Therefore, inside the radius , the ion atmosphere exerts no force. If the force is zero, then the potential is a constant (by definition).

In a combination of the continuously distributed model which gave the Poisson–Boltzmann equation and the model of the point charge, it is assumed that at the radius , there is a continuity of and its first derivative. Thus^{[2]}^{:232}

By the definition of electric potential energy, the potential energy associated with the singled out ion in the ion atmosphere is^{[2]}^{:230 & 232}

Notice that this only requires knowledge of the charge of the singled out ion and the potential of all the other ions.

To calculate the potential energy of the entire electrolyte solution, one must use the multiple-charge generalization for electric potential energy:^{[2]}^{:230 & 232}

### The additional electric term to the thermodynamic potentialEdit

This article is missing information about not yet added. (May 2018) |

## NondimensionalizationEdit

This article's factual accuracy is disputed. (May 2018) (Learn how and when to remove this template message) |

This section was created without reference to the original paper and there are some errors in it (for instance, the ionic strength is off by a factor of two). Once these are rectified, this section should probably be moved to the nondimensionalization article and then be linked from here, since the nondimensional version of the Poisson–Boltzmann equation isn't necessary to understand the D&H theory.

The differential equation is ready for solution (as stated above, the equation only holds for low concentrations):

Using the Buckingham π theorem on this problem results in the following dimensionless groups:

is called the reduced scalar electric potential field. is called the reduced radius. The existing groups may be recombined to form two other dimensionless groups for substitution into the differential equation. The first is what could be called the square of the reduced inverse screening length, . The second could be called the reduced central ion charge, (with a capital Z). Note that, though is already dimensionless, without the substitution given below, the differential equation would still be dimensional.

To obtain the nondimensionalized differential equation and initial conditions, use the groups to eliminate in favor of , then eliminate in favor of while carrying out the chain rule and substituting , then eliminate in favor of (no chain rule needed), then eliminate in favor of , then eliminate in favor of . The resulting equations are as follows:

For table salt in 0.01 M solution at 25 °C, a typical value of is 0.0005636, while a typical value of is 7.017, highlighting the fact that, in low concentrations, is a target for a zero order of magnitude approximation such as perturbation analysis. Unfortunately, because of the boundary condition at infinity, regular perturbation does not work. The same boundary condition prevents us from finding the exact solution to the equations. Singular perturbation may work, however.

## Experimental verification of the theoryEdit

To verify the validity of the Debye–Hückel theory, many experimental ways have been tried, measuring the activity coefficients: the problem is that we need to go towards very high dilutions.
Typical examples are: measurements of vapour pressure, freezing point, osmotic pressure (indirect methods) and measurement of electric potential in cells (direct method).
Going towards high dilutions good results have been found using liquid membrane cells, it has been possible to investigate aqueous media 10^{−4} M and it has been found that for 1:1 electrolytes (as NaCl or KCl) the Debye–Hückel equation is totally correct, but for 2:2 or 3:2 electrolytes it is possible to find negative deviation from the Debye–Hückel limit law: this strange behavior can be observed only in the very dilute area, and in more concentrate regions the deviation becomes positive.
It is possible that Debye–Hückel equation is not able to foresee this behavior because of the linearization of the Poisson–Boltzmann equation, or maybe not: studies about this have been started only during the last years of the 20th century because before it wasn't possible to investigate the 10^{−4} M region, so it is possible that during the next years new theories will be born.

## Extensions of the theoryEdit

Warning: The notation in this section is (currently) different from in the rest of the article.

A number of approaches have been proposed to extend the validity of the law to concentration ranges as commonly encountered in chemistry

One such extended Debye–Hückel equation is given by:

where as its common logarithm is the activity coefficient, is the integer charge of the ion (1 for H^{+}, 2 for Mg^{2+} etc.), is the ionic strength of the aqueous solution, and is the size or **effective diameter** of the ion in angstrom. The effective hydrated radius of the ion, a is the radius of the ion and its closely bound water molecules. Large ions and less highly charged ions bind water less tightly and have smaller hydrated radii than smaller, more highly charged ions. Typical values are 3Å for ions such as H^{+},Cl^{−},CN^{−}, and HCOO^{−}. The effective diameter for the hydronium ion is 9Å. and are constants with values of respectively 0.5085 and 0.3281 at 25 °C in water ^{[2]}.

The extended Debye–Hückel equation provides accurate results for μ ≤ 0.1. For solutions of greater ionic strengths, the Pitzer equations should be used. In these solutions the activity coefficient may actually increase with ionic strength.

The Debye–Hückel equation cannot be used in the solutions of surfactants where the presence of micelles influences on the electrochemical properties of the system (even rough judgement overestimates γ for ~50%).

## See alsoEdit

## NotesEdit

**^**http://homepages.rpi.edu/~keblip/THERMO/chapters/Chapter33.pdf, page 9.- ^
^{a}^{b}^{c}^{d}^{e}^{f}^{g}^{h}^{i}^{j}^{k}^{l}^{m}^{n}^{o}^{p}^{q}^{r}^{s}^{t}^{u}^{v}^{w}^{x}^{y}^{z}^{aa}^{ab}^{ac}^{ad}^{ae}^{af}*The Collected Papers of Peter J. W. Debye*. New York, New York: Interscience Publishers, Inc. 1954. **^**Salzman, William R. (2001-08-21). "Mixtures; Partial Molar Quantities; Ideal Solutions".*Chemical Thermodynamics*. University of Arizona. Archived from the original on 2007-06-07. Retrieved 2007-10-23.**^**http://www.pma.caltech.edu/Courses/ph136/yr2006/text.html, section 19.3**^**http://guava.physics.uiuc.edu/~nigel/courses/569/Essays_2004/files/lu.pdf section 3.1)**^**http://hyperphysics.phy-astr.gsu.edu/hbase/electric/laplace.html**^**http://hyperphysics.phy-astr.gsu.edu/hbase/electric/potsph.html

## ReferencesEdit

- P. Debye; E. Hückel (1923). "Zur Theorie der Elektrolyte. I. Gefrierpunktserniedrigung und verwandte Erscheinungen" [The theory of electrolytes. I. Lowering of freezing point and related phenomena] (PDF).
*Physikalische Zeitschrift*.**24**: 185–206.^{[permanent dead link]} **^**Hamann, Hamnett, and Vielstich (1998).*Electrochemistry*. Weinheim: Wiley-VCH Verlag GmbH. ISBN 3-527-29096-6.CS1 maint: Multiple names: authors list (link)**^**Harris, Daniel C. (2003).*Quantitative Chemical Analysis*(6th ed.). W. H. Freeman & Company. ISBN 0-7167-4464-3.**^**Skoog, Douglas A.*Fundamentals of Analytical Chemistry*. ISBN 0-534-41796-5.- F. Malatesta, R. Zamboni. Activity and osmotic coefficients from the EMF of liquid membrane cells, VI - ZnSO4, MgSO4, CaSO4 and SrSO4 in water at 25 C. Journal of Solution Chemistry 1997, 26, 791–815

## External linksEdit

- For easy calculation of activity coefficients in (non-micellar) solutions, check out the IUPAC open project Aq-solutions
^{[permanent dead link]}(freeware)**.** - Gold Book definition