Dear Thomas and all, I wrote most of my thoughts on the referred page; while sometimes the discrete Lorentzian oscillators suffice to precisely describe the dielectric function in the frequency range of interest, the optical spectra of semiconductors are usually of the harder cases.
You can manually fit the eps' and eps'' spectrum by summing more narrower oscillators. This takes human time, and each oscillator also requires CPU time in each FDTD step, so it is always a tradeoff. The good point is that if you are interested in photovoltaics, your simulations may be perhaps restricted to the wavelengths of 1200--400 nm, where silicon seems to behave in a relatively cultivated way. DominecF 2014-05-19 18:15 GMT+02:00, Thomas LOPEZ <thomas.lo...@isae.fr>: > Dear Meep users, > > I would like to implement in MEEP this model of Silicon permittivity > (see below, from Deinega, A. & John, S. *Effective optical response of > silicon to sunlight in the finite-difference time-domain method*. > /Optics letters/ *37,* 112--114 (2012)): > > This expression is from multiple Debye, Drude and Lorentz terms and is > sufficient to describe the response of bulk silicon to 300nm to 1000nm... > I found some information about modeling Silicon here > http://fzu.cz/~dominecf/misc/eps/ and also in Meep tutorial. > Unfortunately, the imaginary part of epsilon is often far from tabulated > data. > > So do you have any idea how to implement this expression? Are there any > "tricks" using polarization class? > > Many thanks in advance for your assistance. > > > -- > Thomas Lopez > Doctorant/PhD Student > Département Electronique Optronique et Signal (DEOS) > Institut Supérieur de l'Aeronautique et de l'Espace (ISAE) > 10 Avenue Edouard Belin BP 54032 > 31055 Toulouse CEDEX 4 FRANCE > www.isae.fr > > _______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
