Come calcolare il coefficiente di diffusione con la formula di Stokes-Einstein permette di stimare la mobilità di particelle sferiche disciolte o sospese in un fluido viscoso. La formula fu derivata da Albert Einstein nel 1905 (Annus Mirabilis, Bern, Annalen der Physik 17:549) come applicazione della legge dei grandi numeri al moto browniano osservato da Robert Brown nel 1827, combinando il fattore di drag di Stokes 6πηr (George Stokes, Cambridge 1851) con il rapporto di Einstein D = µ·k_B·T tra mobilità µ e diffusività D. La forma standard è D = k_B·T / (6π·η·r), dove k_B = 1,381×10⁻²³ J/K è la costante di Boltzmann, T la temperatura assoluta, η la viscosità dinamica del fluido in Pa·s, e r il raggio idrodinamico della particella sferica. È il fondamento della Dynamic Light Scattering (DLS), dei modelli di drug delivery con nanoparticelle, della microfluidica, dello studio del moto browniano e della biofisica delle proteine — usata da IIT Genova, IFOM Milano, CNR-IBP Genova, IBSA Lugano, e dei laboratori farmaceutici italiani come Chiesi Parma, Bracco Milano, Menarini Firenze.

1. La formula di Stokes-Einstein e il moto browniano

Il moto browniano, scoperto da Robert Brown nel 1827 osservando granelli di polline al microscopio, fu spiegato quantitativamente da Einstein nel 1905 come effetto delle collisioni termiche delle molecole del fluido sulla particella. Lo spostamento quadratico medio in 3D è ⟨r²⟩ = 6·D·t, e la formula di Stokes-Einstein lega D ai parametri microscopici. Il raggio idrodinamico r include la sfera di solvatazione (in acqua, ~3-5 Å aggiunti al raggio molecolare). Per particelle non sferiche, si introduce il fattore di forma f tramite r_eff = r·f, dove f può variare da 1 (sferica) a 1,5-3 per cilindri o aste. La formula vale solo per particelle Brownian (raggio < 1-10 µm in acqua), per regime diluito (concentrazione < 1% volumetrica), e per fluidi newtoniani (η indipendente dal gradiente di velocità).

2. Coefficienti di diffusione tipici

• Glucosio in acqua a 25°C (r ≈ 0,4 nm) — D ≈ 6,7×10⁻¹⁰ m²/s = 670 µm²/s.
• Albumina BSA in acqua a 37°C (r ≈ 3,5 nm) — D ≈ 6,2×10⁻¹¹ m²/s. Misurato in IIT Genova e IFOM Milano per studi proteomici.
• Virus tabacco TMV (r ≈ 9 nm) — D ≈ 4,3×10⁻¹² m²/s. Riferimento storico per moto browniano biologico.
• Nanoparticella d'oro 20 nm in acqua a 25°C — D ≈ 2,1×10⁻¹¹ m²/s. Standard NanoComposix per calibrazione DLS.
• SARS-CoV-2 in saliva (r ≈ 50 nm) — D ≈ 4,3×10⁻¹² m²/s. Studi su trasmissione aerosol di IIT Genova e MPI Polymerforschung Mainz.
• Liposoma 100 nm in plasma (η = 4 mPa·s) — D ≈ 5×10⁻¹³ m²/s. Drug delivery, IBSA Lugano, Chiesi Parma.
• Cellula leucocita in sangue (r ≈ 7 µm) — D ≈ 3×10⁻¹⁵ m²/s. Solo per scale brevi prima di interazioni cellulari.

3. Limiti di validità e correzioni

La formula classica è strettamente valida per: (a) particelle sferiche rigide; per particelle non sferiche servono fattori di forma f (Perrin 1934, anelli/aste 1,5-3); (b) fluido newtoniano omogeneo; per fluidi non newtoniani (sangue a basso shear, polimeri concentrati, gel) η dipende dal gradiente di velocità e si usa la formula generalizzata di Cox-Merz; (c) regime diluito con ϕ_v < 1%; per concentrazioni elevate domina l'effetto di "crowding" descritto da formule semiempiriche (Phillies, Han); (d) solvent slip nullo; per scale nanometriche entra il fattore di Knudsen quando r ≈ λ_libero (~1 nm in acqua); (e) regime termico classico; sotto T ≈ 100 K ed entrano effetti quantistici. Inoltre la formula non descrive diffusione attiva (chemotaxis, motori molecolari) né diffusione anomala (subdiffusione in cellule, superdiffusione in turbolenza). I codici Gromacs, NAMD, LAMMPS implementano simulazioni MD per casi complessi.

4. Esempi numerici verificati

Esempio 1 — Glucosio in acqua a 25°C, r=0,4 nm. η_acqua_25C = 0,89 mPa·s = 8,9×10⁻⁴ Pa·s. D = (1,381×10⁻²³ · 298) / (6π · 8,9×10⁻⁴ · 4×10⁻¹⁰) = 4,116×10⁻²¹ / 6,71×10⁻¹² = 6,13×10⁻¹⁰ m²/s = 6130 µm²/s. Coerente con dati farmacologici per assorbimento del glucosio nei tessuti.

Esempio 2 — Albumina BSA in acqua a 37°C, r=3,5 nm. η_acqua_37C = 0,69 mPa·s. D = (1,381×10⁻²³ · 310) / (6π · 6,9×10⁻⁴ · 3,5×10⁻⁹) = 4,28×10⁻²¹ / 4,55×10⁻¹¹ = 9,4×10⁻¹¹ m²/s ≈ 9 µm²/s. La BSA è uno standard di calibrazione per DLS Malvern Nano-ZS al IIT Genova e IFOM Milano.

Esempio 3 — Nanoparticella d'oro 20 nm in acqua a 25°C. r=10 nm, η=8,9×10⁻⁴ Pa·s. D = (1,381×10⁻²³ · 298) / (6π · 8,9×10⁻⁴ · 10⁻⁸) = 2,46×10⁻¹¹ m²/s = 24,6 µm²/s. Standard per calibrare strumenti DLS — NanoComposix produce nanoparticelle certificate per fitness DLS, e i laboratori italiani Politecnico Milano, IFOM, CNR-IBF Genova ne usano per validare metodi analitici di drug delivery.

5. Applicazioni pratiche

• Dynamic Light Scattering (DLS): Malvern Nano-ZS, Wyatt DynaPro, Brookhaven 90Plus misurano la diffusione di particelle 1-1000 nm via correlation di intensità Stokes-Einstein. Usato in caratterizzazione di proteine, virus, nanoparticelle, esosomi presso IIT Genova, IFOM Milano, IRCCS Humanitas, Politecnico Milano
• Drug delivery con nanoparticelle: liposomi, polymeric nanoparticles, lipid nanoparticles (LNP) per vaccini mRNA (Pfizer, Moderna) sono caratterizzati via Stokes-Einstein. IBSA Lugano, Chiesi Parma, Bracco Milano, Menarini Firenze sviluppano farmaci nanoincapsulati
• Microfluidica e lab-on-chip: D determina i tempi di mixing e le scale di Péclet — IIT Genova, Politecnico Milano (lab Marco Sampietro), Politecnico Torino, ETH Zurich
• Studi sulle batterie: mobilità di Li⁺ in elettroliti, base degli accumulatori automotive — Stellantis, Enel X, IIT Genova-Genova Smart Cities
• Biofisica e diffusione di proteine in cellule: imaging FRAP (Fluorescence Recovery After Photobleaching) usa Stokes-Einstein modificato per crowding cellulare. IFOM Milano, San Raffaele Milano
• Spettroscopia NMR di diffusione (DOSY): Bruker, Varian usano D per separare molecole per dimensione

6. Confronto con altre formule di trasporto

La formula di Knudsen (1909) descrive diffusione di gas a pressioni basse dove il libero cammino medio è grande rispetto a r: D_Kn = (2/3)·v_th·r. La formula di Smoluchowski (1917) modifica Stokes-Einstein per particelle in interazione (electrostatic, hydrodynamic). La formula di Phillies-Han aggiunge correzioni di concentrazione D(c) = D₀·exp(-α·c^β). Per diffusione attiva (Janus particles, microorganismi), serve trattamento Active Brownian Particle (Berg, Purcell). La diffusione di calore di Fourier usa D_T = k/(ρ·c_p). Per diffusione in solidi domina Arrhenius D = D₀·exp(-E_a/k_B·T). Codici COMSOL Multiphysics, Gromacs, NAMD, LAMMPS, Stokes Einstein simulator (NIST) implementano tutti i regimi.

7. Casi speciali e applicazioni avanzate

Per proteine globulari (BSA, lisozima, hemoglobina), la formula richiede correzione per interazioni elettrostatiche (Linderstrom-Lang). Per polimeri in soluzione con scaling di Zimm o Rouse, D ∝ M^(-3/5) o M^(-1) rispettivamente. Per particelle in plasma sanguigno (η ~3-4 mPa·s), serve correzione per non-Newtonianismo a basso shear (Casson, Quemada). Per diffusione attraverso membrane biologiche entrano fattori di partizione e permeabilità (Renkin). Per aerosol di SARS-CoV-2 (r ~50 nm), Stokes-Einstein in aria umida (η_aria=18 µPa·s) dà D ≈ 5×10⁻¹⁰ m²/s — fondamentale per modelli di trasmissione respiratoria studiati a IIT Genova, Politecnico Torino e MPI Polymerforschung Mainz. Per diffusione subcellulare domina il fenomeno di crowding macromolecolare (effetto Minton), che riduce D effettiva di 10-100× rispetto alla formula classica.

Aggiornato a luglio 2026