Introduction & Context

Molecular diffusivity quantifies how fast a solute migrates through a solvent due to random thermal motion. In process engineering it governs the rate of mass-transfer limited steps such as gas absorption, liquid-liquid extraction, crystallization, membrane separation, and heterogeneous catalysis. Accurate estimates allow engineers to size equipment (column heights, residence times, film thicknesses) and to interpret lab-scale kinetic data. The Einstein–Stokes equation links the macroscopic diffusion coefficient to solvent viscosity, temperature, and solute size, providing a quick first-principles estimate when experimental data are unavailable.

Methodology & Formulas

Convert temperature from Celsius to absolute scale \[T\,[\mathrm{K}]=T\,[^\circ\mathrm{C}]+273.15\]
Convert dynamic viscosity from centipoise to SI units \[\mu\,[\mathrm{Pa\,s}]=\mu\,[\mathrm{cP}]\times10^{-3}\]
Convert solute radius from nanometres to metres \[R\,[\mathrm{m}]=R\,[\mathrm{nm}]\times10^{-9}\]
Apply the Einstein–Stokes relation for the diffusion coefficient of a spherical particle in a continuum fluid \[D=\frac{k_{\mathrm{B}}T}{6\pi\mu R}\] where \(k_{\mathrm{B}}\) is the Boltzmann constant, \(1.380649\times10^{-23}\ \mathrm{J\,K^{-1}}\).

Continuum & creeping-flow regime limits
Parameter	Condition	Interpretation
Particle Reynolds number	\(\displaystyle Re_{\mathrm{p}}=\frac{2R\rho v}{\mu}\ll1\)	Inertial effects negligible; Stokes drag valid
Knudsen number	\(\displaystyle Kn=\frac{\lambda}{2R}\ll1\)	Fluid behaves as continuum; no slip at surface
Schmidt number	\(\displaystyle Sc=\frac{\mu}{\rho D}\gg1\)	Momentum diffusivity dominates mass diffusivity

When the above criteria are met, the Einstein–Stokes estimate provides a reliable order-of-magnitude value for design calculations and scale-up analyses.

The relation is reliable when all of the following are true:

Solute diameter > 5× solvent diameter (large spherical or quasi-spherical molecule)
Temperature is well below the solvent normal boiling point (viscosity is Newtonian)
System is dilute (< 1 mol % solute) so coupling effects are negligible
No specific interactions such as hydrogen bonding or micelle formation occur

Outside these limits use a semi-empirical correlation (Wilke-Chang, Hayduk-Minhas, etc.) or obtain experimental data.

Use the hydrodynamic (Stokes) radius, not the van-der-Waals or covalent radius:

Estimate from the molar volume at the normal boiling point: r = (3 V_b / 4πN_A)^1/3
If V_b is unknown, group-contribution methods (e.g., Joback) give ±10 % accuracy
For long-chain or disk-like molecules, take the radius of the equivalent sphere that has the same hydrodynamic drag; aspect-ratio corrections are available in literature

Use the solution viscosity at the composition of interest, not the pure-solvent value:

Measure it with a viscometer or retrieve from DIPPR/Dechema if available
If data are missing, blend pure-component viscosities with a mixing rule (e.g., Grunberg-Nissan) and validate against a single measurement point
Remember that even trace water or cosolvents can change μ by 20–30 %, so book values for “dry” solvent may mislead

For rigid hydrophobic solutes 0.5–2 nm in diameter, expect ±30 % of experimental D_AB; deviations rise to ±50–80 % for:

Small ions (hydration shell changes effective radius)
Flexible polymers (segmental motion not captured)
Supercritical solvents (density, not viscosity, controls friction)

Always benchmark against a lab measurement if the mass-transfer design margin is < 50 %.

Worked Example: Estimating Molecular Diffusivity for a Trace Contaminant in Water

A process engineer needs to estimate how fast a trace pharmaceutical contaminant (approximated as a spherical molecule with radius 0.45 nm) will diffuse through a 25 °C water stream. The stream viscosity is 0.89 cP. Use the Einstein–Stokes equation to obtain the molecular diffusivity.

Knowns

Temperature: 25.0 °C (298.15 K)
Dynamic viscosity of water: 0.89 cP = 0.00089 Pa·s
Molecular radius: 0.45 nm = 4.5 × 10⁻¹⁰ m
Boltzmann constant: 1.381 × 10⁻²³ J·K⁻¹

Step-by-Step Calculation

Convert temperature to kelvin: \( T = 25.0 + 273.15 = 298.15\ \text{K} \)
Convert viscosity to SI units: \( \mu = 0.89\ \text{cP} = 0.00089\ \text{Pa·s} \)
Convert radius to metres: \( r = 0.45\ \text{nm} = 4.5 \times 10^{-10}\ \text{m} \)
Compute the denominator of the Einstein–Stokes equation: \[ 6\pi\mu r = 6 \times 3.142 \times 0.00089\ \text{Pa·s} \times 4.5 \times 10^{-10}\ \text{m} = 7.55 \times 10^{-12}\ \text{kg·m·s}^{-2} \]
Calculate diffusivity: \[ D = \frac{kT}{6\pi\mu r} = \frac{(1.381 \times 10^{-23}\ \text{J·K}^{-1})(298.15\ \text{K})}{7.55 \times 10^{-12}\ \text{kg·m·s}^{-2}} = 5.45 \times 10^{-10}\ \text{m}^2\text{·s}^{-1} \]

Final Answer

The estimated molecular diffusivity of the contaminant in water at 25 °C is 5.45 × 10⁻¹⁰ m²·s⁻¹.

"Un projet n'est jamais trop grand s'il est bien conçu." — André Citroën

"La difficulté attire l'homme de caractère, car c'est en l'étreignant qu'il se réalise." — Charles de Gaulle