pub fn plasma_frequency(density: f32, mass: f32) -> f32
Plasma frequency: omega_p = sqrt(n * q^2 / (eps0 * m))