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:
||charge number of ion i
||= 9.6485·104 Coulomb/mol
||= 8.31446 J/(K mol)
||= 3.7554·106 s · S mol-1
||proportionality constant at 25°C
||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:
||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 Appelo who rearranged Eq.(3) into
where all aspects regarding the ion-ion interaction are put into the correction factor
with I as ionic strength. It is no surprise that this ion-ion correction factor closely resembles the activity model of Debye-Hückel:
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)
Appelo proposed the following parametrization for Eq.(8):
Just this parametrization is used in PhreeqC and aqion. The equation behind the EC calculation is:
This equation is used as the default method in aqion. The diffusion coefficients Di are taken from here.
[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:
We start from Eq.(7b) and get from Kohlrausch’s Square-Root Law:
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. 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
[last modified: 2016-01-03]