Understanding Diffusion and Brownian Motion
Diffusion is the spontaneous spreading of dissolved particles driven by random molecular collisions. Particles move preferentially from high-concentration regions to low-concentration regions until equilibrium is reached. This process, called Brownian motion when observable at the nanoscale, underlies countless natural and industrial phenomena—from drug delivery in the bloodstream to contaminant transport in groundwater.
The kinetics of diffusion are governed by Fick's laws. The first law relates the flux (mass per unit area per unit time) to the concentration gradient. When particles also experience an external force, the equation becomes:
J = −D(dc/dx) + c·f/ξ
where J is the flux, D is the diffusion coefficient, c is concentration, f is an applied force, and ξ is the friction coefficient. From this relationship, Einstein and Smoluchowski independently derived a fundamental connection between diffusion and friction that remains central to modern transport theory.
The Einstein-Smoluchowski Relation
The diffusion coefficient is determined by a remarkably simple formula that balances thermal energy against frictional drag:
D = (kB × T) / ξ
D— Diffusion coefficient (m²/s)k<sub>B</sub>— Boltzmann's constant, 1.380649 × 10⁻²³ J/KT— Absolute temperature in kelvin (K)ξ— Friction coefficient (Pa·s or kg/s)
Friction Coefficients for Common Particle Shapes
The friction coefficient depends critically on particle geometry and solvent viscosity. A sphere experiences less resistance than an elongated ellipsoid of comparable volume. The table below shows expressions for friction coefficient across standard shapes, where η is dynamic viscosity, a is the major dimension, and b is the minor dimension:
- Sphere: ξ = 6π·η·a
- Disk (face-on): ξ = 16·η·a
- Disk (edge-on): ξ = (32/3)·η·a
- Disk (random tumbling): ξ = 12·η·a
- Ellipsoid (lengthwise motion): ξ = 4π·η·a / [ln(2a/b) + 0.5]
- Ellipsoid (sideways motion): ξ = 8π·η·a / [ln(2a/b) + 0.5]
- Ellipsoid (random orientation): ξ = 6π·η·a / ln(2a/b)
Non-spherical particles exhibit orientation-dependent friction. An elongated molecule moving along its long axis encounters less drag than when tumbling randomly—a distinction crucial for predicting behaviour in biological fluids and polymer solutions.
Practical Example: Viral Diffusion in Hot Springs
Consider Sulfolobus ellipsoid virus 1, an unusual rod-shaped virus with major axis a = 115 nm and minor axis b = 78 nm, suspended in a hot spring at 90°C. Water viscosity at this temperature is approximately 0.000315 Pa·s.
For random tumbling motion, the friction coefficient is:
ξ = 6π × 0.000315 × 115 × 10⁻⁹ / ln(2 × 115/78) = 7.23 × 10⁻¹¹ s/kg
The diffusion coefficient then follows:
D = (1.381 × 10⁻²³ × 363.15) / (7.23 × 10⁻¹¹) = 6.6 × 10⁻¹¹ m²/s
This relatively modest diffusion coefficient reflects the large size and irregular shape of the virus, which impedes rapid migration through the viscous aqueous medium.
Common Pitfalls and Considerations
Several practical issues frequently arise when applying the diffusion coefficient formula:
- Unit consistency is essential — Ensure all inputs use SI units: temperature in kelvin (add 273.15 to Celsius), viscosity in Pa·s (not cP), dimensions in metres, and Boltzmann's constant in J/K. Mismatched units are the leading source of errors by orders of magnitude.
- Temperature sensitivity dominates — The diffusion coefficient increases linearly with absolute temperature. A 10°C rise roughly increases diffusion rate by 3%. For temperature-dependent studies, verify your solvent viscosity at the exact working temperature, as viscosity varies sharply and nonlinearly.
- Shape assumptions matter for non-spherical particles — The friction coefficient formulas assume rigid, smooth particles. Real proteins, polymers, and aggregates may be flexible, hydrated, or rough, causing measured friction to exceed theoretical predictions. Empirical validation against experimental data is strongly recommended for novel systems.
- Reverse calculations from measured diffusion — Laboratory instruments (dynamic light scattering, analytical ultracentrifugation) often measure diffusion directly. You can invert the formula to extract friction coefficient: ξ = (k<sub>B</sub> × T) / D_measured. This indirect method sometimes yields cleaner results than direct size measurement.