Transport Phenomena in Electrochemistry

1. Introduction: Ion Transport in Solution

Electrochemistry depends fundamentally on the motion of charged species through electrolyte solutions. Unlike electrons moving through metals, ionic transport involves the coupled effects of diffusion (driven by concentration gradients), migration (driven by electric fields), and convection (driven by bulk fluid motion). Understanding these transport phenomena is essential for designing batteries, fuel cells, electroplating baths, and desalination systems.

When an electrolyte dissolves in a solvent, it dissociates into cations and anions that are solvated by solvent molecules. These solvated ions experience viscous drag as they move through the solvent, characterised by an ionic mobility $u_i$. The interplay between the driving forces and frictional resistance determines the macroscopic transport properties — conductivity, diffusion coefficients, and transference numbers — that govern electrochemical cell performance.

In this chapter we derive the fundamental equations governing ionic transport from first principles: starting from ionic conductivity and Kohlrausch's laws, through the Debye-Hückel theory of ion-ion interactions, to the Nernst-Planck equation that unifies diffusion and migration, and finally transference numbers that describe how current is carried by individual ionic species.

2. Derivation: Ionic Conductivity and Molar Conductivity

Consider an electrolyte solution containing ionic species $i$ with charge number $z_i$, molar concentration $c_i$, and ionic mobility $u_i$ (velocity per unit electric field, in m$^2$ V$^{-1}$ s$^{-1}$). Under an applied electric field $\mathbf{E}$, each ion species drifts with velocity:

The current density contributed by species $i$ is the charge per unit volume times velocity:

where $F$ is Faraday's constant. Summing over all species and comparing with Ohm's law $j = \kappa E$, the electrolytic conductivity is:

Molar Conductivity

The conductivity $\kappa$ depends on concentration. To isolate the intrinsic conducting ability of the electrolyte, we define the molar conductivity:

where $c$ is the overall electrolyte concentration. The units of $\Lambda_m$ are S cm$^2$ mol$^{-1}$ (or S m$^2$ mol$^{-1}$ in SI). For a binary electrolyte $\mathrm{M}_{\nu_+}\mathrm{X}_{\nu_-}$ that dissociates into$\nu_+$ cations and $\nu_-$ anions:

Kohlrausch's Law of Independent Migration

Kohlrausch discovered empirically (1874) that at infinite dilution, each ion contributes independently to the molar conductivity. We define the limiting ionic conductivity $\lambda_i^\circ$ for each ion species. The limiting molar conductivity of the electrolyte is:

This law of independent migration holds because at infinite dilution, ion-ion interactions vanish entirely and each ion moves independently through the solvent. The individual $\lambda_i^\circ$ values are characteristic of each ion and independent of the counter-ion present. For example,$\lambda^\circ(\text{Na}^+) = 50.11$ S cm$^2$ mol$^{-1}$regardless of whether the anion is Cl$^-$, Br$^-$, or NO$_3^-$.

Kohlrausch's Square-Root Law

For strong electrolytes (fully dissociated), Kohlrausch found that the molar conductivity decreases linearly with $\sqrt{c}$:

The constant $K$ depends on the nature of the electrolyte. The physical origin of this behaviour lies in ion-ion interactions: as concentration increases, each ion becomes surrounded by an ionic atmosphere of opposite charge. Under an applied field, two retarding effects emerge:

• • Electrophoretic effect: the ionic atmosphere moves in the opposite direction, creating a viscous drag on the central ion.
• • Relaxation effect: as the ion moves, its atmosphere cannot rearrange instantaneously, creating an asymmetric charge distribution that retards motion.

Onsager (1926) provided a theoretical foundation for the $\sqrt{c}$ dependence by combining Debye-Hückel theory with hydrodynamics, yielding the Debye-Hückel-Onsager equation:

where the first term is the electrophoretic contribution and the second is the relaxation contribution. Since $\kappa_D \propto \sqrt{c}$, this explains the empirical square-root law.

3. Derivation: Debye-Hückel Theory

The Debye-Hückel theory (1923) provides a quantitative model of ion-ion interactions in electrolyte solutions. The central idea is that each ion is surrounded by a diffuse ionic atmosphere carrying a net charge opposite to that of the central ion. We derive the electrostatic potential around an ion and the resulting correction to the chemical potential.

The Poisson-Boltzmann Equation

The electrostatic potential $\phi(\mathbf{r})$ around a central ion obeys Poisson's equation:

where $\varepsilon = \varepsilon_0 \varepsilon_r$ is the permittivity of the solvent and$\rho_{\text{free}}$ is the free charge density due to the mobile ions. The local concentration of ion species $i$ at position $\mathbf{r}$ follows the Boltzmann distribution:

where $c_i^0$ is the bulk concentration. The total charge density is:

Substituting into Poisson's equation gives the full Poisson-Boltzmann equation:

Linearization (Debye-Hückel Approximation)

For dilute solutions where the electrostatic energy is small compared to thermal energy ($|z_i e \phi| \ll k_B T$), we expand the exponential to first order:

The charge density becomes:

The first sum vanishes by electroneutrality ($\sum_i z_i c_i^0 = 0$). Defining the Debye parameter $\kappa$:

where $I = \tfrac{1}{2}\sum_i c_i z_i^2$ is the ionic strength. The linearized Poisson-Boltzmann equation is:

Solution in Spherical Coordinates

For a spherically symmetric potential around a point ion of charge $ze$, the Laplacian in spherical coordinates gives:

