pub fn coulomb_logarithm(n_e_m3: f64, t_e_kev: f64) -> FusionResult<f64>
Coulomb logarithm via NRL formula, clamped to [5, 30].