Electrical Conductivity based on Diffusion Coefficients
The Nernst-Einstein equation establishes the relationship between the molar limiting conductivity and the diffusion coefficient Di for any given ion i:
|zi||charge number of ion i|
|T||in K||absolute temperature|
|F||= 9.6485·104 Coulomb/mol||Faraday’s constant|
|R||= 8.31446 J/(K mol)||gas constant|
|F2/(RT)||= 3.7554·106 s · S mol-1||proportionality constant at 25°C|
|Di||in m2/s||diffusion coefficient of ion i|
|in S m2/mol (= 104 S cm2/mol)||molar limiting conductivity of ion i|
Example. If we enter some well-known values of Di (taken from literature) into Eq.(1) the following molar limiting conductivities are obtained:
|ion||Di [m2 s-1]||[S cm2 mol-1]|
The goal is now to exploit this method for the calculation of electrical conductivities (EC) of aqueous solutions of arbitrary composition.
EC of Ideal Aqueous Solutions
In the limit of infinite dilution (non-interacting ions) we obtain from Eq.(7a) a simple formula that relies on diffusion coefficients:
In the very realistic case of non-ideal solutions, however, the method becomes somewhat more elaborate.
EC of Real (Non-Ideal) Aqueous Solutions
The Nernst-Einstein equation is restricted and valid only for molar limiting conductivities . In contrast, the EC of real aqueous solutions rest upon molar conductivities :
Both quantities are related by Kohlrausch’s Square-Root Law, but this law requires the knowledge of an extra parameter K which depends non-trivially on the type of electrolyte (and which is hardly available in tables or literature). An alternative was proposed by Appelo1 who rearranged Eq.(3) into
where all aspects regarding the ion-ion interaction are put into the correction factor
with A = 0.5085 M-1/2. This strategy looks promising because we replaced the non-trivial Kohlrausch parameter K by the activity constant γ – a quantity that belongs to the standard repertoire of hydrochemical models anyway (available for each ion and aqueous species).
Formally, Eq.(4) can be converted into
The only thing we need is a clever parametrization of α as a fairly constant quantity.
Parametrization à la Appelo (PhreeqC 3)
Appelo1 proposed the following parametrization for Eq.(8):
Just this parametrization is used in PhreeqC and aqion. The equation behind the EC calculation is:
[Remark to γi: While Eq.(6) represents the most simplest form of an activity model (in order to demonstrate the main ideas in the present article), in the program, however, more sophisticated approaches for γi are used.]
Appendix: Rearrangement of the EC Equation
The aim is to convert Eq.(3) into a form similar to the ideal-solution formula:
Using , the last line yields:
Note that c, and not ci, enters the square root. This is because c refers to the ‘medium effect’ of the electrolyte (as a composition of several ion types). In other words, c appears as a sibling of the ionic strength I = ½ Σ zi2 ci.2 Thus, replacing c by I in Eq.(A3) yields the desired equation:
The value of is near to 1, which permits the expansion exp(-a) = 1–a + … :
References & Remarks
The following example for NaCl illustrates the close relationship between the electrolyte concentration c and ionic strength: I = ½ (cNa + cCl) = cNaCl. ↩