Thanks Steven, just to be sure, my system is made of 3 objects, the first is the background material, no polarization. The second and third are made of the same material with polarization parameters. the second is just a cube while the third is a function, described before. As you say, I need one fixed object with the polarization parameters. That would be the second object. But how will meep know that the third object, defined by the dielectric function, is the same? Should I make a fourth object, hidden behind the third and with no size? I have attached a very, very crude figure showing by cell.
As of now, this is how by objects are described:
(define (f p)
(if (> (vector3-y p) (/ (cos (* (vector3-x p) pix2)) (* 2 (/ 1
pitch)))) nbeps xeps))
(set! geometry
(append ; combine lists of
objects:
(list
;-------MEDIUMS----------------------------------------------------------------
(make block (center 0 0 0) (size infinity infinity infinity)
(material (make dielectric (epsilon nbeps)
)))
;-------emitting polymer
layer------------------------------------------------
(make block (center 0 (/ thb -2) 0) (size infinity thb infinity)
(material (make dielectric (epsilon xeps)
(polarizations
(make polarizability
(omega 0.7227) (gamma 0.10) (delta-epsilon De)))
)))
;-------sinusoidal
layer---------------------------------------------------------
(make block (center 0 (/ pitch 2) 0) (size N pitch infinity)
(material (make material-function (epsilon-func f)
(polarizations
(make polarizability
(omega 0.7227) (gamma 0.10) (delta-epsilon De)))
)))
And the output:
Initializing structure...
Working in 2D dimensions.
block, center = (0,0,0)
size (1e+20,1e+20,1e+20)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon = 2.14623
block, center = (0,-0.15,0)
size (1e+20,0.3,1e+20)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon = 3.2761
block, center = (0,0.2,0)
size (1,0.4,1e+20)
axes (1,0,0), (0,1,0), (0,0,1)
time for set_epsilon = 1.4093 s
polarizability: omega=0.7227, gamma=0.1, deps=0.01, esat=0
So it seems the the polarizability parameters don't get attached to any
object in particular? Thanks again.
Nic
<<attachment: cell.png>>
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
