Hello,
Thank you for answering.
I use meep because I am interested of the part of the band diagram
above the light line. I have already try to use mpb to simulate my
slab structure but the part above the light line is too unstable.
Moreover I found on the site of mpb discuss a message of Steven which
told to use meep in this kind of situation ("a program with absorbing
boundary support, like Meep, is much more appropriate for computing
leaky modes" see mpb discuss 13 Aug 2008).
I try to simulate a square slab with 5*5*0.31 dimension with holes. I
have set a geometry to (10 10 10) because first I would like to add
pml layer in the three direction (so I need space).
I didn't know that I can just simulate a (1 1 10) geometry with only
one hole. It's the same as a (5 5 10) geometry with multiple holes?
I will try this, thanks.
It will be great if you can send me a .ctl file example
Regards
Elise
mahmoud hamoud <[email protected]> a écrit :
Hi,
Dr. Steven suggest mpb for computing the band structures, also
your structure consist a hole slab with square lattice, i don't
know why you set your geometry to (10 10 10), you have to set it
only as (1 1 10).
2. in the case when you want to introduce a defect (point defect or
linear defect), you have to set your geometry to (1 6 10).
i can send you a ctl.file example if you want.
hope this help you
with regards
--- On Wed, 7/7/10, [email protected]
<[email protected]> wrote:
From: [email protected] <[email protected]>
Subject: [Meep-discuss] simple square slab
To: [email protected]
Date: Wednesday, July 7, 2010, 1:04 AM
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
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
