Dear Dr. Steven and meep users,
I want to simulate a slab structure (without defect for the moment)
with square lattice. I try to obtain the band diagram (I am
particularly interested of the part above the light line).
;my file
(define-param r 0.3)
(define-param eps 11.36)
(define-param loweps 1.0)
(define-param dpml 1)
(define-param supercell-h 10)
;dimensions of my slab:
(define-param h 0.31)
(define-param sx 5)
(define-param sy 5)
(set! eps-averaging? false)
(set! default-material (make dielectric (epsilon loweps)))
(set! geometry-lattice (make lattice (size 10 10 supercell-h)))
(set! geometry
(list (make block (material (make dielectric (epsilon eps)))
(center 0) (size sx sy h))
(make cylinder
(center 0 0 0) (radius r) (height h)
(material air))))
(set! geometry
(geometric-objects-lattice-duplicates
(list
(make cylinder
(center 0 0 0) (radius r) (height h)
(material air)))
1 1 supercell-h))
;(set! pml-layers (list (make pml (direction Z) (thickness dpml))))
(set! pml-layers (list (make pml (thickness dpml))))
(define-param fcen 0.29)
(define-param df 0.1)
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen)(fwidth df))) (component
Ez) (center 0 0 (+ dpml (* -0.5 supercell-h))) )))
(set! resolution 8)
(define-param k-interp 9)
(run-k-points 150 (interpolate k-interp (list (vector3 0 0 0) (vector3
0.5 0 0))))
When I put pml layers in the three directions, the result is zero (w=0).
When I calculate with only a pml layer in z direction, I obtain some
value of w but if I executed h5topng -0z0 I got only black long
rectangle.
I don't find what is wrong in my file. I don't understand why it's not
working when I put pml layers in x and y direction even far from my
structure.
thanks
Regards,
Elise
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
