pub fn debye_length(temp_ev: f32, density: f32) -> f32
Debye length: lambda_D = sqrt(eps0 * kT / (n * q^2)) In normalized units with eps0=1, kT in eV, n = number density.