Substituting $\phi = \psi/r$ transforms this into $d^2\psi/dr^2 = \kappa^2 \psi$, with solutions $\psi = Ae^{-\kappa r} + Be^{+\kappa r}$. Requiring $\phi \to 0$ as$r \to \infty$ eliminates $B$, and matching to the bare Coulomb potential at short range determines $A$. The result is:

This is a screened Coulomb potential (Yukawa potential). The quantity $\kappa^{-1}$ is the Debye length — the characteristic distance over which the ionic atmosphere screens the central ion:

For a 0.1 mol/L solution of NaCl in water at 25°C, the Debye length is approximately 0.96 nm — just a few molecular diameters.

Activity Coefficient

The excess electrostatic energy of an ion due to its ionic atmosphere gives a correction to the chemical potential. The potential at the surface of the central ion (radius $a \to 0$) due to the atmosphere alone is obtained by subtracting the bare Coulomb potential:

The electrostatic contribution to the chemical potential of the ion is:

Since $\mu_i^{\text{el}} = k_B T \ln \gamma_i$, the mean ionic activity coefficient for a $z_+$:$z_-$ electrolyte is:

Since $\kappa \propto \sqrt{I}$, this yields the Debye-Hückel limiting law: $\ln \gamma_\pm = -A|z_+ z_-|\sqrt{I}$, where $A$ is a solvent-dependent constant ($A \approx 1.172$ (mol/kg)$^{-1/2}$ for water at 25°C).

4. Derivation: Nernst-Planck Equation

The Nernst-Planck equation describes the flux of ionic species under the combined influence of concentration gradients (diffusion) and electric fields (migration). It is the fundamental transport equation in electrochemistry.

Diffusion Flux

Fick's first law gives the diffusive flux of species $i$:

where $D_i$ is the diffusion coefficient of species $i$.

Migration Flux

Under an electric field, ions migrate with velocity $v_i = u_i z_i E = -u_i z_i (d\phi/dx)$. The migration flux is:

Using the relation $u_i = D_i |z_i| e / (k_B T) = D_i z_i F / (RT)$ (to be derived below), the migration flux becomes:

The Combined Nernst-Planck Equation

Adding the diffusion and migration contributions gives the Nernst-Planck equation:

This can be written more compactly using the electrochemical potential$\tilde{\mu}_i = \mu_i^\circ + RT \ln c_i + z_i F \phi$:

This shows that the true driving force for ionic transport is the gradient of the electrochemical potential, not the concentration gradient or electric field alone.

Derivation of the Einstein Relation

The Einstein relation connects the diffusion coefficient to the ionic mobility. Consider an ion in thermodynamic equilibrium under an electric field. At equilibrium, the net flux must vanish:

This gives the Boltzmann distribution $c_i \propto \exp(-z_i F \phi / RT)$. Now, independently, from the definition of mobility, the drift velocity is $v = u_i E$ and the migration flux is $J^{\text{mig}} = u_i c_i E$. Comparing with the Nernst-Planck migration term $J^{\text{mig}} = (z_i F D_i / RT) c_i E$, we identify:

This is the Nernst-Einstein relation. It relates a non-equilibrium transport property (diffusion coefficient) to an equilibrium response property (mobility), and is a manifestation of the fluctuation-dissipation theorem. The relation allows computation of diffusion coefficients from conductivity measurements and vice versa.

5. Derivation: Transference Numbers

In an electrolyte solution, the total current is carried by all ionic species present. The transference number (or transport number) $t_i$ of species $i$ is defined as the fraction of the total current carried by that species:

where $\lambda_i = |z_i| F u_i$ is the ionic conductivity of species $i$. By construction, $\sum_i t_i = 1$.

Derivation from Ion Velocities

Consider a binary electrolyte with cations ($+$) and anions ($-$). Under an applied field $E$, the current densities are:

The total current density is $j = j_+ + j_- = \kappa E$. The transference number of the cation is:

For a symmetric $z$:$z$ electrolyte where $c_+ = c_- = c$ and$|z_+| = |z_-| = z$:

Hittorf Method

The Hittorf method (1853) measures transference numbers by analysing concentration changes in the electrode compartments after electrolysis. When a total charge $Q$ passes through the cell:

• • Cations carry fraction $t_+$ of the current toward the cathode
• • Anions carry fraction $t_-$ of the current toward the anode
• • At each electrode, one equivalent is discharged per Faraday of charge

In the cathode compartment, the number of equivalents lost by deposition is $Q/F$, but $t_+ \cdot Q/F$ equivalents migrate in. The net change is:

Thus $t_-$ is determined from the concentration decrease in the cathode compartment.

Moving Boundary Method

The moving boundary method provides a more direct measurement. A sharp boundary is formed between two electrolyte solutions sharing a common ion. As current flows, the boundary moves with a velocity proportional to the transference number of the leading ion:

where $V$ is the volume swept by the boundary, $c$ is the electrolyte concentration, and $Q$ is the total charge passed. This method is considered more accurate than the Hittorf method because it avoids complications from electrode processes.

6. Applications

7. Historical Context

8. Python Simulation: Debye-Hückel Activity Coefficients

This simulation plots (1) Debye-Hückel limiting law activity coefficients vs ionic strength for different electrolyte types (1:1, 2:1, 2:2, etc.), (2) Debye screening length vs concentration, (3) extended Debye-Hückel vs limiting law comparison, and (4) Kohlrausch square-root law for molar conductivity of several strong electrolytes.

9. Fortran Simulation: Conductivities & Transference Numbers

High-precision computation of molar conductivities, transference numbers, Debye screening lengths, activity coefficients, and diffusion coefficients from the Einstein-Nernst relation for common electrolytes.

Key Equations Summary