polarizability_ellipsoid(wavelength, epsilon, a = 50, b = 30, c = b, medium = 1.33, kuwata = TRUE)The Kuwata version includes semi-empirical terms of radiative correction and dynamic depolarisation to better match the fully retarded dipolar response in a reasonable range of (subwavelength) sizes and aspect ratios.
alpha_kuwata,
inverse_polarizability, La