Dear Rui, I tried to fit the dielectric function of silicon by the Lorentz oscillator model at: http://fzu.cz/~dominecf/misc/eps/
The exact constants of the Lorentz oscillators can be found in http://fzu.cz/~dominecf/misc/meep/scripts/meep_materials.py which is a file I directly import in my Python-MEEP scripts. You may have to modify it to match your length and time scale. For simplicity, I always use the SI units in my scripts. Try to reduce the number of Lorentz-Drude oscillators as possible within your needed accuracy in the required spectral domain - each additional oscillator consumes CPU time! Some remarks on the materials in MEEP and FDTD stability is also here: http://fzu.cz/~dominecf/misc/meep/#mat_mod Regards, Filip Dominec On Fri, 20 Sep 2013 15:55:50 -0700 Rui La <larui...@gmail.com> wrote: > Dear meep users, > > I want to study the short wavelength light(300-700nm) propagation in Si but > found that the real and imaginary part of refractive index of silicon from > has a huge difference at short wavelength. For example, the real part of > refractive index will have maximum of 6.8 at 370nm corresponding to 46.24 > as dielectric constant. But the refractive index will drop to around 3.5 > after 900nm corresponding to 12. The extinction coefficient is also > changing dramatically under 400nm. So how could I define the materials at > short wavelength. Thank you. > > Best, > > Rui _______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss