Diffusion-Controlled Reactions in Solution

Introduction

In solution, two reactants must move through solvent before they can react. If nearly every suitable encounter makes product, the overall rate can be limited by how quickly diffusion brings the molecules together. If most encounters separate without reacting, a chemical activation step remains important. The diffusion-controlled limit is therefore a useful upper benchmark for many bimolecular reactions in liquids, but its numerical value depends on solvent, solute size, charge and the definition of reactive contact.

Core explanation

For two dilute, neutral, approximately spherical solutes A and B, the classical Smoluchowski model treats their relative motion as diffusion toward an absorbing sphere. Let R c be the capture radius, roughly the centre-to-centre separation at which an encounter is counted. Let D A and D B be their diffusion coefficients. If they diffuse independently, relative diffusion is D rel = D A + D B. The encounter-rate constant per particle-number concentration in SI form is k D,particle = 4πR cD rel, with units m³ s⁻¹. To express a second-order constant for concentrations in mol L⁻¹, multiply by N A and 1000 L m⁻³: k D = 4πR cD relN A × 1000, in L mol⁻¹ s⁻¹. An ACS primary discussion of diffusive charge transfer uses the Smoluchowski form with an encounter efficiency factor, while the IUPAC encounter-controlled definition gives an order-of-magnitude benchmark near 10¹⁰ L mol⁻¹ s⁻¹ for bimolecular solutes in water at 25 °C.

The “absorbing sphere” assumption means contact produces reaction with probability effectively one. Real chemistry may require a particular orientation, electronic coupling or passage over an additional activation barrier. Then an encounter pair can separate back into bulk solvent. A simple conceptual scheme is A+B ⇌ {A···B} enc → products. If product formation is fast compared with separation after contact, observed k approaches k D. If product formation is slow, k is below the encounter limit and can be activation-controlled. This is not the same as saying every molecule travels straight toward every other one: Brownian motion is irregular and solvent continually collides with the solutes.

Viscosity affects diffusion. Under a simple Stokes–Einstein approximation for a spherical solute of hydrodynamic radius a in a continuum solvent, D ≈ k B T/(6πηa), with η the dynamic viscosity. Holding temperature and molecular size roughly fixed, larger viscosity lowers D and therefore lowers the Smoluchowski encounter rate. But solvent changes can also alter reaction free energies, dielectric screening, solvation shells and molecular conformations. A correlation between rate and inverse viscosity supports a diffusion contribution but is not by itself proof that every encounter reacts.

Electrostatic interactions modify the neutral-sphere benchmark. Oppositely charged ions attract and can encounter faster than a neutral-particle estimate under some conditions; like-charged ions repel and can encounter more slowly. Ionic strength changes screening and thus both approach probabilities and activities. A finite reaction probability at contact also matters for electron-transfer quenching and other processes, for which diffusion brings partners together but electronic transfer is not guaranteed. State which effect a fitted “diffusion limit” includes before comparing values across solvents.

At very short times or in confined geometries, the steady three-dimensional formula may need modification. Immediately after a photochemical pulse, reactants can be born close together, making the initial pair distribution nonuniform. On a membrane, diffusion is effectively two-dimensional and the simple 4πR cD expression does not apply unchanged. Crowded biological media can have time-dependent or anomalous diffusion. The classic formula remains a baseline to test more detailed models, not a universal ceiling under every geometry.

To decide whether a reaction is diffusion controlled, compare its observed second-order rate with a plausible encounter estimate under matching conditions. A value close to the estimate suggests high reactivity on contact; a much smaller value suggests an additional activation, orientation or gating constraint. Because R c and D values have uncertainties, near-equality is suggestive rather than definitive. Independent viscosity, temperature or spectroscopic measurements can clarify what happens inside the encounter pair.

Step-by-step reasoning

1. Check whether the reaction is bimolecular in solution and whether a three-dimensional dilute model is appropriate. 2. Estimate D A, D B and a physically justified capture radius R c in consistent SI units. 3. Compute k D,particle = 4πR c(D A+D B). 4. Convert from m³ s⁻¹ per particle concentration to L mol⁻¹ s⁻¹ by multiplying by N A and 1000. 5. Compare observed k with the encounter estimate, considering charge and contact efficiency. 6. Use viscosity or solvent changes as tests while accounting for other solvent-dependent chemistry.

Visual explanation

Draw A at the centre of a sphere of radius R c. Place B at several distances outside it, with wandering Brownian paths. One path reaches the sphere and reacts; another reaches it and diffuses away, showing that contact need not be productive. Beside the sphere draw two rate-limiting sketches: slow diffusion with rapid contact reaction, and rapid diffusion with slow chemical conversion inside the pair. Label the first “encounter controlled” and the second “activation controlled.”

Real-world analogy

