PNG 550
Reactive Transport in the Subsurface

4.2 Reaction Thermodynamics and Important Parameters

PrintPrint

Cation exchange reactions

Ion exchange reactions occur when ions exchange their positions at the soid-surface – water interface. These reactions are typically assumed reversible and occur instantaneously at time scales ranging from microseconds to hours. Ion exchange is typically represented in the following form:

\begin{equation}uBX_v(s)+vA^{u+}(aq)\Leftrightarrow vAX_u(s)+uB^{v+}(aq)\end{equation}

Here (s) and (aq) refer to solid and aqueous phases, respectively; X- denotes the negatively charged surface sites that bound cations Au+ and Bv+, with u and v being the charges of A and B, respectively. In this reaction, A and B are cations that compete for sorption sites (Appelo and Postma, 1993; Sposito et al., 1981; Vanselow, 1932). The equilibrium constant Keq of reaction (1) can be expressed as

\begin{equation}K_{eq}=\frac{\left[\mathrm{AX}_{\mathrm{u}}(s)\right]^{\mathrm{y}}\left[\mathrm{B}^{v+}\right]^{\mathrm{u}}}{\left[\mathrm{BX}_{\mathrm{v}}(s)\right]^{\mathrm{u}}\left[\mathrm{A}^{v+}\right]^v}=\frac{\left[\mathrm{AX}_{\mathrm{u}}(s)\right]^{\mathrm{v}}/\left[\mathrm{A}^{\mathrm{u}+}\right]^{\mathrm{v}}}{\left[\mathrm{BX}_{\mathrm{v}}(s)\right]^{\mathrm{u}}/\left[\mathrm{B}^{\mathrm{v}+}\right]^u}\end{equation}

Here the parentheses [ ] represent activities. Aqueous concentrations are easily related to activities through concentration and activity coefficients. Activities of ions on solid phases are typically expressed as a fraction of the total, either as molar fraction or as equivalent fractions. The total number can be based on the number of exchange sites or as the number of exchangeable cations.

The units meq is often used in ion exchange calculations. An meq is the number of ions that sums a specific quantity of electrical charges. For example, an meq of K+ is about 6.02 x 1020 positive charges. On the other hand, an meq of Ca2+ is also 6.02 x 1020 positive charges, however only 3.01 x 1020 ions because Ca2+ has two positive charges. For an ion Ii+ with a charge of i, the equivalent fraction bI is calculated as:

\begin{equation}\beta_I=\dfrac{meq\mathrm{I}-\mathrm{X}_i\text{ per kg sediment }}{CEC}=\dfrac{meq_{I-X_1}}{\sum\limits_{I,J,K..}^{ }meq_{\mathrm{I-X_1}}}\end{equation}

where I, J, K are exchangeable cation with charges i, j, k, respectively. A molar fraction $\beta_{I}^{M}$ is obtained by the following form:

\begin{equation}\beta_I^M=\frac{mmol\mathrm{\ I}-\mathrm{X}_i\text{ per }\mathrm{kg}\text{ sediment }}{TEC}=\frac{\left(meq_{I-X_i}\right)/i}{\sum\limits_{I,J,K...}^{ }\left(meq_{I-X_i}\right)/i}\end{equation}

Here TEC is the total exchangeable cations in mmol/kg sediment, not cation exchange capacity. The use of fractions should give the summation of fractions being 1, that is, β=1 .

Three common conventions used in writing ion exchange equilibrium constants are Gaines-Thomas convention, Gapon convention, and Vanselow convention. If the ion exchange is between cations of the same valence (homovalent exchange), the convention does not make a difference. If the exchange is between cations of different valences (heterovalent), the convention makes a difference. For example, the ion exchange reaction between Na and Ca can be written as follows:

\begin{equation}\mathrm{Na}^{+}+1 / 2 \mathrm{Ca}-X_{2} \leftrightarrow \mathrm{Na}-\mathrm{X}+\mathrm{1} / 2 \mathrm{Ca}^{2+}\end{equation}

With

\begin{equation}K_{\text{Na\Ca}} =\frac{[Na-X]\left[Ca^{2+}\right]^{0.5}}{\left[Ca-X_2\right]^{0.5}\left[Na^+\right]}=\frac{\beta_{Na}\left[Ca^{2+}\right]^{0.5}}{\beta_{Ca}^{0.5}\left[Na^+\right]}\end{equation}

Here if the equivalent fraction of the exchangeable cations is used for $β$ values, the Gaines-Thomas convention is followed (Graines and Thomas, 1953). If we use molar fractions for $β$ values, we follow the Vanselow convention (Vanselow, 1932). If the activities of cations on exchange sites are expressed as a fraction of the number of exchange sites (X-), we follow the Gapon convention, the reaction will be written as follows:

\begin{equation}\mathrm{Na}^{+}+\mathrm{Ca}_{0.5}-X \leftrightarrow \mathrm{Na}-\mathrm{X}+1 / 2 \mathrm{Ca}^{2+}\end{equation}

Where

\begin{equation}K_{N a \mid C a}=\frac{[N a-X]\left[\mathrm{Ca}^{2+}\right]^{0.5}}{\left[\mathrm{Ca}_{0.5}-X\right]\left[\mathrm{Na}^{+}\right]}=\frac{\beta_{\mathrm{Na}}\left[\mathrm{Ca}^{2+}\right]^{0.5}}{\beta_{\mathrm{Ca}}\left[\mathrm{Na}^{+}\right]}\end{equation}

Here the activities are expressed in terms of the mole fraction of the total number of exchangeable sites. To fully understand the difference between different conventions and calculation, please go over Example 6.3 in chapter 6 of Appelo and Postma (2005).

Selectivity coefficient

The capabilities of ions to compete for exchange sites are governed by their affinity to the surface of the exchangers. Selectivity coefficients have been reported in literature for common ions including Na+, K+, Ca2+, and Mg2+, however rarely for trace metals. The larger the selectivity coefficient, the higher affinity of the ion to the solid phase and the more competitive for exchange. In general, ions have high affinity to exchange sites when they have higher valence, are less solvated with water molecules, and react strongly with the surface sites. The following Table 1 lists selectivity coefficients reported in literatures for ion exchange on kaolinite (Appelo and Postma, 2005; Bundschuh and Zilberbrand, 2011).

Table 1. Ion exchange reaction on kaolinite (Bundschuh and Zilberbrand, 2011)
No. Ion exchange reaction log K
1 $2 \mathrm{Na}-X+\mathrm{Mg}^{2+} \Leftrightarrow \mathrm{Na}^{+}+\mathrm{Mg}-X_{2}$ 0.60 (0.44 ~ 0.78)
2 $2 \mathrm{Na}-\mathrm{X}+\mathrm{Ca}^{2+} \Leftrightarrow \mathrm{Na}^{+}+\mathrm{Ca}-X_{2}$ 0.80 (0.44 ~ 0.104)
3 $2 \mathrm{Na}-X+\mathrm{Ba}^{2+} \Leftrightarrow \mathrm{Na}^{+}+\mathrm{Ba}-X_{2}$ 0.91 (0.44 ~ 0.104)
4 $2 \mathrm{Na}-X+\mathrm{Sr}^{2+} \Leftrightarrow \mathrm{Na}^{+}+\mathrm{Sr}-X_{2}$ 0.91 (0.44 ~ 0.104)

For typical ion exchangers, the sequence of affinity is as follows: Ba2+> Pb2+ > Sr2+ > Ca2+ > Ni2+ > Cd2+ > Cu2+ > Co2+> Zn2+ > Mg2+ > Ag+> Cs+ > K+ >NH4+ >Na+ >H+. When there are multiple cations co-exist in the solution, the surface exchange site composition is largely determined by the cation aqueous concentration and their affinity to the exchange sites. To calculate the exchange site composition, that is, the mole fraction of each cation on the exchange site, please follow the example 6.4 in chapter 6 of Appelo and Postma (2005).

Cation exchange capacity

The cation exchange capacity (CEC, meq/kg solid) is a measure of the solid phase capacity for ion exchange reactions (Meunier, 2005). The CEC of different porous media are very much associated with their clay content, organic carbon, and grain size. Different materials have different CEC values. CEC values of clay minerals such as muscovite, illite, kaolinite, and chlorite are high for their grain sizes smaller than 2 μm (Drever, 1982). In general, organic matter has the highest CEC values (1500 - 4000 meq/kg). Iron oxides play a vital role in natural processes and controls nutrient availability and heavy-metal mobility (Houben and Kaufhold, 2011). Iron oxides, such as goethite and hematite, have CEC values from 40 to 1000 meq/kg. Many solid materials have iron oxides or organic matters coated on their surface and therefore have large CEC values. For soils, CEC value is a function of solution pH depending on hydrolysis reactions of surface sites. In general, cation exchange occurs due to the broken bonds around the crystal edges, the substitutions within the lattice, and the hydrogen of exposed surface hydroxyls that may be exchanged. Higher pH values give rise to more negative charges on clay, resulting in higher CEC. CEC values increase as the grain size decreases due to the large surface area associated with smaller grains.

Table 2. Ion exchange properties for clay minerals (Bundschuh and Zilberbrand, 2011)
Minerals Grain size (μm) Surface area (m2/g) CEC (meq/kg)
Kaolinite 0.1-5.0 5-20 30-150
Illite 0.1-2.0 15-40 150-400
Montmorillonite 0.01-1.0 600-800 800-1200