On Sat, 6 Jan 2007, miles wrote:
Thanks for your last reply! however, i am still confused by the
problem that material-function can not return polarizations. My code is
as below.However, whether sigma=0 or sigma=1, the results are the same.
Would you give some instruction, please?
========================================================================================
(define (epsilon-func vCartP)
(if (>= (dIntensity vCartP) 9.65)
(if srr? (make dielectric (epsilon 1)
(polarizations (make polarizability (omega 1e-30)
(gamma 0.0217) (delta-epsilon (* 7.25 7.25 1e+60)) (sigma 1) ))
)
(make dielectric (epsilon 1))
)
(make dielectric (epsilon dEps-lo))
))
(set! default-material (make material-function (material-func epsilon-func)))
=========================================================================================
The point is that you need some ordinary object somewhere, even if it is
outside the computational cell, with that
polarization material. e.g. try adding:
(set! geometry (list (make block (size 0 0 0) (center infinity)
(material (make dielectric (epsilon 1)
(polarizations (make polarizability (omega
1e-30) (gamma 0.0217) (delta-epsilon (* 7.25 7.25 1e+60)) (sigma 1) ))))))
as a "dummy" object to tell Meep to include that polarization.
In the next version of Meep, I should probably just add a new input
variable "extra-materials" or something like that to provide a simpler way
to "declare" dispersive materials that will be used in the computational
cell via material-functions.
What's more, I have another question: how to obtain the complex
transmission and reflection coefficients over a flux region for a
guassian source?
Use the add-flux commands, as explained in the Meep tutorial, to specify
regions where you want to compute flux spectra.
Cordially,
Steven G. Johnson
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