Two people in a large crowd must first meet before exchanging an object. If exchange is immediate on meeting, the meeting rate controls the overall exchange. If they frequently meet but cannot exchange until a complex check is completed, the check controls it. The analogy separates encounter from chemical conversion, but molecular diffusion follows statistical motion and the contact radius is a model parameter rather than a visible handshake distance.

Real-world example

In fluorescence quenching, an excited solute can lose its excitation when it encounters a quencher. Researchers measure a bimolecular quenching constant and compare it with the diffusion estimate. A near-diffusion-limited value suggests many encounters quench successfully. A lower value can indicate that electron transfer or energy transfer requires favorable geometry or energetics. Changing solvent viscosity may alter encounter frequency, while changing solvent polarity can also change the electronic driving force, so both effects must be separated.

Why?

Why does the relative diffusion coefficient equal D A+D B in the simple model? Both particles wander independently. Their relative displacement is the difference of their individual displacements, and independent mean-square displacements add. The pair therefore approaches or separates as if one particle diffused toward a fixed target with diffusion coefficient D A+D B. This simplification fails when hydrodynamic or other correlated motions become important.

Common misconception

“A rate constant near 10¹⁰ L mol⁻¹ s⁻¹ is a universal maximum for all solution reactions.” It is an approximate water benchmark for ordinary three-dimensional encounters at a stated temperature, not a fixed natural constant. Solvent viscosity, charge, capture radius and unusual geometries change the estimate. Another error is using k D = 4πRD without converting the particle-based SI units when observed rates are reported in L mol⁻¹ s⁻¹.

Worked example

Suppose R c = 0.50 nm = 5.0 × 10⁻¹⁰ m and D A+D B = 2.0 × 10⁻⁹ m² s⁻¹. Then 4πR cD rel ≈ 4π(5.0 × 10⁻¹⁰)(2.0 × 10⁻⁹) = 1.26 × 10⁻¹⁷ m³ s⁻¹. Multiply by N A = 6.022 × 10²³ mol⁻¹ and 1000 L m⁻³ to obtain k D ≈ 7.6 × 10⁹ L mol⁻¹ s⁻¹. An observed k of 7 × 10⁹ L mol⁻¹ s⁻¹ would be close to this simple encounter estimate, while 7 × 10⁶ would be much lower and suggest a substantial post-encounter limitation or an incorrect capture model.

Quick check

1. If both D A and D B double while R c stays fixed, what happens to the simple k D? Answer: It doubles because D A+D B doubles.

Exam focus

Write D rel = D A+D B and k D = 4πR cD relN A × 1000 for molar concentrations in mol L⁻¹ when SI R c and D are used. Check dimensions before calculating. Distinguish meeting from reacting after meeting and state when the absorbing-boundary assumption is justified. Explain viscosity trends with Stokes–Einstein only under its continuum assumptions.

Advanced insight

The Smoluchowski model solves a diffusion equation with an absorbing boundary at r = R c. Replacing perfect absorption by a partially reactive boundary leads to radiation-boundary or Collins–Kimball treatments, which combine diffusional supply with finite contact reactivity. In charged systems, the diffusion equation can include a potential of mean force, changing the encounter flux. Time-dependent reaction coefficients can arise before a steady concentration profile around each reactant forms. These extensions preserve the core distinction between transport to an encounter and conversion within it.

Summary

A diffusion-controlled solution reaction is limited mainly by how quickly reactants encounter one another. For dilute neutral spheres in three dimensions, the Smoluchowski estimate is proportional to capture radius and the sum of diffusion coefficients. Conversion to molar units is essential, and a water benchmark is around 10¹⁰ L mol⁻¹ s⁻¹ under common conditions. Rates below the encounter estimate can reflect chemical activation, orientation or electronic constraints; solvent and geometry can modify the limit itself.

Practice questions

1. With D A = 1.0 × 10⁻⁹ and D B = 0.5 × 10⁻⁹ m² s⁻¹, what is D rel in the independent-diffusion approximation? Answer: 1.5 × 10⁻⁹ m² s⁻¹.

2. At fixed T and hydrodynamic radius, what does doubling viscosity do to Stokes–Einstein D? Answer: It approximately halves D and thus tends to lower a diffusion-controlled encounter rate.

3. Why is k D = 4πRD not yet in L mol⁻¹ s⁻¹ when R and D are SI values? Answer: It has m³ s⁻¹ per pair-number concentration; multiply by Avogadro's constant and 1000 L m⁻³ for molar concentration units.

4. Does an observed rate far below k D necessarily identify which activation step is slow? Answer: No. It shows the perfect-contact reaction assumption is inadequate, but orientation, solvation and chemical barriers need further evidence to distinguish.

5. Why might a reaction be slower than its encounter rate? Answer: Many encounter pairs may separate before crossing a chemical barrier, achieving a suitable orientation or undergoing electron transfer.