epsilon_dispersion
epsilon_dispersion(epsilon, wavelength = seq(400, 1000), envir = parent.frame())
list
apply a function to a range of wavelength and return dielectric function