------
I'm confused on what happens when i use the following dielectric
function:
(set! default-material
(make dielectric (epsilon 0.34914)
(polarizations
(make polarizability
(omega 3.30345e-22) (gamma 0.111395) (delta-epsilon
8.72259e+43))
(make polarizability
(omega 3.72683) (gamma 37.3009) (delta-epsilon 66.1831))
)))
which is simply a Drude + Lorentz approximation to actual
experimental data (omega, gamma and delta-epsilon were fitted to
match data).
if i launch a gaussian source (f_min = 0.3, f_max=1.2) into a slab of
such material, i find the total EM energy (via 'meep-fields-field-
energy-in-box') appears to grow without bounds.
even if i change the sample material-dispersion.ctl file to have this
dielectric function, it never finishes on the first k-point.
i could understand the apparent gain in energy if gamma was < 0, but
it isn't. i admit i don't like the value of gamma in the Lorentz
term, but alas, this is what fitting to experimental data gives...
thanks,
gp
------
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
