pub fn einstein_diffusivity(mobility: f64, t_kelvin: f64) -> f64
Einstein relation: D = μ k_B T / q (m²/s).